salim-barzani

1560 Reputation

9 Badges

1 years, 21 days

MaplePrimes Activity


These are replies submitted by salim-barzani

@dharr sorry for that this is same as paper which i have to put in adomian polynomial 

restart

with(inttrans)

with(PDEtools)

with(DEtools)

with(Physics)

declare(u(x, t), quiet); declare(v(x, t), quiet)

undeclare(prime)

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

(1)

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

u(x, t)+I*(diff(diff(u(x, t), x), x))+(2*I)*((diff(u(x, t), x))*conjugate(u(x, t))+u(x, t)*(diff(conjugate(u(x, t)), x)))*u(x, t)+I*u(x, t)^3*conjugate(u(x, t))^2

(2)

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

u(x, t)+I*(diff(diff(u(x, t), x), x)), (2*I)*u(x, t)*(diff(u(x, t), x))*conjugate(u(x, t))+(2*I)*u(x, t)^2*(diff(conjugate(u(x, t)), x))+I*u(x, t)^3*conjugate(u(x, t))^2

(3)

NULL

NULL

NULL

NULL

NULL

B[0] := -I*u[0]^3*conjugate(u[0])^2

-I*u[0]^3*conjugate(u[0])^2

(4)

B1[0] := -(2*I)*u[0]^2*(diff(u[0](x), x))

-(2*I)*u[0]^2*(diff(u[0](x), x))

(5)

T[0] := -(2*I)*u[0]*(diff(u[0](x), x))*conjugate(u[0])

-(2*I)*u[0]*(diff(u[0](x), x))*conjugate(u[0])

(6)

for n to 4 do B[n] := expand(-(2*I)*simplify(diff((sum(u[k]*lambda^k, k = 0 .. n))^3*(sum(v[k]*lambda^k, k = 0 .. n))^2/factorial(n), [`$`(lambda, n)]))); B1[n] := expand(-(2*I)*simplify(diff((sum(u[k]*lambda^k, k = 0 .. n))^2*(sum(v[k, x]*lambda^k, k = 0 .. n))/factorial(n), [`$`(lambda, n)]))); T[n] := expand(-(2*I)*simplify(diff((sum(u[k]*lambda^k, k = 0 .. n))*(sum(v[k]*lambda^k, k = 0 .. n))*(sum(z[k, x]*lambda^k, k = 0 .. n))/factorial(n), [`$`(lambda, n)]))) end do

lambda := 0; for n from 0 to 4 do P[n] := B[n]; Q[n] := B1[n]; R[n] := T[n] end do

0

 

-I*u[0]^3*v[0]^2

 

-(2*I)*u[0]^2*v[0, x]

 

-(2*I)*u[0]*z[0, x]*v[0]

 

-(6*I)*u[0]^2*u[1]*v[0]^2-(4*I)*u[0]^3*v[0]*v[1]

 

-(2*I)*u[0]^2*v[1, x]-(4*I)*u[0]*u[1]*v[0, x]

 

-(2*I)*u[0]*v[0]*z[1, x]-(2*I)*u[0]*v[1]*z[0, x]-(2*I)*u[1]*v[0]*z[0, x]

 

-(2*I)*u[0]^3*v[1]^2-(4*I)*u[0]^3*v[0]*v[2]-(6*I)*u[0]^2*u[2]*v[0]^2-(6*I)*u[0]*u[1]^2*v[0]^2-(12*I)*u[0]^2*u[1]*v[0]*v[1]

 

-(2*I)*u[0]^2*v[2, x]-(4*I)*u[0]*u[1]*v[1, x]-(4*I)*u[0]*u[2]*v[0, x]-(2*I)*u[1]^2*v[0, x]

 

-(2*I)*u[1]*v[0]*z[1, x]-(2*I)*u[1]*v[1]*z[0, x]-(2*I)*u[0]*v[1]*z[1, x]-(2*I)*u[0]*v[2]*z[0, x]-(2*I)*u[2]*v[0]*z[0, x]-(2*I)*u[0]*v[0]*z[2, x]

 

-(12*I)*u[0]^2*u[1]*v[0]*v[2]-(12*I)*u[0]^2*u[2]*v[0]*v[1]-(12*I)*u[0]*u[1]^2*v[0]*v[1]-(12*I)*u[0]*u[1]*u[2]*v[0]^2-(2*I)*u[1]^3*v[0]^2-(4*I)*u[0]^3*v[0]*v[3]-(4*I)*u[0]^3*v[1]*v[2]-(6*I)*u[0]^2*u[1]*v[1]^2-(6*I)*u[0]^2*u[3]*v[0]^2

 

