Maple 18 Questions and Posts

These are Posts and Questions associated with the product, Maple 18

This is a simple encryption method to hide text messages

Mentioned in Arabic manuscrips with more than hundreds years old ...

PRINCIPLE :

Just the place of letters in the sentence rearranged as described below :

For example "ABCDE" we pick up the First letter "A" from the left and write it as the last letter in the Right "......A"

but this time we pick up the letter "E" as the last letter from Right and place it at the Left Side of the previous one  ".....EA"

and this cycle continue until for rest letters ... "CDBEA" .

by this way the text become hard to discover !

It is Amazing that for decoding this message you should repeat the same rearrangment algorithm several times until the readable text appears as the first "ABCDE"

EXample :

"AlbertEinstein"

"iEntsrteebilnA"

"eterbsitlnnEAi"

 "tilsnbnrEeAtie"

"rnEbenAstliiet"

"sAtnleibiEentr"

 "biieElennttArs"

"nenltEteAirisb"

"etAEitrlinsebn"

"lritnisEeAbtne"

"EseiAnbttinrel"

"tbtniAnireeslE"

"inrAeienstlbEt"

"nesitelAbrEnti"

"AlbertEinstein"

the same text appeared after 14 step cycle


 

Arabic Cipher

 

ArabicCipher := proc (x) options operator, arrow; StringTools[Permute](x, [seq(1+iquo(StringTools[Length](x), 2)+((1/2)*i+(1/2)*irem(i, 2))*(-1)^(i+irem(StringTools[Length](x), 2)), i = 0 .. StringTools[Length](x)-1)]) end proc

proc (x) options operator, arrow; StringTools[Permute](x, [seq(1+iquo(StringTools[Length](x), 2)+((1/2)*i+(1/2)*irem(i, 2))*(-1)^(i+irem(StringTools[Length](x), 2)), i = 0 .. StringTools[Length](x)-1)]) end proc

(1.1)

seq((ArabicCipher@@i)("AlbertEinstein"), i = 1 .. 14)

"iEntsrteebilnA", "eterbsitlnnEAi", "tilsnbnrEeAtie", "rnEbenAstliiet", "sAtnleibiEentr", "biieElennttArs", "nenltEteAirisb", "etAEitrlinsebn", "lritnisEeAbtne", "EseiAnbttinrel", "tbtniAnireeslE", "inrAeienstlbEt", "nesitelAbrEnti", "AlbertEinstein"

(1.2)

NULL

seq((ArabicCipher@@i)("FereydoonShekofte"), i = 1 .. 12)

"nSohoedkyoefrteeF", "yokedferotheoeSFn", "otrheefodeeSkFony", "deoefSekeFhorntyo", "eFkheoSrfnetoyeod", "fnreStooeyhekoFde", "eyohoetkSoeFrdnef", "SoketFerodhnoeyfe", "odrhenFoteeykfoeS", "teoeFynkefhoredSo", "efkhnoyrFeedoSeot", "FereydoonShekofte"

(1.3)

``


 

Download Arabic_Cipher.mw

 

 

restart;
T := mu+lambda*H(xi)+(v-1)*H(xi)^2;
                                                2
               mu + lambda H(xi) + (v - 1) H(xi) 
u[0] := a[0]+a[1]*(d+H(xi))+a[2]/(d+H(xi))+a[3]*(d+H(xi))^2+a[4]/(d+H(xi))^2;
                                a[2]                      2
    a[0] + a[1] (d + H(xi)) + --------- + a[3] (d + H(xi)) 
                              d + H(xi)                    

             a[4]    
       + ------------
                    2
         (d + H(xi)) 
diff(u[0], xi);
                            / d        \
                       a[2] |---- H(xi)|
        / d        \        \ dxi      /
   a[1] |---- H(xi)| - -----------------
        \ dxi      /                2   
                         (d + H(xi))    

                                                 / d        \
                                          2 a[4] |---- H(xi)|
                           / d        \          \ dxi      /
      + 2 a[3] (d + H(xi)) |---- H(xi)| - -------------------
                           \ dxi      /                 3    
                                             (d + H(xi))     
