2
$\begingroup$

I want to find the maximum of the following function,

$$ f(x, y) = e^{m e^{-x}+n e^{-y}-x-y}(mrxe^y+nsye^x+mn(r+s)xy), 0 \le x, y \le 1 $$

where $r, s, m,$ and $n$ are positive constants. At the maximum point we have the following equations,

$$ \begin{cases} \frac{\partial{f(x, y)}}{\partial{x}} = 0 \\ \frac{\partial{f(x, y)}}{\partial{y}} = 0 \end{cases} $$

Simplifying the above equations we will get,

$$ \begin{cases} x = -W_0\left(\frac {s(ny+ye^y-e^y)} {m(nry+nsy+rye^y+sye^y-se^y)}\right) \\ y = -W_0\left(\frac {r(mx+xe^x-e^x)} {n(mrx+msx+rxe^x+sxe^x-re^x)}\right) \end{cases} $$

Where $W_0(.)$ is the upper branch of the Lambert W function. We know that $W_0(.)$ cannot be expressed in terms of elementary functions. Does it imply that the system of equations cannot be solved analytically? If not, how can I solve it analytically?

1 Answers 1

1

Hint: One possible solution is using complex analysis, by substituting $z=x+iy$ (and then $x=\frac{z+\overline{z}}{2}$ and $y=\frac{z-\overline{z}}{2}$), you can get a simpler form (using, for example, mathematica). Hopefully, the function would be analytic in $\{z|0\leq\Re(z),\Im(z)\leq 1\}$. Then the maximum would be on the boundary, which is relatively easy to calculate (you get a function of one variable).

  • 0
    @Dennis: I don't know why you have the hope that the function turns out to be holomorphic.2011-05-15
  • 0
    @Fabian: I'm just hoping for the best :)2011-05-15
  • 0
    Hi guys! I edited the question. Could you please let me know your opinion about it now. Thank you.2011-05-16
  • 0
    @Mohsen: Now it's clear that my answer is useless in this situation :(2011-05-16
  • 1
    Nevertheless, the analytic method means that you aren't approximating numerically the extrema points but using the derivative to find the suspected points for extrema. The fact that you got the Lambert W function means that here you will have to solve numerically from now on. Which means that there is no purely analytic solution, but you can use a not purely numeric solution.2011-05-16
  • 0
    All right. It's not clear for me yet, why there is no analytical solution for the system. How can I prove that $x$ and $y$ cannot be shown in terms of $m, n, r$ and $s$ using elementary functions? i.e., one may find another way to drive $x$ and $y$ without the $W$ function. Do you see any contradiction?2011-05-17