Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

# The question is: how to get S knowing E?  All the attemps below after S assigned (suggested by Claude) fail factoring E.  

restart;

interface(showassumed = 0): 
assume(a>=0, b>=0, c>=0);
assume(x::real, y::real,z::real);

S := sqrt(x^2+(a+sqrt(y^2+(b+sqrt(z^2+c^2))^2))^2):
E := expand(S):
simplify(E,sqrt):
combine(E,sqrt):
simplify(combine(E, sqrt), sqrt):
evalindets(E, 'radical', f -> simplify(f, sqrt)):
simplify(radnormal(E), sqrt):

Is there an easier or shorter way to do the following? Any handy package command?

RootOf(-(Int(1/sqrt(2*cos(_a)+RootOf(-(Int(1/sqrt(2*cos(_a)+_Z), _a = 0 .. Pi))+_a+c__2)), _a = 0 .. _Z))+x+c__2)

RootOf(-(Int(1/(2*cos(_a)+RootOf(-(Int(1/(2*cos(_a)+_Z)^(1/2), _a = 0 .. Pi))+_a+c__2))^(1/2), _a = 0 .. _Z))+x+c__2)

(1)

indets(RootOf(-(Int(1/(2*cos(_a)+RootOf(-(Int(1/(2*cos(_a)+_Z)^(1/2), _a = 0 .. Pi))+_a+c__2))^(1/2), _a = 0 .. _Z))+x+c__2), 'specfunc(anything, RootOf)')

{RootOf(-(Int(1/(2*cos(_a)+_Z)^(1/2), _a = 0 .. Pi))+_a+c__2), RootOf(-(Int(1/(2*cos(_a)+RootOf(-(Int(1/(2*cos(_a)+_Z)^(1/2), _a = 0 .. Pi))+_a+c__2))^(1/2), _a = 0 .. _Z))+x+c__2)}

(2)

subs(_a = _b, {RootOf(-(Int(1/(2*cos(_a)+_Z)^(1/2), _a = 0 .. Pi))+_a+c__2), RootOf(-(Int(1/(2*cos(_a)+RootOf(-(Int(1/(2*cos(_a)+_Z)^(1/2), _a = 0 .. Pi))+_a+c__2))^(1/2), _a = 0 .. _Z))+x+c__2)}[1])

RootOf(-(Int(1/(2*cos(_b)+_Z)^(1/2), _b = 0 .. Pi))+_b+c__2)

(3)

subs({RootOf(-(Int(1/(2*cos(_a)+_Z)^(1/2), _a = 0 .. Pi))+_a+c__2), RootOf(-(Int(1/(2*cos(_a)+RootOf(-(Int(1/(2*cos(_a)+_Z)^(1/2), _a = 0 .. Pi))+_a+c__2))^(1/2), _a = 0 .. _Z))+x+c__2)}[1] = RootOf(-(Int(1/(2*cos(_b)+_Z)^(1/2), _b = 0 .. Pi))+_b+c__2), RootOf(-(Int(1/(2*cos(_a)+RootOf(-(Int(1/(2*cos(_a)+_Z)^(1/2), _a = 0 .. Pi))+_a+c__2))^(1/2), _a = 0 .. _Z))+x+c__2))

RootOf(-(Int(1/(2*cos(_a)+RootOf(-(Int(1/(2*cos(_b)+_Z)^(1/2), _b = 0 .. Pi))+_b+c__2))^(1/2), _a = 0 .. _Z))+x+c__2)

(4)

NULL

Related question: Is that substitution mathematically correct or does the _a belong to the outermost RootOf?

Download parameters_in_nested_RootOf.mw

most of time is give me true my substittuetion but  a lot time i saw it is not make my substittuetion true and this time i figure out which author did   and outcome is what i am looking but when i do that is so different where is problem 

f-m.mw

