Maple 2020 Questions and Posts

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

I'm running Maple 2020 on a Mac with OS 14.3.1.  Maple is not seeing current files in folders, and when it saves files, it does so with an older date.  In Settings, Maple has permission to access the Desktop, Documents, and Downloads folders, the only options that my Mac provides.

It's a bit frustrating.  Any help would be appreciated.

Hey!

I am working on a 2-dimensional real curve gamma, that is being deformed by a real 2x2 matrix A and translated by some shift vector s. Both the matrix A and the vector s depend on two real parameters alpha and beta, which vary between 0 and 1. My goal is to create an animated plot that shows the curve depending on the parameters alpha and beta (ideally with some sort of slider, so I can play with different values of alpha and beta). I am having a hard time animating this, also because I am unsure if this is actually possible in Maple. The code I am working with is so far:

with(plots);
with(LinearAlgebra);
f := x -> exp(2*I*Pi*x);
CDFT := Matrix(3, 3, [[1, 1, 1], [1, f(1/3), f(2/3)], [1, f(2/3), f(1/3)]]);

al := 10/100;
be := 10/100;
lam1 := f(al);
lam2 := f(be);
lam3 := f(-al - be);
coef := (1/3*HermitianTranspose(CDFT)) . (Vector[column](3, [lam3, lam2, lam1]));
a := coef[1];
b := coef[2];
c := coef[3];

gam := t -> MatrixVectorMultiply(Matrix(2, 2, [[Re(a^2 + b*c), -Im(a^2 - b*c)], [Im(a^2 + b*c), Re(a^2 - b*c)]]), Vector[column](2, [cos(2*t) - 2*cos(t), sin(2*t) + 2*sin(t)])) + 3*Vector[column](2, [Re(b*c), Im(b*c)]);
P1 := plot([cos(2*t) - 2*cos(t), sin(2*t) + 2*sin(t), t = 0 .. 2*Pi], color = [blue]);
P2 := plot([gam(t)[1], gam(t)[2], t = 0 .. 2*Pi], color = [purple]);
plots:-display([P1, P2]);

Note that you can basically ignore everything up to the definition of the curve gam (short for gamma), apart from the definition of the parameters al and be (alpha and beta). The plot P1 corresponds to the "unperturbed curve", i.e. when we multiply the curve by the identity, which happens for alpha=beta=0. The plot P2 is now the deformed curve. The output is then:

My goal is now to animate this plot such that I can play with different values of alpha and beta, without having to manually insert them. How do I do this?

This worksheet animates the motion of an object (say, a cube which slides frictionlessly) on a rotating carousel. The cube is not self-propelled.

How can the worksheet be modified to handle the combination of the carousel motivated motion and the cube's own generated motion, caused by, say, by a few strategically placed thrusters?

The cube's own generated path could be a straight line, or a curve such as an ellipse. The cube's own motion could have a constant velocity or be accelerating.

Carousel_dynamics.mw

Maple seems to give me very inaccurate results for this computation, I'm wondering if this is a known issue and if there's a way to fix it? Worksheet:

Exact commands I ran:
with(combinat):
B := (n, i, p) -> binomial(n, i)*(p^i)*(1-p)^(n-i)/i;

F:=(n,p)->sum(B(n,i,p),i=1..n);
F(2,1);

This outputs 0, when it should output 1/2. See image attached:

 

Is there a fix for this?

Hello,

How to add the arrow or symbol for the flow function ( psi[1] and psi[2]). I do not know how to use from fieldplot with the option of arrows = SLIM

please see the attached figure.

Thanks so much

h1_h2_ok_ac_2020_8_ac_(1).mw

this is my model. Please give me how to find the DFE and basic reproduction number from maple.

please help me to fixed and find the analytic solution...

restart

with(linalg)

f1 := mu[1]*N-delta*r-beta*N*s(e+i)

mu[1]*N-delta*r-beta*N*s(e+i)

(1)

f2 := beta*N*s(e+i)+omega*v-(mu[1]+pi+gamma)*e

beta*N*s(e+i)+omega*v-(mu[1]+pi+gamma)*e

(2)

``

f3 := pi*e-(mu[1]+mu[2]+sigma)*i

pi*e-(mu[1]+mu[2]+sigma)*i

(3)

f4 := gamma*e+sigma*i-(theta+mu[1]+mu[2])*q

gamma*e+sigma*i-(theta+mu[1]+mu[2])*q

(4)

f5 := theta*q+xi*v-(mu[1]+delta)*r

theta*q+xi*v-(mu[1]+delta)*r

(5)

f6 := eta*N*s-(mu[1]+omega+xi)*v

eta*N*s-(mu[1]+omega+xi)*v

(6)

T := solve({f1, f2, f3, f4, f5, f6}, [s, e, i, q, r, v])

Error, (in solve) cannot solve for an unknown function with other operations in its arguments

 

``

Download cobacoba.mw

How I can solve this, because I can't find the solutions. The display show "Length of output exceeds limit of 1000000"

Please help me.. 

This is my model on picture

The command

plot(x^(1/3), x = -10 .. 10)

plots only the points where x>=0, as for negative values of x the cubic root retuns only the complex root.

Ca can I modify it so that it returns only the real root, so that the whole plot can be viewed?

For a project I need to construct a large symbolic adjoint matrix, hoping it can be factored afterward into nice expressions.
In the worksheet, I present an adjoint matrix using permuted Hadamard products. What puzzles me is that only for even dimensions, I need to multiply it (elementwise) with a parity matrix. Okay, not a specific Maple question, but maybe someone can help me out.

