Maple 2024 Questions and Posts

These are Posts and Questions associated with the product, Maple 2024

I am not sure how to use dsolve for my problem.
CQ_v1.mw

I'm trying to extract the real and imaginary parts of a given PDE. However, during the substitution step, something unexpected occurs. Specifically, my replacement does not yield the expected result, and an extra term appears: D(k)(t*w + x) in the expression P11. I'm unsure why this term arises and would appreciate any insight into what might be going wrong during the substitution process.

restart

with(PDEtools)

with(LinearAlgebra)

NULL

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

declare(u(x, t)); declare(U(xi)); declare(V(xi))

u(x, t)*`will now be displayed as`*u

 

U(xi)*`will now be displayed as`*U

 

V(xi)*`will now be displayed as`*V

(2)

pde := I*(diff(u(x, t), `$`(t, 2))-s^2*(diff(u(x, t), `$`(x, 2))))+(1/24)*c[1]*(diff(u(x, t), t, `$`(x, 4)))-alpha*s*c[1]*(diff(u(x, t), `$`(x, 5)))+diff(c[2]*u(x, t)*U(-t*v+x)^2+c[3]*u(x, t)*U(-t*v+x)^4+c[4]*u(x, t)*(diff(U(-t*v+x)^2, `$`(x, 2))), t)-beta*s*(diff(u(x, t), `$`(x, 5)))+diff(c[2]*u(x, t)*U(-t*v+x)^2+c[3]*u(x, t)*U(-t*v+x)^4+c[4]*u(x, t)*(diff(U(-t*v+x)^2, `$`(x, 2))), x)

I*(diff(diff(u(x, t), t), t)-s^2*(diff(diff(u(x, t), x), x)))+(1/24)*c[1]*(diff(diff(diff(diff(diff(u(x, t), t), x), x), x), x))-alpha*s*c[1]*(diff(diff(diff(diff(diff(u(x, t), x), x), x), x), x))+c[2]*(diff(u(x, t), t))*U(-t*v+x)^2-2*c[2]*u(x, t)*U(-t*v+x)*(D(U))(-t*v+x)*v+c[3]*(diff(u(x, t), t))*U(-t*v+x)^4-4*c[3]*u(x, t)*U(-t*v+x)^3*(D(U))(-t*v+x)*v+c[4]*(diff(u(x, t), t))*(2*(D(U))(-t*v+x)^2+2*U(-t*v+x)*((D@@2)(U))(-t*v+x))+c[4]*u(x, t)*(-6*(D(U))(-t*v+x)*((D@@2)(U))(-t*v+x)*v-2*U(-t*v+x)*((D@@3)(U))(-t*v+x)*v)-beta*s*(diff(diff(diff(diff(diff(u(x, t), x), x), x), x), x))+c[2]*(diff(u(x, t), x))*U(-t*v+x)^2+2*c[2]*u(x, t)*U(-t*v+x)*(D(U))(-t*v+x)+c[3]*(diff(u(x, t), x))*U(-t*v+x)^4+4*c[3]*u(x, t)*U(-t*v+x)^3*(D(U))(-t*v+x)+c[4]*(diff(u(x, t), x))*(2*(D(U))(-t*v+x)^2+2*U(-t*v+x)*((D@@2)(U))(-t*v+x))+c[4]*u(x, t)*(6*(D(U))(-t*v+x)*((D@@2)(U))(-t*v+x)+2*U(-t*v+x)*((D@@3)(U))(-t*v+x))

(3)

``

G1 := U(-t*v+x) = U(xi); G2 := (D(U))(-t*v+x) = diff(U(xi), xi); G3 := ((D@@2)(U))(-t*v+x) = diff(U(xi), `$`(xi, 2)); G4 := ((D@@3)(U))(-t*v+x) = diff(U(xi), `$`(xi, 3)); G5 := ((D@@4)(U))(-t*v+x) = diff(U(xi), `$`(xi, 4)); G6 := ((D@@5)(U))(-t*v+x) = diff(U(xi), `$`(xi, 5))

U(-t*v+x) = U(xi)

 

(D(U))(-t*v+x) = diff(U(xi), xi)

 

((D@@2)(U))(-t*v+x) = diff(diff(U(xi), xi), xi)

 

((D@@3)(U))(-t*v+x) = diff(diff(diff(U(xi), xi), xi), xi)

 

((D@@4)(U))(-t*v+x) = diff(diff(diff(diff(U(xi), xi), xi), xi), xi)

 

((D@@5)(U))(-t*v+x) = diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi)

(4)

NULL

T := xi = -t*v+x; T1 := u(x, t) = U(-t*v+x)*exp(I*k(t*w+x))

xi = -t*v+x

 

u(x, t) = U(-t*v+x)*exp(I*k(t*w+x))

(5)

P1 := I*(diff(u(x, t), `$`(t, 2))-s^2*(diff(u(x, t), `$`(x, 2))))+(1/24)*c[1]*(diff(u(x, t), t, `$`(x, 4)))-alpha*s*c[1]*(diff(u(x, t), `$`(x, 5)))+diff(c[2]*u(x, t)*U(-t*v+x)^2+c[3]*u(x, t)*U(-t*v+x)^4+c[4]*u(x, t)*(diff(U(-t*v+x)^2, `$`(x, 2))), t)-beta*s*(diff(u(x, t), `$`(x, 5)))+diff(c[2]*u(x, t)*U(-t*v+x)^2+c[3]*u(x, t)*U(-t*v+x)^4+c[4]*u(x, t)*(diff(U(-t*v+x)^2, `$`(x, 2))), x)

I*(diff(diff(u(x, t), t), t)-s^2*(diff(diff(u(x, t), x), x)))+(1/24)*c[1]*(diff(diff(diff(diff(diff(u(x, t), t), x), x), x), x))-alpha*s*c[1]*(diff(diff(diff(diff(diff(u(x, t), x), x), x), x), x))+c[2]*(diff(u(x, t), t))*U(-t*v+x)^2-2*c[2]*u(x, t)*U(-t*v+x)*(D(U))(-t*v+x)*v+c[3]*(diff(u(x, t), t))*U(-t*v+x)^4-4*c[3]*u(x, t)*U(-t*v+x)^3*(D(U))(-t*v+x)*v+c[4]*(diff(u(x, t), t))*(2*(D(U))(-t*v+x)^2+2*U(-t*v+x)*((D@@2)(U))(-t*v+x))+c[4]*u(x, t)*(-6*(D(U))(-t*v+x)*((D@@2)(U))(-t*v+x)*v-2*U(-t*v+x)*((D@@3)(U))(-t*v+x)*v)-beta*s*(diff(diff(diff(diff(diff(u(x, t), x), x), x), x), x))+c[2]*(diff(u(x, t), x))*U(-t*v+x)^2+2*c[2]*u(x, t)*U(-t*v+x)*(D(U))(-t*v+x)+c[3]*(diff(u(x, t), x))*U(-t*v+x)^4+4*c[3]*u(x, t)*U(-t*v+x)^3*(D(U))(-t*v+x)+c[4]*(diff(u(x, t), x))*(2*(D(U))(-t*v+x)^2+2*U(-t*v+x)*((D@@2)(U))(-t*v+x))+c[4]*u(x, t)*(6*(D(U))(-t*v+x)*((D@@2)(U))(-t*v+x)+2*U(-t*v+x)*((D@@3)(U))(-t*v+x))

(6)

P11 := eval(P1, {T, T1})

