MaplePrimes Questions

the second ode is giving me zero also when we back to orginal under the condition by using them must the orginal ode be zero but i don't know where is mistake , when Orginal paper use some thing different but i think they must have same results i don't know i use them wrong i am not sure at here just , when U(xi)=y(z) in my mw

restart

with(PDEtools)

NULL

undeclare(prime)

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

(1)

NULL

G := V(xi) = RootOf(3*_Z^2-3*_Z-1)*B[1]+B[1]*(exp(xi)+exp(-xi))/(exp(xi)-exp(-xi))

V(xi) = RootOf(3*_Z^2-3*_Z-1)*B[1]+B[1]*(exp(xi)+exp(-xi))/(exp(xi)-exp(-xi))

(2)

NULL

p := 2*k

2*k

(3)

ode := I*(-(diff(U(xi), xi))*p*exp(I*(k*x-t*w))-I*U(xi)*w*exp(I*(k*x-t*w)))+(diff(diff(U(xi), xi), xi))*exp(I*(k*x-t*w))+(2*I)*(diff(U(xi), xi))*k*exp(I*(k*x-t*w))-U(xi)*k^2*exp(I*(k*x-t*w))+eta*U(xi)*exp(I*(k*x-t*w))+beta*U(xi)^n*U(xi)*exp(I*(k*x-t*w))+gamma*U(xi)^(2*n)*U(xi)*exp(I*(k*x-t*w))+delta*U(xi)^(3*n)*U(xi)*exp(I*(k*x-t*w))+lambda*U(xi)^(4*n)*U(xi)*exp(I*(k*x-t*w)) = 0

I*(-2*(diff(U(xi), xi))*k*exp(I*(k*x-t*w))-I*U(xi)*w*exp(I*(k*x-t*w)))+(diff(diff(U(xi), xi), xi))*exp(I*(k*x-t*w))+(2*I)*(diff(U(xi), xi))*k*exp(I*(k*x-t*w))-U(xi)*k^2*exp(I*(k*x-t*w))+eta*U(xi)*exp(I*(k*x-t*w))+beta*U(xi)^n*U(xi)*exp(I*(k*x-t*w))+gamma*U(xi)^(2*n)*U(xi)*exp(I*(k*x-t*w))+delta*U(xi)^(3*n)*U(xi)*exp(I*(k*x-t*w))+lambda*U(xi)^(4*n)*U(xi)*exp(I*(k*x-t*w)) = 0

(4)

case1 := [beta = 2*RootOf(3*_Z^2-3*_Z-1)*(n+2)/(B[1]*n^2), delta = 2*B[1]*(RootOf(3*_Z^2-3*_Z-1)+1)*(3*n+2)/(3*n^2), eta = (k^2*n^2*B[1]^2-n^2*w*B[1]^2-1)/(n^2*B[1]^2), gamma = -6*RootOf(3*_Z^2-3*_Z-1)*(n+1)/n^2, lambda = B[1]^2*(3*RootOf(3*_Z^2-3*_Z-1)-7)*(2*n+1)/(9*n^2), A[0] = RootOf(3*_Z^2-3*_Z-1)*B[1], A[1] = 0, B[1] = B[1]]

[beta = 2*RootOf(3*_Z^2-3*_Z-1)*(n+2)/(B[1]*n^2), delta = (2/3)*B[1]*(RootOf(3*_Z^2-3*_Z-1)+1)*(3*n+2)/n^2, eta = (k^2*n^2*B[1]^2-n^2*w*B[1]^2-1)/(n^2*B[1]^2), gamma = -6*RootOf(3*_Z^2-3*_Z-1)*(n+1)/n^2, lambda = (1/9)*B[1]^2*(3*RootOf(3*_Z^2-3*_Z-1)-7)*(2*n+1)/n^2, A[0] = RootOf(3*_Z^2-3*_Z-1)*B[1], A[1] = 0, B[1] = B[1]]

(5)

n := 1

1

(6)

G := U(xi) = (B[1]*(RootOf(3*_Z^2-3*_Z-1)+coth(xi)))^(-1/n)

