Maple 18 Questions and Posts

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

I'm trying to do something very simple, but I can not do it. I would like to fill the chart with colors of my choice.

restart: with (plots):
plots [animate] (plot, [[sqrt (x), sqrt (x) -1], x = 0..t, filled = true, view = [0..20, 0..5]], t = 0..20);

The filled = true option fills the graph with random colors.
I tried to use filled = ["Blue", "Red"], but that does not work.

Any tips?

Thank you

 

Please, I need to use Maple to solve Euler-Bernoulli Beam on Pasternak Foundation using Homotopy Perturbation Method.

The governing equation is 

initial conditions are 

the boundary condition is 

The governing equation represents Euler-Bernoulli beam on a generalized Pasternak viscoelastic foundation under an
arbitrary distributed dynamic load. in which E, I , ρ ,A are the parameters of the beam, representing Young’s modulus of elasticity, moment of inertia, density and area of cross section, respectively. K,C and Gp are spring stiffness, damping coefficient, and shear coefficient of the foundation. Moreover, y(x, t) and F(x, t) are defined as the vertical deflection of the beam and the generic arbitrary dynamic loads, respectively, where the loads distribute along the x-axis and t is time.

 

I will appreciate anyone who can help me with a Maple solution.

 

Thank you.

restart;
lambda := 1;
                               1
mu := 1;
                               1
alpha := 1;
                               1
v := 2;
                               2
delta := 1;
                               1
m := 1;
                               1
d := 3;
                               3
l := 1;
                               1
omega := -(1/2)*a*lambda^2-a*m^2+2*a*mu*v-delta*l*m-(1/2)*delta*lambda^2+2*delta*mu*v-l^2*mu-(1/2)*lambda^2*mu+2*mu^2*v-2*a*mu-2*delta*mu-2*mu^2;
                            1      
                            - a + 1
                            2      
a[0] := ((2*d*v-2*d-lambda)*(1/2))*sqrt(2)*sqrt(gamma*(a+delta+mu))/gamma;
                    (1/2)                (1/2)
                 5 2      (gamma (a + 2))     
                 -----------------------------
                            2 gamma           
a[1] := 0;
                               0
a[2] := -(2*(1/2))*sqrt(2)*sqrt(gamma*(a+delta+mu))*(d^2*v-d^2-d*lambda+mu)/gamma;
                     (1/2)                (1/2)
                  7 2      (gamma (a + 2))     
                - -----------------------------
                              gamma            
Omega := lambda^2-4*mu*v+4*mu;
                               -3
H := (-lambda+sqrt(-Omega)*tan((1/2)*sqrt(-Omega)*xi))/(2*(v-1));
                  1   1  (1/2)    /1  (1/2)   \
                - - + - 3      tan|- 3      xi|
                  2   2           \2          /
u := a[0]+a[1]*(d+H)+a[2]/(d+H);
              (1/2)                (1/2)
           5 2      (gamma (a + 2))     
           -----------------------------
                      2 gamma           

                       (1/2)                (1/2)    
                    7 2      (gamma (a + 2))         
              - -------------------------------------
                      /5   1  (1/2)    /1  (1/2)   \\
                gamma |- + - 3      tan|- 3      xi||
                      \2   2           \2          //

eta := -l*y-m*x+omega*t;
                               /1      \  
                      -y - x + |- a + 1| t
                               \2      /  
u := a[0]+a[1]*(d+H)+a[2]/(d+H);
              (1/2)                (1/2)
           5 2      (gamma (a + 2))     
           -----------------------------
                      2 gamma           

                       (1/2)                (1/2)    
                    7 2      (gamma (a + 2))         
              - -------------------------------------
                      /5   1  (1/2)    /1  (1/2)   \\
                gamma |- + - 3      tan|- 3      xi||
                      \2   2           \2          //

f := diff(u, xi);
                                  /                      2\
       (1/2)                (1/2) |3   3    /1  (1/2)   \ |
    7 2      (gamma (a + 2))      |- + - tan|- 3      xi| |
                                  \4   4    \2          / /
    -------------------------------------------------------
                                                 2         
                  /5   1  (1/2)    /1  (1/2)   \\          
            gamma |- + - 3      tan|- 3      xi||          
                  \2   2           \2          //          
S := diff(u, xi);
                                  /                      2\
       (1/2)                (1/2) |3   3    /1  (1/2)   \ |
    7 2      (gamma (a + 2))      |- + - tan|- 3      xi| |
                                  \4   4    \2          / /
    -------------------------------------------------------
                                                 2         
                  /5   1  (1/2)    /1  (1/2)   \\          
            gamma |- + - 3      tan|- 3      xi||          
                  \2   2           \2          //          
xi := x+y-(-2*alpha*m-delta*l-delta*m-2*l*mu)*t;
                          x + y + 6 t