I*(((D@@2)(U))(-t*v+x)*v^2*exp(I*k(t*w+x))-(2*I)*(D(U))(-t*v+x)*v*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(-t*v+x)*((D@@2)(k))(t*w+x)*w^2*exp(I*k(t*w+x))-U(-t*v+x)*(D(k))(t*w+x)^2*w^2*exp(I*k(t*w+x))-s^2*(((D@@2)(U))(-t*v+x)*exp(I*k(t*w+x))+(2*I)*(D(U))(-t*v+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(-t*v+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-U(-t*v+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))))+(1/24)*c[1]*(6*((D@@3)(U))(-t*v+x)*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-4*((D@@3)(U))(-t*v+x)*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))+4*(D(U))(-t*v+x)*(D(k))(t*w+x)^4*w*exp(I*k(t*w+x))-(D(U))(-t*v+x)*v*(D(k))(t*w+x)^4*exp(I*k(t*w+x))+I*U(-t*v+x)*(D(k))(t*w+x)^5*w*exp(I*k(t*w+x))-(24*I)*(D(U))(-t*v+x)*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*U(-t*v+x)*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(-t*v+x)*((D@@2)(k))(t*w+x)^2*w*(D(k))(t*w+x)*exp(I*k(t*w+x))+(6*I)*(D(U))(-t*v+x)*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+12*((D@@2)(U))(-t*v+x)*v*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-18*((D@@2)(U))(-t*v+x)*(D(k))(t*w+x)*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+10*U(-t*v+x)*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+I*((D@@4)(U))(-t*v+x)*(D(k))(t*w+x)*w*exp(I*k(t*w+x))-16*(D(U))(-t*v+x)*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-12*(D(U))(-t*v+x)*((D@@2)(k))(t*w+x)^2*w*exp(I*k(t*w+x))-5*U(-t*v+x)*((D@@4)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-10*U(-t*v+x)*((D@@3)(k))(t*w+x)*w*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(-t*v+x)*((D@@5)(k))(t*w+x)*w*exp(I*k(t*w+x))-I*(D(U))(-t*v+x)*v*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))+3*(D(U))(-t*v+x)*v*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+(4*I)*((D@@3)(U))(-t*v+x)*((D@@2)(k))(t*w+x)*w*exp(I*k(t*w+x))-(6*I)*((D@@2)(U))(-t*v+x)*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))-(6*I)*((D@@3)(U))(-t*v+x)*v*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*((D@@2)(U))(-t*v+x)*v*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+(6*I)*((D@@2)(U))(-t*v+x)*((D@@3)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(D(U))(-t*v+x)*((D@@4)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*((D@@2)(U))(-t*v+x)*v*(D(k))(t*w+x)^3*exp(I*k(t*w+x))-(4*I)*((D@@4)(U))(-t*v+x)*v*(D(k))(t*w+x)*exp(I*k(t*w+x))+4*(D(U))(-t*v+x)*v*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-((D@@5)(U))(-t*v+x)*v*exp(I*k(t*w+x)))-alpha*s*c[1]*(I*U(-t*v+x)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))+(10*I)*((D@@2)(U))(-t*v+x)*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*((D@@2)(U))(-t*v+x)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))-(30*I)*(D(U))(-t*v+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*((D@@4)(U))(-t*v+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-10*U(-t*v+x)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(D(U))(-t*v+x)*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(-t*v+x)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(D(U))(-t*v+x)*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(-t*v+x)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(D(U))(-t*v+x)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*((D@@2)(U))(-t*v+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*U(-t*v+x)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(-t*v+x)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(-t*v+x)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))+(10*I)*((D@@3)(U))(-t*v+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+((D@@5)(U))(-t*v+x)*exp(I*k(t*w+x))+5*(D(U))(-t*v+x)*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*((D@@3)(U))(-t*v+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*(-(D(U))(-t*v+x)*v*exp(I*k(t*w+x))+I*U(-t*v+x)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(-t*v+x)^2-2*c[2]*U(-t*v+x)^2*exp(I*k(t*w+x))*(D(U))(-t*v+x)*v+c[3]*(-(D(U))(-t*v+x)*v*exp(I*k(t*w+x))+I*U(-t*v+x)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(-t*v+x)^4-4*c[3]*U(-t*v+x)^4*exp(I*k(t*w+x))*(D(U))(-t*v+x)*v+c[4]*(-(D(U))(-t*v+x)*v*exp(I*k(t*w+x))+I*U(-t*v+x)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*(2*(D(U))(-t*v+x)^2+2*U(-t*v+x)*((D@@2)(U))(-t*v+x))+c[4]*U(-t*v+x)*exp(I*k(t*w+x))*(-6*(D(U))(-t*v+x)*((D@@2)(U))(-t*v+x)*v-2*U(-t*v+x)*((D@@3)(U))(-t*v+x)*v)-beta*s*(I*U(-t*v+x)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))+(10*I)*((D@@2)(U))(-t*v+x)*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*((D@@2)(U))(-t*v+x)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))-(30*I)*(D(U))(-t*v+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*((D@@4)(U))(-t*v+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-10*U(-t*v+x)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(D(U))(-t*v+x)*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(-t*v+x)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(D(U))(-t*v+x)*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(-t*v+x)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(D(U))(-t*v+x)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*((D@@2)(U))(-t*v+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*U(-t*v+x)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(-t*v+x)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(-t*v+x)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))+(10*I)*((D@@3)(U))(-t*v+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+((D@@5)(U))(-t*v+x)*exp(I*k(t*w+x))+5*(D(U))(-t*v+x)*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*((D@@3)(U))(-t*v+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*((D(U))(-t*v+x)*exp(I*k(t*w+x))+I*U(-t*v+x)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(-t*v+x)^2+2*c[2]*U(-t*v+x)^2*exp(I*k(t*w+x))*(D(U))(-t*v+x)+c[3]*((D(U))(-t*v+x)*exp(I*k(t*w+x))+I*U(-t*v+x)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(-t*v+x)^4+4*c[3]*U(-t*v+x)^4*exp(I*k(t*w+x))*(D(U))(-t*v+x)+c[4]*((D(U))(-t*v+x)*exp(I*k(t*w+x))+I*U(-t*v+x)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*(2*(D(U))(-t*v+x)^2+2*U(-t*v+x)*((D@@2)(U))(-t*v+x))+c[4]*U(-t*v+x)*exp(I*k(t*w+x))*(6*(D(U))(-t*v+x)*((D@@2)(U))(-t*v+x)+2*U(-t*v+x)*((D@@3)(U))(-t*v+x))

(7)

P111 := subs({G1, G2, G3, G4, G5, G6}, P11)NULL

I*((diff(diff(U(xi), xi), xi))*v^2*exp(I*k(t*w+x))-(2*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*w^2*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*w^2*exp(I*k(t*w+x))-s^2*((diff(diff(U(xi), xi), xi))*exp(I*k(t*w+x))+(2*I)*(diff(U(xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))))+(1/24)*c[1]*(6*(diff(diff(diff(U(xi), xi), xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-4*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*(D(k))(t*w+x)^4*w*exp(I*k(t*w+x))-(diff(U(xi), xi))*v*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-12*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*w*exp(I*k(t*w+x))+3*(diff(U(xi), xi))*v*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*w*exp(I*k(t*w+x))+12*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-18*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+10*U(xi)*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+I*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))-16*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*w*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*w*exp(I*k(t*w+x))-I*(diff(U(xi), xi))*v*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))+(4*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*w*exp(I*k(t*w+x))-(6*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))-(6*I)*(diff(diff(diff(U(xi), xi), xi), xi))*v*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*(diff(diff(U(xi), xi), xi))*v*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+(6*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)^3*exp(I*k(t*w+x))-(4*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*v*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*v*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*w*(D(k))(t*w+x)*exp(I*k(t*w+x))+(6*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(24*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x))-alpha*s*c[1]*(-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^2-2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*v*c[2]+c[3]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^4-4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*v*c[3]+c[4]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(-6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))*v-2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))*v)-beta*s*(-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^2+2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*c[2]+c[3]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^4+4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*c[3]+c[4]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))+2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi)))

(8)

``

pde1 := P111 = 0

I*((diff(diff(U(xi), xi), xi))*v^2*exp(I*k(t*w+x))-(2*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*w^2*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*w^2*exp(I*k(t*w+x))-s^2*((diff(diff(U(xi), xi), xi))*exp(I*k(t*w+x))+(2*I)*(diff(U(xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))))+(1/24)*c[1]*(6*(diff(diff(diff(U(xi), xi), xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-4*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*(D(k))(t*w+x)^4*w*exp(I*k(t*w+x))-(diff(U(xi), xi))*v*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-12*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*w*exp(I*k(t*w+x))+3*(diff(U(xi), xi))*v*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*w*exp(I*k(t*w+x))+12*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-18*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+10*U(xi)*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+I*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))-16*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*w*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*w*exp(I*k(t*w+x))-I*(diff(U(xi), xi))*v*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))+(4*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*w*exp(I*k(t*w+x))-(6*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))-(6*I)*(diff(diff(diff(U(xi), xi), xi), xi))*v*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*(diff(diff(U(xi), xi), xi))*v*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+(6*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)^3*exp(I*k(t*w+x))-(4*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*v*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*v*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*w*(D(k))(t*w+x)*exp(I*k(t*w+x))+(6*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(24*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x))-alpha*s*c[1]*(-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^2-2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*v*c[2]+c[3]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^4-4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*v*c[3]+c[4]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(-6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))*v-2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))*v)-beta*s*(-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^2+2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*c[2]+c[3]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^4+4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*c[3]+c[4]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))+2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))) = 0

(9)

numer(lhs(I*((diff(diff(U(xi), xi), xi))*v^2*exp(I*k(t*w+x))-(2*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*w^2*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*w^2*exp(I*k(t*w+x))-s^2*((diff(diff(U(xi), xi), xi))*exp(I*k(t*w+x))+(2*I)*(diff(U(xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))))+(1/24)*c[1]*(-(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*v*exp(I*k(t*w+x))+6*(diff(diff(diff(U(xi), xi), xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-4*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*(D(k))(t*w+x)^4*w*exp(I*k(t*w+x))-12*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*w*exp(I*k(t*w+x))+3*(diff(U(xi), xi))*v*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))-(diff(U(xi), xi))*v*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*w*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*v*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-I*(diff(U(xi), xi))*v*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(6*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))-(6*I)*(diff(diff(diff(U(xi), xi), xi), xi))*v*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*(diff(diff(U(xi), xi), xi))*v*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*w*exp(I*k(t*w+x))+I*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*w*exp(I*k(t*w+x))+(6*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+12*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-18*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+10*U(xi)*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-16*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))+(6*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-(24*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x))-alpha*s*c[1]*(-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^2-2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*v*c[2]+c[3]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^4-4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*v*c[3]+c[4]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(-6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))*v-2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))*v)-beta*s*(-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^2+2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*c[2]+c[3]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^4+4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*c[3]+c[4]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))+2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))) = 0))*denom(rhs(I*((diff(diff(U(xi), xi), xi))*v^2*exp(I*k(t*w+x))-(2*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*w^2*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*w^2*exp(I*k(t*w+x))-s^2*((diff(diff(U(xi), xi), xi))*exp(I*k(t*w+x))+(2*I)*(diff(U(xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))))+(1/24)*c[1]*(-(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*v*exp(I*k(t*w+x))+6*(diff(diff(diff(U(xi), xi), xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-4*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*(D(k))(t*w+x)^4*w*exp(I*k(t*w+x))-12*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*w*exp(I*k(t*w+x))+3*(diff(U(xi), xi))*v*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))-(diff(U(xi), xi))*v*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*w*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*v*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-I*(diff(U(xi), xi))*v*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(6*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))-(6*I)*(diff(diff(diff(U(xi), xi), xi), xi))*v*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*(diff(diff(U(xi), xi), xi))*v*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*w*exp(I*k(t*w+x))+I*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*w*exp(I*k(t*w+x))+(6*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+12*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-18*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+10*U(xi)*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-16*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))+(6*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-(24*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x))-alpha*s*c[1]*(-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^2-2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*v*c[2]+c[3]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^4-4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*v*c[3]+c[4]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(-6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))*v-2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))*v)-beta*s*(-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^2+2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*c[2]+c[3]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^4+4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*c[3]+c[4]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))+2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))) = 0)) = numer(rhs(I*((diff(diff(U(xi), xi), xi))*v^2*exp(I*k(t*w+x))-(2*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*w^2*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*w^2*exp(I*k(t*w+x))-s^2*((diff(diff(U(xi), xi), xi))*exp(I*k(t*w+x))+(2*I)*(diff(U(xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))))+(1/24)*c[1]*(-(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*v*exp(I*k(t*w+x))+6*(diff(diff(diff(U(xi), xi), xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-4*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*(D(k))(t*w+x)^4*w*exp(I*k(t*w+x))-12*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*w*exp(I*k(t*w+x))+3*(diff(U(xi), xi))*v*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))-(diff(U(xi), xi))*v*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*w*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*v*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-I*(diff(U(xi), xi))*v*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(6*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))-(6*I)*(diff(diff(diff(U(xi), xi), xi), xi))*v*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*(diff(diff(U(xi), xi), xi))*v*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*w*exp(I*k(t*w+x))+I*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*w*exp(I*k(t*w+x))+(6*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+12*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-18*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+10*U(xi)*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-16*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))+(6*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-(24*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x))-alpha*s*c[1]*(-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^2-2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*v*c[2]+c[3]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^4-4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*v*c[3]+c[4]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(-6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))*v-2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))*v)-beta*s*(-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^2+2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*c[2]+c[3]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^4+4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*c[3]+c[4]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))+2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))) = 0))*denom(lhs(I*((diff(diff(U(xi), xi), xi))*v^2*exp(I*k(t*w+x))-(2*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*w^2*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*w^2*exp(I*k(t*w+x))-s^2*((diff(diff(U(xi), xi), xi))*exp(I*k(t*w+x))+(2*I)*(diff(U(xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-U(xi)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))))+(1/24)*c[1]*(-(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*v*exp(I*k(t*w+x))+6*(diff(diff(diff(U(xi), xi), xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-4*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*(D(k))(t*w+x)^4*w*exp(I*k(t*w+x))-12*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*w*exp(I*k(t*w+x))+3*(diff(U(xi), xi))*v*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))-(diff(U(xi), xi))*v*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*w*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))+4*(diff(U(xi), xi))*v*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-I*(diff(U(xi), xi))*v*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(6*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))-(6*I)*(diff(diff(diff(U(xi), xi), xi), xi))*v*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*(diff(diff(U(xi), xi), xi))*v*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))-(4*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*w*exp(I*k(t*w+x))+I*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*w*exp(I*k(t*w+x))+(6*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*w*exp(I*k(t*w+x))+(4*I)*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+12*(diff(diff(U(xi), xi), xi))*v*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-18*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+10*U(xi)*(D(k))(t*w+x)^3*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-16*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*w*(D(k))(t*w+x)*exp(I*k(t*w+x))+(6*I)*(diff(U(xi), xi))*v*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*w*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*w*(D(k))(t*w+x)*exp(I*k(t*w+x))-(24*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*w*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x))-alpha*s*c[1]*(-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^2-2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*v*c[2]+c[3]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*U(xi)^4-4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*v*c[3]+c[4]*(-(diff(U(xi), xi))*v*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*w*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(-6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))*v-2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))*v)-beta*s*(-(10*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*exp(I*k(t*w+x))+(5*I)*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))+(10*I)*(diff(diff(U(xi), xi), xi))*((D@@3)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*((D@@5)(k))(t*w+x)*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)^5*exp(I*k(t*w+x))-10*U(xi)*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-15*(diff(U(xi), xi))*((D@@2)(k))(t*w+x)^2*exp(I*k(t*w+x))+10*U(xi)*(D(k))(t*w+x)^3*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(10*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*exp(I*k(t*w+x))-5*U(xi)*((D@@4)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-20*(diff(U(xi), xi))*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*exp(I*k(t*w+x))-30*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(5*I)*(diff(U(xi), xi))*((D@@4)(k))(t*w+x)*exp(I*k(t*w+x))-(15*I)*U(xi)*((D@@2)(k))(t*w+x)^2*(D(k))(t*w+x)*exp(I*k(t*w+x))-(10*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*exp(I*k(t*w+x))-(30*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x))*((D@@2)(k))(t*w+x)+(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))*exp(I*k(t*w+x))+5*(diff(U(xi), xi))*(D(k))(t*w+x)^4*exp(I*k(t*w+x))-10*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2*exp(I*k(t*w+x)))+c[2]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^2+2*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^2*c[2]+c[3]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*U(xi)^4+4*exp(I*k(t*w+x))*(diff(U(xi), xi))*U(xi)^4*c[3]+c[4]*((diff(U(xi), xi))*exp(I*k(t*w+x))+I*U(xi)*(D(k))(t*w+x)*exp(I*k(t*w+x)))*(2*(diff(U(xi), xi))^2+2*U(xi)*(diff(diff(U(xi), xi), xi)))+c[4]*U(xi)*exp(I*k(t*w+x))*(6*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))+2*U(xi)*(diff(diff(diff(U(xi), xi), xi), xi))) = 0))

-exp(I*k(t*w+x))*(-I*U(xi)*(D(k))(t*w+x)^5*c[1]*w-(4*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*c[1]*v+(6*I)*c[1]*w*(D(k))(t*w+x)^3*(diff(diff(U(xi), xi), xi))+(4*I)*c[1]*v*((D@@3)(k))(t*w+x)*(diff(diff(U(xi), xi), xi))-(6*I)*((D@@3)(k))(t*w+x)*(diff(diff(U(xi), xi), xi))*c[1]*w-(4*I)*((D@@4)(k))(t*w+x)*(diff(U(xi), xi))*c[1]*w+(4*I)*c[1]*v*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)-I*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)*c[1]*w-I*U(xi)*((D@@5)(k))(t*w+x)*c[1]*w+(6*I)*c[1]*v*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)-(4*I)*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)*c[1]*w+(24*I)*s*beta*U(xi)*(D(k))(t*w+x)^5-(240*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*beta*s+(240*I)*s*beta*((D@@3)(k))(t*w+x)*(diff(diff(U(xi), xi), xi))+(120*I)*s*beta*((D@@4)(k))(t*w+x)*(diff(U(xi), xi))+(120*I)*s*beta*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)+(24*I)*s*beta*((D@@5)(k))(t*w+x)*U(xi)+(240*I)*s*beta*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)-(48*I)*U(xi)^2*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*c[4]-(48*I)*U(xi)*(diff(U(xi), xi))^2*(D(k))(t*w+x)*c[4]-(24*I)*c[2]*w*U(xi)^3*(D(k))(t*w+x)-(24*I)*c[3]*w*U(xi)^5*(D(k))(t*w+x)+24*c[1]*s*alpha*(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))+48*c[4]*v*(diff(diff(diff(U(xi), xi), xi), xi))*U(xi)^2+c[1]*v*(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))+24*s*beta*(diff(diff(diff(diff(diff(U(xi), xi), xi), xi), xi), xi))-48*c[4]*(diff(diff(diff(U(xi), xi), xi), xi))*U(xi)^2+I*c[1]*v*((D@@4)(k))(t*w+x)*(diff(U(xi), xi))-10*c[1]*w*((D@@2)(k))(t*w+x)*U(xi)*(D(k))(t*w+x)^3-12*c[1]*v*((D@@2)(k))(t*w+x)*(D(k))(t*w+x)*(diff(diff(U(xi), xi), xi))+18*c[1]*w*((D@@2)(k))(t*w+x)*(D(k))(t*w+x)*(diff(diff(U(xi), xi), xi))-4*c[1]*v*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*(diff(U(xi), xi))+16*c[1]*w*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*(diff(U(xi), xi))+5*c[1]*w*((D@@4)(k))(t*w+x)*U(xi)*(D(k))(t*w+x)+10*c[1]*w*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*U(xi)+120*c[1]*s*alpha*(D(k))(t*w+x)^4*(diff(U(xi), xi))-360*c[1]*s*alpha*((D@@2)(k))(t*w+x)^2*(diff(U(xi), xi))-240*c[1]*s*alpha*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2+192*c[4]*v*U(xi)*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))+240*s*beta*((D@@2)(k))(t*w+x)*U(xi)*(D(k))(t*w+x)^3-720*s*beta*((D@@2)(k))(t*w+x)*(D(k))(t*w+x)*(diff(diff(U(xi), xi), xi))-480*s*beta*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*(diff(U(xi), xi))-120*s*beta*((D@@4)(k))(t*w+x)*U(xi)*(D(k))(t*w+x)-240*s*beta*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*U(xi)-(240*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*c[1]*alpha*s-(360*I)*U(xi)*(D(k))(t*w+x)*((D@@2)(k))(t*w+x)^2*c[1]*alpha*s-(720*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*((D@@2)(k))(t*w+x)*c[1]*alpha*s+240*c[1]*s*alpha*((D@@2)(k))(t*w+x)*U(xi)*(D(k))(t*w+x)^3-720*c[1]*s*alpha*((D@@2)(k))(t*w+x)*(D(k))(t*w+x)*(diff(diff(U(xi), xi), xi))-480*c[1]*s*alpha*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)*(diff(U(xi), xi))-120*c[1]*s*alpha*((D@@4)(k))(t*w+x)*U(xi)*(D(k))(t*w+x)-240*c[1]*s*alpha*((D@@3)(k))(t*w+x)*((D@@2)(k))(t*w+x)*U(xi)+(15*I)*c[1]*w*((D@@2)(k))(t*w+x)^2*U(xi)*(D(k))(t*w+x)+(24*I)*c[1]*s*alpha*U(xi)*(D(k))(t*w+x)^5-(240*I)*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)^3*c[1]*alpha*s+(240*I)*c[1]*s*alpha*((D@@3)(k))(t*w+x)*(diff(diff(U(xi), xi), xi))+(120*I)*c[1]*s*alpha*((D@@4)(k))(t*w+x)*(diff(U(xi), xi))+(120*I)*c[1]*s*alpha*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*(D(k))(t*w+x)+(24*I)*c[1]*s*alpha*((D@@5)(k))(t*w+x)*U(xi)+(240*I)*c[1]*s*alpha*(diff(diff(diff(U(xi), xi), xi), xi))*((D@@2)(k))(t*w+x)-(48*I)*U(xi)^2*(diff(diff(U(xi), xi), xi))*(D(k))(t*w+x)*c[4]*w-48*c[4]*(diff(U(xi), xi))^3-(48*I)*U(xi)*(diff(U(xi), xi))^2*(D(k))(t*w+x)*c[4]*w-(720*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*((D@@2)(k))(t*w+x)*beta*s-(240*I)*U(xi)*((D@@3)(k))(t*w+x)*(D(k))(t*w+x)^2*beta*s-(360*I)*U(xi)*(D(k))(t*w+x)*((D@@2)(k))(t*w+x)^2*beta*s-(6*I)*(diff(U(xi), xi))*(D(k))(t*w+x)^2*((D@@2)(k))(t*w+x)*c[1]*v+(24*I)*c[1]*w*((D@@2)(k))(t*w+x)*(D(k))(t*w+x)^2*(diff(U(xi), xi))+(10*I)*c[1]*w*((D@@3)(k))(t*w+x)*U(xi)*(D(k))(t*w+x)^2-(24*I)*s^2*U(xi)*(D(k))(t*w+x)^2-(24*I)*c[2]*U(xi)^3*(D(k))(t*w+x)-(24*I)*c[3]*U(xi)^5*(D(k))(t*w+x)+(24*I)*w^2*U(xi)*(D(k))(t*w+x)^2+24*w^2*((D@@2)(k))(t*w+x)*U(xi)+c[1]*v*(D(k))(t*w+x)^4*(diff(U(xi), xi))-4*c[1]*w*(D(k))(t*w+x)^4*(diff(U(xi), xi))-3*c[1]*v*((D@@2)(k))(t*w+x)^2*(diff(U(xi), xi))+12*c[1]*w*((D@@2)(k))(t*w+x)^2*(diff(U(xi), xi))-6*c[1]*v*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2+4*c[1]*w*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2+72*c[2]*v*U(xi)^2*(diff(U(xi), xi))+120*c[3]*v*U(xi)^4*(diff(U(xi), xi))+120*s*beta*(D(k))(t*w+x)^4*(diff(U(xi), xi))-360*s*beta*((D@@2)(k))(t*w+x)^2*(diff(U(xi), xi))-240*s*beta*(diff(diff(diff(U(xi), xi), xi), xi))*(D(k))(t*w+x)^2-192*c[4]*U(xi)*(diff(U(xi), xi))*(diff(diff(U(xi), xi), xi))-48*w*(D(k))(t*w+x)*v*(diff(U(xi), xi))-(24*I)*v^2*(diff(diff(U(xi), xi), xi))+(24*I)*s^2*(diff(diff(U(xi), xi), xi))-48*s^2*(D(k))(t*w+x)*(diff(U(xi), xi))-24*s^2*((D@@2)(k))(t*w+x)*U(xi)+48*c[4]*v*(diff(U(xi), xi))^3-72*c[2]*U(xi)^2*(diff(U(xi), xi))-120*c[3]*U(xi)^4*(diff(U(xi), xi))) = 0

