We could use the Lagrange multiplier L :
f(x,y,L)=yexp(x)−xexp(y)+L(xy−4)
dfdx=yexp(x)−exp(y)+Ly=0
dfdy=exp(x)−xexp(y)+Lx=0
dfdL=xy−4=0⇒y=4x
⇒(4x)exp(x)−exp(4x)+4Lx=0
⇒exp(x)−xexp(4x)+Lx=0
Multiply the last equation by (4/x) :
⇒(4x)exp(x)−4exp(4x)+4L=0
Subtract this equation from the first :
⇒3exp(4x)+4Lx−4L=0
⇒4L(1−1x)=3exp(4x)
⇒L=(34)exp(4x)1−1x
Fill in this value for L in the second equation :
⇒exp(x)−xexp(4x)+(34)exp(4x)x2x−1=0
⇒exp(x)+exp(4x)[(34)x2x−1−x]=0
⇒exp(x)=exp(4x)x(x−1)−(34)x2x−1
⇒exp(x)=exp(4x)x(14)x−1x−1
⇒exp(x−4x)=x(14)x−1x−1
This looks like a transcendental equation.
I am stopping here, but if one solves this numerically, one gets
x and then L and also y.