Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

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

In the Maple 2025 release, the "old/former" user interface was provided under a menu entry called "Maple 2025 for Screen Readers". I used this version because of shortcommings of the new ribbon interface. Some of the shortcommings have been fixed in the 2026 release but it is still slow to use (I am missing functions in the quick access toolbar that I frequently use and no customization of the bar seems possible).

I could not find a similar menu entry in the 2026 release under Windows 11. Is the old unser interface still available?

I am solving a hybrid nanofluid flow problem in a bifurcated artery using Maple. The governing equations for velocity and temperature are solved using dsolve(..., numeric, method=bvp[midrich]).

My Maple code successfully produces for both the artery  parentartery_and_daughter_artery_error.mw.

The velocity profiles are obtained correctly using odeplot.

However, I want to compute additional physical quantities and generate plots similar to the velocity profiles.

Specifically I want to plot:

  1. Flow rate Q versus axial distance z

  2. Impedance (flow resistance) λ versus z

  3. Wall shear stress τ versus z

for different values of Hartmann number Ha.

The formulas I am using are

Flow rate:

Q=2π(R∫01ηw(η) dη+R2∫01w(η) dη)Q = 2\pi \left( R \int_0^1 \eta w(\eta)\,d\eta + R_2 \int_0^1 w(\eta)\,d\eta \right)Q=2π(R∫01​ηw(η)dη+R2​∫01​w(η)dη)

Wall shear stress:

τ=μ∣dwdr∣\tau = \mu \left|\frac{dw}{dr}\right|τ=μ​drdw​​

Impedance:

λ=∣dp/dz∣Q\lambda = \frac{|dp/dz|}{Q}λ=Q∣dp/dz∣​
Please help me to solve this question.

For decades, Maple has been built around one of the world’s most powerful mathematics engines—helping students, educators, engineers, and researchers explore ideas, solve complex problems, and communicate mathematics clearly.

Maple 2026 builds on that foundation with major advances in the math engine, expanding the kinds of problems Maple can solve while improving reliability and performance.

At the same time, Maple 2026 introduces new AI-powered tools that help you work faster—finding commands, generating visualizations, explaining concepts, and helping you explore ideas. The key difference is that these tools sit on top of Maple’s math engine, so the results are grounded in real computation rather than guesswork.

If you’ve been following along with our recent Mathy teaser videos and sneak peek posts, you may already have seen hints of some of these features. Now I’m excited to finally share them in full.

One of the most exciting additions in Maple 2026 is the new AI Assistant.

AI tools are incredibly useful for exploring ideas, writing code, and learning new topics. But when the mathematics becomes more involved, relying on AI alone can be risky. The Maple AI Assistant brings those productivity benefits into Maple while keeping the mathematics grounded in Maple’s trusted computation engine.

You can ask the AI Assistant questions in natural language and have it help you:

  • find Maple commands or formulas
  • generate Maple code
  • create visualizations
  • explain mathematical concepts
  • draft examples, worksheets, or reports

Because Maple performs the underlying computations where appropriate, the results are grounded in Maple’s powerful math engine. The AI Assistant becomes a productivity partner that helps you accomplish tasks in Maple faster and more easily, combining the flexibility of AI with mathematics you can trust.

Watch the AI Assistant in action.

 
Turn Documents into Live Mathematics

Another feature I’m particularly excited about is Document Import.

Many of us have years of mathematical content stored in PDFs, lecture notes, journal articles, slides, or even handwritten pages. Traditionally these documents are static—you can read them, but you can’t interact with the mathematics inside them.

With Maple 2026, that changes.

Document Import allows Maple to convert many document formats—including PDFs, DOCX files, and presentations—into Maple worksheets where the mathematics becomes live and executable. 

The image below illustrates the transformation.

On the left (“Before”), scribbled handwritten notes from a Calculus III lecture were saved in a Word document. The notes include hand-drawn sketches, formulas, and written explanations.