-(2*I)*u[0]^2*v[3, x]-(2*I)*u[1]^2*v[1, x]-(4*I)*u[0]*u[1]*v[2, x]-(4*I)*u[0]*u[2]*v[1, x]-(4*I)*u[0]*u[3]*v[0, x]-(4*I)*u[1]*u[2]*v[0, x]

 

-(2*I)*u[0]*v[0]*z[3, x]-(2*I)*u[0]*v[1]*z[2, x]-(2*I)*u[0]*v[2]*z[1, x]-(2*I)*u[0]*v[3]*z[0, x]-(2*I)*u[1]*v[0]*z[2, x]-(2*I)*u[1]*v[1]*z[1, x]-(2*I)*u[1]*v[2]*z[0, x]-(2*I)*u[2]*v[0]*z[1, x]-(2*I)*u[2]*v[1]*z[0, x]-(2*I)*u[3]*v[0]*z[0, x]

 

-(4*I)*u[0]^3*v[1]*v[3]-(12*I)*u[0]*u[1]^2*v[0]*v[2]-(6*I)*u[0]*u[1]^2*v[1]^2-(24*I)*u[0]*u[1]*u[2]*v[0]*v[1]-(12*I)*u[0]^2*u[2]*v[0]*v[2]-(12*I)*u[0]^2*u[1]*v[0]*v[3]-(6*I)*u[0]*u[2]^2*v[0]^2-(4*I)*u[1]^3*v[0]*v[1]-(2*I)*u[0]^3*v[2]^2-(6*I)*u[0]^2*u[4]*v[0]^2-(12*I)*u[0]*u[1]*u[3]*v[0]^2-(12*I)*u[0]^2*u[3]*v[0]*v[1]-(6*I)*u[1]^2*u[2]*v[0]^2-(6*I)*u[0]^2*u[2]*v[1]^2-(12*I)*u[0]^2*u[1]*v[1]*v[2]-(4*I)*u[0]^3*v[0]*v[4]

 

-(2*I)*u[0]^2*v[4, x]-(2*I)*u[1]^2*v[2, x]-(2*I)*u[2]^2*v[0, x]-(4*I)*u[0]*u[1]*v[3, x]-(4*I)*u[0]*u[2]*v[2, x]-(4*I)*u[0]*u[3]*v[1, x]-(4*I)*u[0]*u[4]*v[0, x]-(4*I)*u[1]*u[2]*v[1, x]-(4*I)*u[1]*u[3]*v[0, x]

 

-(2*I)*u[0]*v[0]*z[4, x]-(2*I)*u[0]*v[1]*z[3, x]-(2*I)*u[0]*v[2]*z[2, x]-(2*I)*u[0]*v[4]*z[0, x]-(2*I)*u[1]*v[2]*z[1, x]-(2*I)*u[0]*v[3]*z[1, x]-(2*I)*u[1]*v[0]*z[3, x]-(2*I)*u[1]*v[1]*z[2, x]-(2*I)*u[1]*v[3]*z[0, x]-(2*I)*u[2]*v[0]*z[2, x]-(2*I)*u[2]*v[1]*z[1, x]-(2*I)*u[2]*v[2]*z[0, x]-(2*I)*u[4]*v[0]*z[0, x]-(2*I)*u[3]*v[0]*z[1, x]-(2*I)*u[3]*v[1]*z[0, x]

(7)
 

NULL

Download Dr.D.mw

@dharr this is just one term of nonlinear i have to  put this function in adomian polynomial for calculatin P[0],p[1],P[2] and also for Q and R then add them to get A[0],A[1] which A[n] is all term of nonlinearity, we don't want calculate  one by one of them we calculate all term together one time , 
how you define in this adomian polynomian 

if you watch the paper you will see the directly  i did the same 

Download Dr.D.mw

@dharr i send you in mail i hope you got the mean, if you undrestand the problem it will be solved automatically

@dharr is just about aranging function with letter all thus function i used are u and when we use that laplace we have just one function which is u[0] when we substitute we get u[1] and for u[2] we use u[1] for get u[3] we have to use u[2] but in adomian polynomial i can't use conjugate(u) or diff(u,x) is not work like that becuase of that i use v for conjugate and z for derivative of u . 

@janhardo  the code are too comlicate , if you watch the accer answer is so short and  mine is true too but i add equal sign sometime code not work with = sign .

