Question: there is any way to use explicate for remove the name of expression in equations?

I want to remove the Lambert function (LambertW) from my equation, but I don't know how. I tried using the explicit option, but it didn't work. How can I express the equation without LambertW?

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, y, z, t))

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

(2)

declare(f(x, y, z, t))

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

(3)

pde := diff(diff(u(x, y, z, t), t)+6*u(x, y, z, t)*(diff(u(x, y, z, t), x))+diff(u(x, y, z, t), `$`(x, 3)), x)-lambda*(diff(u(x, y, z, t), `$`(y, 2)))+diff(alpha*(diff(u(x, y, z, t), x))+beta*(diff(u(x, y, z, t), y))+gamma*(diff(u(x, y, z, t), z)), x)

diff(diff(u(x, y, z, t), t), x)+6*(diff(u(x, y, z, t), x))^2+6*u(x, y, z, t)*(diff(diff(u(x, y, z, t), x), x))+diff(diff(diff(diff(u(x, y, z, t), x), x), x), x)-lambda*(diff(diff(u(x, y, z, t), y), y))+alpha*(diff(diff(u(x, y, z, t), x), x))+beta*(diff(diff(u(x, y, z, t), x), y))+gamma*(diff(diff(u(x, y, z, t), x), z))

(4)

pde_nonlinear, pde_linear := selectremove(proc (term) options operator, arrow; not has((eval(term, u(x, y, t) = a*u(x, y, t)))/a, a) end proc, expand(pde))

0, diff(diff(u(x, y, z, t), t), x)+6*(diff(u(x, y, z, t), x))^2+6*u(x, y, z, t)*(diff(diff(u(x, y, z, t), x), x))+diff(diff(diff(diff(u(x, y, z, t), x), x), x), x)-lambda*(diff(diff(u(x, y, z, t), y), y))+alpha*(diff(diff(u(x, y, z, t), x), x))+beta*(diff(diff(u(x, y, z, t), x), y))+gamma*(diff(diff(u(x, y, z, t), x), z))

(5)

thetai := t*w[i]+x*k[i]+y*l[i]+z*r[i]+eta[i]; eval(pde_linear, u(x, y, z, t) = exp(thetai)); eq15 := isolate(%, w[i])

t*w[i]+x*k[i]+y*l[i]+z*r[i]+eta[i]

 

w[i]*k[i]*exp(t*w[i]+x*k[i]+y*l[i]+z*r[i]+eta[i])+12*k[i]^2*(exp(t*w[i]+x*k[i]+y*l[i]+z*r[i]+eta[i]))^2+k[i]^4*exp(t*w[i]+x*k[i]+y*l[i]+z*r[i]+eta[i])-lambda*l[i]^2*exp(t*w[i]+x*k[i]+y*l[i]+z*r[i]+eta[i])+alpha*k[i]^2*exp(t*w[i]+x*k[i]+y*l[i]+z*r[i]+eta[i])+beta*k[i]*l[i]*exp(t*w[i]+x*k[i]+y*l[i]+z*r[i]+eta[i])+gamma*k[i]*r[i]*exp(t*w[i]+x*k[i]+y*l[i]+z*r[i]+eta[i])

 

w[i] = -(t*k[i]^4+gamma*t*k[i]*r[i]+alpha*t*k[i]^2+beta*t*k[i]*l[i]-lambda*t*l[i]^2+LambertW(12*t*k[i]*exp(-(t*k[i]^4+alpha*t*k[i]^2+beta*t*k[i]*l[i]+gamma*t*k[i]*r[i]-lambda*t*l[i]^2-x*k[i]^2-y*k[i]*l[i]-z*k[i]*r[i]-eta[i]*k[i])/k[i]))*k[i])/(t*k[i])

(6)

sol := solve(eq15, w[i], explicit)

-(t*k[i]^4+gamma*t*k[i]*r[i]+alpha*t*k[i]^2+beta*t*k[i]*l[i]-lambda*t*l[i]^2+LambertW(12*t*k[i]*exp(-(t*k[i]^4+alpha*t*k[i]^2+beta*t*k[i]*l[i]+gamma*t*k[i]*r[i]-lambda*t*l[i]^2-x*k[i]^2-y*k[i]*l[i]-z*k[i]*r[i]-eta[i]*k[i])/k[i]))*k[i])/(t*k[i])

(7)
 

NULL

Download remove.mw

Please Wait...