Maple 2024 Questions and Posts

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

i need at least 10 iteration but i can go 4 of them untill then need a lot of time to waiting i saw the number are so big in fractional manner i want to reduce them and do something which be shorter there is any way for doing that?

ite1.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

i got some issue regarding to solving numerical example, he got combine function of infinite term i don't know how when he give alpha a number ? and when didn't give number  to order of fractional laplace and stay alpha letter he got outcome i want apply laplace while is for general alpha betwen n-1<alpha<n when n is number of initial condition, i want to get same result as he did steps are not wrong but outcome have problem .

i have same problem in three example in one paper, is related to update maple to 2025, i saw talking about calculating inverse laplace at new version  any one  know how fix this?

Ex1.mw

EX2.mw

Ex3.mw

in thus example all of them are write in shape of matrix but all of them are linear differential equation which have critical point zero ,i want to plot thus example and decided the kind of critical point by eagenvalue and by eagenvector i can find the trajectory how i can plot by matrix of each example and show that that critical point is which type as mention in picture?

system-phase-examples.mw

in my iteration is appear Error, (in assuming) when calling 'evala/preproc4'. Received: 'floats not handled' which i know is  becuase assuming but how fix it?

and how i can get that table ?

error-problem.mw

i want plot like that but i can't  and there is anyway for finding the equalibriom point of system? 

restart

with(PDEtools)

with(LinearAlgebra)

with(DEtools)

with(DynamicSystems)

sys := {diff(x(t), t) = 2*x(t)+3*y(t), diff(y(t), t) = 2*x(t)+y(t)}

{diff(x(t), t) = 2*x(t)+3*y(t), diff(y(t), t) = 2*x(t)+y(t)}

(1)

fns := {x(t), y(t)}

{x(t), y(t)}

(2)

sol := dsolve(sys, fns)

{x(t) = c__1*exp(4*t)+c__2*exp(-t), y(t) = (2/3)*c__1*exp(4*t)-c__2*exp(-t)}

(3)

ode := [diff(x(t), t) = 2*x(t)+3*y(t), diff(y(t), t) = 2*x(t)+y(t)]; S := dsolve(ode)

[diff(x(t), t) = 2*x(t)+3*y(t), diff(y(t), t) = 2*x(t)+y(t)]

 

{x(t) = c__1*exp(4*t)+c__2*exp(-t), y(t) = (2/3)*c__1*exp(4*t)-c__2*exp(-t)}

(4)

Student:-ODEs:-ODESteps(ode, {x(t), y(t)})

"[[,,"Let's solve"],[,,[(&DifferentialD;)/(&DifferentialD;t) x(t)=2 x(t)+3 y(t),(&DifferentialD;)/(&DifferentialD;t) y(t)=2 x(t)+y(t)]],["&bullet;",,"Define vector"],[,,x(t)=[?]],["&bullet;",,"Convert system into a vector equation"],[,,(&DifferentialD;)/(&DifferentialD;t) x(t)=[?]*x(t)+[?]],["&bullet;",,"System to solve"],[,,(&DifferentialD;)/(&DifferentialD;t) x(t)=[?]*x(t)],["&bullet;",,"Define the coefficient matrix"],[,,A=[?]],["&bullet;",,"Rewrite the system as"],[,,(&DifferentialD;)/(&DifferentialD;t) x(t)=A*x(t)],["&bullet;",,"To solve the system, find the eigenvalues and eigenvectors of" A],["&bullet;",,"Eigenpairs of" A],[,,[[-1,[?]],[4,[?]]]],["&bullet;",,"Consider eigenpair"],[,,[-1,RTABLE(18446744074191517278,MATRIX([[-1], [1]]),Vector[column])]],["&bullet;",,"Solution to homogeneous system from eigenpair"],[,,(x)[1]=[]],["&bullet;",,"Consider eigenpair"],[,,[4,RTABLE(18446744074192645174,MATRIX([[3/2], [1]]),Vector[column])]],["&bullet;",,"Solution to homogeneous system from eigenpair"],[,,(x)[2]=[]],["&bullet;",,"General solution to the system of ODEs"],[,,x=`c__1` (x)[1]+`c__2` (x)[2]],["&bullet;",,"Substitute solutions into the general solution"],[,,x=[]+[]],["&bullet;",,"Substitute in vector of dependent variables"],[,,[?]=[?]],["&bullet;",,"Solution to the system of ODEs"],[,,{x(t)=-`c__1` (e)^(-t)+(3 `c__2` (e)^(4 t))/2,y(t)=`c__1` (e)^(-t)+`c__2` (e)^(4 t)}]]"

(5)
 

NULL

Download Plot-1.mw

Hi,

I’m not very familiar with how to handle complex numbers in Maple, but I’d like to solve a final exam exercise (June 2025) on complex numbers in this environment. I’ve already solved it using pen and paper. My goal is also to illustrate the different results graphically, for pedagogical purposes, so I’d like some guidance to translate my pen-and-paper solution.

Bac_25_SM_Complexes.mw

Epreuve_Bac_SM_2025_Ss1_Sol-Exo-02.pdf

we have a lot term which we add thus term we can get something near to our solution but i don't know how it is work , did any one did pade approximate of series solution?

restart

with(inttrans)

with(PDEtools)

with(DEtools)

with(Physics)

declare(u(x, t), quiet); declare(v(x, t), quiet); declare(U(x, t), quiet)

undeclare(prime)

approx_sol := -8*x + 3 - 208/3*x^3 + 24*x^2 + 2279/15*x^5 + 200*x^4 + 443/63*x^7 + 676/15*x^6 + 62/2835*x^9 + 62/105*x^8