In here i did try my best and my equation outcome are ok but is not same as author did i dont know why, beside this he try to use two ode for constructing  a new ode which find of one solution of this can be the third solution of ode!
 i have to use eq(5) for my orginal ode  but eq(5) contain f(x) and g(x) which by some assumption  and taking two other ode eq(6) and eq(7)  they construct new one which is eq(14) and by f(xi) and g(xi) have corelation with W(xi) which is third ode  as in eq(9) and eq(11) mentioned, i try to use the solution which author mentioned but is not give me solution of third ode by using corelation what is problem here? also in eq(21)  and eq(25) when thus parameter are satisfy must our odetest be zero

i will update two maple file which realted separatly for constructing  equations and other is for apply and satisfy the solution for ode!

F-p.mw

ode-17.mw

 a while ago there is a code for changing function from trig to hyperbolic and viceversa  but i can't find that code except changind xi=I*xi there is another one?

restart

S3 := G(xi) = -(sqrt(Omega)*(tanh(sqrt(Omega)*xi)+I*sech(sqrt(Omega)*xi))+B)/(2*C); S4 := G(xi) = -(sqrt(Omega)*(tanh(sqrt(Omega)*xi)-I*sech(sqrt(Omega)*xi))+B)/(2*C)

G(xi) = -(1/2)*(Omega^(1/2)*(tanh(Omega^(1/2)*xi)-I*sech(Omega^(1/2)*xi))+B)/C

(1)

convert(rhs(S3), trig)

-(1/2)*(Omega^(1/2)*(tanh(Omega^(1/2)*xi)+I*sech(Omega^(1/2)*xi))+B)/C

(2)

NULL

Download convert.mw

Maple 2026 and Maple 2025.2

Is this a bug in limit? or as designed?

Doing 

limit(sol,[_C3 = 0, _C4 = 0])

Gives internal error. But

limit(sol,_C3 = 0);
limit(%,_C4 = 0);

works and no error.

Worksheet below. I've had problems before with multilimit. I think I need to change my code to do limit one by one from now on.

interface(version);

`Standard Worksheet Interface, Maple 2026.0, Windows 10, March 05 2026 Build ID 2001916`

restart;

sol:=(piecewise(t <= 0,0,t <= Pi,(arctan(tan(t))*cos(t)*_C3*_C4+arctan(tan(t))*sin(t)*_C4^2-cos(t)*_C3*_C4*t+sin(t)*_C3^2*t)/(_C3*cos(t)+_C4*sin(t)),Pi < t,Pi*_C3*(tan(t)*_C3-_C4)/(tan(t)*_C4+_C3))+2*_C3^2+2*_C4^2)*(_C3*cos(t)+_C4*sin(t))/(2*_C3^2+2*_C4^2)

sol := (piecewise(t <= 0, 0, t <= Pi, (arctan(tan(t))*cos(t)*_C3*_C4+arctan(tan(t))*sin(t)*_C4^2-cos(t)*_C3*_C4*t+sin(t)*_C3^2*t)/(cos(t)*_C3+_C4*sin(t)), Pi < t, Pi*_C3*(tan(t)*_C3-_C4)/(tan(t)*_C4+_C3))+2*_C3^2+2*_C4^2)*(cos(t)*_C3+_C4*sin(t))/(2*_C3^2+2*_C4^2)

limit(sol,[_C3 = 0, _C4 = 0])

Error, (in limit/multi/ldegree1) invalid input: limit/multi/ReIm expects its 1st argument, f, to be of type polynom, but received _DIR1*cos(t)+_DIR2*sin(t)

limit(sol,_C3 = 0)

piecewise(t <= 0, _C4*sin(t), t <= Pi, (1/2)*arctan(tan(t))*sin(t)+_C4*sin(t), _C4*sin(t))

limit(%,_C4 = 0)

piecewise(t <= 0, 0, t <= Pi, (1/2)*arctan(tan(t))*sin(t), 0)

 

 

Download limit_problem_april_2_2026.mw

Dear MaplePrimes Community,

With the rapid advancement of technology, it has become increasingly important to improve our workflows and reduce repetitive manual tasks. In my research, I frequently use Maple to solve  PDEs and obtain analytical solutions. However, one major difficulty I face is the time-consuming process of manually converting each result into LaTeX format.