collect(%, diff(H(xi), xi));
/           a[2]                               2 a[4]   \ / d       
|a[1] - ------------ + 2 a[3] (d + H(xi)) - ------------| |---- H(xi
|                  2                                   3| \ dxi     
\       (d + H(xi))                         (d + H(xi)) /           

   \
  )|
   /
d[1] := (a[1]-a[2]/(d+H(xi))^2+2*a[3]*(d+H(xi))-2*a[4]/(d+H(xi))^3)*T;
 /           a[2]                               2 a[4]   \ /  
 |a[1] - ------------ + 2 a[3] (d + H(xi)) - ------------| \mu
 |                  2                                   3|    
 \       (d + H(xi))                         (d + H(xi)) /    

                                  2\
    + lambda H(xi) + (v - 1) H(xi) /
diff(d[1], xi);
/       / d        \                                / d        \\ 
|2 a[2] |---- H(xi)|                         6 a[4] |---- H(xi)|| 
|       \ dxi      /          / d        \          \ dxi      /| 
|------------------- + 2 a[3] |---- H(xi)| + -------------------| 
|              3              \ dxi      /                 4    | 
\   (d + H(xi))                                 (d + H(xi))     / 

  /                                 2\   /           a[2]    
  \mu + lambda H(xi) + (v - 1) H(xi) / + |a[1] - ------------
                                         |                  2
                                         \       (d + H(xi)) 

                             2 a[4]   \ /       / d        \
   + 2 a[3] (d + H(xi)) - ------------| |lambda |---- H(xi)|
                                     3| \       \ dxi      /
                          (d + H(xi)) /                     

                     / d        \\
   + 2 (v - 1) H(xi) |---- H(xi)||
                     \ dxi      //
collect(%, diff(H(xi), xi));
//   2 a[2]                  6 a[4]   \ /                 
||------------ + 2 a[3] + ------------| \mu + lambda H(xi)
||           3                       4|                   
\\(d + H(xi))             (d + H(xi)) /                   

                  2\   /           a[2]                         
   + (v - 1) H(xi) / + |a[1] - ------------ + 2 a[3] (d + H(xi))
                       |                  2                     
                       \       (d + H(xi))                      

        2 a[4]   \                           \ / d        \
   - ------------| (lambda + 2 (v - 1) H(xi))| |---- H(xi)|
                3|                           | \ dxi      /
     (d + H(xi)) /                           /             
d[2] := ((2*a[2]/(d+H(xi))^3+2*a[3]+6*a[4]/(d+H(xi))^4)*(mu+lambda*H(xi)+(v-1)*H(xi)^2)+(a[1]-a[2]/(d+H(xi))^2+2*a[3]*(d+H(xi))-2*a[4]/(d+H(xi))^3)*(lambda+(2*(v-1))*H(xi)))*T;
//   2 a[2]                  6 a[4]   \ /                 
||------------ + 2 a[3] + ------------| \mu + lambda H(xi)
||           3                       4|                   
\\(d + H(xi))             (d + H(xi)) /                   

                  2\   /           a[2]                         
   + (v - 1) H(xi) / + |a[1] - ------------ + 2 a[3] (d + H(xi))
                       |                  2                     
                       \       (d + H(xi))                      

        2 a[4]   \                           \ /                 
   - ------------| (lambda + 2 (v - 1) H(xi))| \mu + lambda H(xi)
                3|                           |                   
     (d + H(xi)) /                           /                   

                  2\
   + (v - 1) H(xi) /

eq := (2*k*k)*w*beta*d[2]-(2*alpha*k*k)*d[1]-2*w*u[0]+k*u[0]*u[0];
   2        //   2 a[2]                  6 a[4]   \ /  
