Maple 2024 Questions and Posts

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

in this equation i solved without polar which is not shown well without polar specially when a=0 without polar if we solve this  and the other point are complex which not shown like a point in plot, except point (0,0) other not fine, but i want to solve by polar for finding  limit cycle the polar way are so usefull and also for some nonlinear system how i can  change the system to polar and find r* and theta* and plot them ?

polar.mw

i want to plot the function appear in P but there is a variable t which appear from variable which my system is depending on u,v and must v appear in my function instead of exponential so exponential contain the t variable  there is any way for fixing this?

restart;

with(plots):

 

f := (u,v) -> -u+u^3;
g := (u,v) -> -2*v;

proc (u, v) options operator, arrow; -u+u^3 end proc

 

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

(1)

 

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

{u = 0, v = 0}, {u = 1, v = 0}, {u = -1, v = 0}

(2)

 

de1 := diff(u(t),t) = f(u(t),v(t));
de2 := diff(v(t),t) = g(u(t),v(t));

diff(u(t), t) = -u(t)+u(t)^3

 

diff(v(t), t) = -2*v(t)

(3)
   

PDEtools:-ConservedCurrents({de1, de2}, [u(t), v(t)]); P1 := -(1/2)*op(1, rhs(op(%)))

[_J[t](t, u(t), v(t)) = f__1((-u(t)^2+1)*exp(-2*t)/u(t)^2, v(t)*exp(2*t))]

 

-(1/2)*(-u(t)^2+1)*exp(-2*t)/u(t)^2

(4)

P := -(-u^2+1)*exp(-2*t)/(2*u^2)

-(1/2)*(-u^2+1)*exp(-2*t)/u^2

(5)

 

equilibria;

{u = 0, v = 0}, {u = 1, v = 0}, {u = -1, v = 0}

(6)

 

p1 := contourplot(P, u=-1.5..2, v=-1.5..1.5, scaling=constrained,
    colorscheme="DivergeRainbow", contours=[seq](x, x=-0.4..0.4,0.1)):

Error, (in plot/iplot2d:-Levels) could not evaluate expression

 

 

p2 := pointplot([[0,0],[1,0],[-1,0]], symbol=solidcircle, symbolsize=15, color=black):

 

p3 := fieldplot([f(u,v), g(u,v)], u=-1.5..2, v=-1.5..1.5,
        arrows=medium, fieldstrength=fixed(0.4), grid=[10,10], labels=["",""]):

 

display(p1,p2,p3, axes=normal);

 

 

Download p1-1.mw

it is apear by sign of taking derivative from equation and substitute in it the sign less zero is stable otherwise is unstable but i want to shown by graph like shown in graph

plot-stablity-1equation.mw

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

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.

i need a way for making my loop working fast, i need something for calculating such term fast without showing result but thus term make a lot of time for coming out ?

fasting.mw

Hi,

I’m trying to remove the parentheses (for the dᵢ values) to achieve optimal display. Ideas ? Thanks

S4_Droites_Implicite.mw

i want to do export  the table to latex there is command to do?

table.mw

in my iteration when i got answer the error is so big and the method is powerfull i don't know i did mistake or my exact function i have to change it, becuase error is so big , how i can make error be smaller i did any mistake ?
and how i can plot them together exact and approximate in one graph also 3D and 2D

Num-e1.mw

5 6 7 8 9 10 11 Last Page 7 of 48