eq := (-omega-a*m*m-delta*l*m-mu*l*l)*u-(gamma*u*u)*u+(a+delta+mu)*S;
            /   (1/2)                (1/2)
/  3      \ |5 2      (gamma (a + 2))     
|- - a - 3| |-----------------------------
\  2      / |           2 gamma           
            |                             
            \                             

                 (1/2)                (1/2)          \         /
              7 2      (gamma (a + 2))               |         |
   - ------------------------------------------------| - gamma |
           /5   1  (1/2)    /1  (1/2)              \\|         |
     gamma |- + - 3      tan|- 3      (x + y + 6 t)|||         |
           \2   2           \2                     ///         \

     (1/2)                (1/2)
  5 2      (gamma (a + 2))     
  -----------------------------
             2 gamma           

                 (1/2)                (1/2)          \     
              7 2      (gamma (a + 2))               |     
   - ------------------------------------------------|^3 + 
           /5   1  (1/2)    /1  (1/2)              \\|     
     gamma |- + - 3      tan|- 3      (x + y + 6 t)|||     
           \2   2           \2                     ///     

                                                    /          
                          1                         |          
  ------------------------------------------------- |7 (a + 2) 
                                                  2 \          
        /5   1  (1/2)    /1  (1/2)              \\             
  gamma |- + - 3      tan|- 3      (x + y + 6 t)||             
        \2   2           \2                     //             

                              /                                 2
   (1/2)                (1/2) |3   3    /1  (1/2)              \ 
  2      (gamma (a + 2))      |- + - tan|- 3      (x + y + 6 t)| 
                              \4   4    \2                     / 

  \\
  ||
  ||
  //
value(%);
            /   (1/2)                (1/2)
/  3      \ |5 2      (gamma (a + 2))     
|- - a - 3| |-----------------------------
\  2      / |           2 gamma           
            |                             
            \                             

                 (1/2)                (1/2)          \         /
              7 2      (gamma (a + 2))               |         |
   - ------------------------------------------------| - gamma |
           /5   1  (1/2)    /1  (1/2)              \\|         |
     gamma |- + - 3      tan|- 3      (x + y + 6 t)|||         |
           \2   2           \2                     ///         \

     (1/2)                (1/2)
  5 2      (gamma (a + 2))     
  -----------------------------
             2 gamma           

                 (1/2)                (1/2)          \     
              7 2      (gamma (a + 2))               |     
   - ------------------------------------------------|^3 + 
           /5   1  (1/2)    /1  (1/2)              \\|     
     gamma |- + - 3      tan|- 3      (x + y + 6 t)|||     
           \2   2           \2                     ///     

                                                    /          
                          1                         |          
  ------------------------------------------------- |7 (a + 2) 
                                                  2 \          
        /5   1  (1/2)    /1  (1/2)              \\             
  gamma |- + - 3      tan|- 3      (x + y + 6 t)||             
        \2   2           \2                     //             

                              /                                 2
   (1/2)                (1/2) |3   3    /1  (1/2)              \ 
  2      (gamma (a + 2))      |- + - tan|- 3      (x + y + 6 t)| 
                              \4   4    \2                     / 

  \\
  ||
  ||
  //

simplify(%);
Error, (in simplify/tools/_zn) too many levels of recursion

Q := exp(I*(-y-x+((1/2)*a+1)*t))*((5/2)*sqrt(2)*sqrt(gamma*(a+2))/gamma-7*sqrt(2)*sqrt(gamma*(a+2))/(gamma*(5/2+(1/2)*sqrt(3)*tan((1/2)*sqrt(3)*(x+y+6*t)))));
                                /   (1/2)                (1/2)
     /  /         /1      \  \\ |5 2      (gamma (a + 2))     
  exp|I |-y - x + |- a + 1| t|| |-----------------------------
     \  \         \2      /  // |           2 gamma           
                                |                             
                                \                             

                   (1/2)                (1/2)          \
                7 2      (gamma (a + 2))               |
     - ------------------------------------------------|
             /5   1  (1/2)    /1  (1/2)              \\|
       gamma |- + - 3      tan|- 3      (x + y + 6 t)|||
             \2   2           \2                     ///
                                /   (1/2)                (1/2)
     /  /         /1      \  \\ |5 2      (gamma (a + 2))     
  exp|I |-y - x + |- a + 1| t|| |-----------------------------
     \  \         \2      /  // |           2 gamma           
                                |                             
                                \                             

                   (1/2)                (1/2)          \
                7 2      (gamma (a + 2))               |
     - ------------------------------------------------|
             /5   1  (1/2)    /1  (1/2)              \\|
       gamma |- + - 3      tan|- 3      (x + y + 6 t)|||
             \2   2           \2                     ///