2 k  w beta ||------------ + 2 a[3] + ------------| \mu
            ||           3                       4|    
            \\(d + H(xi))             (d + H(xi)) /    

                                 2\   /           a[2]    
   + lambda H(xi) + (v - 1) H(xi) / + |a[1] - ------------
                                      |                  2
                                      \       (d + H(xi)) 

                             2 a[4]   \                          
   + 2 a[3] (d + H(xi)) - ------------| (lambda + 2 (v - 1) H(xi)
                                     3|                          
                          (d + H(xi)) /                          

   \ /                                 2\            2 /    
  )| \mu + lambda H(xi) + (v - 1) H(xi) / - 2 alpha k  |a[1]
   |                                                   |    
   /                                                   \    

         a[2]                               2 a[4]   \ /  
   - ------------ + 2 a[3] (d + H(xi)) - ------------| \mu
                2                                   3|    
     (d + H(xi))                         (d + H(xi)) /    

                                 2\       /    
   + lambda H(xi) + (v - 1) H(xi) / - 2 w |a[0]
                                          |    
                                          \    

                          a[2]                      2
   + a[1] (d + H(xi)) + --------- + a[3] (d + H(xi)) 
                        d + H(xi)                    

         a[4]    \     /                            a[2]   
   + ------------| + k |a[0] + a[1] (d + H(xi)) + ---------
                2|     |                          d + H(xi)
     (d + H(xi)) /     \                                   

                     2       a[4]    \  
   + a[3] (d + H(xi))  + ------------|^2
                                    2|  
                         (d + H(xi)) /  
value(%);
   2        //   2 a[2]                  6 a[4]   \ /  
2 k  w beta ||------------ + 2 a[3] + ------------| \mu
            ||           3                       4|    
            \\(d + H(xi))             (d + H(xi)) /    

                                 2\   /           a[2]    
   + lambda H(xi) + (v - 1) H(xi) / + |a[1] - ------------
                                      |                  2
                                      \       (d + H(xi)) 

                             2 a[4]   \                          
   + 2 a[3] (d + H(xi)) - ------------| (lambda + 2 (v - 1) H(xi)
                                     3|                          
                          (d + H(xi)) /                          

   \ /                                 2\            2 /    
  )| \mu + lambda H(xi) + (v - 1) H(xi) / - 2 alpha k  |a[1]
   |                                                   |    
   /                                                   \    

         a[2]                               2 a[4]   \ /  
   - ------------ + 2 a[3] (d + H(xi)) - ------------| \mu
                2                                   3|    
     (d + H(xi))                         (d + H(xi)) /    

                                 2\       /    
   + lambda H(xi) + (v - 1) H(xi) / - 2 w |a[0]
                                          |    
                                          \    

                          a[2]                      2
   + a[1] (d + H(xi)) + --------- + a[3] (d + H(xi)) 
                        d + H(xi)                    

         a[4]    \     /                            a[2]   
   + ------------| + k |a[0] + a[1] (d + H(xi)) + ---------
                2|     |                          d + H(xi)
     (d + H(xi)) /     \                                   

                     2       a[4]    \  
   + a[3] (d + H(xi))  + ------------|^2
                                    2|  
                         (d + H(xi)) /  
expr := simplify(%);
Error, (in simplify) too many levels of recursion
temp := algsubs(d+H(xi) = freeze(d+H(xi)), numer(expr));
                              expr
thaw(collect(temp, freeze(d+H(xi)))/denom(expr));
                              expr
collect(%, H(xi));
 

Hello everyone,

I'm struggling to solve inequalities with conditions.

I have this inequality with 4 variables, which I have some conditions. However, I can't implement this conditions to the inequality and solve using the 'solve' command.

Can anybody help me?

inequality.mw

restart;
lambda := 1;
                               1
mu := 1;
                               1