U(xi) = 1/(B[1]*(RootOf(3*_Z^2-3*_Z-1)+coth(xi)))

(7)

pde3 := eval(ode, case1)

I*(-2*(diff(U(xi), xi))*k*exp(I*(k*x-t*w))-I*U(xi)*w*exp(I*(k*x-t*w)))+(diff(diff(U(xi), xi), xi))*exp(I*(k*x-t*w))+(2*I)*(diff(U(xi), xi))*k*exp(I*(k*x-t*w))-U(xi)*k^2*exp(I*(k*x-t*w))+(k^2*B[1]^2-w*B[1]^2-1)*U(xi)*exp(I*(k*x-t*w))/B[1]^2+6*RootOf(3*_Z^2-3*_Z-1)*U(xi)^2*exp(I*(k*x-t*w))/B[1]-12*RootOf(3*_Z^2-3*_Z-1)*U(xi)^3*exp(I*(k*x-t*w))+(10/3)*B[1]*(RootOf(3*_Z^2-3*_Z-1)+1)*U(xi)^4*exp(I*(k*x-t*w))+(1/3)*B[1]^2*(3*RootOf(3*_Z^2-3*_Z-1)-7)*U(xi)^5*exp(I*(k*x-t*w)) = 0

(8)

odetest(eval(G, case1), pde3)