plot3d(Re(Q), x = -40 .. -37, y = -40 .. -37, t = -10 .. 10);
 

Hy everyone,

i'm writing a code wich ends with an Explore including multiple parameters (8) ; my problem is that if i want to change every parameter, elaborating time becames too long, because Maple evaluates every variation of every parameter every time.
I'm looking for something like a button, which works like this: when i change parameters in explore maple doesn't evaluate the function, and when i click the button maple evaluates one time all the changed parameters.

Any suggestion? 

Thanks in advance

restart;
Solve({(6*(-8*L*alpha^2*a[1]^3+2*L*a[1]^3-2*a[1]^3))*d-24*L*alpha^2*a[0]*a[1]^2+6*L*a[0]*a[1]^2-6*a[0]*a[1]^2+12*alpha^2*k^2*lambda*v*a[1]-12*alpha^2*k^2*lambda*a[1]-3*k^2*lambda*v*a[1]+3*k^2*lambda*a[1]+(3*(8*alpha^2*k^2*v^2*a[1]-16*alpha^2*k^2*v*a[1]+8*alpha^2*k^2*a[1]-2*k^2*v^2*a[1]+4*k^2*v*a[1]-2*k^2*a[1]))*d = 0, (-8*L*alpha^2*a[1]^3+2*L*a[1]^3-2*a[1]^3)*d^6+(-24*L*alpha^2*a[0]*a[1]^2+6*L*a[0]*a[1]^2-6*a[0]*a[1]^2)*d^5+(-24*L*alpha^2*a[0]^2*a[1]-24*L*alpha^2*a[1]^2*a[2]-4*alpha^4*a[1]+6*L*a[0]^2*a[1]+6*L*a[1]^2*a[2]-4*alpha^2*beta*a[1]+alpha^2*a[1]-6*a[0]^2*a[1]-6*a[1]^2*a[2]+beta*a[1])*d^4+(4*alpha^2*k^2*lambda*mu*a[1]-8*L*alpha^2*a[0]^3-48*L*alpha^2*a[0]*a[1]*a[2]-4*alpha^4*a[0]-k^2*lambda*mu*a[1]+2*L*a[0]^3+12*L*a[0]*a[1]*a[2]-4*alpha^2*beta*a[0]+alpha^2*a[0]-2*a[0]^3-12*a[0]*a[1]*a[2]+beta*a[0])*d^3+(-24*L*alpha^2*a[0]^2*a[2]-24*L*alpha^2*a[1]*a[2]^2-4*alpha^4*a[2]+6*L*a[0]^2*a[2]+6*L*a[1]*a[2]^2-4*alpha^2*beta*a[2]+alpha^2*a[2]-6*a[0]^2*a[2]-6*a[1]*a[2]^2+beta*a[2])*d^2+(-4*alpha^2*k^2*lambda*mu*a[2]-24*L*alpha^2*a[0]*a[2]^2+k^2*lambda*mu*a[2]+6*L*a[0]*a[2]^2-6*a[0]*a[2]^2)*d+8*alpha^2*k^2*mu^2*a[2]-2*a[2]^3-2*k^2*mu^2*a[2]-8*L*alpha^2*a[2]^3+2*L*a[2]^3 = 0, (6*(-8*L*alpha^2*a[1]^3+2*L*a[1]^3-2*a[1]^3))*d^5+(5*(-24*L*alpha^2*a[0]*a[1]^2+6*L*a[0]*a[1]^2-6*a[0]*a[1]^2))*d^4+(4*(-24*L*alpha^2*a[0]^2*a[1]-24*L*alpha^2*a[1]^2*a[2]-4*alpha^4*a[1]+6*L*a[0]^2*a[1]+6*L*a[1]^2*a[2]-4*alpha^2*beta*a[1]+alpha^2*a[1]-6*a[0]^2*a[1]-6*a[1]^2*a[2]+beta*a[1]))*d^3+(3*(4*alpha^2*k^2*lambda*mu*a[1]-8*L*alpha^2*a[0]^3-48*L*alpha^2*a[0]*a[1]*a[2]-4*alpha^4*a[0]-k^2*lambda*mu*a[1]+2*L*a[0]^3+12*L*a[0]*a[1]*a[2]-4*alpha^2*beta*a[0]+alpha^2*a[0]-2*a[0]^3-12*a[0]*a[1]*a[2]+beta*a[0]))*d^2+(4*alpha^2*k^2*lambda^2*a[1]+8*alpha^2*k^2*mu*v*a[1]-8*alpha^2*k^2*mu*a[1]-k^2*lambda^2*a[1]-2*k^2*mu*v*a[1]+2*k^2*mu*a[1])*d^3+(2*(-24*L*alpha^2*a[0]^2*a[2]-24*L*alpha^2*a[1]*a[2]^2-4*alpha^4*a[2]+6*L*a[0]^2*a[2]+6*L*a[1]*a[2]^2-4*alpha^2*beta*a[2]+alpha^2*a[2]-6*a[0]^2*a[2]-6*a[1]*a[2]^2+beta*a[2]))*d+12*alpha^2*k^2*lambda*mu*a[2]-24*L*alpha^2*a[0]*a[2]^2-3*k^2*lambda*mu*a[2]+6*L*a[0]*a[2]^2-6*a[0]*a[2]^2+(-4*alpha^2*k^2*lambda^2*a[2]-8*alpha^2*k^2*mu*v*a[2]+8*alpha^2*k^2*mu*a[2]+k^2*lambda^2*a[2]+2*k^2*mu*v*a[2]-2*k^2*mu*a[2])*d = 0, (15*(-8*L*alpha^2*a[1]^3+2*L*a[1]^3-2*a[1]^3))*d^2+(5*(-24*L*alpha^2*a[0]*a[1]^2+6*L*a[0]*a[1]^2-6*a[0]*a[1]^2))*d+(3*(12*alpha^2*k^2*lambda*v*a[1]-12*alpha^2*k^2*lambda*a[1]-3*k^2*lambda*v*a[1]+3*k^2*lambda*a[1]))*d+(3*(8*alpha^2*k^2*v^2*a[1]-16*alpha^2*k^2*v*a[1]+8*alpha^2*k^2*a[1]-2*k^2*v^2*a[1]+4*k^2*v*a[1]-2*k^2*a[1]))*d^2+8*alpha^2*k^2*mu*v*a[1]-6*a[0]^2*a[1]-6*a[1]^2*a[2]-4*alpha^4*a[1]+4*alpha^2*k^2*lambda^2*a[1]-8*alpha^2*k^2*mu*a[1]-2*k^2*mu*v*a[1]-24*L*alpha^2*a[0]^2*a[1]-24*L*alpha^2*a[1]^2*a[2]-k^2*lambda^2*a[1]+2*k^2*mu*a[1]+6*L*a[0]^2*a[1]+6*L*a[1]^2*a[2]-4*alpha^2*beta*a[1]+alpha^2*a[1]+beta*a[1] = 0, (15*(-8*L*alpha^2*a[1]^3+2*L*a[1]^3-2*a[1]^3))*d^4+(10*(-24*L*alpha^2*a[0]*a[1]^2+6*L*a[0]*a[1]^2-6*a[0]*a[1]^2))*d^3+(6*(-24*L*alpha^2*a[0]^2*a[1]-24*L*alpha^2*a[1]^2*a[2]-4*alpha^4*a[1]+6*L*a[0]^2*a[1]+6*L*a[1]^2*a[2]-4*alpha^2*beta*a[1]+alpha^2*a[1]-6*a[0]^2*a[1]-6*a[1]^2*a[2]+beta*a[1]))*d^2+(12*alpha^2*k^2*lambda*v*a[1]-12*alpha^2*k^2*lambda*a[1]-3*k^2*lambda*v*a[1]+3*k^2*lambda*a[1])*d^3+(3*(4*alpha^2*k^2*lambda*mu*a[1]-8*L*alpha^2*a[0]^3-48*L*alpha^2*a[0]*a[1]*a[2]-4*alpha^4*a[0]-k^2*lambda*mu*a[1]+2*L*a[0]^3+12*L*a[0]*a[1]*a[2]-4*alpha^2*beta*a[0]+alpha^2*a[0]-2*a[0]^3-12*a[0]*a[1]*a[2]+beta*a[0]))*d+(3*(4*alpha^2*k^2*lambda^2*a[1]+8*alpha^2*k^2*mu*v*a[1]-8*alpha^2*k^2*mu*a[1]-k^2*lambda^2*a[1]-2*k^2*mu*v*a[1]+2*k^2*mu*a[1]))*d^2+(-12*alpha^2*k^2*lambda*v*a[2]+12*alpha^2*k^2*lambda*a[2]+3*k^2*lambda*v*a[2]-3*k^2*lambda*a[2])*d+8*alpha^2*k^2*mu*v*a[2]-6*a[0]^2*a[2]-6*a[1]*a[2]^2-4*alpha^4*a[2]+4*alpha^2*k^2*lambda^2*a[2]-8*alpha^2*k^2*mu*a[2]-2*k^2*mu*v*a[2]-24*L*alpha^2*a[0]^2*a[2]-24*L*alpha^2*a[1]*a[2]^2-k^2*lambda^2*a[2]+2*k^2*mu*a[2]+6*L*a[0]^2*a[2]+6*L*a[1]*a[2]^2-4*alpha^2*beta*a[2]+alpha^2*a[2]+beta*a[2] = 0, (20*(-8*L*alpha^2*a[1]^3+2*L*a[1]^3-2*a[1]^3))*d^3+(10*(-24*L*alpha^2*a[0]*a[1]^2+6*L*a[0]*a[1]^2-6*a[0]*a[1]^2))*d^2+(4*(-24*L*alpha^2*a[0]^2*a[1]-24*L*alpha^2*a[1]^2*a[2]-4*alpha^4*a[1]+6*L*a[0]^2*a[1]+6*L*a[1]^2*a[2]-4*alpha^2*beta*a[1]+alpha^2*a[1]-6*a[0]^2*a[1]-6*a[1]^2*a[2]+beta*a[1]))*d+(8*alpha^2*k^2*v^2*a[1]-16*alpha^2*k^2*v*a[1]+8*alpha^2*k^2*a[1]-2*k^2*v^2*a[1]+4*k^2*v*a[1]-2*k^2*a[1])*d^3+(3*(4*alpha^2*k^2*lambda^2*a[1]+8*alpha^2*k^2*mu*v*a[1]-8*alpha^2*k^2*mu*a[1]-k^2*lambda^2*a[1]-2*k^2*mu*v*a[1]+2*k^2*mu*a[1]))*d+(3*(12*alpha^2*k^2*lambda*v*a[1]-12*alpha^2*k^2*lambda*a[1]-3*k^2*lambda*v*a[1]+3*k^2*lambda*a[1]))*d^2+(-8*alpha^2*k^2*v^2*a[2]+16*alpha^2*k^2*v*a[2]-8*alpha^2*k^2*a[2]+2*k^2*v^2*a[2]-4*k^2*v*a[2]+2*k^2*a[2])*d+4*alpha^2*k^2*lambda*v*a[2]+2*L*a[0]^3-4*alpha^4*a[0]-48*L*alpha^2*a[0]*a[1]*a[2]+4*alpha^2*k^2*lambda*mu*a[1]-4*alpha^2*k^2*lambda*a[2]-k^2*lambda*v*a[2]-k^2*lambda*mu*a[1]+12*L*a[0]*a[1]*a[2]+k^2*lambda*a[2]-8*L*alpha^2*a[0]^3-12*a[0]*a[1]*a[2]-4*alpha^2*beta*a[0]+alpha^2*a[0]+beta*a[0]-2*a[0]^3 = 0, 8*alpha^2*k^2*v^2*a[1]-8*L*alpha^2*a[1]^3-16*alpha^2*k^2*v*a[1]+8*alpha^2*k^2*a[1]-2*k^2*v^2*a[1]+2*L*a[1]^3+4*k^2*v*a[1]-2*k^2*a[1]-2*a[1]^3 = 0}, {alpha, a[0], a[1], a[2]});
     / /  /          2     3           3         3\  