(10)

%/(-exp(I*k(t*w+x)))

-()/exp(k(t*w+x)*I)

(11)
 

``

Download tr.mw

I want to convert my code output into LaTeX format, but the current formatting isn't suitable for presentation. For example, when I use simplify, it sometimes introduces unnecessary fractions, making the expression look cluttered and less elegant on paper. I'm looking for a way to simplify expressions preferably by factoring terms, without introducing extra fractions, so the final LaTeX result appears clean and well-structured.

restart

with(PDEtools)

with(LinearAlgebra)

NULL

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

declare(u(x, t)); declare(U(xi)); declare(V(xi))

u(x, t)*`will now be displayed as`*u

 

U(xi)*`will now be displayed as`*U

 

V(xi)*`will now be displayed as`*V

(2)

``

RR := (2*v*(1/3)+2*alpha*beta*(1/3))*U(xi)^3+(-3*alpha*k^2*lambda-4*alpha^2*k-k^2*v+2*k^2*w-4*k*v*w)*U(xi)+(alpha*lambda+v)*(diff(diff(U(xi), xi), xi)) = 0

((2/3)*v+(2/3)*alpha*beta)*U(xi)^3+(-3*alpha*k^2*lambda-4*alpha^2*k-k^2*v+2*k^2*w-4*k*v*w)*U(xi)+(alpha*lambda+v)*(diff(diff(U(xi), xi), xi)) = 0

(3)

IM := -2*(diff(diff(U(xi), xi), xi))*v*k-2*U(xi)*k^2*w^2-2*(diff(diff(U(xi), xi), xi))*alpha^2+2*(diff(diff(U(xi), xi), xi))*v^2+(diff(diff(U(xi), xi), xi))*k*w+2*U(xi)^3*k*w-2*U(xi)^3*alpha*beta*k-3*(diff(diff(U(xi), xi), xi))*alpha*k*lambda-U(xi)*k^3*w+2*U(xi)*alpha^2*k^2+U(xi)*alpha*k^3*lambda = 0

-2*(diff(diff(U(xi), xi), xi))*v*k-2*U(xi)*k^2*w^2-2*(diff(diff(U(xi), xi), xi))*alpha^2+2*(diff(diff(U(xi), xi), xi))*v^2+(diff(diff(U(xi), xi), xi))*k*w+2*U(xi)^3*k*w-2*U(xi)^3*alpha*beta*k-3*(diff(diff(U(xi), xi), xi))*alpha*k*lambda-U(xi)*k^3*w+2*U(xi)*alpha^2*k^2+U(xi)*alpha*k^3*lambda = 0

(4)

collect(%, {U(xi), diff(diff(U(xi), xi), xi)})

(-2*alpha*beta*k+2*k*w)*U(xi)^3+(alpha*k^3*lambda+2*alpha^2*k^2-k^3*w-2*k^2*w^2)*U(xi)+(-3*alpha*k*lambda-2*alpha^2-2*k*v+k*w+2*v^2)*(diff(diff(U(xi), xi), xi)) = 0

(5)

P := %

(-2*alpha*beta*k+2*k*w)*U(xi)^3+(alpha*k^3*lambda+2*alpha^2*k^2-k^3*w-2*k^2*w^2)*U(xi)+(-3*alpha*k*lambda-2*alpha^2-2*k*v+k*w+2*v^2)*(diff(diff(U(xi), xi), xi)) = 0

(6)

NULL

NULL

C1 := v = solve(2*v*(1/3)+2*alpha*beta*(1/3) = 0, v)

v = -alpha*beta

(7)

C2 := k = solve(-3*alpha*k^2*lambda-4*alpha^2*k-k^2*v+2*k^2*w-4*k*v*w = 0, k)

k = (0, -4*(alpha^2+v*w)/(3*alpha*lambda+v-2*w))

(8)

C22 := subs(C1, C2)

k = (0, -4*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w))

(9)

C222 := k = -(4*(-alpha*beta*w+alpha^2))/(-alpha*beta+3*alpha*lambda-2*w)

k = -4*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)

(10)

ode := subs({C1, C222}, P)

(8*alpha*beta*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)-8*(-alpha*beta*w+alpha^2)*w/(-alpha*beta+3*alpha*lambda-2*w))*U(xi)^3+(-64*alpha*(-alpha*beta*w+alpha^2)^3*lambda/(-alpha*beta+3*alpha*lambda-2*w)^3+32*alpha^2*(-alpha*beta*w+alpha^2)^2/(-alpha*beta+3*alpha*lambda-2*w)^2+64*(-alpha*beta*w+alpha^2)^3*w/(-alpha*beta+3*alpha*lambda-2*w)^3-32*(-alpha*beta*w+alpha^2)^2*w^2/(-alpha*beta+3*alpha*lambda-2*w)^2)*U(xi)+(12*alpha*(-alpha*beta*w+alpha^2)*lambda/(-alpha*beta+3*alpha*lambda-2*w)-2*alpha^2-8*alpha*beta*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)-4*(-alpha*beta*w+alpha^2)*w/(-alpha*beta+3*alpha*lambda-2*w)+2*alpha^2*beta^2)*(diff(diff(U(xi), xi), xi)) = 0

(11)

numer(lhs((8*alpha*beta*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)-8*(-alpha*beta*w+alpha^2)*w/(-alpha*beta+3*alpha*lambda-2*w))*U(xi)^3+(-64*alpha*(-alpha*beta*w+alpha^2)^3*lambda/(-alpha*beta+3*alpha*lambda-2*w)^3+32*alpha^2*(-alpha*beta*w+alpha^2)^2/(-alpha*beta+3*alpha*lambda-2*w)^2+64*(-alpha*beta*w+alpha^2)^3*w/(-alpha*beta+3*alpha*lambda-2*w)^3-32*(-alpha*beta*w+alpha^2)^2*w^2/(-alpha*beta+3*alpha*lambda-2*w)^2)*U(xi)+(12*alpha*(-alpha*beta*w+alpha^2)*lambda/(-alpha*beta+3*alpha*lambda-2*w)-2*alpha^2-8*alpha*beta*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)-4*(-alpha*beta*w+alpha^2)*w/(-alpha*beta+3*alpha*lambda-2*w)+2*alpha^2*beta^2)*(diff(diff(U(xi), xi), xi)) = 0))*denom(rhs((8*alpha*beta*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)-8*(-alpha*beta*w+alpha^2)*w/(-alpha*beta+3*alpha*lambda-2*w))*U(xi)^3+(-64*alpha*(-alpha*beta*w+alpha^2)^3*lambda/(-alpha*beta+3*alpha*lambda-2*w)^3+32*alpha^2*(-alpha*beta*w+alpha^2)^2/(-alpha*beta+3*alpha*lambda-2*w)^2+64*(-alpha*beta*w+alpha^2)^3*w/(-alpha*beta+3*alpha*lambda-2*w)^3-32*(-alpha*beta*w+alpha^2)^2*w^2/(-alpha*beta+3*alpha*lambda-2*w)^2)*U(xi)+(12*alpha*(-alpha*beta*w+alpha^2)*lambda/(-alpha*beta+3*alpha*lambda-2*w)-2*alpha^2-8*alpha*beta*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)-4*(-alpha*beta*w+alpha^2)*w/(-alpha*beta+3*alpha*lambda-2*w)+2*alpha^2*beta^2)*(diff(diff(U(xi), xi), xi)) = 0)) = numer(rhs((8*alpha*beta*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)-8*(-alpha*beta*w+alpha^2)*w/(-alpha*beta+3*alpha*lambda-2*w))*U(xi)^3+(-64*alpha*(-alpha*beta*w+alpha^2)^3*lambda/(-alpha*beta+3*alpha*lambda-2*w)^3+32*alpha^2*(-alpha*beta*w+alpha^2)^2/(-alpha*beta+3*alpha*lambda-2*w)^2+64*(-alpha*beta*w+alpha^2)^3*w/(-alpha*beta+3*alpha*lambda-2*w)^3-32*(-alpha*beta*w+alpha^2)^2*w^2/(-alpha*beta+3*alpha*lambda-2*w)^2)*U(xi)+(12*alpha*(-alpha*beta*w+alpha^2)*lambda/(-alpha*beta+3*alpha*lambda-2*w)-2*alpha^2-8*alpha*beta*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)-4*(-alpha*beta*w+alpha^2)*w/(-alpha*beta+3*alpha*lambda-2*w)+2*alpha^2*beta^2)*(diff(diff(U(xi), xi), xi)) = 0))*denom(lhs((8*alpha*beta*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)-8*(-alpha*beta*w+alpha^2)*w/(-alpha*beta+3*alpha*lambda-2*w))*U(xi)^3+(-64*alpha*(-alpha*beta*w+alpha^2)^3*lambda/(-alpha*beta+3*alpha*lambda-2*w)^3+32*alpha^2*(-alpha*beta*w+alpha^2)^2/(-alpha*beta+3*alpha*lambda-2*w)^2+64*(-alpha*beta*w+alpha^2)^3*w/(-alpha*beta+3*alpha*lambda-2*w)^3-32*(-alpha*beta*w+alpha^2)^2*w^2/(-alpha*beta+3*alpha*lambda-2*w)^2)*U(xi)+(12*alpha*(-alpha*beta*w+alpha^2)*lambda/(-alpha*beta+3*alpha*lambda-2*w)-2*alpha^2-8*alpha*beta*(-alpha*beta*w+alpha^2)/(-alpha*beta+3*alpha*lambda-2*w)-4*(-alpha*beta*w+alpha^2)*w/(-alpha*beta+3*alpha*lambda-2*w)+2*alpha^2*beta^2)*(diff(diff(U(xi), xi), xi)) = 0))

