Maple 18 Questions and Posts

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

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/

 

 

Hello;

 I am trying to verify the analytic solution of a electric and magnetic fields created by a small dipole antenna (also called "Hertzian dipole"). The study of a small dipole is ground zero of anyone learning about antennas as calculations are "relatively" easy if a mathematical software is used. As the title suggests, the fieldplot3d returns an empty box.

Here is some introduction for the problem in question:

 The procedure is relatively straightforward, first, current density vector is defined, from there, magnetic vector potential (named "vector A") is calculated. The curl of vector A gives the magnetic field (named "B-field") produced by the antenna. From B-field, H field is deduced as it is only a multiplication of the B-field by a constant.

 At this point, I try to plot the H-field, and it works like a charm. No problem at all.

The electric field (named "E-field") then, may be calculated by taking the curl of the H-field and multiplying by a constant.

 At this second point, I try to plot the E-field, however, Maple returns an empty box. 

First, I thaught, maybe it was a problem of division by 0, however, after redefining the axis ranges, the problem still persists. I am attaching the code and the images. Any help will be greatly appreciated.

PS: This is my first post and I am very new to maple, please indulgde me if I make some formatting and/or post mistakes

 

KB

First step: Verify calculations for Hertzian dipole

 

restart;

with(plots):

with(LinearAlgebra):

with(VectorCalculus):

#IMPORTANT: R is constant for the calculation of A

 

 

 

J:=Vector[column]([ 0 ,
                 0 ,
                 I_0/s ]);

J := I_0*e[z]/s

(1)

 

A := VectorCalculus:-`*`(VectorCalculus:-`*`(mu_0, 1/VectorCalculus:-`*`(4, Pi)), int(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(J, exp(VectorCalculus:-`-`(VectorCalculus:-`*`(VectorCalculus:-`*`(I, k), R)))), s), 1/R), z = VectorCalculus:-`-`(VectorCalculus:-`*`(l, 1/2)) .. VectorCalculus:-`*`(l, 1/2)))

A := (1/4)*mu_0*exp(-I*k*R)*I_0*l*e[z]/(Pi*R)

(2)

A[1];

0

(3)

A[2];

0

(4)

A[3];

(1/4)*mu_0*exp(-I*k*R)*I_0*l/(Pi*R)

(5)

#Taking the curl of A:

#IMPORTANT: R is a function of x,y,z:

R:=sqrt(x^2+y^2+z^2);

(x^2+y^2+z^2)^(1/2)

(6)

 

 

 

#Defining B by taking the curl

B[1] := VectorCalculus:-`+`(diff(A[3], y), VectorCalculus:-`-`(diff(A[2], z))):

 

 

B[2] := VectorCalculus:-`-`(VectorCalculus:-`+`(diff(A[3], x), VectorCalculus:-`-`(diff(A[1], z)))):

 

 

B[3] := VectorCalculus:-`+`(diff(A[2], x), VectorCalculus:-`-`(diff(A[1], y))):

B:=Vector[column]([ B[1] ,
                 B[2] ,
                 B[3] ]):

 

 

mu_0:=4*Pi*10^(-7):

I_0:=2400:

f:=2500:

omega:=2*Pi*f:

c:=3*10^8:

k:=omega/c:

l:=3*10^(-2):

epsilon_0:=1/(mu_0*c^2):

 

 

 

B_plot:=fieldplot3d([B[1],B[2],B[3]], x=-1..1,y=-1..1,z=-1..1,fieldstrength=log,arrows=SLIM):

 

 

H:=(1/mu_0)*B:

 

H_plot:=fieldplot3d([H[1],H[2],H[3]], x=-1..1,y=-1..1,z=-1..1,fieldstrength=log,arrows=SLIM):

 

# Taking curl of H to find the E field:

 

E[1] := VectorCalculus:-`*`(1/VectorCalculus:-`*`(VectorCalculus:-`*`(I, omega), epsilon_0), VectorCalculus:-`+`(diff(H[3], y), VectorCalculus:-`-`(diff(H[2], z)))):