Solve|{ 6 \-8 L alpha  a[1]  + 2 L a[1]  - 2 a[1] / d
     \ \                                             

               2          2                2              2
   - 24 L alpha  a[0] a[1]  + 6 L a[0] a[1]  - 6 a[0] a[1] 

             2  2                         2  2            
   + 12 alpha  k  lambda v a[1] - 12 alpha  k  lambda a[1]

        2                    2                 /       2  2  2   
   - 3 k  lambda v a[1] + 3 k  lambda a[1] + 3 \8 alpha  k  v  a[

               2  2                 2  2           2  2     
  1] - 16 alpha  k  v a[1] + 8 alpha  k  a[1] - 2 k  v  a[1]

        2             2     \        /          2     3
   + 4 k  v a[1] - 2 k  a[1]/ d = 0, \-8 L alpha  a[1] 

             3         3\  6
   + 2 L a[1]  - 2 a[1] / d 

     /           2          2                2              2\  5   
   + \-24 L alpha  a[0] a[1]  + 6 L a[0] a[1]  - 6 a[0] a[1] / d  + 

  /           2     2                  2     2     
  \-24 L alpha  a[0]  a[1] - 24 L alpha  a[1]  a[2]

            4                2                2     
   - 4 alpha  a[1] + 6 L a[0]  a[1] + 6 L a[1]  a[2]

            2                  2              2     
   - 4 alpha  beta a[1] + alpha  a[1] - 6 a[0]  a[1]

           2                 \  4   /       2  2               
   - 6 a[1]  a[2] + beta a[1]/ d  + \4 alpha  k  lambda mu a[1]

              2     3             2               
   - 8 L alpha  a[0]  - 48 L alpha  a[0] a[1] a[2]

            4         2                          3
   - 4 alpha  a[0] - k  lambda mu a[1] + 2 L a[0] 

                                  2                  2     
   + 12 L a[0] a[1] a[2] - 4 alpha  beta a[0] + alpha  a[0]

           3                                \  3   /
   - 2 a[0]  - 12 a[0] a[1] a[2] + beta a[0]/ d  + \
           2     2                  2          2          4     
-24 L alpha  a[0]  a[2] - 24 L alpha  a[1] a[2]  - 4 alpha  a[2]

             2                     2          2          
   + 6 L a[0]  a[2] + 6 L a[1] a[2]  - 4 alpha  beta a[2]

          2              2                   2            \  2   
   + alpha  a[2] - 6 a[0]  a[2] - 6 a[1] a[2]  + beta a[2]/ d  + 

  /        2  2                            2          2
  \-4 alpha  k  lambda mu a[2] - 24 L alpha  a[0] a[2] 

      2                               2              2\  
   + k  lambda mu a[2] + 6 L a[0] a[2]  - 6 a[0] a[2] / d

            2  2   2              3      2   2     
   + 8 alpha  k  mu  a[2] - 2 a[2]  - 2 k  mu  a[2]

              2     3           3        /          2     3
   - 8 L alpha  a[2]  + 2 L a[2]  = 0, 6 \-8 L alpha  a[1] 

             3         3\  5
   + 2 L a[1]  - 2 a[1] / d 

       /           2          2                2              2\ 
   + 5 \-24 L alpha  a[0] a[1]  + 6 L a[0] a[1]  - 6 a[0] a[1] / 

   4     /           2     2                  2     2     
  d  + 4 \-24 L alpha  a[0]  a[1] - 24 L alpha  a[1]  a[2]

            4                2                2     
   - 4 alpha  a[1] + 6 L a[0]  a[1] + 6 L a[1]  a[2]

            2                  2              2     
   - 4 alpha  beta a[1] + alpha  a[1] - 6 a[0]  a[1]

           2                 \  3     /       2  2               
   - 6 a[1]  a[2] + beta a[1]/ d  + 3 \4 alpha  k  lambda mu a[1]

              2     3             2               
   - 8 L alpha  a[0]  - 48 L alpha  a[0] a[1] a[2]

            4         2                          3
   - 4 alpha  a[0] - k  lambda mu a[1] + 2 L a[0] 

                                  2                  2     
   + 12 L a[0] a[1] a[2] - 4 alpha  beta a[0] + alpha  a[0]

           3                                \  2   /       2  2 
   - 2 a[0]  - 12 a[0] a[1] a[2] + beta a[0]/ d  + \4 alpha  k  

        2               2  2                    2  2        
  lambda  a[1] + 8 alpha  k  mu v a[1] - 8 alpha  k  mu a[1]

      2       2           2                2        \  3     /
   - k  lambda  a[1] - 2 k  mu v a[1] + 2 k  mu a[1]/ d  + 2 \
           2     2                  2          2          4     
-24 L alpha  a[0]  a[2] - 24 L alpha  a[1] a[2]  - 4 alpha  a[2]

             2                     2          2          
   + 6 L a[0]  a[2] + 6 L a[1] a[2]  - 4 alpha  beta a[2]

          2              2                   2            \  
   + alpha  a[2] - 6 a[0]  a[2] - 6 a[1] a[2]  + beta a[2]/ d

             2  2                            2          2
   + 12 alpha  k  lambda mu a[2] - 24 L alpha  a[0] a[2] 

        2                               2              2   /
   - 3 k  lambda mu a[2] + 6 L a[0] a[2]  - 6 a[0] a[2]  + \
        2  2       2               2  2          
-4 alpha  k  lambda  a[2] - 8 alpha  k  mu v a[2]

            2  2            2       2           2          
   + 8 alpha  k  mu a[2] + k  lambda  a[2] + 2 k  mu v a[2]

        2        \           /          2     3           3
   - 2 k  mu a[2]/ d = 0, 15 \-8 L alpha  a[1]  + 2 L a[1] 

           3\  2
   - 2 a[1] / d 

       /           2          2                2              2\       
   + 5 \-24 L alpha  a[0] a[1]  + 6 L a[0] a[1]  - 6 a[0] a[1] / d + 3 

  /        2  2                         2  2            
  \12 alpha  k  lambda v a[1] - 12 alpha  k  lambda a[1]

        2                    2            \       /       2  2  2 
   - 3 k  lambda v a[1] + 3 k  lambda a[1]/ d + 3 \8 alpha  k  v  

                 2  2                 2  2           2  2     
  a[1] - 16 alpha  k  v a[1] + 8 alpha  k  a[1] - 2 k  v  a[1]

        2             2     \  2          2  2          
   + 4 k  v a[1] - 2 k  a[1]/ d  + 8 alpha  k  mu v a[1]

           2              2               4     
   - 6 a[0]  a[1] - 6 a[1]  a[2] - 4 alpha  a[1]

            2  2       2               2  2        
   + 4 alpha  k  lambda  a[1] - 8 alpha  k  mu a[1]

        2                       2     2     
   - 2 k  mu v a[1] - 24 L alpha  a[0]  a[1]

               2     2         2       2           2        
   - 24 L alpha  a[1]  a[2] - k  lambda  a[1] + 2 k  mu a[1]

             2                2               2          
   + 6 L a[0]  a[1] + 6 L a[1]  a[2] - 4 alpha  beta a[1]

          2                          /          2     3
   + alpha  a[1] + beta a[1] = 0, 20 \-8 L alpha  a[1] 

             3         3\  3
   + 2 L a[1]  - 2 a[1] / d 

        /           2          2                2              2\ 
   + 10 \-24 L alpha  a[0] a[1]  + 6 L a[0] a[1]  - 6 a[0] a[1] / 

   2     /           2     2                  2     2     
  d  + 4 \-24 L alpha  a[0]  a[1] - 24 L alpha  a[1]  a[2]

            4                2                2     
   - 4 alpha  a[1] + 6 L a[0]  a[1] + 6 L a[1]  a[2]

            2                  2              2     
   - 4 alpha  beta a[1] + alpha  a[1] - 6 a[0]  a[1]

           2                 \     /       2  2  2     
   - 6 a[1]  a[2] + beta a[1]/ d + \8 alpha  k  v  a[1]

             2  2                 2  2           2  2     
   - 16 alpha  k  v a[1] + 8 alpha  k  a[1] - 2 k  v  a[1]

        2             2     \  3     /       2  2       2     
   + 4 k  v a[1] - 2 k  a[1]/ d  + 3 \4 alpha  k  lambda  a[1]

            2  2                    2  2        
   + 8 alpha  k  mu v a[1] - 8 alpha  k  mu a[1]

      2       2           2                2        \       /   
   - k  lambda  a[1] - 2 k  mu v a[1] + 2 k  mu a[1]/ d + 3 \12 

       2  2                         2  2            
  alpha  k  lambda v a[1] - 12 alpha  k  lambda a[1]

        2                    2            \  2   /
   - 3 k  lambda v a[1] + 3 k  lambda a[1]/ d  + \
        2  2  2                2  2                 2  2     
-8 alpha  k  v  a[2] + 16 alpha  k  v a[2] - 8 alpha  k  a[2]

        2  2           2             2     \  
   + 2 k  v  a[2] - 4 k  v a[2] + 2 k  a[2]/ d

            2  2                         3          4     
   + 4 alpha  k  lambda v a[2] + 2 L a[0]  - 4 alpha  a[0]

               2                         2  2               
   - 48 L alpha  a[0] a[1] a[2] + 4 alpha  k  lambda mu a[1]

            2  2                2              
   - 4 alpha  k  lambda a[2] - k  lambda v a[2]

      2                                         2            
   - k  lambda mu a[1] + 12 L a[0] a[1] a[2] + k  lambda a[2]

              2     3                              2          
   - 8 L alpha  a[0]  - 12 a[0] a[1] a[2] - 4 alpha  beta a[0]

          2                          3         /          2     3
   + alpha  a[0] + beta a[0] - 2 a[0]  = 0, 15 \-8 L alpha  a[1] 

             3         3\  4
   + 2 L a[1]  - 2 a[1] / d 

        /           2          2                2              2\ 
   + 10 \-24 L alpha  a[0] a[1]  + 6 L a[0] a[1]  - 6 a[0] a[1] / 

   3     /           2     2                  2     2     
  d  + 6 \-24 L alpha  a[0]  a[1] - 24 L alpha  a[1]  a[2]

            4                2                2     
   - 4 alpha  a[1] + 6 L a[0]  a[1] + 6 L a[1]  a[2]

            2                  2              2     
   - 4 alpha  beta a[1] + alpha  a[1] - 6 a[0]  a[1]

           2                 \  2   /        2  2              
   - 6 a[1]  a[2] + beta a[1]/ d  + \12 alpha  k  lambda v a[1]

             2  2                  2              
   - 12 alpha  k  lambda a[1] - 3 k  lambda v a[1]

        2            \  3     /       2  2               
   + 3 k  lambda a[1]/ d  + 3 \4 alpha  k  lambda mu a[1]

              2     3             2               
   - 8 L alpha  a[0]  - 48 L alpha  a[0] a[1] a[2]

            4         2                          3
   - 4 alpha  a[0] - k  lambda mu a[1] + 2 L a[0] 

                                  2                  2     
   + 12 L a[0] a[1] a[2] - 4 alpha  beta a[0] + alpha  a[0]

           3                                \       /       2  2 
   - 2 a[0]  - 12 a[0] a[1] a[2] + beta a[0]/ d + 3 \4 alpha  k  

        2               2  2                    2  2        
  lambda  a[1] + 8 alpha  k  mu v a[1] - 8 alpha  k  mu a[1]

      2       2           2                2        \  2   /
   - k  lambda  a[1] - 2 k  mu v a[1] + 2 k  mu a[1]/ d  + \
         2  2                         2  2            
-12 alpha  k  lambda v a[2] + 12 alpha  k  lambda a[2]

        2                    2            \  
   + 3 k  lambda v a[2] - 3 k  lambda a[2]/ d

            2  2                   2                   2
   + 8 alpha  k  mu v a[2] - 6 a[0]  a[2] - 6 a[1] a[2] 

            4               2  2       2     
   - 4 alpha  a[2] + 4 alpha  k  lambda  a[2]

            2  2              2          
   - 8 alpha  k  mu a[2] - 2 k  mu v a[2]

               2     2                  2          2
   - 24 L alpha  a[0]  a[2] - 24 L alpha  a[1] a[2] 

      2       2           2                   2     
   - k  lambda  a[2] + 2 k  mu a[2] + 6 L a[0]  a[2]

                  2          2                  2     
   + 6 L a[1] a[2]  - 4 alpha  beta a[2] + alpha  a[2]

                           2  2  2                 2     3
   + beta a[2] = 0, 8 alpha  k  v  a[1] - 8 L alpha  a[1] 

             2  2                 2  2           2  2     
   - 16 alpha  k  v a[1] + 8 alpha  k  a[1] - 2 k  v  a[1]

             3      2             2              3    \   
   + 2 L a[1]  + 4 k  v a[1] - 2 k  a[1] - 2 a[1]  = 0 }, 
                                                      /   

                           \
  {alpha, a[0], a[1], a[2]}|
                           /