@dharr  let me give you the more detail i know each steps but coding for me is prblem and i did my trail, thus B[0] and B1[0] and T[o] they non linear term in my pdf which i changed to P[i] and Q[i] and R[i] as paper did that  and when i add three term of non linear part it give  A[0]=P[0]+Q[0]+R[0] and so on for other A[i], 
the second part i have to define thus nonlinear term by adomian polynomial i have two definition i try to apply more easier but i fail so at end i did somehow which make my term true  but  in function susbtitution i got problem  becuase in definition of adomian polynomial i change each cojugate(u(x,t)) to v(x,t)  and for diff(u(x,t),x) i try use z[i,x] which this i is change and x is derivative and for diff(conjugate(u(x,t)),x) i use v[i,x] but all of thus z and v are u function, really is make a lot hard for me to arange them , and my invistagation is stuck becuase of that there is 10 equation i am stoped here for this step 

i will update the figure here also adomian definition for nonlinear term  for more undrestanding i hope you find out really i needed that

if you watch my step are not wrong and outcome are completly true  
if you need other information please just mention it 

note: |u|^2=u*conjugate(u)  by complex property

and this is adomian polynomial  for non linear term

@janhardo  don't be angry i didn't mean that, is perfect for some one who know how work with it and i am know you are expert with it and 100 year i work on that ode which you give me the idea i could not found that idea is depend on the question sometime give us answer sometime not 

@acer Ai is trash in 100 time just one time work , when i am stuck i will post here and my work is rare i cant find it any where i have to do step by step,  many thanks for you

@acer Dr. david and maccdara write thus code each one write one of them i try to used but i am fail , i have a lot question they easy but i  don't know how apply code when i apply not work 

@janhardo  i have to change my topic, i will go to another topic and i will ask a lot question in future

@janhardo  i put condition on it but didn't work it and i don't know how he reach this goal.

@janhardo  but again i didn't get result

ode1 := diff(G(xi), xi)^2 = A^2 + 2*A*B*G(xi) + 2*A*C*G(xi)^2 + 2*B*C*G(xi)^3 + C^2*G(xi)^4;
                              /         
ode1 := Typesetting:-mcomplete|G[ξ], 
                              \         

                    /[ d        ]\\      2
  Typesetting:-_Hold|[---- G(xi)]||^2 = A 
                    \[ dxi      ]//       

   + 2 A B Typesetting:-mcomplete(G, Typesetting:-_Hold([G(xi)]))

   + 2 A C 

                                                        2
  Typesetting:-mcomplete(G, Typesetting:-_Hold([G(xi)])) 

   + 2 B C 

                                                        3
  Typesetting:-mcomplete(G, Typesetting:-_Hold([G(xi)])) 

      2                                                       4
   + C  Typesetting:-mcomplete(G, Typesetting:-_Hold([G(xi)])) 


W := B = -2*A*C;
                        W := B = -2 A C

ode2 := subs(W, ode1);
                              /         
ode2 := Typesetting:-mcomplete|G[ξ], 
                              \         

                    /[ d        ]\\      2
  Typesetting:-_Hold|[---- G(xi)]||^2 = A 
                    \[ dxi      ]//       

        2                                                         
   - 4 A  C Typesetting:-mcomplete(G, Typesetting:-_Hold([G(xi)]))

   + 2 A C 

                                                        2
  Typesetting:-mcomplete(G, Typesetting:-_Hold([G(xi)])) 

          2 
   - 4 A C  

                                                        3
  Typesetting:-mcomplete(G, Typesetting:-_Hold([G(xi)])) 

      2                                                       4
   + C  Typesetting:-mcomplete(G, Typesetting:-_Hold([G(xi)])) 


;
G25 := G(xi) = -1/(2*C)*(sqrt(-2*A*C) + sqrt(-6*A*C)*tanh(1/2*sqrt(-6*A*C)*xi));
G25 := Typesetting:-mcomplete(G, Typesetting:-_Hold([G(xi)])) = - 

          (1/2)           (1/2)     /1         (1/2)   \
  (-2 A C)      + (-6 A C)      tanh|- (-6 A C)      xi|
                                    \2                 /
  ------------------------------------------------------
                           2 C                          


(simplify(odetest(G25, ode2)) assuming (A < 0, 0 < C));

@janhardo  we have to change the ode when we change the ode we have to replace A,B,C changed to thus in eq4 and 5 but how satisfy i don't have idia

@janhardo  after equation 24 the reds one is not satisfy the equation i have to use equation 4-5 condition but i don't know how they use and satisfy the ode 

@janhardo  in equation 4 to 5 there is some substitution maybe they use that can you arrange the condition

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