v := 2;
                               2
r := lambda*(v-1);
                               1
g := mu*(v-1);
                               1
a[2] := 0;
                               0

omega := -(1/2)*alpha*l[1]*l[2]*lambda^2+2*alpha*l[1]*l[2]*mu*v-2*alpha*l[1]*l[2]*mu-alpha*h[1]*h[2];
              3                                  
              - alpha l[1] l[2] - alpha h[1] h[2]
              2                                  
a[0] := -(1/2)*(2*d*v-2*d-lambda)*alpha*l[1]*l[2]/(h[1]*beta*(sqrt(h[1]*beta*alpha*l[1]*l[2])/(h[1]*beta)));
                    (2 d - 1) alpha l[1] l[2]     
              - ----------------------------------
                                             (1/2)
                2 (h[1] beta alpha l[1] l[2])     
a[1] := sqrt(h[1]*beta*alpha*l[1]*l[2])*(v-1)/(beta*h[1]);
                                           (1/2)
                (h[1] beta alpha l[1] l[2])     
                --------------------------------
                           beta h[1]            

Omega := lambda^2-4*mu*v+4*mu;
                               -3
H := (-lambda+sqrt(-Omega)*{tan(sqrt(-Omega)*xi)+sec(sqrt(-Omega)*xi)})/(2*(v-1));
         1   1  (1/2)  /   / (1/2)   \      / (1/2)   \\ 
       - - + - 3      { tan\3      xi/ + sec\3      xi/ }
         2   2         \                               / 

u := a[0]+a[1]*(d+H)+a[2]/(d+H);
            (2 d - 1) alpha l[1] l[2]            1     /
      - ---------------------------------- + --------- |
                                     (1/2)   beta h[1] \
        2 (h[1] beta alpha l[1] l[2])                   

                                   (1/2) /    1
        (h[1] beta alpha l[1] l[2])      |d - -
                                         \    2

           1  (1/2)  /   / (1/2)   \      / (1/2)   \\ \\
         + - 3      { tan\3      xi/ + sec\3      xi/ }||
           2         \                               / //
f := diff(u, xi);
               /                                         // 
        1      |                           (1/2)  (1/2) { | 
   ----------- \(h[1] beta alpha l[1] l[2])      3       \\1
   2 beta h[1]                                              

                      2\       
           / (1/2)   \ |  (1/2)
      + tan\3      xi/ / 3     

                                            \ \
           / (1/2)   \    / (1/2)   \  (1/2) }|
      + sec\3      xi/ tan\3      xi/ 3     / /
S := diff(f, xi);
             /                                         /      
      1      |                           (1/2)  (1/2) {      /
 ----------- \(h[1] beta alpha l[1] l[2])      3       \6 tan\
 2 beta h[1]                                                  

              /                  2\
    (1/2)   \ |       / (1/2)   \ |
   3      xi/ \1 + tan\3      xi/ /

                                     2
           / (1/2)   \    / (1/2)   \ 
    + 3 sec\3      xi/ tan\3      xi/ 

                       /                  2\\ \
           / (1/2)   \ |       / (1/2)   \ | }|
    + 3 sec\3      xi/ \1 + tan\3      xi/ // /

eq := -(alpha*h[1]*h[2]+omega)*u-(2*beta*h[1]*u*u)*u+alpha*l[1]*l[2]*S;
  3                 /      (2 d - 1) alpha l[1] l[2]        
