Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

what is problem in here

restart;

with(plots):

with(LinearAlgebra):

with(DEtools):

diff(u(x), x) = f(u(x), v(x)), diff(v(x), x) = g(u(x), v(x)) for the two differential equations.

f := (u,v) -> u+cos(v);
g := (u,v) -> u*v-v+sin(v);

proc (u, v) options operator, arrow; u+cos(v) end proc

 

proc (u, v) options operator, arrow; v*u-v+sin(v) end proc

(1)

The equilibria:

equilibria := solve({f(u,v)=0, g(u,v)=0}, {u,v},explicit);

{u = 1, v = Pi}, {u = -cos(RootOf(sin(_Z)*_Z^2+sin(_Z)-2*_Z)), v = RootOf(sin(_Z)*_Z^2+sin(_Z)-2*_Z)}

(2)

allvalues(RootOf(sin(_Z)*_Z^2+sin(_Z)-2*_Z))

RootOf(sin(_Z)*_Z^2+sin(_Z)-2*_Z, 1.306542374), RootOf(sin(_Z)*_Z^2+sin(_Z)-2*_Z, -1.306542374), RootOf(sin(_Z)*_Z^2+sin(_Z)-2*_Z, -2.331122370), RootOf(sin(_Z)*_Z^2+sin(_Z)-2*_Z, 2.331122370), 0

(3)
 

NULL

Download remove.mw

when we do the linearization we do transfer the critical point to the orgin, i want see the plot of nonlinear system before the linearization and after linearization which is critical point is tranfered to orgin, also i try to find `conserved quantity` and make a good plot with the critical point shown in graph how i can do the linearization by maple, also  we  by hand when we solved we show that in new system when we give him variable we have to write for each new system  u=x-x[0] when is u is new system for linearization we need that or not? becuase result is not change when we take jacobian result is same !

 

Linearization.mw

i want plot a system of differential equation and do phaseportrait i did but when i want make it a little bit more clear and colorfull like rainbow when i find the C.Q i don't know how set the option for ploting?

e1.mw

i have a system but i don't know how i can plot chaotic or is need time for ploting or not ?

Has anyone else run into the following bug in the PolynomialTools package?

> PolynomialTools[PolynomialToPDE]([z*f], [z], [f]);

                            d
                           [-- f(z)]
                            dz

> PolynomialTools[PolynomialToPDE]([z[1]*f], [z[1]], [f]);

                           [f(z[1])]

I believe the first behavior is correct, but when I use subscripted variables, I get the second behavior for linear terms.  (Nonlinear terms appear to work correctly:


>  PolynomialTools[PolynomialToPDE]([z[1]^2*f], [z[1]], [f]);
                                2
                               d
                             [------ f(z[1])]
                                  2
                              dz[1]

I have observed this in Maple 2024.2 and previous versions.  Does anyone know if it's fixed in Maple 2025?  (I don't have access yet.)  Does anyone know a good workaround for polynomials with both linear and nonlinear terms?

Thanks!

----Josh

On my journey of discovery through the world of Maple, I would like to ask for help again:
How are the properties of variables and the indexing of sequences handled/determined? For this purpose, I chose an old problem from a challenging MO as an exercise and tried it in the attached file, which, of course, failed.

restart

a[i] := 2*cos(t/2^i)-1

2*cos(t/2^i)-1

(1)

b[n] := product(a[i], i = 1 .. n)

product(2*cos(t/2^i)-1, i = 1 .. n)

(2)

limit(b[n], n = infinity)

limit(product(2*cos(t/2^i)-1, i = 1 .. n), n = infinity)

(3)

simplify(limit(b[n], n = infinity))

limit(product(2*cos(t*2^(-i))-1, i = 1 .. n), n = infinity)

(4)

NULL

Download test.mw

We are excited to announce that the Maple Conference will be held Novemeber 5-7, 2025!