After importing the document into Maple (“After”), the mathematical expressions were recognized and converted into live, editable Maple mathematics. The text was preserved, and the hand-drawn sketches were retained as images. The resulting worksheet supports evaluation, editing, and further computation.

Once imported, you can:

  • evaluate expressions
  • modify formulas
  • extend derivations
  • add visualizations
  • explore variations of the mathematics

Instead of recreating examples from scratch, you can bring existing material directly into Maple and start exploring.

While the new AI features are exciting, the heart of Maple has always been its mathematics engine—and Maple 2026 delivers significant advances here.

One particularly notable improvement is Maple’s expanded ability to solve linear recurrence equations. Through improvements to the rsolve command and major extensions to the LREtools package, Maple can now solve dramatically more recurrence relations than before, including many third- and fourth-order cases that were previously beyond reach.

In fact, Maple can now fully solve over 94% of the 55,979 entries in the Online Encyclopedia of Integer Sequences (OEIS) that that can be shown to satisfy a linear recurrence relation. These advances reflect ongoing research into linear difference equations and their algorithmic implementation in Maple, continuing Maple’s long tradition of advancing the state of computer algebra.

Beyond recurrence solving, Maple 2026 includes many improvements across its core symbolic and numeric algorithms. Maple’s assumption system has been strengthened to improve reasoning under mathematical assumptions, and enhancements to the simplify, combine, and evalc commands allow Maple to produce more compact and mathematically natural forms for a wider range of expressions.

There are also improvements to Maple’s differential equation solvers, polynomial system solving, and numerical solving routines such as fsolve, along with updates to other foundational parts of the math library used throughout the system.

Taken together, these improvements expand the range of problems Maple can solve and improve the robustness, correctness, and efficiency of the results.

Maple has always offered extensive control over plotting options, but achieving consistent visual styling across multiple plots could require specifying many settings each time.

Maple 2026 introduces Plotting Themes, which allow you to define a plotting style once and apply it across many plots with a single option.

Themes make it easy to maintain consistent visual styles in worksheets, teaching materials, reports, and publications, while still allowing individual plots to override specific options when needed.

The image below shows an example of creating and applying a custom plotting theme. 

 

Maple continues to be widely used in classrooms around the world, and Maple 2026 includes several improvements designed to support teaching and learning.

The Check My Work system has been enhanced so Maple can recognize a wider variety of valid student solution steps and provide more accurate feedback.

Maple 2026 also improves the generation of similar practice problems, making it easier to create variations of a problem while preserving its mathematical structure.

In addition, Maple’s step-by-step solutions have been expanded to support more types of expressions, helping students better understand the reasoning behind the mathematics they’re learning.

Maple 2026 also introduces improvements for developers building advanced applications, along with performance enhancements across the system.

One particularly interesting addition is the new VectorSearch package, which implements a vector database directly inside Maple.

If you’re not familiar with vector databases, one way to think about them is through recommendation systems like Netflix or Spotify. Each movie or song can be represented by a vector containing thousands of numbers describing its characteristics—things like genre, pacing, or mood. When you watch something, the system finds other items whose vectors are closest to it, which is how recommendations are generated.

With the new VectorSearch package, Maple can store thousands (or more) of vectors and efficiently find the ones most similar to a given vector. This makes it easier to build applications involving machine learning, data analysis, and modern AI workflows directly in Maple.

Maple 2026 also delivers significant performance improvements. For example, operations involving quantities with units have been greatly optimized—some computations now run over 90 times faster, making Maple even more efficient for engineering and scientific workflows.

Maple 2026 also expands the benefits available through the Maplesoft Elite Maintenance Program (EMP). The new benefits include access to additional Maplesoft products and services:

  • Maple Learn, the online environment for teaching and learning mathematics
  • Maple Calculator Premium, bringing the power of Maple to your phone with full access to features like Solution Steps and Check My Work
  • Maple MCP, which allows you to connect Maple’s math engine to external AI tools so they can produce mathematical results you can trust

