Monday, 12 May 2014

Solving an equation with Lambert's W function? Or by any other means?




I am trying to solve the following equation for x in terms of y and c (with x,y[0,1])
log(x1xy)+x1xy+y1xy=c



I can solve this easier equation
log(x1xy)+x1xy=c




Let
z=x1xy



Then I can solve for x using Lambert's W function
log(z)+z=cz=exp(c)exp(z)zexp(z)=exp(c)z=W(exp(c))x=(1y)W(exp(c))1+W(exp(c))



Can anyone help me solve the harder equation? Is Lambert's W function helpful here?



Thanks!


Answer



ln(x1xy)+x1xy+y1xy=c



ln(x1xy)=c+1xy11xy,ln(1y1xy1)=c+111xy,ln(1y1xy1)=c(11xy1),ln(1y1xy1)=c11y(1y1xy(1y)),ln(1y1xy1)=c11y(1y1xy(1y)),ln(1y1xy1)=cy1y11y(1y1xy1).



Let



1y1xy1=z,cy1y=u,11y=v



and we have an equation



lnz=u+vz,



which has a standard solution for z in terms Lambert W function



z=W(vexp(u))v,x=(1y)(111+z).


No comments:

Post a Comment

real analysis - How to find limhrightarrow0fracsin(ha)h

How to find limh0sin(ha)h without lhopital rule? I know when I use lhopital I easy get $$ \lim_{h\rightarrow 0}...