- - alpha l[1] l[2] |- ---------------------------------- + 
  2                 |                               (1/2)   
                    \  2 (h[1] beta alpha l[1] l[2])        

      1     /                           (1/2) /    1
  --------- |(h[1] beta alpha l[1] l[2])      |d - -
  beta h[1] \                                 \    2

     1  (1/2)  /   / (1/2)   \      / (1/2)   \\ \\\             
   + - 3      { tan\3      xi/ + sec\3      xi/ }||| - 2 beta h[1
     2         \                               / //|             
                                                   /             

    /      (2 d - 1) alpha l[1] l[2]            1     /
  ] |- ---------------------------------- + --------- |
    |                               (1/2)   beta h[1] \
    \  2 (h[1] beta alpha l[1] l[2])                   

                             (1/2) /    1
  (h[1] beta alpha l[1] l[2])      |d - -
                                   \    2

     1  (1/2)  /   / (1/2)   \      / (1/2)   \\ \\\     
   + - 3      { tan\3      xi/ + sec\3      xi/ }|||^3 + 
     2         \                               / //|     
                                                   /     

              /                                                 
       1      |                                           (1/2) 
  ----------- \alpha l[1] l[2] (h[1] beta alpha l[1] l[2])      
  2 beta h[1]                                                   

          /                 /                  2\
   (1/2) {      / (1/2)   \ |       / (1/2)   \ |
  3       \6 tan\3      xi/ \1 + tan\3      xi/ /

                                    2
          / (1/2)   \    / (1/2)   \ 
   + 3 sec\3      xi/ tan\3      xi/ 

                      /                  2\\ \
          / (1/2)   \ |       / (1/2)   \ | }|
   + 3 sec\3      xi/ \1 + tan\3      xi/ // /
value(%);
  3                 /      (2 d - 1) alpha l[1] l[2]        
- - alpha l[1] l[2] |- ---------------------------------- + 
  2                 |                               (1/2)   
                    \  2 (h[1] beta alpha l[1] l[2])        

      1     /                           (1/2) /    1
  --------- |(h[1] beta alpha l[1] l[2])      |d - -
  beta h[1] \                                 \    2

     1  (1/2)  /   / (1/2)   \      / (1/2)   \\ \\\             
   + - 3      { tan\3      xi/ + sec\3      xi/ }||| - 2 beta h[1
     2         \                               / //|             
                                                   /             

    /      (2 d - 1) alpha l[1] l[2]            1     /
  ] |- ---------------------------------- + --------- |
    |                               (1/2)   beta h[1] \
    \  2 (h[1] beta alpha l[1] l[2])                   

                             (1/2) /    1
  (h[1] beta alpha l[1] l[2])      |d - -
                                   \    2

     1  (1/2)  /   / (1/2)   \      / (1/2)   \\ \\\     
   + - 3      { tan\3      xi/ + sec\3      xi/ }|||^3 + 
     2         \                               / //|     
                                                   /     

              /                                                 
       1      |                                           (1/2) 
  ----------- \alpha l[1] l[2] (h[1] beta alpha l[1] l[2])      
  2 beta h[1]                                                   

          /                 /                  2\
   (1/2) {      / (1/2)   \ |       / (1/2)   \ |
  3       \6 tan\3      xi/ \1 + tan\3      xi/ /

                                    2
          / (1/2)   \    / (1/2)   \ 
   + 3 sec\3      xi/ tan\3      xi/ 

                      /                  2\\ \
          / (1/2)   \ |       / (1/2)   \ | }|
   + 3 sec\3      xi/ \1 + tan\3      xi/ // /
simplify(%);
Error, (in simplify/power) invalid input: ln expects its 1st argument, x, to be of type algebraic, but received {(sin(3^(1/2)*xi)+1)/cos(3^(1/2)*xi)}
 

Hi guys

I want to solve the following differential equation but I can not. please help me in this way

diff(phi((8*R^(3/2)-W)*sqrt(2)/(24*sqrt(M))), W$2)=lambda*phi((8*R^(3/2)-W)*sqrt(2)/(24*sqrt(M)))

 

with the best regard

restart;


K := -3;
                               -3
m := 1;
                               1
w := -4*K;
                               12
alpha[0] := -2;
                               -2
alpha[1] := 0;
                               0
a := 2;
                               2
b := 3;
                               3
                               1