These additions extend Maple beyond the desktop, giving users powerful tools for learning, teaching, and exploring mathematics across web and mobile platforms, as well as through integrations with external AI tools.

This post only scratches the surface of what’s new in Maple 2026. There are many more improvements across the math library, programming tools, and performance.

To learn more about all the new features and enhancements in Maple 2026, visit the What’s New in Maple page on our website.

 

 

I am studying a nonlinear wave equation and trying to reproduce the energy balance method shown in a research paper. First, the original partial differential equation is reduced to an ordinary differential equation using a traveling wave transformation. After obtaining the reduced equation, the paper rewrites it in a form suitable for the energy balance method and derives the corresponding variational principle and Hamiltonian invariant. Then a trial periodic solution in cosine form is assumed. Using the Hamiltonian invariant and some initial conditions, the parameters of the trial function are determined and a periodic solution is obtained.

I would like to know how to implement this procedure in Maple. Specifically,  compute the Hamiltonian invariant from the equation, substitute the cosine trial function, and determine the unknown parameter in the trial solution using the energy balance method. I will attach images from the paper that show the derivation steps I am trying to reproduce. Any guidance on how to perform these symbolic steps in Maple would be very helpful.

f-s.mw

Hi,

For a pedagogical purpose, I am trying to illustrate the orthogonal projection H of a point A on to a plane P1

I constructed the line l1​ passing through A and perpendicular to the plane P1
However, in the graphical visualization, the line does not appear to be perpendicular to the plane. Visually, it gives the impression that the line is not orthogonal to the plane.

Do you have any idea what might cause this effect?

Thank you for your help.

Q_Espace.mw

While teaching a linear programming course I put together a worksheet to illustrate finding the largest disk inside a convex polygon as in section 2.6 of Understanding and Using Linear Programming by Jiří Matoušek and Bernd Gärtner.  I used both the Optimization[LPSolve] and the simplex[maximize] routines with the same objective and the same constraints.  Optimization[LPSolve] gives the correct answer but simplex[maximize] does not.  Is this a bug or did I do something wrong?

Below is the worksheet.

This is a worksheet I put together to illustrate finding the largest disk inside a convex polygon

as in section 2.6 of Understanding and Using Linear Programming by Jiří Matoušek and
Bernd Gärtner.

 

Load some packages

restart;
with(plots):
with(LinearAlgebra):
with(simplex);
with(Optimization);
with(plottools):

[basis, convexhull, cterm, define_zero, display, dual, feasible, maximize, minimize, pivot, pivoteqn, pivotvar, ratio, setup, standardize]

 

[ImportMPS, Interactive, LPSolve, LSSolve, Maximize, Minimize, NLPSolve, QPSolve]

(1)

This is the definition of the polygon with a plot.

constr := [x>=0,y>=0,y-x<=1,x+y<=5,2*x-y<=6];
Vector(%);
P1 := inequal(constr,x=-1..5,y=-1..5,optionsfeasible=[color=yellow]);

These are the outward normals for the polygon constraints and the "b" values.

C := [<-1,0>,<0,-1>,<-1,1>,<1,1>,<2,-1>];
b:= [0,0,1,5,6];

This is a check that a point is inside the polygon with plots of the circle that show the distance

from the pont to each defining line.

Pt := <1,.5>;
Pta := convert(Pt,list);
map(is,subs({x=2,y=1},constr));

DistV := [seq((b[i]-DotProduct(C[i],Pt))/Norm(C[i],2),i=1..5)]

P2 := NULL:
for i from 1 to 5 do
   P2 := P2,plottools:-circle(Pta,DistV[i],color=cat("Bright ",i),thickness=2);
   end do:
display([P2,P1],scaling=constrained);
display(map(display,[seq([P1,[P2][i]],i=1..5)]),insequence=true,scaling=constrained);

Now the distance from a point in the polygon to the boundary of the polygon is defined

and plotted over the polygon.

f := proc(x,y)
     local i;
      uses LinearAlgebra;
     min([seq(((b[i]-DotProduct(C[i],<x,y>))/Norm(C[i],2)),i=1..5)]);
     end proc;