In many cases, Maple produces multiple solutions, and I must copy and convert each one individually into LaTeX. This process is tedious, inefficient, and prone to formatting errors.

Therefore, I would like to ask:

Is it possible to establish an automated link or workflow between Maple and LaTeX such that, after solving a problem in Maple, all resulting expressions are automatically exported and inserted into a LaTeX document in the desired format and location?

Ideally, I am looking for a system where:

  • Maple solutions are automatically converted into LaTeX code,

  • The equations are directly placed into a specified LaTeX file or section,

  • The process works for multiple solutions without manual intervention.

Furthermore, if such integration is possible, can this idea be extended to connect Maple with multiple tools (for example, Maple → LaTeX → plotting or document-generation systems), creating a more advanced automated workflow?

I would greatly appreciate any guidance, examples, or best practices for implementing such a system.

Thank you very much for your help.

Hello :)

I am totally lost with ( I guess) the LREtools package.

See the ws.

The two integer sequences come from the OEIS.

I don't understand why i have a FAILED result.

intseq := [0, 1, 4, 2, 131, 129, 3, 5, 16, 14, 12, 10, 8, 6, 31, 29, 27, 25, 23, 99734, 7, 9, 11, 13, 15, 17, 64, 62, 60, 58, 56, 54, 52, 50, 48, 46, 44, 42, 40, 38, 111, 22, 20, 18, 28, 30, 32, 222, 220, 218, 216, 214, 212, 210, 208, 206, 204, 202, 200, 198, 196]

[0, 1, 4, 2, 131, 129, 3, 5, 16, 14, 12, 10, 8, 6, 31, 29, 27, 25, 23, 99734, 7, 9, 11, 13, 15, 17, 64, 62, 60, 58, 56, 54, 52, 50, 48, 46, 44, 42, 40, 38, 111, 22, 20, 18, 28, 30, 32, 222, 220, 218, 216, 214, 212, 210, 208, 206, 204, 202, 200, 198, 196]

(1)

rec := LREtools:-GuessRecurrence(intseq, q(n))

FAIL

(2)

soln := rsolve(rec, q)

Error, (in rsolve/single) Equations do not involve function q()

 

restart

NULL

intseq := [0, 1, 3, 6, 2, 7, 13, 20, 12, 21, 11, 22, 10, 23, 9, 24, 8, 25, 43, 62, 42, 63, 41, 18, 42, 17, 43, 16, 44, 15, 45, 14, 46, 79, 113, 78, 114, 77, 39, 78, 38, 79, 37, 80, 36, 81, 35, 82, 34, 83, 33, 84, 32, 85, 31, 86, 30, 87, 29, 88, 28, 89, 27, 90, 26, 91, 157, 224, 156, 225, 155]

[0, 1, 3, 6, 2, 7, 13, 20, 12, 21, 11, 22, 10, 23, 9, 24, 8, 25, 43, 62, 42, 63, 41, 18, 42, 17, 43, 16, 44, 15, 45, 14, 46, 79, 113, 78, 114, 77, 39, 78, 38, 79, 37, 80, 36, 81, 35, 82, 34, 83, 33, 84, 32, 85, 31, 86, 30, 87, 29, 88, 28, 89, 27, 90, 26, 91, 157, 224, 156, 225, 155]

(3)

rec := LREtools:-GuessRecurrence(intseq, q(n))

FAIL

(4)

NULL

restart

with(LREtools)

intseq := [0, 1, 3, 6, 2, 7, 13, 20, 12, 21, 11, 22, 10, 23, 9, 24, 8, 25, 43, 62, 42, 63, 41, 18, 42, 17, 43, 16, 44, 15, 45, 14, 46, 79, 113, 78, 114, 77, 39, 78, 38, 79, 37, 80, 36, 81, 35, 82, 34, 83, 33, 84, 32, 85, 31, 86, 30, 87, 29, 88, 28, 89, 27, 90, 26, 91, 157, 224, 156, 225, 155]