beta[1] := (12*(m^2+K))/(a+b);
                              -24
                              ---
                               5 
xi := -t*w+x+y;
                         -12 t + x + y

F := -sqrt(-K)*tanh(sqrt(-K)*xi);
                (1/2)     / (1/2)                \
              -3      tanh\3      (-12 t + x + y)/
U := alpha[0]+alpha[1]*(m+F)+beta[1]/(m+F);
                                 24                     
        -2 - -------------------------------------------
               /     (1/2)     / (1/2)                \\
             5 \1 - 3      tanh\3      (-12 t + x + y)//
                               1
y := 0;
                               0

plot3d(U, x = -10 .. -10, t = -10 .. 10);

t := 0;
                               0
plot(U, x = -10 .. 10);

RootOf(_Z^2*beta*h[1]-alpha*l[1]*l[2], label = _L2)

Hi I have the following ln function that I want to differentiate wrt variable c:

 

I2 := -sqrt(-c^2+1)*ln(abs((.9*sqrt(-c^2+1)-c*sqrt(1-.9^2))/(-.9*sqrt(-c^2+1)-c*sqrt(1-.9^2))))

 

When I differentiate I obtain an expression that involves 

...abs(1, (.9*sqrt(-c^2+1)-.4358898944*c)/(-.9*sqrt(-c^2+1)-.4358898944*c))...

Why does it give a comma at the abs expression? how to get rid of that.

Hello everyone,

I have a function of 5 variables, A, P, N, k. I want to solve and express as a power series of 'k'

w=F*k+G+H/k+...

and gather the coefficients of each power.

However, the result I obtain in my code differs from the analytical value I found. What am I doing wrong?

Thanks!

collect_mp.mw

How to get a plot for different values of Mh.

like Mh=[1 2 3 4]

Code:

restart;
with(DEtools,odeadvisor);

m:=10;H:=1;Mh:=1;b:=0.02; a:=0.05;V:=array(0..m); V[0]:=1-exp(-t);

for k from 1 to m do

if k=1 then chi:=0;

 chi:=1;

 fi;

 p:=0;

 for j from 0 to k-1 do

   p:=p+(V[k-1-j]*diff(V[j],t$2)-diff(V[k-1-j],t)*diff(V[j],t)-a*(2*diff(V[k-1-j],t)*diff(V[j],t$3)-diff(V[k-1-j],t$2)*diff(V[j],t$2)-V[k-1-j]*diff(V[j],t$4)));  od;

p:=(p+diff(V[k-1],t$3)-b*(diff(V[k-1],t$2)+t*diff(V[k-1],t$3))-Mh*diff(V[k-1],t))*h*H;

p:=factor(p);

V[k]:=(-int(p,t)+0.5*exp(t)*int(exp(-t)*p,t)+0.5*exp(-t)*int(exp(t)*p,t)+chi*V[k-1]+C1+C3*exp(-t));

v:=unapply(V[k],t);

V[k]:=frontend(expand,[V[k]]);  V[k]:=subs(C3=solve(eval(subs(t=0,diff(V[k],t))),C3),V[k]); V[k]:=frontend(expand,[V[k]]);

V[k]:=subs(C1=solve(eval(subs(t=0,-V[k]-diff(V[k],t))),C1),V[k]);

od:

appr:=0;

for k from 0 to m do

 appr:=appr+V[k];

od:

u_appr:=unapply(appr,(h,t)):

u_appr_1:=unapply(diff(u_appr(h,t),t),(h,t)):

evalf(u_appr_1(-0.4,t)):

with(plots);

plot([u_appr_1(-0.4,t)],t=0..4,0..1.2,color=[black],axes=frame):

 

 

this plot for Mh=1:

 

How to apply two for loops to solve ode problem.

code:

restart; with(plots); fcns := {T(eta), f(eta)};
m := .5; bet := 1; na := 1/6; N := 5;
eq1 := (diff(f(eta), `$`(eta, 3)))*pr+m-m*(diff(f(eta), `$`(eta, 1)))+((m+1)*(1/2))*(diff(f(eta), `$`(eta, 2)))*f(eta) = 0;
eq2 := diff(T(eta), `$`(eta, 2))+((m+1)*(1/2))*(diff(T(eta), `$`(eta, 1)))*f(eta) = 0;
bc := f(0) = 0, (D(f))(0) = 0, (D(f))(N) = 1, (D(T))(0) = -bi*(1-T(0)), T(N) = 0;
bi:= [seq(1..4,0.1)];  NN := nops(bi);  
pr:=[seq(1..2,0.1)];  NN1 := nops(pr);
for i  from 1 to NN do    
for j from 1 to NN1 do  

R := dsolve(eval({bc, eq1,eq2}, bi[i],pr[j]), fcns, type = numeric, method = bvp[midrich], maxmesh=2400):  
X1||[i,j]:=rhs(-R(0)[3]):
end do:  
end do:  

Have a good day.
 

I am trying to expand a multivariable (more specifically 4 variables) function in powers of one of its variables when it goes to infinity.

However, the result I get is always zero, even if I input (or not) values for some of the other variables.

Can anybody help?

series_expansion.mw

P.s.: I want to do the same for the other two functions I defined in the worksheet as well.

Dear friends,

Greetings.

How to get the second solution.

how to change the guess value in maple.

figure 1 plot in Matlab with two different initial guesses.

 

TWOSOLUTION.mw

 



 

Hi,

Due to an unexpected maintenance operation, I had to uninstall Maple 18 from my Windows 7-64 bit PC.

Later on I installed it without any problems but, to my surprise now I can't configure it properly. This is:

* Enabling Maple Text as Input (classical input method)
* Removing numbers from equations
* Setting Maple language to English (it took Spanish by default because of Windows)
* Hiding left panel

As usual, this is performed under Tools/Options/Interface and so on...

After I modify my preferred settings, apply globally and close-open again, the program is again in its original form.

How can I do in this case? I have installed it several times. It is also worth noting that the maintenance that I performed was related to deep Windows registry modifications. 

Thanks and regards.

 

restart;
T := -S(xi)*S(xi)+mu*R(xi)-lambda;
                        2                    
                  -S(xi)  + mu R(xi) - lambda
Q := -S(xi)*R(xi);
                          -S(xi) R(xi)
u := a[0]+a[1]*S(xi)+b[1]*R(xi);
                 a[0] + a[1] S(xi) + b[1] R(xi)
diff(u, xi);
                  / d        \        / d        \
             a[1] |---- S(xi)| + b[1] |---- R(xi)|
                  \ dxi      /        \ dxi      /
Fr := Q*b[1]+T*a[1];
                         /      2                    \     
     -S(xi) R(xi) b[1] + \-S(xi)  + mu R(xi) - lambda/ a[1]
diff(Fr, xi);
       / d        \                    / d        \     
      -|---- S(xi)| R(xi) b[1] - S(xi) |---- R(xi)| b[1]
       \ dxi      /                    \ dxi      /     

           /         / d        \      / d        \\     
         + |-2 S(xi) |---- S(xi)| + mu |---- R(xi)|| a[1]
           \         \ dxi      /      \ dxi      //     
d := -T*R(xi)*b[1]-S(xi)*Q*b[1]+(-2*S(xi)*T+mu*Q)*a[1];
 /      2                    \                   2           
-\-S(xi)  + mu R(xi) - lambda/ R(xi) b[1] + S(xi)  R(xi) b[1]

     /         /      2                    \                 \   
   + \-2 S(xi) \-S(xi)  + mu R(xi) - lambda/ - mu S(xi) R(xi)/ a[

  1]
diff(d, xi);
 /         / d        \      / d        \\           