-2*alpha*(-24*U(xi)^3*alpha^4*beta^2*lambda+36*U(xi)^3*alpha^4*beta*lambda^2-12*U(xi)^3*alpha^2*beta^3*w^2-16*U(xi)*alpha^4*beta^3*w^2-16*U(xi)*alpha^2*beta^3*w^4+9*(diff(diff(U(xi), xi), xi))*alpha^4*beta^4*lambda-27*(diff(diff(U(xi), xi), xi))*alpha^4*beta^3*lambda^2+27*(diff(diff(U(xi), xi), xi))*alpha^4*beta^2*lambda^3+12*U(xi)^3*alpha^3*beta^2*w-36*U(xi)^3*alpha^3*lambda^2*w+32*U(xi)*alpha^5*beta^2*w+32*U(xi)*alpha^3*beta^2*w^3+32*U(xi)*alpha*beta^2*w^5-2*(diff(diff(U(xi), xi), xi))*alpha^3*beta^4*w+48*U(xi)^3*alpha^2*lambda*w^2-16*U(xi)*alpha^4*beta*w^2-48*U(xi)*alpha^4*lambda*w^2-64*U(xi)*alpha^2*beta*w^4+21*(diff(diff(U(xi), xi), xi))*alpha^4*beta^2*lambda-45*(diff(diff(U(xi), xi), xi))*alpha^4*beta*lambda^2+6*(diff(diff(U(xi), xi), xi))*alpha^2*beta^3*w^2-12*(diff(diff(U(xi), xi), xi))*alpha^3*beta^2*w-36*(diff(diff(U(xi), xi), xi))*alpha^3*lambda^2*w+16*(diff(diff(U(xi), xi), xi))*alpha*beta^2*w^3-12*(diff(diff(U(xi), xi), xi))*alpha^2*beta*w^2+12*(diff(diff(U(xi), xi), xi))*alpha^2*lambda*w^2-4*U(xi)^3*alpha^3*beta^4*w+24*U(xi)^3*alpha^3*beta^3*lambda*w-36*U(xi)^3*alpha^3*beta^2*lambda^2*w+32*U(xi)*alpha^3*beta^3*lambda*w^3+24*U(xi)^3*alpha^2*beta^2*lambda*w^2+36*U(xi)^3*alpha^2*beta*lambda^2*w^2-48*U(xi)*alpha^4*beta^2*lambda*w^2-48*U(xi)*alpha^2*beta^2*lambda*w^4-24*U(xi)^3*alpha^3*beta*lambda*w-48*U(xi)^3*alpha*beta*lambda*w^3+96*U(xi)*alpha^3*beta*lambda*w^3+6*(diff(diff(U(xi), xi), xi))*alpha^3*beta^3*lambda*w+18*(diff(diff(U(xi), xi), xi))*alpha^3*beta^2*lambda^2*w-54*(diff(diff(U(xi), xi), xi))*alpha^3*beta*lambda^3*w-48*(diff(diff(U(xi), xi), xi))*alpha^2*beta^2*lambda*w^2+90*(diff(diff(U(xi), xi), xi))*alpha^2*beta*lambda^2*w^2+48*(diff(diff(U(xi), xi), xi))*alpha^3*beta*lambda*w-48*(diff(diff(U(xi), xi), xi))*alpha*beta*lambda*w^3+4*U(xi)^3*alpha^4*beta^3-(diff(diff(U(xi), xi), xi))*alpha^4*beta^5+16*U(xi)^3*beta*w^4-16*U(xi)*alpha^6*beta+16*U(xi)*alpha^6*lambda-3*(diff(diff(U(xi), xi), xi))*alpha^4*beta^3+27*(diff(diff(U(xi), xi), xi))*alpha^4*lambda^3-16*U(xi)^3*alpha*w^3+32*U(xi)*alpha^3*w^3+8*(diff(diff(U(xi), xi), xi))*beta*w^4) = 0

(12)

simplify(-2*alpha*(24*U(xi)^3*alpha^3*beta^3*lambda*w-36*U(xi)^3*alpha^3*beta^2*lambda^2*w+32*U(xi)*alpha^3*beta^3*lambda*w^3+24*U(xi)^3*alpha^2*beta^2*lambda*w^2+36*U(xi)^3*alpha^2*beta*lambda^2*w^2-48*U(xi)*alpha^4*beta^2*lambda*w^2-48*U(xi)*alpha^2*beta^2*lambda*w^4-24*U(xi)^3*alpha^3*beta*lambda*w-48*U(xi)^3*alpha*beta*lambda*w^3+96*U(xi)*alpha^3*beta*lambda*w^3+6*(diff(diff(U(xi), xi), xi))*alpha^3*beta^3*lambda*w+18*(diff(diff(U(xi), xi), xi))*alpha^3*beta^2*lambda^2*w-54*(diff(diff(U(xi), xi), xi))*alpha^3*beta*lambda^3*w-48*(diff(diff(U(xi), xi), xi))*alpha^2*beta^2*lambda*w^2+90*(diff(diff(U(xi), xi), xi))*alpha^2*beta*lambda^2*w^2+48*(diff(diff(U(xi), xi), xi))*alpha^3*beta*lambda*w-48*(diff(diff(U(xi), xi), xi))*alpha*beta*lambda*w^3-24*U(xi)^3*alpha^4*beta^2*lambda+36*U(xi)^3*alpha^4*beta*lambda^2-12*U(xi)^3*alpha^2*beta^3*w^2-16*U(xi)*alpha^4*beta^3*w^2-16*U(xi)*alpha^2*beta^3*w^4+9*(diff(diff(U(xi), xi), xi))*alpha^4*beta^4*lambda-27*(diff(diff(U(xi), xi), xi))*alpha^4*beta^3*lambda^2+27*(diff(diff(U(xi), xi), xi))*alpha^4*beta^2*lambda^3+12*U(xi)^3*alpha^3*beta^2*w-36*U(xi)^3*alpha^3*lambda^2*w+32*U(xi)*alpha^5*beta^2*w+32*U(xi)*alpha^3*beta^2*w^3+32*U(xi)*alpha*beta^2*w^5-2*(diff(diff(U(xi), xi), xi))*alpha^3*beta^4*w+48*U(xi)^3*alpha^2*lambda*w^2-16*U(xi)*alpha^4*beta*w^2-48*U(xi)*alpha^4*lambda*w^2-64*U(xi)*alpha^2*beta*w^4+21*(diff(diff(U(xi), xi), xi))*alpha^4*beta^2*lambda-45*(diff(diff(U(xi), xi), xi))*alpha^4*beta*lambda^2+6*(diff(diff(U(xi), xi), xi))*alpha^2*beta^3*w^2-12*(diff(diff(U(xi), xi), xi))*alpha^3*beta^2*w-36*(diff(diff(U(xi), xi), xi))*alpha^3*lambda^2*w+16*(diff(diff(U(xi), xi), xi))*alpha*beta^2*w^3-12*(diff(diff(U(xi), xi), xi))*alpha^2*beta*w^2+12*(diff(diff(U(xi), xi), xi))*alpha^2*lambda*w^2-4*U(xi)^3*alpha^3*beta^4*w-(diff(diff(U(xi), xi), xi))*alpha^4*beta^5+4*U(xi)^3*alpha^4*beta^3+16*U(xi)^3*beta*w^4-16*U(xi)*alpha^6*beta+16*U(xi)*alpha^6*lambda-3*(diff(diff(U(xi), xi), xi))*alpha^4*beta^3+27*(diff(diff(U(xi), xi), xi))*alpha^4*lambda^3-16*U(xi)^3*alpha*w^3+32*U(xi)*alpha^3*w^3+8*(diff(diff(U(xi), xi), xi))*beta*w^4) = 0)

-64*alpha*((1/4)*((1/2)*(-beta^3+3*beta^2*lambda-3*beta+3*lambda)*alpha^2+beta*w*(beta-3*lambda)*alpha+beta*w^2)*((1/2)*(beta-3*lambda)*alpha+w)^2*(diff(diff(U(xi), xi), xi))+((1/2)*(-alpha*beta+w)*((1/2)*(beta-3*lambda)*alpha+w)^2*U(xi)^2+alpha*((1/2)*(-beta+lambda)*alpha^3+w*beta*alpha^2*lambda-(1/2)*w^2*(beta+3*lambda)*alpha+w^3)*(beta*w-alpha))*(beta*w-alpha)*U(xi)) = 0

(13)

%/(-64*alpha)

(1/4)*((1/2)*(-beta^3+3*beta^2*lambda-3*beta+3*lambda)*alpha^2+beta*w*(beta-3*lambda)*alpha+beta*w^2)*((1/2)*(beta-3*lambda)*alpha+w)^2*(diff(diff(U(xi), xi), xi))+((1/2)*(-alpha*beta+w)*((1/2)*(beta-3*lambda)*alpha+w)^2*U(xi)^2+alpha*((1/2)*(-beta+lambda)*alpha^3+w*beta*alpha^2*lambda-(1/2)*w^2*(beta+3*lambda)*alpha+w^3)*(beta*w-alpha))*(beta*w-alpha)*U(xi) = 0

(14)

PP := numer(lhs((1/4)*((1/2)*(-beta^3+3*beta^2*lambda-3*beta+3*lambda)*alpha^2+beta*w*(beta-3*lambda)*alpha+beta*w^2)*((1/2)*(beta-3*lambda)*alpha+w)^2*(diff(diff(U(xi), xi), xi))+((1/2)*(-alpha*beta+w)*((1/2)*(beta-3*lambda)*alpha+w)^2*U(xi)^2+alpha*((1/2)*(-beta+lambda)*alpha^3+w*beta*alpha^2*lambda-(1/2)*w^2*(beta+3*lambda)*alpha+w^3)*(beta*w-alpha))*(beta*w-alpha)*U(xi) = 0))*denom(rhs((1/4)*((1/2)*(-beta^3+3*beta^2*lambda-3*beta+3*lambda)*alpha^2+beta*w*(beta-3*lambda)*alpha+beta*w^2)*((1/2)*(beta-3*lambda)*alpha+w)^2*(diff(diff(U(xi), xi), xi))+((1/2)*(-alpha*beta+w)*((1/2)*(beta-3*lambda)*alpha+w)^2*U(xi)^2+alpha*((1/2)*(-beta+lambda)*alpha^3+w*beta*alpha^2*lambda-(1/2)*w^2*(beta+3*lambda)*alpha+w^3)*(beta*w-alpha))*(beta*w-alpha)*U(xi) = 0)) = numer(rhs((1/4)*((1/2)*(-beta^3+3*beta^2*lambda-3*beta+3*lambda)*alpha^2+beta*w*(beta-3*lambda)*alpha+beta*w^2)*((1/2)*(beta-3*lambda)*alpha+w)^2*(diff(diff(U(xi), xi), xi))+((1/2)*(-alpha*beta+w)*((1/2)*(beta-3*lambda)*alpha+w)^2*U(xi)^2+alpha*((1/2)*(-beta+lambda)*alpha^3+w*beta*alpha^2*lambda-(1/2)*w^2*(beta+3*lambda)*alpha+w^3)*(beta*w-alpha))*(beta*w-alpha)*U(xi) = 0))*denom(lhs((1/4)*((1/2)*(-beta^3+3*beta^2*lambda-3*beta+3*lambda)*alpha^2+beta*w*(beta-3*lambda)*alpha+beta*w^2)*((1/2)*(beta-3*lambda)*alpha+w)^2*(diff(diff(U(xi), xi), xi))+((1/2)*(-alpha*beta+w)*((1/2)*(beta-3*lambda)*alpha+w)^2*U(xi)^2+alpha*((1/2)*(-beta+lambda)*alpha^3+w*beta*alpha^2*lambda-(1/2)*w^2*(beta+3*lambda)*alpha+w^3)*(beta*w-alpha))*(beta*w-alpha)*U(xi) = 0))

-24*U(xi)^3*alpha^4*beta^2*lambda+36*U(xi)^3*alpha^4*beta*lambda^2-12*U(xi)^3*alpha^2*beta^3*w^2-16*U(xi)*alpha^4*beta^3*w^2-16*U(xi)*alpha^2*beta^3*w^4+9*(diff(diff(U(xi), xi), xi))*alpha^4*beta^4*lambda-27*(diff(diff(U(xi), xi), xi))*alpha^4*beta^3*lambda^2+27*(diff(diff(U(xi), xi), xi))*alpha^4*beta^2*lambda^3+12*U(xi)^3*alpha^3*beta^2*w-36*U(xi)^3*alpha^3*lambda^2*w+32*U(xi)*alpha^5*beta^2*w+32*U(xi)*alpha^3*beta^2*w^3+32*U(xi)*alpha*beta^2*w^5-2*(diff(diff(U(xi), xi), xi))*alpha^3*beta^4*w+48*U(xi)^3*alpha^2*lambda*w^2-16*U(xi)*alpha^4*beta*w^2-48*U(xi)*alpha^4*lambda*w^2-64*U(xi)*alpha^2*beta*w^4+21*(diff(diff(U(xi), xi), xi))*alpha^4*beta^2*lambda-45*(diff(diff(U(xi), xi), xi))*alpha^4*beta*lambda^2+6*(diff(diff(U(xi), xi), xi))*alpha^2*beta^3*w^2-12*(diff(diff(U(xi), xi), xi))*alpha^3*beta^2*w-36*(diff(diff(U(xi), xi), xi))*alpha^3*lambda^2*w+16*(diff(diff(U(xi), xi), xi))*alpha*beta^2*w^3-12*(diff(diff(U(xi), xi), xi))*alpha^2*beta*w^2+12*(diff(diff(U(xi), xi), xi))*alpha^2*lambda*w^2-4*U(xi)^3*alpha^3*beta^4*w+24*U(xi)^3*alpha^3*beta^3*lambda*w-36*U(xi)^3*alpha^3*beta^2*lambda^2*w+32*U(xi)*alpha^3*beta^3*lambda*w^3+24*U(xi)^3*alpha^2*beta^2*lambda*w^2+36*U(xi)^3*alpha^2*beta*lambda^2*w^2-48*U(xi)*alpha^4*beta^2*lambda*w^2-48*U(xi)*alpha^2*beta^2*lambda*w^4-24*U(xi)^3*alpha^3*beta*lambda*w-48*U(xi)^3*alpha*beta*lambda*w^3+96*U(xi)*alpha^3*beta*lambda*w^3+6*(diff(diff(U(xi), xi), xi))*alpha^3*beta^3*lambda*w+18*(diff(diff(U(xi), xi), xi))*alpha^3*beta^2*lambda^2*w-54*(diff(diff(U(xi), xi), xi))*alpha^3*beta*lambda^3*w-48*(diff(diff(U(xi), xi), xi))*alpha^2*beta^2*lambda*w^2+90*(diff(diff(U(xi), xi), xi))*alpha^2*beta*lambda^2*w^2+48*(diff(diff(U(xi), xi), xi))*alpha^3*beta*lambda*w-48*(diff(diff(U(xi), xi), xi))*alpha*beta*lambda*w^3+4*U(xi)^3*alpha^4*beta^3-(diff(diff(U(xi), xi), xi))*alpha^4*beta^5+16*U(xi)^3*beta*w^4-16*U(xi)*alpha^6*beta+16*U(xi)*alpha^6*lambda-3*(diff(diff(U(xi), xi), xi))*alpha^4*beta^3+27*(diff(diff(U(xi), xi), xi))*alpha^4*lambda^3-16*U(xi)^3*alpha*w^3+32*U(xi)*alpha^3*w^3+8*(diff(diff(U(xi), xi), xi))*beta*w^4 = 0

(15)

NULL

collect(PP, {U(xi), diff(U(xi), xi), diff(diff(U(xi), xi), xi)})