E[2] := VectorCalculus:-`*`(1/VectorCalculus:-`*`(VectorCalculus:-`*`(I, omega), epsilon_0), VectorCalculus:-`-`(VectorCalculus:-`+`(diff(H[3], x), VectorCalculus:-`-`(diff(H[1], z))))):

E[3] := VectorCalculus:-`*`(1/VectorCalculus:-`*`(VectorCalculus:-`*`(I, omega), epsilon_0), VectorCalculus:-`+`(diff(H[2], x), VectorCalculus:-`-`(diff(H[1], y)))):

E:=Vector[column]([ E[1] ,
                 E[2] ,
                 E[3] ]):

E_plot:=fieldplot3d([E[1],E[2],E[3]], x=1..500,y=1..500,z=1..500,fieldstrength=log,arrows=SLIM):

 

 

subs(x=1,y=1,z=1,H):

 

H_plot;

 

 

 

E_plot;

 

 

``


 

Download short_dipole_matrix_way_old_school.mw 

 

Foucault’s Pendulum Exploration Using MAPLE18

https://www.ias.ac.in/describe/article/reso/024/06/0653-0659

In this article, we develop the traditional differential equation for Foucault’s pendulum from physical situation and solve it from
standard form. The sublimation of boundary condition eliminates the constants and choice of the local parameters (latitude, pendulum specifications) offers an equation that can be used for a plot followed by animation using MAPLE. The fundamental conceptual components involved in preparing differential equation viz; (i) rotating coordinate system, (ii) rotation of the plane of oscillation and its dependence on the latitude, (iii) effective gravity with latitude, etc., are discussed in detail. The accurate calculations offer quantities up to the sixth decimal point which are used for plotting and animation. This study offers a hands-on experience. Present article offers a know-how to devise a Foucault’s pendulum just by plugging in the latitude of reader’s choice. Students can develop a miniature working model/project of the pendulum.

Hi! my name is Euge, i live in Argentina!
I´m solving a problem that ask me to find an ODE from a family of curves, 
I know this is pretty simple, I have to derive the equation two times, but i can´t find a comand in maple that helps me do that. Is there any?
Thank you! 

hi, I am trying to evaluate the following integrations (Note that; the values of alpha, lambda_1 and lambda_2) are random and they change with every iteration 

`λ1_B`[so] := 10.94:

``

W[1] := evalf[4](int(`α_B`[so]^2/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])), z = 1 .. infinity)); W[2] := evalf[4](int(z^(`λ2_B`[so]/`λ1_B`[so])*ln(z)/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^2), z = 1 .. infinity)); W[3] := evalf[4](int(1/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^2), z = 1 .. infinity)); W[4] := evalf[4](int(z^(`λ2_B`[so]/`λ1_B`[so])*ln(z)^2/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^2), z = 1 .. infinity)); W[5] := evalf[4](int(1/((z-1+`α_B`[so])^3*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])), z = 1 .. infinity)); W[6] := evalf[4](int((z^(`λ2_B`[so]/`λ1_B`[so]))^2*ln(z)^2/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^3), z = 1 .. infinity)); W[7] := evalf[4](int(z^(`λ2_B`[so]/`λ1_B`[so])*ln(z)/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^3), z = 1 .. infinity)); W[8] := evalf[4](int(1/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^3), z = 1 .. infinity)); W[9] := evalf[4](int(z^(`λ2_B`[so]/`λ1_B`[so])*ln(z)/((z-1+`α_B`[so])^3*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^2), z = 1 .. infinity)); W[10] := evalf[4](int(1/((z-1+`α_B`[so])^4*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])), z = 1 .. infinity)); W[11] := evalf[4](int(1/((z-1+`α_B`[so])^3*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^2), z = 1 .. infinity))

int(79.977249/((z+7.943)^2*(z^.1408592322+7.943)), z = 1 .. infinity)

 

Warning,  computation interrupted

 

Warning,  computation interrupted

 

0.1744e-1

 

Warning,  computation interrupted

 

0.3348e-2

 

Warning,  computation interrupted

 

Warning,  computation interrupted

 

Warning,  computation interrupted

 

Warning,  computation interrupted

 

Warning,  computation interrupted

 

``

``


 

Download ask_maple.mw

 

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.

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