salim-barzani

1715 Reputation

9 Badges

1 years, 184 days

MaplePrimes Activity


These are questions asked by salim-barzani

i did for two of them base on the information but one of them is not make my odetest to be zero? where is problem

test.mw

i try to use explor option but is not work like i want and i try to use other way but they so classical i want to ask there is any other way for do plotting of this kind of function?
in second part even plot not working?

 

Bgraph1.mw

i did substitution but my result is so different from the author i think he just take the linear term of theta but i didn't do that so how take just linear term of that function and find unknwon , and how afeter replacing eq(12) inside eq(11) we can remove thus exponential and find w? also i think author did a mistake which the equation 12 is theta(x,t) not Q(x,t)

restart

with(PDEtools)

undeclare(prime, quiet)

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

pde := diff(u(x, t), `$`(t, 2))-s^2*(diff(u(x, t), `$`(x, 2)))+(2*I)*(diff(u(x, t)*U^2, t))-(2*I)*alpha*s*(diff(u(x, t)*U^2, t))+I*(diff(u(x, t), `$`(x, 2), t))-I*beta*s*(diff(u(x, t), `$`(x, 3)))

diff(diff(u(x, t), t), t)-s^2*(diff(diff(u(x, t), x), x))+(2*I)*(diff(u(x, t), t))*U^2-(2*I)*alpha*s*(diff(u(x, t), t))*U^2+I*(diff(diff(diff(u(x, t), t), x), x))-I*beta*s*(diff(diff(diff(u(x, t), x), x), x))

(1)

T := u(x, t) = (sqrt(Q)+theta(x, t))*exp(I*(Q^2*epsilon*gamma+Q*q)*t); T1 := U = sqrt(Q)+theta(x, t)

u(x, t) = (Q^(1/2)+theta(x, t))*exp(I*(Q^2*epsilon*gamma+Q*q)*t)

 

U = Q^(1/2)+theta(x, t)

(2)

P := collect(eval(subs({T, T1}, pde)), exp)/exp(I*(Q^2*gamma*`ε`+Q*q)*t)

diff(diff(theta(x, t), t), t)+(2*I)*(diff(theta(x, t), t))*(Q^2*epsilon*gamma+Q*q)-(Q^(1/2)+theta(x, t))*(Q^2*epsilon*gamma+Q*q)^2-s^2*(diff(diff(theta(x, t), x), x))+(2*I)*(diff(theta(x, t), t)+I*(Q^(1/2)+theta(x, t))*(Q^2*epsilon*gamma+Q*q))*(Q^(1/2)+theta(x, t))^2-(2*I)*alpha*s*(diff(theta(x, t), t)+I*(Q^(1/2)+theta(x, t))*(Q^2*epsilon*gamma+Q*q))*(Q^(1/2)+theta(x, t))^2+I*(diff(diff(diff(theta(x, t), t), x), x)+I*(diff(diff(theta(x, t), x), x))*(Q^2*epsilon*gamma+Q*q))-I*beta*s*(diff(diff(diff(theta(x, t), x), x), x))

(3)

 

TT := Q = alpha[1]*exp(I*(k*x-t*w))+alpha[2]*exp(-I*(k*x-t*w))

Q = alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w))

(4)

S := eval(subs(TT, P))

diff(diff(theta(x, t), t), t)+(2*I)*(diff(theta(x, t), t))*(gamma*epsilon*(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^2+(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))*q)-((alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^(1/2)+theta(x, t))*(gamma*epsilon*(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^2+(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))*q)^2-s^2*(diff(diff(theta(x, t), x), x))+(2*I)*(diff(theta(x, t), t)+I*((alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^(1/2)+theta(x, t))*(gamma*epsilon*(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^2+(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))*q))*((alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^(1/2)+theta(x, t))^2-(2*I)*alpha*s*(diff(theta(x, t), t)+I*((alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^(1/2)+theta(x, t))*(gamma*epsilon*(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^2+(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))*q))*((alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^(1/2)+theta(x, t))^2+I*(diff(diff(diff(theta(x, t), t), x), x)+I*(diff(diff(theta(x, t), x), x))*(gamma*epsilon*(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^2+(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))*q))-I*beta*s*(diff(diff(diff(theta(x, t), x), x), x))

(5)

Download steps.mw

or for this equation 

steps-2.mw

I test a lot of them but some of them make a problem i  don't know i am do it in wrong way or the author did wrong i need verifying thus solution of odes specially in case 4 when we have not equal sign how use that?

and case 5 is Weierstrass elliptic function which i don't know how set up and use i think is a on kinf of odes but why they use that sign for this function?

ode-17.mw

Already  by Help of my favorite Dr.David he did find the thus three step for non schrodinger equation but in here i got some issue of coding which is so different from before, is about transformation of pdes to two parts od real and imaginary part and then substitution our function the functions is clear but combining them and findinf leading exponent and resonance point and finding function in step 3 in different and jsut the function is different with eperate the real and imaginary part for finding step one ...

note:=q=u*exp(#) then |q|=u

schrodinger-test.mw

paper-1

paper-2

4 5 6 7 8 9 10 Last Page 6 of 39