I've looked up previous answers to this but __ double underscore only works in math input, and ctrl + shift + ^ or ctrl + shift + _ does not seem to work. I'm using Maple 18, windows computer.

Hi

I have a first oder PDE, I use pdsolve I obtained a solution depend on function F

condition_unique_solution.mw

 

My question: The boundary condition  f(x,y) = 1 is supplied on the line y = k x, where k is a constant. For which k
does there exist a unique solution for f(x, y)?

 

Many thanks for your help

 

Hi 

I solve the laplace equation written in polar coordinates in annular domain.
The code run without any error 

But there is no solution displayed after running the code, note that I use Maple 18

Laplace_annulardomain.mw

Many thinks for your help

i want to design a packaging container to hold 320 sphere-shaped chocolates that each has a diameter 1.8 cm and weights about 3.2g each. i hope can get all posible shape using maple18 .

Hi

I try to solve the laplace equation with some special boundary conditions.

But, i get the follwoing error

Error, (in pdsolve/sys) the given system is not polynomial in the variables {f}

 

 laplace_equation.mw

Thank you for any help

 


Hi, I am struggling a bit animating a result that I have.

Let me explain.

I solved numerically some equations for me movement of a rolling disc and a bar bolted to the center of the disc.

So I have a long list of values for:

x[h] for the horizontal position of the disc
beta[h] for the rotation of the disc
theta[h] for the rotation of the bar