Download Adjoint.mw

How to fix the error?

How are they (tanh(a+b) , tanh(a-b)) defined in Maple?

Ger.mw

How I can solve a PDE on two regions with matching conditions at the common boundary?  

T1.mw

The uploaded worksheet begins to uniformly tile the Poincare disk with pentagons using hyperbolic reflection .

Although relatively easy to create the central pentagon and the first adjacent pentagon, it becomes increasingly difficult to determine which lines to reflect to create the remaining pentagons in the first tier adjacent to the central pentagon and more so to create the pentagons of the second tier adjacent to those in the first tier and so on.

Is there a better technique for accomplishing this?

In particular can Mobius tranformations be employed to do this? If so, please replay with or point to a working example of this for me to follow.

 Tile_Poincare_disk.mw

Sorry, I forgot that respondents to this question must establich their own link to the DirectSearch package.

Hello
I have a problem in writing the Maple code of the image below, I don't know why the 3.5 answers are not available?

which one is better?

123.mw

0123.mw

Why aren't all the variables in fin 1 equation?

And the answers are different from the solutions?

 

restart

with(student)

eq1 := 12*gamma^3*rho[3]^2*(diff(w(psi), `$`(psi, 2)))+(-3*gamma*rho[2]^2+4*omega*rho[3]^2)*w(psi)+gamma*rho[3]^2*(rho[1]+2*rho[3])*w(psi)^3

12*gamma^3*rho[3]^2*(diff(diff(w(psi), psi), psi))+(-3*gamma*rho[2]^2+4*omega*rho[3]^2)*w(psi)+gamma*rho[3]^2*(rho[1]+2*rho[3])*w(psi)^3

(1)

NULL

"w(psi):=kappa[0]+sum(kappa[i]*((diff(E,psi))^(i))/((E(psi))^(i)),i=1..1)+sum(h[i]*(((diff(E,psi))^())/((E(psi))^()))^(-i),i=1..1);"

proc (psi) options operator, arrow, function_assign; kappa[0]+sum(kappa[i]*(diff(E, psi))^i/E(psi)^i, i = 1 .. 1)+sum(h[i]*((diff(E, psi))/E(psi))^(-i), i = 1 .. 1) end proc

(2)

"E(psi):=((epsilon[1]*jacobiCN(Zeta[1]*psi))+(epsilon[2]*jacobiSN(Zeta[2]*psi)))/((epsilon[3]*jacobiCN(Zeta[3]*psi))+(epsilon[4]*jacobiSN(Zeta[4]*psi))) ;"

proc (psi) options operator, arrow, function_assign; (varepsilon[1]*jacobiCN(Zeta[1]*psi)+varepsilon[2]*jacobiSN(Zeta[2]*psi))/(varepsilon[3]*jacobiCN(Zeta[3]*psi)+varepsilon[4]*jacobiSN(Zeta[4]*psi)) end proc

(3)

 

NULL

fin1 := simplify(eq1)

kappa[0]*(gamma*rho[3]^2*(rho[1]+2*rho[3])*kappa[0]^2-3*gamma*rho[2]^2+4*omega*rho[3]^2)

(4)

Sol := solve(fin1, {omega, Zeta[1], Zeta[2], Zeta[3], Zeta[4], epsilon[1], epsilon[2], epsilon[3], epsilon[4], h[1], kappa[0], kappa[1]})

{omega = omega, Zeta[1] = Zeta[1], Zeta[2] = Zeta[2], Zeta[3] = Zeta[3], Zeta[4] = Zeta[4], h[1] = h[1], kappa[0] = 0, kappa[1] = kappa[1], varepsilon[1] = varepsilon[1], varepsilon[2] = varepsilon[2], varepsilon[3] = varepsilon[3], varepsilon[4] = varepsilon[4]}, {omega = -(1/4)*gamma*(kappa[0]^2*rho[1]*rho[3]^2+2*kappa[0]^2*rho[3]^3-3*rho[2]^2)/rho[3]^2, Zeta[1] = Zeta[1], Zeta[2] = Zeta[2], Zeta[3] = Zeta[3], Zeta[4] = Zeta[4], h[1] = h[1], kappa[0] = kappa[0], kappa[1] = kappa[1], varepsilon[1] = varepsilon[1], varepsilon[2] = varepsilon[2], varepsilon[3] = varepsilon[3], varepsilon[4] = varepsilon[4]}

(5)

for i to 2 do Case[i] := allvalues(Sol[i]) end do

{omega = omega, Zeta[1] = Zeta[1], Zeta[2] = Zeta[2], Zeta[3] = Zeta[3], Zeta[4] = Zeta[4], h[1] = h[1], kappa[0] = 0, kappa[1] = kappa[1], varepsilon[1] = varepsilon[1], varepsilon[2] = varepsilon[2], varepsilon[3] = varepsilon[3], varepsilon[4] = varepsilon[4]}

 

{omega = -(1/4)*gamma*(kappa[0]^2*rho[1]*rho[3]^2+2*kappa[0]^2*rho[3]^3-3*rho[2]^2)/rho[3]^2, Zeta[1] = Zeta[1], Zeta[2] = Zeta[2], Zeta[3] = Zeta[3], Zeta[4] = Zeta[4], h[1] = h[1], kappa[0] = kappa[0], kappa[1] = kappa[1], varepsilon[1] = varepsilon[1], varepsilon[2] = varepsilon[2], varepsilon[3] = varepsilon[3], varepsilon[4] = varepsilon[4]}

(6)

NULL

NULL

Download 0123.mw

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