Please join us at this free virtual event as it will be an excellent opportunity to meet other members of the Maple community, get the latest news about our products, and hear from the experts about the challenges and opportunities that our technology brings to teaching, learning, and research. More importantly, it's a chance for you to share the work you've been doing with Maple and Maple Learn. 

The Call for Participation is now open. We are inviting submissions of presentation proposals on a range of topics related to Maple, including Maple in education, algorithms and software, and applications. We also encourage submission of proposals related to Maple Learn. 

You can find more information about the themes of the conference and how to submit a presentation proposal at the Call for Participation page. Applications are due July 25, 2025.

After the conference, all accepted presenters and invited speakers will be invited to submit related content to the Maple Transactions journal for consideration.

Registration for attending the conference will open in July.  Watch for further announcements in the coming weeks.

We hope all of you in the Maple Primes community will join us for this event!

Kaska Kowalska
Contributed Program Co-Chair

Regarding to this equation i want to do the same as paper did, but i am not sure how i can determine the both  point in plot and how thus point show effect of system of equation and  how find jacobian in the two different point  i think he find it in general but i think we have to find it for each one  of equalibriom point  or we just find in general .
thus plot also are emazing did he use special code in matlab or mathematica? there is some app for plot thus kind of phase portrait but they not like this this is must have a special code for plotting in matlab or mathematica

restart

with(PDEtools)

undeclare(prime, quiet)

with(LinearAlgebra)

declare(u(x, t), quiet); declare(U(xi), quiet); declare(V(xi), quiet)

NULL

DU := V

V

(1)

DV := U^2*beta/(delta*k^2)+(alpha*w/(delta*k^3)+1/(delta*k^2)-w^2/(delta*k^4))*U

U^2*beta/(delta*k^2)+(alpha*w/(delta*k^3)+1/(delta*k^2)-w^2/(delta*k^4))*U

(2)

var := U, V

U, V

(3)

J := Student:-MultivariateCalculus:-Jacobian([DU, DV], [var], 'output' = 'matrix')

Matrix(%id = 36893490902008029780)

(4)

Eigenvalues(J, implicit, output = 'list')

[RootOf(_Z^2*delta*k^4-2*U*beta*k^2-alpha*k*w-k^2+w^2, index = 1), RootOf(_Z^2*delta*k^4-2*U*beta*k^2-alpha*k*w-k^2+w^2, index = 2)]

(5)

Eigenvalues(J)

Vector[column](%id = 36893490902007993756)

(6)

NULL

NULL

DU1 := V

V

(7)

V = solve(DU1, V)

V = 0

(8)

DV1 := U^2*F[2]+U*F[1]

U^2*F[2]+U*F[1]

(9)

U = solve(DV1, U)

U = (0, -F[1]/F[2])

(10)

so*we*have*three*equalibriom*point*(0, 0)

so*we*have*three*equalibriom*point*(0, 0)

(11)

J1 := Student:-MultivariateCalculus:-Jacobian([DU1, DV1], [var], 'output' = 'matrix')

Matrix(%id = 36893490902007256836)

(12)

Eigenvalues(J1, implicit, output = 'list')

[RootOf(-2*U*F[2]+_Z^2-F[1], index = 1), RootOf(-2*U*F[2]+_Z^2-F[1], index = 2)]

(13)

Eigenvalues(J1)

Vector[column](%id = 36893490901938203092)

(14)
 

NULL

Download bi-1.mw

The attached problem is from a 1988 MO. It can be solved using complete induction, paper, and pencil, and with some effort, yields a simple answer. It's quite challenging to do by hand, but with "derive", it only takes three lines and a fraction of a second. Mow test.mw

restart

"puzzle(n):=(∑)(cos(k*Pi/(2*n+1))^())^(4);"

Download test.mw

I can't do it with Maple because I'm doing something wrong again. Therefore, I'm asking for help.

Recently, I was surprised by the fact that it is not possible to create a table, where the values of entries are tables.