(-4*alpha^3*beta^4*w+24*alpha^3*beta^3*lambda*w-36*alpha^3*beta^2*lambda^2*w+4*alpha^4*beta^3-24*alpha^4*beta^2*lambda+36*alpha^4*beta*lambda^2-12*alpha^2*beta^3*w^2+24*alpha^2*beta^2*lambda*w^2+36*alpha^2*beta*lambda^2*w^2+12*alpha^3*beta^2*w-24*alpha^3*beta*lambda*w-36*alpha^3*lambda^2*w-48*alpha*beta*lambda*w^3+48*alpha^2*lambda*w^2+16*beta*w^4-16*alpha*w^3)*U(xi)^3+(32*alpha^3*beta^3*lambda*w^3-16*alpha^4*beta^3*w^2-48*alpha^4*beta^2*lambda*w^2-16*alpha^2*beta^3*w^4-48*alpha^2*beta^2*lambda*w^4+32*alpha^5*beta^2*w+32*alpha^3*beta^2*w^3+96*alpha^3*beta*lambda*w^3+32*alpha*beta^2*w^5-16*alpha^6*beta+16*alpha^6*lambda-16*alpha^4*beta*w^2-48*alpha^4*lambda*w^2-64*alpha^2*beta*w^4+32*alpha^3*w^3)*U(xi)+(-alpha^4*beta^5+9*alpha^4*beta^4*lambda-27*alpha^4*beta^3*lambda^2+27*alpha^4*beta^2*lambda^3-2*alpha^3*beta^4*w+6*alpha^3*beta^3*lambda*w+18*alpha^3*beta^2*lambda^2*w-54*alpha^3*beta*lambda^3*w-3*alpha^4*beta^3+21*alpha^4*beta^2*lambda-45*alpha^4*beta*lambda^2+27*alpha^4*lambda^3+6*alpha^2*beta^3*w^2-48*alpha^2*beta^2*lambda*w^2+90*alpha^2*beta*lambda^2*w^2-12*alpha^3*beta^2*w+48*alpha^3*beta*lambda*w-36*alpha^3*lambda^2*w+16*alpha*beta^2*w^3-48*alpha*beta*lambda*w^3-12*alpha^2*beta*w^2+12*alpha^2*lambda*w^2+8*beta*w^4)*(diff(diff(U(xi), xi), xi)) = 0

(16)
 

NULL

Download B-R.mw

How can we eliminate nonlinear terms involving two functions in a differential equation?

restart

with(PDEtools)

with(LinearAlgebra)

NULL

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

declare(G(xi)); declare(Q(x, t)); declare(Q1(x, t))

G(xi)*`will now be displayed as`*G

 

Q(x, t)*`will now be displayed as`*Q

 

Q1(x, t)*`will now be displayed as`*Q1

(2)

NULL

q := (sqrt(P)+Q(x, t))*exp(I*gamma*P*t); B := (sqrt(P)+Q(x, t))*exp(-I*gamma*P*t); B1 := sqrt(P)+Q(x, t); P+sqrt(P)*(Q1(x, t)+Q(x, t))

GeF := I*(diff(q, t))+alpha[1]*(diff(q, x, x))+alpha[2]*(P+sqrt(P)*(Q1(x, t)+Q(x, t)))*q+alpha[3]*(P+sqrt(P)*(Q1(x, t)+Q(x, t)))^2*q+alpha[4]*(P+sqrt(P)*(Q1(x, t)+Q(x, t)))^3*q+alpha[5]*(diff(P+sqrt(P)*(Q1(x, t)+Q(x, t)), x, x))*q = 0

K := simplify(GeF*exp(-I*gamma*P*t))

(Q(x, t)^4*alpha[4]+3*Q(x, t)^3*Q1(x, t)*alpha[4]+3*(Q1(x, t)^2*alpha[4]+alpha[3])*Q(x, t)^2+Q1(x, t)*(Q1(x, t)^2*alpha[4]+4*alpha[3])*Q(x, t)+Q1(x, t)^2*alpha[3]+alpha[2]-gamma)*P^(3/2)+(6*Q(x, t)^2*alpha[4]+9*Q(x, t)*Q1(x, t)*alpha[4]+3*Q1(x, t)^2*alpha[4]+alpha[3])*P^(5/2)+(P^(1/2)*Q(x, t)*alpha[5]+P*alpha[5]+alpha[1])*(diff(diff(Q(x, t), x), x))+alpha[5]*(P^(1/2)*Q(x, t)+P)*(diff(diff(Q1(x, t), x), x))+P^(7/2)*alpha[4]+I*(diff(Q(x, t), t))+alpha[2]*(Q1(x, t)*Q(x, t)+Q(x, t)^2)*P^(1/2)+(4*P^2*alpha[4]+P*alpha[3])*Q(x, t)^3+(9*P^2*alpha[4]+2*P*alpha[3])*Q1(x, t)*Q(x, t)^2+((6*P^2*alpha[4]+P*alpha[3])*Q1(x, t)^2+4*alpha[4]*P^3+3*alpha[3]*P^2+(2*alpha[2]-gamma)*P)*Q(x, t)+P^2*Q1(x, t)^3*alpha[4]+(3*P^3*alpha[4]+2*P^2*alpha[3]+P*alpha[2])*Q1(x, t) = 0

(3)

remove(has, K, {Q(x, t)^2, Q(x, t)^3, Q(x, t)^4, Q1(x, t)^2, Q1(x, t)^3, Q1(x, t)^4})

() = 0

(4)

NULL

NULL

AA := (alpha[2]-gamma)*P^(3/2)+alpha[3]*P^(5/2)+(P*alpha[5]+alpha[1])*(diff(Q(x, t), x, x))+alpha[5]*P*(diff(Q1(x, t), x, x))+P^(7/2)*alpha[4]+I*(diff(Q(x, t), t))+(4*alpha[4]*P^3+3*alpha[3]*P^2+(2*alpha[2]-gamma)*P)*Q(x, t)+(3*P^3*alpha[4]+2*P^2*alpha[3]+P*alpha[2])*Q1(x, t) = 0

(alpha[2]-gamma)*P^(3/2)+P^(5/2)*alpha[3]+(P*alpha[5]+alpha[1])*(diff(diff(Q(x, t), x), x))+P*(diff(diff(Q1(x, t), x), x))*alpha[5]+P^(7/2)*alpha[4]+I*(diff(Q(x, t), t))+(4*alpha[4]*P^3+3*alpha[3]*P^2+(2*alpha[2]-gamma)*P)*Q(x, t)+(3*P^3*alpha[4]+2*P^2*alpha[3]+P*alpha[2])*Q1(x, t) = 0

(5)
 

test := (alpha[2]-gamma)*P^(3/2)+alpha[3]*P^(5/2)+(P*alpha[5]+alpha[1])*0+P*alpha[5]*0+P^(7/2)*alpha[4]+I*0 = 0

(alpha[2]-gamma)*P^(3/2)+P^(5/2)*alpha[3]+P^(7/2)*alpha[4] = 0

(6)
 

NULL

Download remove.mw

I’ve successfully plotted this exercise, but I’m exploring ways to improve the visualization. Specifically, I’d like to know if it's possible to combine a 2D plot as a smaller inset within a 3D plot, positioned in a corner of the 3D graph.

Additionally, I’ve noticed that although the same data is plotted in different examples, the design and style of the plots can vary significantly, some look much more polished or professional. Are there recommended techniques, functions, or toolboxes in Maple that can help improve the visual design or aesthetic of the plots

Download plot.mw

in here How we can seperate the coefficent of conjugate this conjugate sign how remove from my equation ?

restart

with(PDEtools)

with(LinearAlgebra)

NULL

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

declare(u(x, t)); declare(U(xi)); declare(V(xi)); declare(P(x, t)); declare(q(x, t))

u(x, t)*`will now be displayed as`*u

 

U(xi)*`will now be displayed as`*U

 

V(xi)*`will now be displayed as`*V

 

P(x, t)*`will now be displayed as`*P

 

q(x, t)*`will now be displayed as`*q

(2)

pde := I*(diff(u(x, t), t))+diff(u(x, t), `$`(x, 2))+abs(u(x, t))^2*u(x, t) = 0

I*(diff(u(x, t), t))+diff(diff(u(x, t), x), x)+abs(u(x, t))^2*u(x, t) = 0

(3)

S := u(x, t) = (sqrt(a)+P(x, t))*exp(I*a*t)

u(x, t) = (a^(1/2)+P(x, t))*exp(I*a*t)

(4)

S1 := conjugate(u(x, t)) = (sqrt(a)+conjugate(P(x, t)))*exp(-I*a*t)

conjugate(u(x, t)) = (a^(1/2)+conjugate(P(x, t)))*exp(-I*a*t)

(5)

Q := abs(u(x, t))^2 = u(x, t)*conjugate(u(x, t))

abs(u(x, t))^2 = u(x, t)*conjugate(u(x, t))

(6)

F1 := expand(simplify(subs({S, S1}, rhs(Q))))

a+a^(1/2)*P(x, t)+a^(1/2)*conjugate(P(x, t))+abs(P(x, t))^2

(7)

F2 := abs(u(x, t))^2 = remove(has, F1, abs(P(x, t))^2)

abs(u(x, t))^2 = a+a^(1/2)*P(x, t)+a^(1/2)*conjugate(P(x, t))

(8)

FF := collect(F2, sqrt(a))

abs(u(x, t))^2 = a+(P(x, t)+conjugate(P(x, t)))*a^(1/2)

(9)

F3 := abs(u(x, t))^2*u(x, t) = (a+(P(x, t)+conjugate(P(x, t)))*sqrt(a))*rhs(S)

abs(u(x, t))^2*u(x, t) = (a+(P(x, t)+conjugate(P(x, t)))*a^(1/2))*(a^(1/2)+P(x, t))*exp(I*a*t)

(10)

F4 := remove(has, F3, P(x, t)*conjugate(P(x, t)))

abs(u(x, t))^2*u(x, t) = (a+(P(x, t)+conjugate(P(x, t)))*a^(1/2))*(a^(1/2)+P(x, t))*exp(I*a*t)

(11)

expand(%)

abs(u(x, t))^2*u(x, t) = exp(I*a*t)*a^(3/2)+2*exp(I*a*t)*a*P(x, t)+exp(I*a*t)*a^(1/2)*P(x, t)^2+exp(I*a*t)*a*conjugate(P(x, t))+exp(I*a*t)*a^(1/2)*conjugate(P(x, t))*P(x, t)

(12)

pde_linear, pde_nonlinear := selectremove(proc (term) options operator, arrow; not has((eval(term, P(x, t) = T*P(x, t)))/T, T) end proc, expand(%))

() = (), abs(u(x, t))^2*u(x, t) = exp(I*a*t)*a^(3/2)+2*exp(I*a*t)*a*P(x, t)+exp(I*a*t)*a^(1/2)*P(x, t)^2+exp(I*a*t)*a*conjugate(P(x, t))+exp(I*a*t)*a^(1/2)*conjugate(P(x, t))*P(x, t)

(13)

F6 := abs(u(x, t))^2*u(x, t) = exp(I*a*t)*a^(3/2)+2*exp(I*a*t)*a*P(x, t)+exp(I*a*t)*a*conjugate(P(x, t))

abs(u(x, t))^2*u(x, t) = exp(a*t*I)*a^(3/2)+2*exp(a*t*I)*a*P(x, t)+exp(a*t*I)*a*conjugate(P(x, t))

(14)

subs({F6, S}, pde)

I*(diff((a^(1/2)+P(x, t))*exp(a*t*I), t))+diff(diff((a^(1/2)+P(x, t))*exp(a*t*I), x), x)+exp(a*t*I)*a^(3/2)+2*exp(a*t*I)*a*P(x, t)+exp(a*t*I)*a*conjugate(P(x, t)) = 0

(15)

eval(%)

I*((diff(P(x, t), t))*exp(a*t*I)+I*(a^(1/2)+P(x, t))*a*exp(a*t*I))+(diff(diff(P(x, t), x), x))*exp(a*t*I)+exp(a*t*I)*a^(3/2)+2*exp(a*t*I)*a*P(x, t)+exp(a*t*I)*a*conjugate(P(x, t)) = 0

(16)

expand(%)

I*(diff(P(x, t), t))*exp(a*t*I)+exp(a*t*I)*a*P(x, t)+(diff(diff(P(x, t), x), x))*exp(a*t*I)+exp(a*t*I)*a*conjugate(P(x, t)) = 0

(17)

expand(%/exp(I*a*t))

I*(diff(P(x, t), t))+a*P(x, t)+diff(diff(P(x, t), x), x)+a*conjugate(P(x, t)) = 0

(18)

PP := collect(%, a)

(P(x, t)+conjugate(P(x, t)))*a+I*(diff(P(x, t), t))+diff(diff(P(x, t), x), x) = 0

(19)

U1 := P(x, t) = r[1]*exp(I*(l*x-m*t))+r[2]*exp(-I*(l*x-m*t))

P(x, t) = r[1]*exp(I*(l*x-m*t))+r[2]*exp(-I*(l*x-m*t))

(20)

eval(subs(U1, PP))

(r[1]*exp(I*(l*x-m*t))+r[2]*exp(-I*(l*x-m*t))+conjugate(r[1]*exp(I*(l*x-m*t))+r[2]*exp(-I*(l*x-m*t))))*a+I*(-I*r[1]*m*exp(I*(l*x-m*t))+I*r[2]*m*exp(-I*(l*x-m*t)))-r[1]*l^2*exp(I*(l*x-m*t))-r[2]*l^2*exp(-I*(l*x-m*t)) = 0

(21)

simplify((r[1]*exp(I*(l*x-m*t))+r[2]*exp(-I*(l*x-m*t))+conjugate(r[1]*exp(I*(l*x-m*t))+r[2]*exp(-I*(l*x-m*t))))*a+I*(-I*r[1]*m*exp(I*(l*x-m*t))+I*r[2]*m*exp(-I*(l*x-m*t)))-r[1]*l^2*exp(I*(l*x-m*t))-r[2]*l^2*exp(-I*(l*x-m*t)) = 0)

conjugate(r[1]*exp(I*(l*x-m*t))+r[2]*exp(-I*(l*x-m*t)))*a+r[2]*(-l^2+a-m)*exp(-I*(l*x-m*t))+r[1]*exp(I*(l*x-m*t))*(-l^2+a+m) = 0

(22)

J := eval(%)

conjugate(r[1]*exp(I*(l*x-m*t))+r[2]*exp(-I*(l*x-m*t)))*a+r[2]*(-l^2+a-m)*exp(-I*(l*x-m*t))+r[1]*exp(I*(l*x-m*t))*(-l^2+a+m) = 0

(23)

expand(%)

a*conjugate(r[1])*exp(I*conjugate(m)*conjugate(t))/exp(I*conjugate(l)*conjugate(x))+a*conjugate(r[2])*exp(I*conjugate(l)*conjugate(x))/exp(I*conjugate(m)*conjugate(t))-r[2]*exp(I*m*t)*l^2/exp(I*l*x)+r[2]*exp(I*m*t)*a/exp(I*l*x)-r[2]*exp(I*m*t)*m/exp(I*l*x)-r[1]*exp(I*l*x)*l^2/exp(I*m*t)+r[1]*exp(I*l*x)*a/exp(I*m*t)+r[1]*exp(I*l*x)*m/exp(I*m*t) = 0

(24)

indets(%)

{a, l, m, t, x, r[1], r[2], exp(I*l*x), exp(I*m*t), exp(I*conjugate(l)*conjugate(x)), exp(I*conjugate(m)*conjugate(t)), conjugate(l), conjugate(m), conjugate(t), conjugate(x), conjugate(r[1]), conjugate(r[2])}