[0, 1, 3, 6, 2, 7, 13, 20, 12, 21, 11, 22, 10, 23, 9, 24, 8, 25, 43, 62, 42, 63, 41, 18, 42, 17, 43, 16, 44, 15, 45, 14, 46, 79, 113, 78, 114, 77, 39, 78, 38, 79, 37, 80, 36, 81, 35, 82, 34, 83, 33, 84, 32, 85, 31, 86, 30, 87, 29, 88, 28, 89, 27, 90, 26, 91, 157, 224, 156, 225, 155]

(5)

NULL

rec := GuessRecurrence(intseq, q(n))

FAIL

(6)

NULL

NULL Thank you everyone.

Jean-Michel

Download rsolve_and_LREtools.mw

There seems to be a regression in Maple 2026 in the XMLTools:-ParseFile function.

As Maple2026 is not yet in the list of products to be chosen, I have added it in the subject.

Error, (in XMLTools:-ParseFile) invalid input: too many and/or wrong type of arguments passed to XMLTools:-NSXML:-Parser:-ParseFile; first unused argument is prolog = true

The test file is right from the help related to ParseFile.
Test_XML.mw

Objective: Solve a system of two equations.

Obstacle: Generating these two equations depends on millions of previous combinations as well as derivatives.

In other words, we've reached the maximum limit that Maple on my computer can handle.

What would be better, to leave the equations aside or to upgrade my computer?

restart

with(plots)

with(linalg)

H01 := -gamma11*S11-gamma12*S12-gamma13*S13-gamma14*S14-gamma15*S15-gamma16*S16-gamma17*S17-gamma18*S18-gamma19*S19-gamma110*S110-gamma111*S111-gamma112*S112-eta1*(S11^2+S110^2+S111^2+S112^2+S12^2+S13^2+S14^2+S15^2+S16^2+S17^2+S18^2+S19^2)-J1*(S11*S12+S12*S13+S13*S14+S14*S18+S18*S112+S112*S111+S111*S110+S110*S19+S19*S15+S15*S11+S16*(S12+S110+S15+S17)+S17*(S113+S111+S18))

H02 := -gamma21*S21-gamma22*S22-gamma23*S23-gamma24*S24-gamma25*S25-gamma26*S26-gamma27*S27-gamma28*S28-gamma29*S29-gamma210*S210-gamma211*S211-gamma212*S212-eta2*(S21^2+S210^2+S211^2+S212^2+S22^2+S23^2+S24^2+S25^2+S26^2+S27^2+S28^2+S29^2)-J1*(S21*S22+S22*S23+S23*S24+S24*S28+S28*S212+S212*S211+S211*S210+S210*S29+S29*S25+S25*S21+S26*(S22+S210+S25+S27)+S27*(S213+S211+S28))

Z01 := exp(-beta*H01)

Z01 := add(add(add(add(add(add(add(add(add(add(add(add(Z01, S11 = [-2, -1, 0, 1, 2]), S12 = [-2, -1, 0, 1, 2]), S13 = [-2, -1, 0, 1, 2]), S14 = [-2, -1, 0, 1, 2]), S15 = [-2, -1, 0, 1, 2]), S16 = [-2, -1, 0, 1, 2]), S17 = [-2, -1, 0, 1, 2]), S18 = [-2, -1, 0, 1, 2]), S19 = [-2, -1, 0, 1, 2]), S110 = [-2, -1, 0, 1, 2]), S111 = [-2, -1, 0, 1, 2]), S112 = [-2, -1, 0, 1, 2])

NULL

Z02 := exp(-beta*H02)