-8*x+3-(208/3)*x^3+24*x^2+(2279/15)*x^5+200*x^4+(443/63)*x^7+(676/15)*x^6+(62/2835)*x^9+(62/105)*x^8

(1)
 

NULL

Download pade-approximate.mw

in this substitution really is hard to substitute which i don't have idea how find u[1]  u[2] so on becuase  in that term some function like u[i,x] appear and i don't have idea how define them to remove it in my iteration ? and how remove the conjugate of function u[0] i did try evalc(u[0]) and make it worst than this 

which i did some special substitution which i don't know true or not but result is true i didn't have other idea for that
which  instead of any conjugate i use v and istead of any diff(u,x) i use z[i,x] and if conjugate have derivative v[i,x]

k4.mw

is easy to determine the linear and non linear part but how we can do it by command specially when contain conjugate part of function even i try to use another function instead of conjugate but stil i didn't got  the result?

restart

with(inttrans)

with(PDEtools)

with(DEtools)

with(DifferentialAlgebra)

"with(Student[ODEs][Solve]): "

with(IntegrationTools)

with(inttrans)

with(PDEtools)

with(Physics)

with(PolynomialTools)

with(RootFinding)

with(SolveTools)

with(LinearAlgebra)

with(sumtools)

declare(u(x, t))*conjugate(u(x, t))*declare(v(x, t))

u(x, t)*`will now be displayed as`*u

 

v(x, t)*`will now be displayed as`*v

 

conjugate(u(x, t))

(1)

pde := u(x, t)+I*(diff(u(x, t), `$`(x, 2)))+2*(diff(u(x, t)*conjugate(u(x, t)), x))*u(x, t)+u(x, t)^2*conjugate(u(x, t))^2*u(x, t) = 0

u(x, t)+I*(diff(diff(u(x, t), x), x))+2*((diff(u(x, t), x))*conjugate(u(x, t))+u(x, t)*(diff(conjugate(u(x, t)), x)))*u(x, t)+u(x, t)^3*conjugate(u(x, t))^2 = 0

(2)

pde_linear, pde_nonlinear := selectremove(proc (term) options operator, arrow; not has((eval(term, u(x, t) = T*u(x, t)))/T, T) end proc, expand(pde))

() = 0, u(x, t)+I*(diff(diff(u(x, t), x), x))+2*u(x, t)*(diff(u(x, t), x))*conjugate(u(x, t))+2*u(x, t)^2*(diff(conjugate(u(x, t)), x))+u(x, t)^3*conjugate(u(x, t))^2 = ()

(3)

u_occurrences := map(proc (i) options operator, arrow; numelems(select(has, [op([op(i)])], u)) end proc, oppde); linear_op_indices := ListTools:-SearchAll(1, u_occurrences); pde_linear := add(oppde[[linear_op_indices]]); pde_nonlinear := expand(simplify(expand(pde)-pde_linear))

Error, invalid input: ListTools:-SearchAll expects its 2nd argument, L, to be of type {list, rtable}, but received 0

 

[linear_op_indices]

 

Error, (in simpl/relopsum) invalid terms in sum: u(x,t)+I*diff(diff(u(x,t),x),x)+2*u(x,t)*diff(u(x,t),x)*conjugate(u(x,t))+2*u(x,t)^2*diff(conjugate(u(x,t)),x)+u(x,t)^3*conjugate(u(x,t))^2 = 0

 
 

NULL

Download solving.mw

I am trying to check if a matrix quailfies as being skew-symmetric. So I added it to it's transpose to get a matrix of zeroes. 

How do I check it  equals the ZeroMatrix?

restart

NULL

M := Matrix(3, 3, {(1, 1) = 0, (1, 2) = (40*sqrt(7)+140)/(560+sqrt(7)), (1, 3) = (80*sqrt(7)+35)/(560+sqrt(7)), (2, 1) = (-40*sqrt(7)-140)/(560+sqrt(7)), (2, 2) = 0, (2, 3) = (280-32*sqrt(7))/(560+sqrt(7)), (3, 1) = (-80*sqrt(7)-35)/(560+sqrt(7)), (3, 2) = (-280+32*sqrt(7))/(560+sqrt(7)), (3, 3) = 0})

Matrix(%id = 36893490078408578100)

(1)

mtest := simplify(M+M^%T)

Matrix(%id = 36893490078408571348)

(2)

LinearAlgebra:-ZeroMatrix(3)

Matrix(%id = 36893490078408566412)

(3)

is(mtest = LinearAlgebra:-ZeroMatrix(3))

false

(4)

NULL

Download 2025-05-30_is_Matrix_skew-symmetric.mw

I used Maple’s odetest to check the validity of 36 exact solutions.

About half of them return zero, meaning they satisfy the ODE correctly. But the other half don’t — odetest gives nonzero expressions.

My question is:

Could there be a mistake in my assumptions, or are some of these solutions only valid under specific conditions?

I’d appreciate any clarification or suggestions.

ode-test-36.mw

Hello!

I use Maple to calculate a variety of physiology questions. But i keeps using the standardized SI unites and not the ones i plug in. Is there a way to get maple to behave differently? Ive added a screenshot of an example. It comes out in m3/s when i just want it in ml/min. 

GFR = (18.2*Unit(('mg')/('ml')))*1.5*Unit(('ml')/('min'))/(0.201*Unit(('mg')/('ml')));
 = 
                                            / 3\
                                            |m |
                GFR = 0.000002263681592 Unit|--|
                                            \s /

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