Maple 18 Questions and Posts

how to choose the colors of the "filled" option...

Asked by: Paulo Baumbach 50

June 07 2019

1 5

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

How do I use Homotopy Perturbation Method to find ...

Asked by: kennywhyte 0

May 23 2019

0 0

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.

fix simplify problem...

Asked by: bashar27 55

May 22 2019

0 2

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

fix the problem pls...

Asked by: bashar27 55

May 22 2019

0 1

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);

How to not evaluate worksheet until i want...

Asked by: federicie 10

May 22 2019

2 3

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

fix the solve problem pls...

Asked by: bashar27 55

May 21 2019

0 1

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]}|
/

How do I write a subscript/superscript into text s...

Asked by: John Morrell 10

May 05 2019

1 6

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.

Unique solution using pdsolve...

Asked by: Zeineb 520

May 03 2019

0 0

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

Maple 18: laplace in annular domain...

Asked by: Zeineb 520

May 02 2019

0 6

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

How to write codes using maple to stimulate all de...

Asked by: Awesome123 15

May 02 2019

1 3

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 .

Laplace equation pdsolve...

Asked by: Zeineb 520

May 01 2019

1 11

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

Animate Numerical set of data...

Asked by: Don_Caraota 150

April 20 2019

1 4

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

How do I plot each point of a parametrized curve w...

Asked by: sbleher 5

April 19 2019

1 3

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.

Drawing together a function of two variable and it...

Asked by: Vahid_math 15

April 17 2019

0 1

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

How do I solve system of non linear equations? ...

Asked by: umairaslam 15

April 15 2019

2 1

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

x+y=1

y=x^2-5

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

Ask a Question

Create a Post