plot3d(f(x,y),x=-1..5,y=-1..5,view=[0..4,0..4,0..2],grid=[151,151]);

 

This is the definition of the LP to solve the largest disk problem.

 

We maximize z, a variable less than or equal to the distance from the point x, y to each of the

constraining lines where x, yis in the polygon.

 

C := [<-1,0>,<0,-1>,<-1,1>,<1,1>,<2,-1>];
b:= [0,0,1,5,6];

Here are the constraints.

obj := z;
constr1 := [x>=0,y>=0,y-x<=1,x+y<=5,2*x-y<=6];
constr2 := [seq(z<=(b[i]-DotProduct(C[i],<x,y>))/Norm(C[i],2),i=1..5)];

This is the solution using LPSolve.

soln1 := LPSolve(z,[op(constr1),op(constr2)],maximize);

This is the same problem solved with simplex[maximize].

soln2 := simplex:-maximize(z,[op(constr1),op(constr2)]);
evalf(%);

The answer from the simplex package is wrong as the following plots demonstrate.

 

What, if anything, did I do wrong?

constr := [x>=0,y>=0,y-x<=1,x+y<=5,2*x-y<=6]:
P1 := inequal(constr,x=-1..5,y=-1..5,optionsfeasible=[color=yellow]):
Pt1 := subs(soln1[2],[x,y]);
C1 := disk(Pt1,soln1[1],color=red):
display(C1,P1);
Pt2 := subs(soln2,[x,y]);
r2 := subs(soln2,z);
C2 := disk(Pt2,r2,color=red):
display(C2,P1);

 

 

Download inscribe_circle.mw

I have generated the regional plot, but the boundaries between the regions are not clearly marked. How can I highlight the boundaries using a black line or a black dashed line, similar to the example image? What syntax should I use? I have attached my generated image file.

Regional_new.mw


how we solve x^(4)-12*x-12=0, 

mapple sent after solving: RootOf(_Z^4 - 12*_Z - 12, index = 1), RootOf(_Z^4 - 12*_Z - 12, index = 2), RootOf(_Z^4 - 12*_Z - 12, index = 3), RootOf(_Z^4 - 12*_Z - 12, index = 4)

sol:=ln( (y-1)^(1/3)* (y^2+y+1)^(1/3) ) - ln(y) = 2/5* ln(t^2+1)+_C1;
solve( sol,y);

in real domain is fine also. But all my attempts failed. I waited 3-4 minutes each time and stopped it.

Any one can find a trick? Below worksheet showing my attempts and also solution by Mathematica which took 0.3 seconds

Make sure to save all your work first. This problem is known to crash Maple !

restart;

sol:=ln( (y-1)^(1/3)* (y^2+y+1)^(1/3) ) - ln(y) = 2/5* ln(t^2+1)+_C1;
solve( sol,y) assuming real;

ln((y-1)^(1/3)*(y^2+y+1)^(1/3))-ln(y) = (2/5)*ln(t^2+1)+_C1

Warning,  computation interrupted

restart;

sol:=ln( (y-1)^(1/3)* (y^2+y+1)^(1/3) ) - ln(y) = 2/5* ln(t^2+1)+_C1;
RealDomain:-solve( sol,y);

ln((y-1)^(1/3)*(y^2+y+1)^(1/3))-ln(y) = (2/5)*ln(t^2+1)+_C1

Warning,  computation interrupted

restart;
sol:=ln( (y-1)^(1/3)* (y^2+y+1)^(1/3) ) - ln(y) = 2/5* ln(t^2+1)+_C1;
solve( sol,y,real);

ln((y-1)^(1/3)*(y^2+y+1)^(1/3))-ln(y) = (2/5)*ln(t^2+1)+_C1

Warning,  computation interrupted

 


 

Download solve_problem_march_7_2026.mw

 

2 3 4 5 6 7 8 Last Page 4 of 2247