and so on...

Then I created a procedure that, for a given instant "h", maple plots the disc and the bar correctly.

Now I want to anumate this procedure that should depend only on the instant h and I cannot perform it.
Maybe the error comes from the fact that there are not functions involved but numerical data, or maybe it is the fact that disc radius and bar lenght are assigned outside the procedure.

Any help would be appreciated.

This is the file:

Mechanics.mw
 

I am using Maple 18 to plot iterates of a discrete map.  A vector R1 is parametrized by the label b that runs from 0 - 2*Pi.  The following initial iterate gives me a blue circle:  

>a := 1.0
>q := a*cos(b)
>p := a*sin(b)
>R1 := PositionVector([q, p], cartesian[x, y])
>PlotPositionVector(R1, b = 0 .. 2*Pi);

Is there a way for me to assign a gradient colorscheme to each point of the circle, using b as parameter from 0 - 2*Pi? I have tried using the following, which does not give an error message, but it also just gives me a blue circle:

>PlotPositionVector(R1, b = 0 .. 2*Pi, curveoptions = [colorscheme = ["valuesplit"]]);

I am not sure how to use the vector and mapping options in the command >plotcommand(plotargs, colorscheme=["valuesplit", V, mapping]) to obtain the color gradient along the circle.

How can I draw together the function f(x,y)=(x.y^2)+1 and its tangent plane at the point (2,1,f(1,2)) in each of the following rectangles? 

(Note that the equation of the tangent plane at the given point is z=x+4y-3)

a. 2<=x<=3 and 1<=y<=3

b. 2<=x<=2.1 and 1<=y<=1.2

c. 2<=x<=2.01 and 1<=y<=1.02

Sir please help me to solve following system of non linear equation, 

x+y=1

y=x^2-5

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