(25)

subs({exp(-I*(l*x-m*t)) = Y, exp(I*(l*x-m*t)) = X}, J)

conjugate(X*r[1]+Y*r[2])*a+r[2]*(-l^2+a-m)*Y+r[1]*X*(-l^2+a+m) = 0

(26)

collect(%, {X, Y})

conjugate(X*r[1]+Y*r[2])*a+r[2]*(-l^2+a-m)*Y+r[1]*X*(-l^2+a+m) = 0

(27)

Download conjugate.mw

I’m currently trying to collect terms in an expression f^G(xi), but the result is not behaving as expected. I attempted two different coding approaches, but both resulted in errors. This particular case of collecting terms seems to be different from what I’ve encountered before, and I’m unsure how to resolve it.

Could you please advise on how to properly collect the terms in this situation and avoid the errors? Any insight into why this case behaves differently would also be appreciated.

Thank you for your help.

collect_term.mw

I am currently working with an ordinary differential equation (ODE) that I find difficult to express and solve accurately. In this ODE, the symbol f represents an exponential function rather than a typical variable, which adds to the confusion. Although I have followed the format used in related research papers, the results I obtain are not satisfactory.

Since this type of ODE is new and somewhat unfamiliar to me, I would greatly appreciate your guidance in:

  1. Properly formulating the ODE.

  2. Understanding the role of f in the context of exponential functions.

  3. Finding the correct and complete solutions.

  4. Learning how to clearly present each solution step by step.

Thank you in advance for your support.

AA.mw

Manually factoring each equation in this system one by one is time-consuming and inefficient. Is there a way to automate the factoring of expressions into two multiplicative terms—some of which may be single-term factors—using code?

restart

with(PDEtools)

NULL

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