Toy example:

T1:=table(["1"=table(["first"=12, "second"=15]), "2"= table(["first"=9, "second"=7])]);

The code

T1["1"];

doesn't return table(["first"=12, "second"=15])

If table has lists in right-hand side

T1:=table(["1"=[12, 15], "2"= [9, 7]]);

all works properly.

 Based on the type of coefficients in the linear equations, SolveTools[Linear] provide several method including method = Rational, Polynomial etc. 
The Polynomial method of SolveTools[Linear] however, cannot be directly called by SolveTools[Linear]. The interpreter complains that no such methd called "Polynomial".

This could only be solved if I import the SolveTools package in advance and call Linear.

I hope Maple could solve this issue. 

I'm trying to reproduce a manual asymptotic analysis (see the attached pdf file) in Maple for a two-soliton solution. Specifically, I want to evaluate the limit of a function (e.g., r[2]r[2]r[2] or ∂xq[2]\partial_x q[2]∂x​q[2]). How can I properly perform the limit a2→+−∞ as t​→+−∞ in Maple, either by substitution or by reparametrization, in order to study the asymptotic behavior of a multi-variable expression symbolically? 

restart

with(Student[Calculus1])

lambda1 := I*mu1; lambda2 := I*mu2; a1 := -2*x/mu1+mu1*t; a2 := -2*x/mu2+mu2*t

numer_r := lambda2*cosh(a1)-lambda1*cosh(a2)

denom_r := (lambda1^2+lambda2^2)*cosh(a1)*cosh(a2)-2*lambda1*lambda2*(1+sinh(a1)*sinh(a2))

r2 := I*(-lambda1^2+lambda2^2)*numer_r/denom_r

numer_dq := (sinh(a1)-sinh(a2))^2; denom_dq := ((lambda1^2+lambda2^2)*cosh(a1)*cosh(a2)-2*lambda1*lambda2*(1+sinh(a1)*sinh(a2)))^2

dq2 := 1-2*(lambda1^2-lambda2^2)^2*numer_dq/denom_dq

limit_r2 := limit(r2, a2 = -infinity); simplify(limit_r2)

Error, invalid input: limit expects its 2nd argument, p, to be of type Or(name = algebraic,set(name = algebraic),list(name = algebraic)), but received -2*x/mu2+mu2*t = -infinity

 

limit_r2

(1)

limit_r2_pos := limit(r2, a2 = infinity); simplify(limit_r2_pos)

Error, invalid input: limit expects its 2nd argument, p, to be of type Or(name = algebraic,set(name = algebraic),list(name = algebraic)), but received -2*x/mu2+mu2*t = infinity

 

limit_r2_pos

(2)

NULL

Download asymptotic.mw cooocp_(2).pdf

Any suggestions for reformulating the limit or change of variables would be appreciated. 

I have solution of this problem but how I can plot simulation on Maple?Thanks, In attach my problem

Maple 2025.1

We have just released an update to Maple. Maple 2025.1 includes several enhancements to the new interface, as well as various small corrections throughout the product. As always, we recommend that all Maple 2025 users install this update.

In particular, please note that this update includes a fix to the problem where new documents were opening in a new window instead of a new tab.  Thanks for helping us, and other users, by letting us know!

This update is available through Tools>Check for Updates in Maple, and is also available from the Maple 2025.1 download page on web site, where you can also find more details.

MapleSim 2025

We are happy to announce that we just released MapleSim 2025. This release includes a new component library to support the modeling of motor drives and updates to several in-product apps that make it even easier to perform optimization and analysis.

See What’s New in MapleSim for details.

Good day, everyone

I am trying to code HPM but it's giving me the error code "Error, invalid subscript selector
" once I increase the number of iterations above 2. 

Attached below is the code. 

HPM.mw

Thanks 

4 5 6 7 8 9 10 Last Page 6 of 2217