Z02 := add(add(add(add(add(add(add(add(add(add(add(add(Z02, S21 = [-2, -1, 0, 1, 2]), S22 = [-2, -1, 0, 1, 2]), S23 = [-2, -1, 0, 1, 2]), S24 = [-2, -1, 0, 1, 2]), S25 = [-2, -1, 0, 1, 2]), S26 = [-2, -1, 0, 1, 2]), S27 = [-2, -1, 0, 1, 2]), S28 = [-2, -1, 0, 1, 2]), S29 = [-2, -1, 0, 1, 2]), S210 = [-2, -1, 0, 1, 2]), S211 = [-2, -1, 0, 1, 2]), S212 = [-2, -1, 0, 1, 2])

NULL

G0 := -(1/24)*N*ln(Z01*Z02)/beta

NULL

m01 := -24*(diff(G0, gamma11))/N

m02 := -24*(diff(G0, gamma21))/N

NULL

beta := 11.605/T; gamma11 := 2*J1*m1+2*J2*m2; eta1 := Delta; gamma21 := 2*J1*m2+2*J2*m1; eta2 := Delta; gamma12 := J1*m1+2*J2*m2; gamma22 := J1*m2+2*J2*m1; gamma13 := J1*m1+2*J2*m2; gamma23 := J1*m2+2*J2*m1; gamma14 := 2*J1*m1+2*J2*m2; gamma24 := 2*J1*m2+2*J2*m1; gamma15 := J1*m1+2*J2*m2; gamma25 := J1*m2+2*J2*m1; gamma16 := 2*J2*m2; gamma26 := 2*J2*m1; gamma17 := 2*J2*m2; gamma27 := 2*J2*m1; gamma18 := J1*m1+2*J2*m2; gamma28 := J1*m2+2*J2*m1; gamma19 := 2*J1*m1+2*J2*m2; gamma29 := 2*J1*m2+2*J2*m1; gamma110 := J1*m1+2*J2*m2; gamma210 := J1*m2+2*J2*m1; gamma111 := J1*m1+2*J2*m2; gamma211 := J1*m2+2*J2*m1; gamma112 := 2*J1*m1+2*J2*m2; gamma212 := 2*J1*m2+2*J2*m1

NULL

eq1 := m1 = m01

eq2 := m2 = m02

fsolve(subs(J1 = 2*.83, N = 1, J2 = -2*.58, Delta = 0, m1 = 0.1e-1, {eq1, eq2}), {T, m2}, T = 0 .. 220, m2 = -.1 .. .1)

fsolve(subs(J1 = 2*.83, N = 1, J2 = -2*.58, T = 1, Delta = 0, {eq1, eq2}), {m1, m2}, m1 = -5 .. 5, m2 = -5 .. 5)

{m1 = 2.000000000, m2 = -2.000000000}

(1)

NULL

Download Maple_forum_test.mw

After system restart on Windows 11:

Maple 2026 was installed as usual with administrator rights and with import of preferences from Maple 2025. Other Maple versions are working on the same computer.

Does this never-seen-before output gives a hint what is wrong with my installation. Any suggestions what to do next?

That's the installed version:

On the same Windows 11 installation, Maple 2025 displays equation labels but Maple 2026 not. Maple 2025 settings were imported for the 2026 installation. (Crtl-l + number does return "invalid label". I assume that not labels have been generated)

Has anbody observed the same? Any suggestions what I could check/do?

Other observation: The output font does not look the same

 

Is there an easy way to read jld files in Maple?  Apparently JLD (Julia Data) and JLD2 files are binary formats primarily designed for saving and loading Julia variables, preserving types.

That's how it looks like in Maple 2026.0 for the Examples on the Help-page for topic solve, on Windows 11

For a fraction of a second I see output rendered in blue and Math-2D. Then it turns to the above.

Is this a regression or a new feature? How to get the output back to Math-2D and blue?

Glad that 2026 has been released now, and I will certainly use it as the default version in the future.

While there are a lot of new features, one thing that keeps annoying me is the inferior font quality, compared to other software.

Here's a screenshot of 4 different programs, all with font Arial 11pt and 100% zoom factor. Font AntiAliasing is set to enabled in Maple.

Judge it for yourself, but in my opinion it gives a clear picture that Maple is much worse to read than any of the other software packages (Word, LibreOffice, pdfXChange).

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