Question

请教各位大虾用MATLAB解方程 dy/dx=[(d-x)/y-ya/(x^2+y^2-xd)]/[1+(d-x)b/(x^2+y^2-xd)] 其中a b d都是常

不好意思，小弟没分，请各位大虾帮忙

Answer 1

这样肯定解不出来，就是解出来 也看不懂

Answer 2

>> dsolve('Dy=((d-x)/y-y*a/(x^2+y^2-x*d))/(1+(d-x)*b/(x^2+y^2-x*d))')

ans =

1/(-d+x+a)*((2*d*x^2-x*d^2-x^3+exp((lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*b*d^2-2*lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*d*b*x+lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*b*x^2-lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*a*x^2+lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*a*x*d-lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*a*b*d+lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*a*b*x-2*d*x^2+x*d^2+x^3-2*t*d^2+4*t*x*d+4*t*a*d-2*t*x^2-4*t*a*x-2*t*a^2-2*C1*d^2+4*C1*x*d+4*C1*a*d-2*C1*x^2-4*C1*a*x-2*C1*a^2)/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)))*(-d+x+a))^(1/2)
-1/(-d+x+a)*((d-x-a)*(-2*d*x^2+x*d^2+x^3-exp((lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*b*d^2-2*lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*d*b*x+lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*b*x^2-lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*a*x^2+lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*a*x*d-lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*a*b*d+lambertw(1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*exp(-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*d*x^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x*d^2+1/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*x^3-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*t*a^2-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*d^2+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x*d+4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*d-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*x^2-4/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a*x-2/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x)*C1*a^2))*a*b*x-2*d*x^2+x*d^2+x^3-2*t*d^2+4*t*x*d+4*t*a*d-2*t*x^2-4*t*a*x-2*t*a^2-2*C1*d^2+4*C1*x*d+4*C1*a*d-2*C1*x^2-4*C1*a*x-2*C1*a^2)/(-b*d^2+2*d*b*x-b*x^2+a*x^2-a*x*d+a*b*d-a*b*x))))^(1/2)

>>