79584*exp(I*k*x-I*t*w+12*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-127440*exp(I*k*x-I*t*w+10*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+8352*exp(I*k*x-I*t*w+6*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-27792*exp(I*k*x-I*t*w+8*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-24*exp(2*xi-I*t*w+I*k*x)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+4752*exp(I*k*x-I*t*w+4*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+288*exp(2*xi-I*t*w+I*k*x)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-479376*exp(I*k*x-I*t*w+12*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+138240*exp(I*k*x-I*t*w+10*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+70560*exp(18*xi+I*k*x-I*t*w)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-492912*exp(I*k*x-I*t*w+16*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+777888*exp(I*k*x-I*t*w+14*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+16608*exp(I*k*x-I*t*w+8*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-1056*exp(I*k*x-I*t*w+6*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+288*exp(I*k*x-I*t*w+4*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-39000*exp(18*xi+I*k*x-I*t*w)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-50400*exp(I*k*x-I*t*w+16*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+121440*exp(I*k*x-I*t*w+14*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+27000*RootOf(3*_Z^2-3*_Z-1)*exp(18*xi+I*k*x-I*t*w)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-55080*RootOf(3*_Z^2-3*_Z-1)*exp(18*xi+I*k*x-I*t*w)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+7200*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+16*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+394416*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+16*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-205920*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+14*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-609984*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+14*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-244512*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+12*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+366768*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+12*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-42480*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+10*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-144720*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+10*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-48672*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+8*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+8208*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+8*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+9504*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+6*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+20736*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+6*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+288*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+4*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+18576*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+4*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-72*RootOf(3*_Z^2-3*_Z-1)*exp(2*xi-I*t*w+I*k*x)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+1080*RootOf(3*_Z^2-3*_Z-1)*exp(2*xi-I*t*w+I*k*x)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))

(9)

simplify(-492912*exp(I*k*x-I*t*w+16*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+777888*exp(I*k*x-I*t*w+14*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+16608*exp(I*k*x-I*t*w+8*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-1056*exp(I*k*x-I*t*w+6*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+288*exp(I*k*x-I*t*w+4*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-39000*exp(18*xi+I*k*x-I*t*w)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-50400*exp(I*k*x-I*t*w+16*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+121440*exp(I*k*x-I*t*w+14*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+79584*exp(I*k*x-I*t*w+12*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-127440*exp(I*k*x-I*t*w+10*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+8352*exp(I*k*x-I*t*w+6*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-27792*exp(I*k*x-I*t*w+8*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-24*exp(2*xi-I*t*w+I*k*x)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+4752*exp(I*k*x-I*t*w+4*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+288*exp(2*xi-I*t*w+I*k*x)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-479376*exp(I*k*x-I*t*w+12*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+138240*exp(I*k*x-I*t*w+10*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+70560*exp(18*xi+I*k*x-I*t*w)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-55080*RootOf(3*_Z^2-3*_Z-1)*exp(18*xi+I*k*x-I*t*w)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+7200*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+16*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+394416*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+16*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-205920*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+14*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-609984*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+14*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-244512*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+12*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+366768*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+12*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-42480*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+10*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-144720*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+10*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-48672*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+8*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+8208*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+8*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+9504*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+6*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+20736*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+6*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+288*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+4*xi)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+18576*RootOf(3*_Z^2-3*_Z-1)*exp(I*k*x-I*t*w+4*xi)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))-72*RootOf(3*_Z^2-3*_Z-1)*exp(2*xi-I*t*w+I*k*x)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+1080*RootOf(3*_Z^2-3*_Z-1)*exp(2*xi-I*t*w+I*k*x)/(B[1]^3*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1))+27000*RootOf(3*_Z^2-3*_Z-1)*exp(18*xi+I*k*x-I*t*w)/(B[1]*(3125*exp(20*xi)+25000*exp(18*xi)+76875*exp(16*xi)+108000*exp(14*xi)+55650*exp(12*xi)-12432*exp(10*xi)-11130*exp(8*xi)+4320*exp(6*xi)-615*exp(4*xi)+40*exp(2*xi)-1)))

(((244512*B[1]^2-366768)*exp(10*xi)+(205920*B[1]^2+609984)*exp(12*xi)+(-7200*B[1]^2-394416)*exp(14*xi)+42480*exp(8*xi)*B[1]^2-27000*exp(16*xi)*B[1]^2-288*exp(2*xi)*B[1]^2-9504*exp(4*xi)*B[1]^2+48672*exp(6*xi)*B[1]^2+72*B[1]^2+144720*exp(8*xi)+55080*exp(16*xi)-18576*exp(2*xi)-20736*exp(4*xi)-8208*exp(6*xi)-1080)*RootOf(3*_Z^2-3*_Z-1)+(-79584*B[1]^2+479376)*exp(10*xi)+(-121440*B[1]^2-777888)*exp(12*xi)+(50400*B[1]^2+492912)*exp(14*xi)+127440*exp(8*xi)*B[1]^2+39000*exp(16*xi)*B[1]^2-288*exp(2*xi)*B[1]^2+1056*exp(4*xi)*B[1]^2-16608*exp(6*xi)*B[1]^2+24*B[1]^2-138240*exp(8*xi)-70560*exp(16*xi)-4752*exp(2*xi)-8352*exp(4*xi)+27792*exp(6*xi)-288)*exp(2*xi-I*t*w+I*k*x)/(B[1]^3*(-3125*exp(20*xi)-25000*exp(18*xi)-76875*exp(16*xi)-108000*exp(14*xi)-55650*exp(12*xi)+12432*exp(10*xi)+11130*exp(8*xi)-4320*exp(6*xi)+615*exp(4*xi)-40*exp(2*xi)+1))

(10)

Download ode.mw

this is my first time something like that   coming up my equation after taking integral exponential coming up why?

g1.mw

the series is so complicated but have a strange pattern if you watch the index of parameter  they are not repeated 

 

 

Hi,

I am looking for a simpler way to find the equation of a parabola passing through 3 points. I see that using the Geometry package requires defining the parabola with 5 distinct points. Thank you for your guidance.QuestionParabole.mw

Do you think the result of String(0.016)  should be "0.016"  instead of ".16e-1" ?

Any reason why it gives the second form and not the first?  Now have to keep using sprintf to force formating as decimal point. Is this documented somewhere? quick search did not find anything do far.

Maple 2024.2 on windows.

s:="0.016";

"0.016"

z:= :-parse(s);

0.16e-1

String(z);

".16e-1"

sprintf("%0.3f",z);

"0.016"

 

 

Download string_of_decimal_number.mw

I'd like to know the details of the method Statistics:-Mean uses to numerically estimate the expectation of a random variable.

showstat seems of no use and neither seems to be LibraryTools[Browse]();

Here are two examples: the first one (1D) suggests  Statistics:-Mean could use some evalf/Int method, but the conclusion to draw from the second example (R2 --> R) is less clear.

How_does_Mean_proceed.mw

Thanks in advance

PS: I already asked a similar question months ago but didn't get any reply.
       Even answers such as “We don't know” or “We don't care” would suit me better than their absence.

restart;
with(plots);
with(geometry);
_EnvHorizontalName := 'x';
_EnvVerticalName := 'y';
HC := HorizontalCoord;
VC := VerticalCoord;
a := 11;
b := 7;
t := (3*Pi)/8;
c := sqrt(a^2 - b^2);
ellipse(e1, x^2/a^2 + y^2/b^2 = 1);
point(Oo, 0, 0);
point(A, a*cos(t), b*sin(t));
point(B, a*cos(t + 2/3*Pi), b*sin(t + 2/3*Pi));
point(C, a*cos(t + 4/3*Pi), b*sin(t + 4/3*Pi));
point(G, (A[1] + B[1] + C[1])/3, (A[2] + B[2] + C[2])/3);
eval(coordinates(G));
line(NorA, y - A[2] = a^2*A[2]*(x - A[1])/b^2, [x, y]);
line(NorB, y - B[2] = a^2*B[2]*(x - B[1])/b^2, [x, y]);
line(NorC, y - C[2] = a^2*C[2]*(x - C[1])/b^2, [x, y]);
lieu := a^2*x^2 + b^2*y^2 - c^4/4 = 0;
Lieu := implicitplot(lieu, x = -a .. a, y = -b .. b, color = green);
tx := textplot([[coordinates(A)[], "A"], [coordinates(B)[], "B"], [coordinates(C)[], "C"], [coordinates(Oo)[], "O"]], font = [times, bold, 16], align = [above, left]);
dr := draw([e1(color = blue), NorA(color = red), NorB(color = red), NorC(color = red), A(color = red, symbol = solidcircle, symbolsize = 12), B(color = red, symbol = solidcircle, symbolsize = 12), C(color = red, symbol = solidcircle, symbolsize = 12), Oo(color = red, symbol = solidcircle, symbolsize = 12)]);
display([dr, tx, Lieu], scaling = constrained, axes = normal, title = "Ellipse et normales ", titlefont = [HELVETICA, 14]);
      [1        1        1       1        1        1     ]
      [- A[1] + - B[1] + - C[1], - A[2] + - B[2] + - C[2]]
      [3        3        3       3        3        3     ]

                              NorA

                              NorB

                              NorC

Warning, data could not be converted to float Matrix
Warning, data could not be converted to float Matrix
Warning, data could not be converted to float Matrix

I have a simple nested for loop in a worksheet:

[> for x in [ -1, 1 ] do for y in [ -1, 1 ] do x*y end do end do

When I press Enter, nothing happens. What am I (as a novice) missing?

Scott

I am calling a proc inside a module, where this module had local module variable initialized to false;

This works fine when calling the proc normally.

When calling it from Grid:-Run() it says it is not seeing the module local variable at all.

Why does this happen and is there a workaround this? Here a simple worksheet showing this.

I hope there is a way to make the proc inside the module see the local module variables when using it from Grid node. Otherwise, this whole thing will not work for me, as I have few variables initialized at modules level in number of places.

restart;

A := module()
  local DEBUG_MSG::truefalse:=false;
  export work_proc:=proc(n::integer)
     print("My flag is ",DEBUG_MSG);
  end proc;
end module;

 

_m1690230418240

#THis works as expected
A:-work_proc(0)

"My flag is ", false

#this does not work. THe call works OK, but inside A:-workproc() it does not see module local variables.

Grid:-Set(A:-work_proc):
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

"My flag is ", DEBUG_MSG

 

 

Download grid_with_module_local_variable_feb_7_2025.mw


Update 2/7/2025

I tried to change the access to the module local variables from the module local procs, by adding explicit A:- to each variable name where A here is the module name.

This works for normal calls, but not when using Grid to call the proc.

This workaround fail, it gives error that module does not export `%1`

May be I have to redo all my code not to use local module variables. But in some places I have to do this, in order to detect if something has happened before or not. I use module local variables to store state the presist after call is completed.

It looks like Grid is not meant to be used for calling proc() that uses/lives inside modules. But this makes Grid not very useful then for large application.

restart;

A := module()
  local DEBUG_MSG::truefalse:=false;
  export work_proc:=proc(n::integer)
     print("My flag is ",A:-DEBUG_MSG);
  end proc;
end module;

 

_m1690230418240

#THis works as expected
A:-work_proc(0)

"My flag is ", false

#this does not work. THe call works OK, but inside A:-workproc() it does not see module local variables.

Grid:-Set(A:-work_proc):
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

Error, (in work_proc) module does not export `%1`

 

 

Download grid_with_module_local_variable_feb_7_2025_V2.mw

I also tried to make the module variable as export instead of local, and that also did not work. I am starting to run out of ideas what else to try...

restart;

A := module()
  export DEBUG_MSG::truefalse:=false;
  export work_proc:=proc(n::integer)
     print("My flag is ",A:-DEBUG_MSG);
  end proc;
end module;

 

_m1690230418208

#THis works as expected
A:-work_proc(0)

"My flag is ", false

#this does not work. THe call works OK, but inside A:-workproc() it does not see module local variables.

Grid:-Set(A:-work_proc):
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

Error, (in work_proc) `%1` does not evaluate to a module

 

 

Download grid_with_module_local_variable_feb_7_2025_V3.mw

Update

I think now that using module level local variables will not work with parallel processing anyway. It is like using global variable in parallel processing. Does not work without synchorization of access to this shared variable. 

So this means I have to change my code and not use any module local variables, and pass any information between functions via arguments only.

Using module local variable is more convenient, but I now realize this design is not good for parallel processing.  This means some code changes I have to do. 

On the positive side, I find using Grid can really speed things up. On some tests, up to 10 times faster. The larger the number of problems to process, the more speed up is gained.

So I think it is worth to rewrite my code and remove any use of local level module variables.

Update Feb 11, 2024

Found the solution!  With grid, one must pass in all variables that to be used on a node. So the above now becomes this

restart;

A := module()
  export DEBUG_MSG::truefalse:=false;
  export work_proc:=proc(n::integer)
     print("My flag is ",DEBUG_MSG);
  end proc;
end module;
 

_m1739526597408

#THis works as expected
A:-work_proc(0)

"My flag is ", false

#this does not work. Because we did not pass all variables to be used at node
Grid:-Set(0,A:-work_proc);
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

"My flag is ", DEBUG_MSG

#this now worksheetdir
Grid:-Set(0,A:-work_proc,'A:-DEBUG_MSG');
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

"My flag is ", false

 

 

Download grid_with_module_V4_feb_11_2025.mw

Actually the above could be done easier this way. By passing in the name of the whole module A, now nodes are able to see all module A local vaiables. Like this

restart;

A := module()
  export DEBUG_MSG::truefalse:=false;
  export work_proc:=proc(n::integer)
     print("My flag is ",A:-DEBUG_MSG);
  end proc;
end module;
 

_m1739526597408

#THis works as expected
A:-work_proc(0)

"My flag is ", false

#this does not work. Because we did not pass the name of the A
Grid:-Set(0,A:-work_proc);
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

Error, (in work_proc) `%1` does not evaluate to a module

#this now works, becuase we pass the module A
Grid:-Set(0,'A:-work_proc','A');
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

"My flag is ", false

 

 

Download grid_with_module_V6_feb_11_2025.mw

And to be able to access global variables in another module from inside the worker module, we must pass the other module name. Here is an example

restart;

B := module()
  export some_value::integer:=10;
end module;

A := module()
  export work_proc:=proc(n::integer)
     print("the value is ",B:-some_value);
  end proc;
end module;
 

_m1739526597280

_m1739594742272

#THis works as expected
A:-work_proc(0)

"the value is ", 10

#this does not work. Because we did not pass the name of the other module B in this case
Grid:-Set(0,A:-work_proc);
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

Error, (in work_proc) `%1` does not evaluate to a module

#this now works, becuase we pass the module name B to the node, so it sees it
Grid:-Set(0,'A:-work_proc','B');
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

"the value is ", 10

 

 

Download grid_with_module_V5_feb_11_2025.mw

The main point, if one is using modules with Grid, names of all modules to be accessed from a node need to be passed in in the Set call. Otherwise the node will not see them.

Update For the above method to work, all module variables to be accessed at node, have to export and not local.

Update 2/11/2025  4 AM

I think I'll give up on Grid. It is completely useless. I found out that I have to make all my modules proc's export now, even though they are only used in that one module, just in order to call them from same module when running in Grid node.

Otherwise, I have to list each proc name in each module and sub module in the Set command. This is ridiculous design. I have 5,000 proc's. I am not going to list them all each time.

Here is an example where a boo() is local proc to a module. But it is not called, because the proc is local.

restart;

A := module()

  export foo:=proc()
       print("entering A:-foo(), before calling local proc boo()");
       boo(); #this call will not work in Grid setting
  end proc;

  local boo:=proc()
     print("inside local A:-foo()");
  end proc;
end module:

A:-foo();

"entering A:-foo(), before calling local proc boo()"

"inside local A:-foo()"

Grid:-Set(0,A:-foo,'A');
Grid:-Run(0,A:-foo,[0]);
Grid:-Wait();

"entering A:-foo(), before calling local proc boo()"

 


 

Download calling_proc_in_module_V1.mw

To make it work, I have to make all my local proc export, and not only that, I have to call it using A:-boo() as well. So the above becomes


 

restart;

A := module()

  export foo:=proc()
       print("entering A:-foo(), before calling local proc boo()");
       A:-boo(); #this call will not work in Grid setting
  end proc;

  export boo:=proc()
     print("inside local A:-foo()");
  end proc;
end module:

A:-foo();

"entering A:-foo(), before calling local proc boo()"

"inside local A:-foo()"

Grid:-Set(0,A:-foo,'A');
Grid:-Run(0,A:-foo,[0]);
Grid:-Wait();

"entering A:-foo(), before calling local proc boo()"

"inside local A:-foo()"

 

 

Download calling_proc_in_module_V2.mw

I am not going to edit 100,000 lines of Maple code and change all the procs to export and make sure each call have the fully qualified module name attched to it even though it is local to the module.

I will not use Grid. It simply does not work with modules. Maplesoft should fix Grid so code works as with minimal changes.

How i can get FN from N-soliton to FN for M-lump by applying long wave method and using limit How we can get the series of lump , there is some example for nowing how lump work like f[2] for 1 lump and f[4] for two lump and f[6] for 3 lump i, and also in last of work i showed the series of for n soliton which we can change it to m-lump  but i don't know how work and a[12] will change to b[12] by applying long wave method

m-lump.mw

After my maple (2023) had a crash (froze) and i forc closed it all my files are opened and then changed. All equation fields are changed into "Maple Input" type equations. This completely ruins all my "document" type files where the appearance of the equations was important. More than that most of the changed equations give rise to errors. So the file is basically destroyed. The odd thing is that the file looks fine the first few seconds, but then it is changed completely.
I have uninstalled 2023, and installed Maple 2024 without success. 

Plz help!

Hello

I encountered a few problems. One is that in the first section, I wanted to use the definition above instead of f (s ) and g (s ), which means that when the variable changes under the integral sign, it should detect and replace it.

And the next is that in the Equality section, I should sort by p and set the coefficients to zero. And then, for example, solve for the zero power of p and get the value of f0 and use it in subsequent solutions. Can you help me?

restart;
EQUATIONS

equ1:=diff(f(t),t)-1-t-t^2-g(t)-int(f(s)+g(s),s=0..t)=0;

equ2:=diff(g(t),t)+1+t-f(t)+int(f(s)-g(s),s=0..t)=0;
 

diff(f(t), t)-1-t-t^2-g(t)-(int(f(s)+g(s), s = 0 .. t)) = 0

 

diff(g(t), t)+1+t-f(t)+int(f(s)-g(s), s = 0 .. t) = 0

(1)

f(t):=sum(f[i](t)*p^i,i=0..1);

f[0](t)+f[1](t)*p

(2)

g(t):=sum(g[i](t)*p^i,i=0..1);

g[0](t)+g[1](t)*p

(3)


HPMs

hpm1:=(1-p)*(diff(f(t),t)-1-t-t^2)+p*(-diff(f(t),t)+1+t+t^2-g(t)-int(f(s)+g(s),s=0..t))=0;

hpm2:=(1-p)*(diff(g(t),t)+1+t)+p*(diff(g(t),t)-1-t+f(t)-int(f(s)-g(s),s=0..t))=0;

(1-p)*(diff(f[0](t), t)+(diff(f[1](t), t))*p-1-t-t^2)+p*(-(diff(f[0](t), t))-(diff(f[1](t), t))*p+1+t+t^2-g[0](t)-g[1](t)*p-(int(f(s)+g(s), s = 0 .. t))) = 0

 

(1-p)*(diff(g[0](t), t)+(diff(g[1](t), t))*p+1+t)+p*(diff(g[0](t), t)+(diff(g[1](t), t))*p-1-t+f[0](t)+f[1](t)*p-(int(f(s)-g(s), s = 0 .. t))) = 0

(4)

``

Collect

A:=collect(hpm1,p);

(-2*(diff(f[1](t), t))-g[1](t))*p^2+(2*t^2-2*(diff(f[0](t), t))+diff(f[1](t), t)-g[0](t)-(int(f(s)+g(s), s = 0 .. t))+2*t+2)*p-t^2+diff(f[0](t), t)-t-1 = 0

(5)

EqualityNULL

for i from 0 to degree(A,p) do EQ[i]:=simplify(coeff(A,p,i)); end do;

Error, final value in for loop must be numeric or character

 
   

Download HPMsystem.mw

This is a error what all the time, shows up with the ai programs.

U:=U(xi);
Error, recursive assignment

U and U(xi)  seems to be related ? 

How do I get Maple 2023 to simplify/combine units in results?  For example,

Quantity(3.3897859*Unit(km^2)/Unit(m^2), 0.2) ;

should simplify to

Quantity(3389785.9, 0.2) ;

The full example where the problem occurs is given below.

alias(l = log10, l100 = log[100], pi[0] = Pi, r = sqrt, S_ellipsoid = ellipsoid) :
with(ScientificErrorAnalysis) :
with(Units) :

_km := Unit(km) :
_lm := Unit(lm) :
_lx := Unit(lx) :
_m := Unit(m) :
_rev := Unit(rev) :
alias(Q = Quantity) :
AddUnit(astronomical_unit, context = astronomy, default = true, conversion = 149597870700*m) :
pi := pi[0] :
S_spheroid := (a, b) -> S_ellipsoid(a, a, b) :

_AU := Unit(AU) :
a_Mars := Q(227939366., 1.)*_km : ##
a_Terra := Q(149598023., 1.)*_km : ##
alpha_Phobos := Q(9517.58, 0.01)*_km : ##
B_Phobos := Q(0.071, 0.012) :
d_x_Phobos := Q(25.90, 0.08)*_km :
d_y_Phobos := Q(22.60, 0.08)*_km :
d_z_Phobos := Q(18.32, 0.06)*_km :
E_I := _lx :
H_Mars := -Q(1.5, 0.1) : ##
L_Sol := Q(3.75E28, 0.01E28)*_lm : ##
m_E := -Q(14.18, 0.01) : ##
pi_Phobos := Q(9234.42, 0.01)*_km : ##
r_e_Mars := Q(3396.2, 0.1)*_km :
r_e_Terra := Q(6378137.0, 0.1)*_m : ##
S_sphere := r -> S_spheroid(r, r) :
theta_rev := _rev :

a_Phobos := (alpha_Phobos + pi_Phobos)/2 :
Delta_Sol := a_Terra - r_e_Terra : 
l_A := _AU :
m := E -> m_E - 5*l100(E/E_I) :
r_x_Phobos := d_x_Phobos/2 :
r_y_Phobos := d_y_Phobos/2 :
r_z_Phobos := d_z_Phobos/2 :
theta_rev2 := theta_rev/2 :

Delta_Mars := r(a_Mars^2 + Delta_Sol^2) : 
q := theta -> 2*((sin(theta)/pi) + (1 - (theta/theta_rev2))*cos(theta))/3 :
r_e_Phobos := (r_x_Phobos + r_y_Phobos)/2 :
S_Phobos := S_ellipsoid(r_x_Phobos, r_y_Phobos, r_z_Phobos) :

Delta_Phobos := Delta_Mars : 
L_Phobos := B_Phobos*L_Sol*S_Phobos/(2*S_sphere(a_Mars)) :
mu_Mars := 5*l(a_Mars*Delta_Mars/l_A^2) : 
rho_e_Mars := arcsin(r_e_Mars/Delta_Mars) :
theta_Mars := arccos((a_Mars^2 + Delta_Mars^2 - Delta_Sol^2)/(2*a_Mars*Delta_Mars)) :

E_Phobos := L_Phobos/S_sphere(Delta_Phobos) :
m_Mars := H_Mars + mu_Mars - 5*l100(q(theta_Mars)) :
rho_e_Phobos := arcsin(r_e_Phobos/Delta_Phobos) :
rho_o_Phobos := arctan(a_Phobos/Delta_Phobos) :

Deltarho_Phobos := rho_o_Phobos - rho_e_Mars - rho_e_Phobos :
m_Phobos := m(E_Phobos) :

Deltam_Phobos := m_Phobos - m_Mars :

Deltarhostar_Phobos := Deltarho_Phobos/Deltam_Phobos :

"Phobos apparent logarithmic brightness in astronomical magnitudes" = combine(m_Phobos, errors) ;
"Mars-Phobos apparent angular separation in arcseconds" = combine(convert(Deltarhostar_Phobos, units, arcsec), errors) ;

(*
https://www.nicolesharp.net/wiki/Solar_System_data_for_Maple
https://en.wikipedia.org/wiki/astronomical_unit
https://en.wikipedia.org/wiki/Star_Sol
https://en.wikipedia.org/wiki/Planet_Terra
https://en.wikipedia.org/wiki/Terran_radius
https://en.wikipedia.org/wiki/WGS84
https://en.wikipedia.org/wiki/Planet_Mars
https://en.wikipedia.org/wiki/Satellite_Phobos
https://en.wikipedia.org/wiki/Solar_System_by_size
https://en.wikipedia.org/wiki/illuminance
*)

 

Why does the code below work when I use a standalone "assume" statement but not "assuming"?

That is, why don't the first two attempts at calculating the limit use the assumption contained in those statements, ie why don't those statements return infinity and not a signum like the last attempt at the limit?

`ω__b` := proc (alpha) options operator, arrow; `ω__0`*sqrt(1+alpha+sqrt(alpha+alpha^2)) end proc

proc (alpha) options operator, arrow; omega__0*sqrt(1+alpha+sqrt(alpha+alpha^2)) end proc

(1)

`assuming`([limit(`ω__b`(alpha), alpha = infinity)], [`ω__0`::positive]) = signum(omega__0)*infinityNULL

`assuming`([limit(`ω__b`(alpha), alpha = infinity)], [`ω__0` > 0]) = signum(omega__0)*infinityNULL

NULL

assume(`ω__0` > 0)

limit(`ω__b`(alpha), alpha = infinity) = infinityNULL

NULL

Download assuming.mw

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