G1 := 5*lambda^2*alpha[1]^4*alpha[0]*a[4]+lambda^2*alpha[1]^4*a[3]-10*lambda*alpha[1]^2*alpha[0]^3*a[4]+lambda*k^2*a[1]*alpha[1]^2-6*lambda*alpha[1]^2*alpha[0]^2*a[3]+alpha[0]^5*a[4]-k^2*a[1]*alpha[0]^2-3*lambda*alpha[1]^2*alpha[0]*a[2]+alpha[0]^4*a[3]+lambda*w*alpha[1]^2+alpha[0]^3*a[2]-w*alpha[0]^2+((lambda^2*a[4]*alpha[1]^5-10*lambda*a[4]*alpha[0]^2*alpha[1]^3-4*lambda*a[3]*alpha[0]*alpha[1]^3+5*a[4]*alpha[0]^4*alpha[1]-2*k^2*a[1]*alpha[0]*alpha[1]-lambda*a[2]*alpha[1]^3+4*a[3]*alpha[0]^3*alpha[1]+3*a[2]*alpha[0]^2*alpha[1]-2*w*alpha[0]*alpha[1])*(diff(G(xi), xi))+lambda^2*beta[0]*a[5]*alpha[1]^2-4*mu*lambda*alpha[1]^4*a[3]+5*lambda^2*beta[0]*alpha[1]^4*a[4]-3*lambda*beta[0]*alpha[1]^2*a[2]-lambda*beta[0]*a[5]*alpha[0]^2-(1/2)*lambda*a[1]*alpha[0]*beta[0]-2*k^2*a[1]*alpha[0]*beta[0]+12*mu*alpha[1]^2*alpha[0]^2*a[3]+6*mu*alpha[1]^2*alpha[0]*a[2]-2*mu*k^2*a[1]*alpha[1]^2-(1/2)*mu*lambda*alpha[1]^2*a[1]+20*mu*alpha[1]^2*alpha[0]^3*a[4]-20*mu*lambda*alpha[1]^4*alpha[0]*a[4]-2*mu*lambda*alpha[1]^2*a[5]*alpha[0]-30*lambda*beta[0]*alpha[1]^2*alpha[0]^2*a[4]-12*lambda*beta[0]*alpha[1]^2*alpha[0]*a[3]-2*w*alpha[0]*beta[0]+5*beta[0]*alpha[0]^4*a[4]+4*beta[0]*alpha[0]^3*a[3]+3*beta[0]*alpha[0]^2*a[2]-2*mu*w*alpha[1]^2)/G(xi)+((1/4)*(3*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*alpha[1]^2*a[1]+6*mu*beta[0]*alpha[1]^2*a[2]+3*mu*beta[0]*a[5]*alpha[0]^2-6*lambda*beta[0]^2*alpha[1]^2*a[3]-2*lambda*beta[0]^2*a[5]*alpha[0]+(6*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^2*alpha[0]^2*a[3]+(3*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^2*alpha[0]*a[2]-12*mu^2*alpha[1]^2*a[5]*alpha[0]+3*mu*a[1]*alpha[0]*beta[0]*(1/2)+10*beta[0]^2*alpha[0]^3*a[4]+6*beta[0]^2*alpha[0]^2*a[3]+3*beta[0]^2*alpha[0]*a[2]-k^2*a[1]*beta[0]^2+(10*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^2*alpha[0]^3*a[4]-(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*k^2*a[1]*alpha[1]^2+(5*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*alpha[1]^4*alpha[0]*a[4]+(4*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*alpha[1]^2*a[5]*alpha[0]+(1/2)*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*alpha[1]^2*lambda*a[1]-9*mu^2*alpha[1]^2*a[1]*(1/4)-(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*w*alpha[1]^2+(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2)*alpha[1]^4*a[3]-(1/4)*lambda*beta[0]^2*a[1]-30*lambda*beta[0]^2*alpha[1]^2*alpha[0]*a[4]+24*mu*beta[0]*alpha[1]^2*alpha[0]*a[3]+60*mu*beta[0]*alpha[1]^2*alpha[0]^2*a[4]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^2*lambda*a[5]*alpha[0]-20*mu*lambda*beta[0]*alpha[1]^4*a[4]-7*mu*lambda*beta[0]*a[5]*alpha[1]^2+(2*mu*alpha[1]^3*a[2]-2*w*alpha[1]*beta[0]-4*lambda*beta[0]*alpha[1]^3*a[3]+8*mu*alpha[1]^3*alpha[0]*a[3]+mu*alpha[1]*a[5]*alpha[0]^2+(1/2)*mu*alpha[1]*alpha[0]*a[1]+20*mu*alpha[1]^3*alpha[0]^2*a[4]-4*mu*lambda*alpha[1]^5*a[4]-mu*lambda*alpha[1]^3*a[5]+20*beta[0]*alpha[1]*alpha[0]^3*a[4]+12*beta[0]*alpha[1]*alpha[0]^2*a[3]+6*beta[0]*alpha[1]*alpha[0]*a[2]-2*k^2*a[1]*alpha[1]*beta[0]-(1/2)*lambda*beta[0]*alpha[1]*a[1]-20*lambda*beta[0]*alpha[1]^3*alpha[0]*a[4]-2*lambda*beta[0]*a[5]*alpha[1]*alpha[0])*(diff(G(xi), xi))-w*beta[0]^2)/G(xi)^2+(((lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*alpha[1]^3*a[2]+(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2)*alpha[1]^5*a[4]+(2*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*alpha[1]^3*a[5]+3*beta[0]^2*alpha[1]*a[2]+3*mu*beta[0]*alpha[1]*a[1]*(1/2)+8*mu*beta[0]*alpha[1]^3*a[3]-2*lambda*beta[0]^2*a[5]*alpha[1]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^3*alpha[0]*a[3]+(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]*a[5]*alpha[0]^2+(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*alpha[1]*alpha[0]*a[1]+(10*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^3*alpha[0]^2*a[4]+(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^3*lambda*a[5]+30*beta[0]^2*alpha[1]*alpha[0]^2*a[4]+12*beta[0]^2*alpha[1]*alpha[0]*a[3]-6*mu^2*alpha[1]^3*a[5]-10*lambda*beta[0]^2*alpha[1]^3*a[4]+40*mu*beta[0]*alpha[1]^3*alpha[0]*a[4]+8*mu*beta[0]*a[5]*alpha[1]*alpha[0])*(diff(G(xi), xi))+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^4*a[3]+(5*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*beta[0]*alpha[1]^4*a[4]+(6*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*beta[0]*a[5]*alpha[1]^2-10*lambda*beta[0]^3*alpha[1]^2*a[4]+(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^2*a[1]+(3*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]^2*a[2]+(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*a[5]*alpha[0]^2+(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*a[1]*alpha[0]*beta[0]+12*mu*beta[0]^2*alpha[1]^2*a[3]+6*mu*beta[0]^2*a[5]*alpha[0]+(20*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^4*alpha[0]*a[4]+(10*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^2*a[5]*alpha[0]+beta[0]^3*a[2]-14*mu^2*beta[0]*a[5]*alpha[1]^2+(30*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]^2*alpha[0]^2*a[4]+(5*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*lambda*a[5]*alpha[1]^2+(12*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]^2*alpha[0]*a[3]+60*mu*beta[0]^2*alpha[1]^2*alpha[0]*a[4]+mu*beta[0]^2*a[1]-lambda*beta[0]^3*a[5]+10*beta[0]^3*alpha[0]^2*a[4]+4*beta[0]^3*alpha[0]*a[3])/G(xi)^3+((4*beta[0]^3*alpha[1]*a[3]+(1/2)*(3*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]*a[1]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]^3*a[3]+7*mu*beta[0]^2*a[5]*alpha[1]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^5*a[4]+(5*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^3*a[5]+20*beta[0]^3*alpha[1]*alpha[0]*a[4]+20*mu*beta[0]^2*alpha[1]^3*a[4]+(20*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]^3*alpha[0]*a[4]+(8*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*a[5]*alpha[1]*alpha[0])*(diff(G(xi), xi))+20*mu*beta[0]^3*alpha[1]^2*a[4]+(6*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*alpha[1]^2*a[3]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*a[5]*alpha[0]+5*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^4*alpha[0]*a[4]+4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^2*a[5]*alpha[0]+(17*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*beta[0]*a[5]*alpha[1]^2+(20*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*beta[0]*alpha[1]^4*a[4]+beta[0]^4*a[3]+(30*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*alpha[1]^2*alpha[0]*a[4]+(1/4)*(3*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*a[1]+3*mu*beta[0]^3*a[5]+5*beta[0]^4*alpha[0]*a[4]+(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^4*a[3]+3*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^2*a[1]*(1/4))/G(xi)^4+(((lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^5*a[4]+2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^3*a[5]+5*beta[0]^4*alpha[1]*a[4]+(6*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*a[5]*alpha[1]+(10*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*alpha[1]^3*a[4])*(diff(G(xi), xi))+(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^3*a[5]+(10*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^3*alpha[1]^2*a[4]+5*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*beta[0]*alpha[1]^4*a[4]+6*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*beta[0]*a[5]*alpha[1]^2+beta[0]^5*a[4])/G(xi)^5 = 0

indets(G1)

{k, lambda, mu, w, xi, B[1], B[2], a[1], a[2], a[3], a[4], a[5], alpha[0], alpha[1], beta[0], G(xi), diff(G(xi), xi)}

(2)

``

(3)

eq0 := 5*lambda^2*a[4]*alpha[0]*alpha[1]^4+lambda^2*a[3]*alpha[1]^4-10*lambda*a[4]*alpha[0]^3*alpha[1]^2+k^2*lambda*a[1]*alpha[1]^2-6*lambda*a[3]*alpha[0]^2*alpha[1]^2+a[4]*alpha[0]^5-k^2*a[1]*alpha[0]^2-3*lambda*a[2]*alpha[0]*alpha[1]^2+a[3]*alpha[0]^4+lambda*w*alpha[1]^2+a[2]*alpha[0]^3-w*alpha[0]^2 = 0

``

eq1 := lambda^2*a[4]*alpha[1]^5-10*lambda*a[4]*alpha[0]^2*alpha[1]^3-4*lambda*a[3]*alpha[0]*alpha[1]^3+5*a[4]*alpha[0]^4*alpha[1]-2*k^2*a[1]*alpha[0]*alpha[1]-lambda*a[2]*alpha[1]^3+4*a[3]*alpha[0]^3*alpha[1]+3*a[2]*alpha[0]^2*alpha[1]-2*w*alpha[0]*alpha[1] = 0

eq2 := lambda^2*beta[0]*a[5]*alpha[1]^2+6*mu*alpha[1]^2*alpha[0]*a[2]-2*mu*k^2*a[1]*alpha[1]^2-(1/2)*mu*alpha[1]^2*lambda*a[1]+20*mu*alpha[1]^2*alpha[0]^3*a[4]+12*mu*alpha[1]^2*alpha[0]^2*a[3]-(1/2)*lambda*a[1]*alpha[0]*beta[0]-2*k^2*a[1]*alpha[0]*beta[0]-3*lambda*beta[0]*alpha[1]^2*a[2]-lambda*beta[0]*a[5]*alpha[0]^2+5*lambda^2*beta[0]*alpha[1]^4*a[4]-4*mu*lambda*alpha[1]^4*a[3]-2*mu*w*alpha[1]^2+5*beta[0]*alpha[0]^4*a[4]+4*beta[0]*alpha[0]^3*a[3]+3*beta[0]*alpha[0]^2*a[2]-2*w*alpha[0]*beta[0]-20*mu*lambda*alpha[1]^4*alpha[0]*a[4]-2*mu*alpha[1]^2*lambda*a[5]*alpha[0]-30*lambda*beta[0]*alpha[1]^2*alpha[0]^2*a[4]-12*lambda*beta[0]*alpha[1]^2*alpha[0]*a[3] = 0

NULL

eq3 := (1/4)*(3*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*alpha[1]^2*a[1]-(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*k^2*a[1]*alpha[1]^2+(1/2)*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*alpha[1]^2*lambda*a[1]+(5*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*alpha[1]^4*alpha[0]*a[4]+(10*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^2*alpha[0]^3*a[4]+(6*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^2*alpha[0]^2*a[3]-30*lambda*beta[0]^2*alpha[1]^2*alpha[0]*a[4]-20*mu*beta[0]*lambda*alpha[1]^4*a[4]+(4*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*alpha[1]^2*a[5]*alpha[0]-12*mu^2*alpha[1]^2*a[5]*alpha[0]+(3*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^2*alpha[0]*a[2]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^2*lambda*a[5]*alpha[0]-7*mu*beta[0]*lambda*a[5]*alpha[1]^2+24*mu*beta[0]*alpha[1]^2*alpha[0]*a[3]-9*mu^2*alpha[1]^2*a[1]*(1/4)-w*beta[0]^2+3*beta[0]^2*alpha[0]*a[2]-(1/4)*lambda*beta[0]^2*a[1]-k^2*a[1]*beta[0]^2+10*beta[0]^2*alpha[0]^3*a[4]+6*beta[0]^2*alpha[0]^2*a[3]-(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*w*alpha[1]^2+3*mu*a[1]*alpha[0]*beta[0]*(1/2)+(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2)*alpha[1]^4*a[3]+3*mu*beta[0]*a[5]*alpha[0]^2-6*lambda*beta[0]^2*alpha[1]^2*a[3]-2*lambda*beta[0]^2*a[5]*alpha[0]+6*mu*beta[0]*alpha[1]^2*a[2]+60*mu*beta[0]*alpha[1]^2*alpha[0]^2*a[4] = 0

eq4 := 2*mu*alpha[1]^3*a[2]-2*w*alpha[1]*beta[0]-20*lambda*beta[0]*alpha[1]^3*alpha[0]*a[4]-2*lambda*beta[0]*a[5]*alpha[1]*alpha[0]-2*k^2*a[1]*alpha[1]*beta[0]+20*beta[0]*alpha[1]*alpha[0]^3*a[4]+12*beta[0]*alpha[1]*alpha[0]^2*a[3]+6*beta[0]*alpha[1]*alpha[0]*a[2]+8*mu*alpha[1]^3*alpha[0]*a[3]+mu*alpha[1]*a[5]*alpha[0]^2+(1/2)*mu*alpha[1]*alpha[0]*a[1]-4*lambda*beta[0]*alpha[1]^3*a[3]-lambda*alpha[1]^3*mu*a[5]-(1/2)*lambda*beta[0]*alpha[1]*a[1]+20*mu*alpha[1]^3*alpha[0]^2*a[4]-4*mu*lambda*alpha[1]^5*a[4] = 0

eq5 := -6*mu^2*alpha[1]^3*a[5]+(2*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*alpha[1]^3*a[5]+(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*alpha[1]^3*a[2]+(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2)*alpha[1]^5*a[4]+3*beta[0]^2*alpha[1]*a[2]+40*mu*beta[0]*alpha[1]^3*alpha[0]*a[4]+8*mu*beta[0]*a[5]*alpha[1]*alpha[0]+30*beta[0]^2*alpha[1]*alpha[0]^2*a[4]+12*beta[0]^2*alpha[1]*alpha[0]*a[3]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^3*alpha[0]*a[3]+(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]*a[5]*alpha[0]^2+(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*alpha[1]*alpha[0]*a[1]+8*mu*beta[0]*alpha[1]^3*a[3]+3*mu*beta[0]*alpha[1]*a[1]*(1/2)-10*lambda*beta[0]^2*alpha[1]^3*a[4]-2*lambda*beta[0]^2*a[5]*alpha[1]+(10*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^3*alpha[0]^2*a[4]+(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^3*lambda*a[5] = 0

eq6 := -14*mu^2*beta[0]*a[5]*alpha[1]^2+beta[0]^3*a[2]+(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)*a[1]*alpha[0]*beta[0]+12*mu*beta[0]^2*alpha[1]^2*a[3]+6*mu*beta[0]^2*a[5]*alpha[0]-10*lambda*beta[0]^3*alpha[1]^2*a[4]+(6*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*beta[0]*a[5]*alpha[1]^2+(3*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]^2*a[2]+(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*a[5]*alpha[0]^2+(5*(-(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*lambda+4*mu^2))*beta[0]*alpha[1]^4*a[4]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^4*a[3]+(2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^2*a[1]+10*beta[0]^3*alpha[0]^2*a[4]+4*beta[0]^3*alpha[0]*a[3]-lambda*beta[0]^3*a[5]+mu*beta[0]^2*a[1]+(20*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^4*alpha[0]*a[4]+(10*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^2*a[5]*alpha[0]+(30*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]^2*alpha[0]^2*a[4]+(5*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*lambda*a[5]*alpha[1]^2+(12*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]^2*alpha[0]*a[3]+60*mu*beta[0]^2*alpha[1]^2*alpha[0]*a[4] = 0

eq7 := 4*beta[0]^3*alpha[1]*a[3]+(20*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]^3*alpha[0]*a[4]+(8*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*a[5]*alpha[1]*alpha[0]+20*beta[0]^3*alpha[1]*alpha[0]*a[4]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]^3*a[3]+(5*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*alpha[1]^3*mu*a[5]+(1/2)*(3*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*alpha[1]*a[1]+20*mu*beta[0]^2*alpha[1]^3*a[4]+7*mu*beta[0]^2*a[5]*alpha[1]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*mu*alpha[1]^5*a[4] = 0

eq8 := 4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^2*a[5]*alpha[0]+5*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^4*alpha[0]*a[4]+beta[0]^4*a[3]+(6*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*alpha[1]^2*a[3]+(4*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*a[5]*alpha[0]+20*mu*beta[0]^3*alpha[1]^2*a[4]+(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^4*a[3]+3*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^2*a[1]*(1/4)+5*beta[0]^4*alpha[0]*a[4]+3*mu*beta[0]^3*a[5]+(1/4)*(3*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*a[1]+(30*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*alpha[1]^2*alpha[0]*a[4]+(17*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*mu*a[5]*alpha[1]^2+(20*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]*mu*alpha[1]^4*a[4] = 0

eq9 := (10*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*alpha[1]^3*a[4]+(6*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^2*a[5]*alpha[1]+5*beta[0]^4*alpha[1]*a[4]+(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^5*a[4]+2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*alpha[1]^3*a[5] = 0

eq10 := (2*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^3*a[5]+beta[0]^5*a[4]+5*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*beta[0]*alpha[1]^4*a[4]+6*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda)^2*beta[0]*a[5]*alpha[1]^2+(10*(lambda*B[1]^2-lambda*B[2]^2-mu^2/lambda))*beta[0]^3*alpha[1]^2*a[4] = 0

 

with(LargeExpressions)

COEFFS := solve({eq0, eq1, eq10, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9}, {w, a[1], a[2], alpha[0], alpha[1], beta[0]})

Download by_hand!.mw

Substituting the solutions into the ODE doesn't yield zero, despite the code appearing correct—suggesting either complexity, symbolic limits, or an implementation issue.

 

 

17-ode.mw

 

also in this ode why solution is like this how i can fixed this too

restart

with(PDEtools)

with(LinearAlgebra)

with(Physics)

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

ode := diff(G(xi), xi) = sqrt(C*G(xi)^4+B*G(xi)^2+A)

diff(G(xi), xi) = (C*G(xi)^4+B*G(xi)^2+A)^(1/2)

(2)

dsolve(ode, G(xi))

xi-Intat(1/(C*_a^4+B*_a^2+A)^(1/2), _a = G(xi))+c__1 = 0

(3)
 

NULL

Download v1.mw

I'm trying to collect all terms involving the expression diff(G(xi), xi)/G(xi) in a symbolic equation using . While it's straightforward to do this by hand, I want to automate it in code — ideally by extracting the coefficient of this entire expression directly. However, when I use collect, Maple treats diff(G(xi), xi) and G(xi) separately, and I can't seem to group terms properly by the full ratio diff(G(xi), xi)/G(xi).

Is there a clean way or built-in Maple function to automatically collect or isolate the coefficient of diff(G(xi), xi)/G(xi) as a whole, without having to manually substitute or restructure the expression?

restart

with(PDEtools)

with(LinearAlgebra)

NULL

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

L := (-6*f^3*g*a[2]+6*f*g^3*a[2]+6*f*g*a[2]^2)*(diff(G(xi), xi))^4/G(xi)^4+(f*g^3*(10*lambda*a[2]+2*a[1])-f^3*g*(10*lambda*a[2]+2*a[1])+12*f*g*a[1]*a[2])*(diff(G(xi), xi))^3/G(xi)^3+(f*g^3*(2*a[2]*lambda^2+3*lambda*a[1]+10*mu*a[2]+2*a[2]*(lambda^2-mu))-4*f*p*a[2]-6*k*l*a[2]-f^3*g*(2*a[2]*lambda^2+3*lambda*a[1]+10*mu*a[2]+2*a[2]*(lambda^2-mu))+6*f*g*(2*a[0]*a[2]+a[1]^2))*(diff(G(xi), xi))^2/G(xi)^2+(f*g^3*(a[1]*(lambda^2-mu)+3*a[1]*mu+6*lambda*a[2]*mu)-4*f*p*a[1]-6*k*l*a[1]-f^3*g*(a[1]*(lambda^2-mu)+3*a[1]*mu+6*lambda*a[2]*mu)+12*f*g*a[0]*a[1])*(diff(G(xi), xi))/G(xi)+f*g^3*(lambda*mu*a[1]+2*mu^2*a[2])-4*f*p*a[0]-6*k*l*a[0]-f^3*g*(lambda*mu*a[1]+2*mu^2*a[2])+6*f*g*a[0]^2 = 0

(-6*f^3*g*a[2]+6*f*g^3*a[2]+6*f*g*a[2]^2)*(diff(G(xi), xi))^4/G(xi)^4+(f*g^3*(10*lambda*a[2]+2*a[1])-f^3*g*(10*lambda*a[2]+2*a[1])+12*f*g*a[1]*a[2])*(diff(G(xi), xi))^3/G(xi)^3+(f*g^3*(3*lambda*a[1]+2*a[2]*lambda^2+10*mu*a[2]+2*a[2]*(lambda^2-mu))-4*f*p*a[2]-6*k*l*a[2]-f^3*g*(3*lambda*a[1]+2*a[2]*lambda^2+10*mu*a[2]+2*a[2]*(lambda^2-mu))+6*f*g*(2*a[0]*a[2]+a[1]^2))*(diff(G(xi), xi))^2/G(xi)^2+(f*g^3*(a[1]*(lambda^2-mu)+3*a[1]*mu+6*a[2]*lambda*mu)-4*f*p*a[1]-6*k*l*a[1]-f^3*g*(a[1]*(lambda^2-mu)+3*a[1]*mu+6*a[2]*lambda*mu)+12*f*g*a[0]*a[1])*(diff(G(xi), xi))/G(xi)+f*g^3*(lambda*mu*a[1]+2*mu^2*a[2])-4*f*p*a[0]-6*k*l*a[0]-f^3*g*(lambda*mu*a[1]+2*mu^2*a[2])+6*f*g*a[0]^2 = 0

(2)

``

(3)

collect(%, {1/(diff(G(xi), xi)), G(xi)})

(-6*f^3*g*a[2]+6*f*g^3*a[2]+6*f*g*a[2]^2)*(diff(G(xi), xi))^4/G(xi)^4+(f*g^3*(10*lambda*a[2]+2*a[1])-f^3*g*(10*lambda*a[2]+2*a[1])+12*f*g*a[1]*a[2])*(diff(G(xi), xi))^3/G(xi)^3+(f*g^3*(3*lambda*a[1]+2*a[2]*lambda^2+10*mu*a[2]+2*a[2]*(lambda^2-mu))-4*f*p*a[2]-6*k*l*a[2]-f^3*g*(3*lambda*a[1]+2*a[2]*lambda^2+10*mu*a[2]+2*a[2]*(lambda^2-mu))+6*f*g*(2*a[0]*a[2]+a[1]^2))*(diff(G(xi), xi))^2/G(xi)^2+(f*g^3*(a[1]*(lambda^2-mu)+3*a[1]*mu+6*a[2]*lambda*mu)-4*f*p*a[1]-6*k*l*a[1]-f^3*g*(a[1]*(lambda^2-mu)+3*a[1]*mu+6*a[2]*lambda*mu)+12*f*g*a[0]*a[1])*(diff(G(xi), xi))/G(xi)+f*g^3*(lambda*mu*a[1]+2*mu^2*a[2])-4*f*p*a[0]-6*k*l*a[0]-f^3*g*(lambda*mu*a[1]+2*mu^2*a[2])+6*f*g*a[0]^2 = 0

(4)

L1 := %

num := numer(lhs(L)); num := expand(num); num_collected := collect(num, [1/(diff(G(xi), xi)), G(xi)]); eqs := [seq(coeff(num_collected, {1/(diff(G(xi), xi)), G(xi)}, i) = 0, i = 0 .. 8)]

[(-6*f^3*g*a[2]+6*f*g^3*a[2]+6*f*g*a[2]^2)*(diff(G(xi), xi))^4+(-10*f^3*g*lambda*a[2]+10*f*g^3*lambda*a[2]-2*f^3*g*a[1]+2*f*g^3*a[1]+12*f*g*a[1]*a[2])*G(xi)*(diff(G(xi), xi))^3+(-4*f^3*g*lambda^2*a[2]+4*f*g^3*lambda^2*a[2]-3*f^3*g*lambda*a[1]-8*f^3*g*mu*a[2]+3*f*g^3*lambda*a[1]+8*f*g^3*mu*a[2]+12*f*g*a[0]*a[2]+6*f*g*a[1]^2-4*f*p*a[2]-6*k*l*a[2])*G(xi)^2*(diff(G(xi), xi))^2+(-f^3*g*lambda^2*a[1]-6*f^3*g*lambda*mu*a[2]+f*g^3*lambda^2*a[1]+6*f*g^3*lambda*mu*a[2]-2*f^3*g*mu*a[1]+2*f*g^3*mu*a[1]+12*f*g*a[0]*a[1]-4*f*p*a[1]-6*k*l*a[1])*G(xi)^3*(diff(G(xi), xi))+(-f^3*g*lambda*mu*a[1]-2*f^3*g*mu^2*a[2]+f*g^3*lambda*mu*a[1]+2*f*g^3*mu^2*a[2]+6*f*g*a[0]^2-4*f*p*a[0]-6*k*l*a[0])*G(xi)^4 = 0, 0 = 0, 0 = 0, 0 = 0, 0 = 0, 0 = 0, 0 = 0, 0 = 0, 0 = 0]

(5)

Download collect-coe.mw

I’m trying to verify a solution given in the form from  using Maple's odeTest. Even though the paper claims the solution satisfies the ODE, Maple does not simplify the result to zero. Could someone explain why the test fails or suggest the correct way to verify it in Maple?

restart

with(PDEtools)

with(LinearAlgebra)

with(Physics)

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

_local(gamma)

Warning, A new binding for the name `gamma` has been created. The global instance of this name is still accessible using the :- prefix, :-`gamma`.  See ?protect for details.

 

declare(u(x, t)); declare(U(xi)); declare(u(x, y, z, t)); declare(Q(xi)); declare(V(xi))

u(x, t)*`will now be displayed as`*u

 

U(xi)*`will now be displayed as`*U

 

u(x, y, z, t)*`will now be displayed as`*u

 

Q(xi)*`will now be displayed as`*Q

 

V(xi)*`will now be displayed as`*V

(2)

NULL

ode := diff(G(xi), xi) = G(xi)*sqrt(A+B*G(xi)^2)

diff(G(xi), xi) = G(xi)*(A+B*G(xi)^2)^(1/2)

(3)

S1 := G(xi) = -sqrt(A/B)*csch(sqrt(A)*(xi+xi[0]))

G(xi) = -(A/B)^(1/2)*csch(A^(1/2)*(xi+xi[0]))

(4)

res := simplify(odetest(S1, ode))

(A/B)^(1/2)*csch(A^(1/2)*(xi+xi[0]))*(A^(1/2)*coth(A^(1/2)*(xi+xi[0]))+(A*coth(A^(1/2)*(xi+xi[0]))^2)^(1/2))

(5)

S2 := G(xi) = sqrt(-A/B)*sec(sqrt(-A)*(xi+xi[0]))

G(xi) = (-A/B)^(1/2)*sec((-A)^(1/2)*(xi+xi[0]))

(6)

res := simplify(odetest(S2, ode))

(-A/B)^(1/2)*sec((-A)^(1/2)*(xi+xi[0]))*((-A)^(1/2)*tan((-A)^(1/2)*(xi+xi[0]))-(-A*tan((-A)^(1/2)*(xi+xi[0]))^2)^(1/2))

(7)

S3 := G(xi) = sqrt(-A/B)*sech(sqrt(A)*(xi+xi[0]))

G(xi) = (-A/B)^(1/2)*sech(A^(1/2)*(xi+xi[0]))

(8)

res := simplify(odetest(S3, ode))

(-A^(1/2)*tanh(A^(1/2)*(xi+xi[0]))-(A*tanh(A^(1/2)*(xi+xi[0]))^2)^(1/2))*(-A/B)^(1/2)*sech(A^(1/2)*(xi+xi[0]))

(9)

S4 := G(xi) = sqrt(-A/B)*csc(sqrt(-A)*(xi+xi[0]))

G(xi) = (-A/B)^(1/2)*csc((-A)^(1/2)*(xi+xi[0]))

(10)

res := simplify(odetest(S4, ode))

(-(-A)^(1/2)*cot((-A)^(1/2)*(xi+xi[0]))-(-A*cot((-A)^(1/2)*(xi+xi[0]))^2)^(1/2))*(-A/B)^(1/2)*csc((-A)^(1/2)*(xi+xi[0]))

(11)

S5 := G(xi) = cos(sqrt(-A)*(xi+xi[0]))+sin(sqrt(-A)*(xi+xi[0]))

G(xi) = cos((-A)^(1/2)*(xi+xi[0]))+sin((-A)^(1/2)*(xi+xi[0]))

(12)

res := simplify(odetest(S5, ode))

(cos((-A)^(1/2)*(xi+xi[0]))-sin((-A)^(1/2)*(xi+xi[0])))*(-A)^(1/2)+(B*sin(2*(-A)^(1/2)*(xi+xi[0]))+A+B)^(1/2)*(-cos((-A)^(1/2)*(xi+xi[0]))-sin((-A)^(1/2)*(xi+xi[0])))

(13)

S6 := G(xi) = 1/(sqrt(B)*(xi+xi[0]))

G(xi) = 1/(B^(1/2)*(xi+xi[0]))

(14)

odetest(S6, subs(A = 0, ode))

-csgn(1/(xi+xi[0]))/(B^(1/2)*(xi+xi[0])^2)-1/(B^(1/2)*(xi+xi[0])^2)

(15)

S7 := G(xi) = 1/(sqrt(-B)*(xi+xi[0]))

G(xi) = 1/((-B)^(1/2)*(xi+xi[0]))

(16)

odetest(S7, subs(A = 0, ode))

-(-1/(xi+xi[0])^2)^(1/2)*xi[0]/((-B)^(1/2)*(xi+xi[0])^2)-(-1/(xi+xi[0])^2)^(1/2)*xi/((-B)^(1/2)*(xi+xi[0])^2)-1/((-B)^(1/2)*(xi+xi[0])^2)

(17)

ode2 := diff(G(xi), xi) = A+B*G(xi)^2

diff(G(xi), xi) = A+B*G(xi)^2

(18)

S8 := G(xi) = sgn(A)*sqrt(A/B)*tan(sqrt(A*B)*(xi+xi[0]))

G(xi) = sgn(A)*(A/B)^(1/2)*tan((A*B)^(1/2)*(xi+xi[0]))

(19)

res := simplify(odetest(S8, ode2))

(sgn(A)*(A*B)^(1/2)*(tan((A*B)^(1/2)*xi[0])^2+1)*(tan((A*B)^(1/2)*xi)^2+1)*(A/B)^(1/2)-A*((tan((A*B)^(1/2)*xi)+tan((A*B)^(1/2)*xi[0]))^2*sgn(A)^2+(tan((A*B)^(1/2)*xi[0])*tan((A*B)^(1/2)*xi)-1)^2))/(tan((A*B)^(1/2)*xi[0])*tan((A*B)^(1/2)*xi)-1)^2

(20)

NULL

S9 := G(xi) = -sgn(A)*sqrt(A/B)*cot(sqrt(A*B)*(xi+xi[0]))

G(xi) = -sgn(A)*(A/B)^(1/2)*cot((A*B)^(1/2)*(xi+xi[0]))

(21)

res := simplify(odetest(S9, ode2))

(sgn(A)*(A*B)^(1/2)*(cot((A*B)^(1/2)*xi[0])^2+1)*(cot((A*B)^(1/2)*xi)^2+1)*(A/B)^(1/2)-A*((cot((A*B)^(1/2)*xi[0])*cot((A*B)^(1/2)*xi)-1)^2*sgn(A)^2+(cot((A*B)^(1/2)*xi[0])+cot((A*B)^(1/2)*xi))^2))/(cot((A*B)^(1/2)*xi[0])+cot((A*B)^(1/2)*xi))^2

(22)

NULL

S9 := G(xi) = sgn(A)*sqrt(-A/B)*tanh(sqrt(-A*B)*(xi+xi[0]))

G(xi) = sgn(A)*(-A/B)^(1/2)*tanh((-A*B)^(1/2)*(xi+xi[0]))

(23)

res := simplify(odetest(S9, ode2))

(A*(sgn(A)^2-1)*cosh(2*(-A*B)^(1/2)*(xi+xi[0]))-sgn(A)^2*A+2*sgn(A)*(-A/B)^(1/2)*(-A*B)^(1/2)-A)/(1+cosh(2*(-A*B)^(1/2)*(xi+xi[0])))

(24)

NULL

S10 := G(xi) = sgn(A)*sqrt(-A/B)*coth(sqrt(-A*B)*(xi+xi[0]))

G(xi) = sgn(A)*(-A/B)^(1/2)*coth((-A*B)^(1/2)*(xi+xi[0]))

(25)

odetest(S10, ode2)

(sgn(A)^2*A*cosh(2*(-A*B)^(1/2)*(xi+xi[0]))+sgn(A)^2*A-2*sgn(A)*(-A/B)^(1/2)*(-A*B)^(1/2)-A*cosh(2*(-A*B)^(1/2)*(xi+xi[0]))+A)/(-1+cosh(2*(-A*B)^(1/2)*(xi+xi[0])))

(26)

NULL

S11 := G(xi) = -1/(B*(xi+xi[0]))

G(xi) = -1/(B*(xi+xi[0]))

(27)

odetest(S11, subs(A = 0, ode2))

0

(28)

S12 := G(xi) = A*(xi+xi[0])

G(xi) = A*(xi+xi[0])

(29)

odetest(S12, subs(B = 0, ode2))

0

(30)

Download Z1.mw

Why doesn't this ODE return zero when using odetest? Did I do something wrong?

Would you like help checking the equation or debugging the issue?

restart

with(PDEtools)

declare(P(mu))

P(mu)*`will now be displayed as`*P

(1)

assume(A::real, r::real, rho::real, lambda::real)

Psol := P(mu) = 2*A*lambda/((A^2+r)*exp(lambda*rho*mu)+r*exp(-lambda*rho*mu))

P(mu) = 2*A*lambda/((A^2+r)*exp(lambda*rho*mu)+r*exp(-lambda*rho*mu))

(2)

ode := (diff(P(mu), mu))^2-rho^2*P(mu)^2*(1+r*P(mu)^2) = 0

(diff(P(mu), mu))^2-rho^2*P(mu)^2*(1+r*P(mu)^2) = 0

(3)

res := odetest(Psol, ode)

4*A^6*rho^2*lambda^4*exp(6*lambda*rho*mu)/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4+8*r*A^4*rho^2*lambda^4*exp(6*lambda*rho*mu)/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4-24*A^4*lambda^4*rho^2*exp(4*lambda*rho*mu)*r/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4-4*A^6*rho^2*lambda^2*exp(6*lambda*rho*mu)/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4+4*A^2*r^2*rho^2*lambda^4*exp(6*lambda*rho*mu)/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4-8*A^2*lambda^4*rho^2*exp(4*lambda*rho*mu)*r^2/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4+4*A^2*r^2*rho^2*lambda^4*exp(2*lambda*rho*mu)/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4-8*r*A^4*rho^2*lambda^2*exp(6*lambda*rho*mu)/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4-8*A^4*lambda^2*rho^2*exp(4*lambda*rho*mu)*r/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4-4*A^2*r^2*rho^2*lambda^2*exp(6*lambda*rho*mu)/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4-8*A^2*lambda^2*rho^2*exp(4*lambda*rho*mu)*r^2/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4-4*A^2*r^2*rho^2*lambda^2*exp(2*lambda*rho*mu)/(A^2*exp(2*lambda*rho*mu)+exp(2*lambda*rho*mu)*r+r)^4

(4)

simplify(res)

4*lambda^2*exp(2*lambda*rho*mu)*(((-2*lambda^2-2)*r^2+(-6*lambda^2-2)*A^2*r)*exp(2*lambda*rho*mu)+(lambda+1)*((A^2+r)^2*exp(4*lambda*rho*mu)+r^2)*(lambda-1))*A^2*rho^2/((A^2+r)*exp(2*lambda*rho*mu)+r)^4

(5)

P_hyper := P(mu) = 2*A*lambda/((A^2+r)*cosh(rho*mu)+(A^2-r)*sinh(rho*mu))

P(mu) = 2*A*lambda/((A^2+r)*cosh(rho*mu)+(A^2-r)*sinh(rho*mu))

(6)

res_hyper := simplify(odetest(P_hyper, ode), symbolic)

-16*A^4*lambda^2*rho^2*r*(lambda^2+1)/((A^2+r)*cosh(rho*mu)+(A^2-r)*sinh(rho*mu))^4

(7)
 

NULL

Download ode.mw

I’m trying to test a specific function as a solution to a nonlinear ODE in Maple. The equation is of the Riccati type, and my candidate solution involves parameters A, B, and C.

I've used assuming to specify the condition (4AC−B2)>0 and (4AC - B^2) <0, but when I use odetest to verify the solution, I still get a nonzero result. Additionally, when I apply the assumption, Maple sometimes introduces a negation sign in the output (e.g., changing sqrt(...) into -sqrt(...)), which wasn't part of the original solution.

restart

with(PDEtools)

with(LinearAlgebra)

with(Physics)

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

E := diff(G(xi), xi) = A+B*G(xi)+C*G(xi)^2

diff(G(xi), xi) = A+B*G(xi)+C*G(xi)^2

(2)

S1 := G(xi) = (sqrt(4*A*C-B^2)*tan((1/2)*sqrt(4*A*C-B^2)*(d[0]+xi))-B)/(2*C)

G(xi) = (1/2)*((4*A*C-B^2)^(1/2)*tan((1/2)*(4*A*C-B^2)^(1/2)*(d[0]+xi))-B)/C

(3)

odetest(S1, E)

0

(4)

S2 := G(xi) = -(sqrt(4*A*C-B^2)*cot((1/2)*sqrt(4*A*C-B^2)*(d[0]+xi))+B)/(2*C)

G(xi) = -(1/2)*((4*A*C-B^2)^(1/2)*cot((1/2)*(4*A*C-B^2)^(1/2)*(d[0]+xi))+B)/C

(5)

odetest(S2, E)

0

(6)

assume(4*A*C-B^2 < 0)

S3 := G(xi) = -(sqrt(4*A*C-B^2)*tanh((1/2)*sqrt(4*A*C-B^2)*(d[0]+xi))+B)/(2*C)

G(xi) = -(1/2)*((4*A*C-B^2)^(1/2)*tanh((1/2)*(4*A*C-B^2)^(1/2)*(d[0]+xi))+B)/C

(7)

odetest(S3, E)

-2*A+(1/2)*B^2/C

(8)

Download A2.mw

i did every thing coreectly but nothing happen not apply where is my mistake?

``

restart

with(PDEtools)

with(LinearAlgebra)

with(Physics)

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

NULL

S := (diff(G(xi), xi))^2-r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2) = 0

(diff(G(xi), xi))^2-r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2) = 0

(2)

SS := diff(G(xi), xi) = sqrt(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))

diff(G(xi), xi) = (r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2)

(3)

Se := sqrt(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2)) = diff(G(xi), xi)

(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2) = diff(G(xi), xi)

(4)

dub := diff(SS, xi)

diff(diff(G(xi), xi), xi) = (1/2)*(2*r^2*G(xi)*(a+b*G(xi)+l*G(xi)^2)*(diff(G(xi), xi))+r^2*G(xi)^2*(b*(diff(G(xi), xi))+2*l*G(xi)*(diff(G(xi), xi))))/(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2)

(5)

Dubl2 := simplify(diff(diff(G(xi), xi), xi) = (1/2)*(2*r^2*G(xi)*(a+b*G(xi)+l*G(xi)^2)*(diff(G(xi), xi))+r^2*G(xi)^2*(b*(diff(G(xi), xi))+2*l*G(xi)*(diff(G(xi), xi))))/(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2))

diff(diff(G(xi), xi), xi) = (1/2)*r^2*G(xi)*(diff(G(xi), xi))*(4*l*G(xi)^2+3*b*G(xi)+2*a)/(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2)

(6)

subs(SA, Dubl2)

diff((r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2), xi) = (1/2)*r^2*G(xi)*(4*l*G(xi)^2+3*b*G(xi)+2*a)

(7)

subs(Se, Dubl2)

diff(diff(G(xi), xi), xi) = (1/2)*r^2*G(xi)*(diff(G(xi), xi))*(4*l*G(xi)^2+3*b*G(xi)+2*a)/(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2)

(8)

subs(lhs(Se) = rhs(Se), Dubl2)

diff(diff(G(xi), xi), xi) = (1/2)*r^2*G(xi)*(diff(G(xi), xi))*(4*l*G(xi)^2+3*b*G(xi)+2*a)/(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2)

(9)
 

NULL

Download subs.mw

1 2 3 4 5 6 7 Last Page 3 of 42