-|-2 S(xi) |---- S(xi)| + mu |---- R(xi)|| R(xi) b[1]
 \         \ dxi      /      \ dxi      //           

     /      2                    \ / d        \     
   - \-S(xi)  + mu R(xi) - lambda/ |---- R(xi)| b[1]
                                   \ dxi      /     

                        / d        \        2 / d        \        /
   + 2 S(xi) R(xi) b[1] |---- S(xi)| + S(xi)  |---- R(xi)| b[1] + |
                        \ dxi      /          \ dxi      /        \
   / d        \ /      2                    \
-2 |---- S(xi)| \-S(xi)  + mu R(xi) - lambda/
   \ dxi      /                              

             /         / d        \      / d        \\
   - 2 S(xi) |-2 S(xi) |---- S(xi)| + mu |---- R(xi)||
             \         \ dxi      /      \ dxi      //

        / d        \                  / d        \\     
   - mu |---- S(xi)| R(xi) - mu S(xi) |---- R(xi)|| a[1]
        \ dxi      /                  \ dxi      //     
h := -(-2*S(xi)*T+mu*Q)*R(xi)*b[1]-(-S(xi)^2+mu*R(xi)-lambda)*Q*b[1]+2*S(xi)*R(xi)*b[1]*T+S(xi)^2*Q*b[1]+(-2*T*(-S(xi)^2+mu*R(xi)-lambda)-2*S(xi)*(-2*S(xi)*T+mu*Q)-mu*T*R(xi)-mu*S(xi)*Q)*a[1];
 /         /      2                    \                 \       
-\-2 S(xi) \-S(xi)  + mu R(xi) - lambda/ - mu S(xi) R(xi)/ R(xi) 

           /      2                    \                 
  b[1] + 3 \-S(xi)  + mu R(xi) - lambda/ S(xi) R(xi) b[1]

                         /                                2         
          3              |   /      2                    \          
   - S(xi)  R(xi) b[1] + \-2 \-S(xi)  + mu R(xi) - lambda/  - 2 S(xi

    /         /      2                    \                 \
  ) \-2 S(xi) \-S(xi)  + mu R(xi) - lambda/ - mu S(xi) R(xi)/

                                                             \   
        /      2                    \                 2      |   
   - mu \-S(xi)  + mu R(xi) - lambda/ R(xi) + mu S(xi)  R(xi)/ a[

  1]
collect(expand(h+3*Fr*Fr+(4*omega+3)*Fr), S(xi), R(xi));
     /      2         \      4
R(xi)\3 a[1]  - 6 a[1]/ S(xi) 

                                                  3        /  
   + R(xi)(6 R(xi) a[1] b[1] - 6 b[1] R(xi)) S(xi)  + R(xi)\3 

       2     2                  2                   
  R(xi)  b[1]  - 6 R(xi) mu a[1]  + 12 a[1] mu R(xi)

                  2                                        \ 
   + 6 lambda a[1]  - 8 a[1] lambda - 4 omega a[1] - 3 a[1]/ 

       2        /        2                       2        
  S(xi)  + R(xi)\-6 R(xi)  mu a[1] b[1] + 6 R(xi)  mu b[1]

   + 6 R(xi) lambda a[1] b[1] - 5 R(xi) lambda b[1]

                                      \              /      2   2 
   - 4 R(xi) omega b[1] - 3 b[1] R(xi)/ S(xi) + R(xi)\3 a[1]  mu  

       2            2      2         2                
  R(xi)  - 3 a[1] mu  R(xi)  - 6 a[1]  mu R(xi) lambda

   + 5 a[1] mu R(xi) lambda + 4 omega a[1] mu R(xi)

           2       2                                  2
   + 3 a[1]  lambda  + 3 a[1] mu R(xi) - 2 a[1] lambda 

                                        \
   - 4 omega a[1] lambda - 3 a[1] lambda/

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