salim-barzani

265 Reputation

3 Badges

0 years, 222 days

MaplePrimes Activity


These are questions asked by salim-barzani

how fixed this for ode test

restart

with(PDEtools)

with(Physics)

undeclare(prime)

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

(1)

``

pde := -I*(diff(U(xi), xi))*gamma*k*mu+I*gamma*(diff(U(xi), xi))*sigma*w+(diff(diff(U(xi), xi), xi))*gamma*k*w+U(xi)*gamma*mu*sigma+(2*I)*(diff(U(xi), xi))*k*sigma+2*alpha*U(xi)^3+(diff(diff(U(xi), xi), xi))*k^2-I*(diff(U(xi), xi))*w-U(xi)*sigma^2-U(xi)*mu

-I*gamma*(diff(U(xi), xi))*k*mu+I*gamma*(diff(U(xi), xi))*sigma*w+gamma*(diff(diff(U(xi), xi), xi))*k*w+gamma*U(xi)*mu*sigma+(2*I)*(diff(U(xi), xi))*k*sigma+2*alpha*U(xi)^3+(diff(diff(U(xi), xi), xi))*k^2-I*(diff(U(xi), xi))*w-U(xi)*sigma^2-U(xi)*mu

(2)

case1 := [mu = -(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1), A[0] = 0, A[1] = -RootOf(_Z^2*alpha+gamma*k*w+k^2), B[1] = RootOf(_Z^2*alpha+gamma*k*w+k^2), w = (gamma*k*mu-2*k*sigma)/(gamma*sigma-1)]

[mu = -(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1), A[0] = 0, A[1] = -RootOf(_Z^2*alpha+gamma*k*w+k^2), B[1] = RootOf(_Z^2*alpha+gamma*k*w+k^2), w = (gamma*k*mu-2*k*sigma)/(gamma*sigma-1)]

(3)

G1 := U(xi) = 2*RootOf(_Z^2*alpha+gamma*k*w+k^2)/sinh(2*xi)

U(xi) = 2*RootOf(_Z^2*alpha+gamma*k*w+k^2)/sinh(2*xi)

(4)

pde1 := subs(case1, pde)

I*gamma*(diff(U(xi), xi))*k*(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1)+I*gamma*(diff(U(xi), xi))*sigma*(gamma*k*mu-2*k*sigma)/(gamma*sigma-1)+gamma*(diff(diff(U(xi), xi), xi))*k*(gamma*k*mu-2*k*sigma)/(gamma*sigma-1)-gamma*U(xi)*(4*gamma*k*w+4*k^2-sigma^2)*sigma/(gamma*sigma-1)+(2*I)*(diff(U(xi), xi))*k*sigma+2*alpha*U(xi)^3+(diff(diff(U(xi), xi), xi))*k^2-I*(diff(U(xi), xi))*(gamma*k*mu-2*k*sigma)/(gamma*sigma-1)-U(xi)*sigma^2+U(xi)*(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1)

(5)

pde2 := subs(case1, pde1)

I*gamma*(diff(U(xi), xi))*k*(4*gamma*(gamma*k*mu-2*k*sigma)*k/(gamma*sigma-1)+4*k^2-sigma^2)/(gamma*sigma-1)+I*gamma*(diff(U(xi), xi))*sigma*(-gamma*k*(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1)-2*k*sigma)/(gamma*sigma-1)+gamma*(diff(diff(U(xi), xi), xi))*k*(-gamma*k*(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1)-2*k*sigma)/(gamma*sigma-1)-gamma*U(xi)*(4*gamma*(gamma*k*mu-2*k*sigma)*k/(gamma*sigma-1)+4*k^2-sigma^2)*sigma/(gamma*sigma-1)+(2*I)*(diff(U(xi), xi))*k*sigma+2*alpha*U(xi)^3+(diff(diff(U(xi), xi), xi))*k^2-I*(diff(U(xi), xi))*(-gamma*k*(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1)-2*k*sigma)/(gamma*sigma-1)-U(xi)*sigma^2+U(xi)*(4*gamma*(gamma*k*mu-2*k*sigma)*k/(gamma*sigma-1)+4*k^2-sigma^2)/(gamma*sigma-1)

(6)

odetest(G1, pde2)

 

NULL

Download test_sol_for_PDE1.mw

I don't know how make my graph be beter for real part and imaginary part and abs part which part how work with parameter can any one explain on this example?

G.mw

i did two case of this equation and odetest is worked good but in this case the odetest is not worked well anyone can determine what is mistake ?

F_P_Correct_case_three.mw

I get my on results but the results are not the same please help me if i did any mistake in my code

 

symmetry_PDESYS_3_time_fraction[1].mw

how we can make pretty of changing equation or case in maple to latex but be more readble and be simplify and the scientific referee don't reject the writting for example in this picture when i do convert i get something like that

 

if you watch the result 1 it is better denomenator which is 2 come out will be better and then times multiply by radical so how we do something like that?

latex2.mw

First 12 13 14 15 16 17 18 Page 14 of 20