Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I need to compute complicated numerical integral, and my Maple integration does not produce very stable results. Is the a way to do that on line? 

In the below I had to add the assumption x>=0 to get simplifications. Am I wrong with my interpretation that the other assumptions should have been sufficient?

We are looking for the smallest natural number n with the property that both the digit sum Q(n) of the number n and the digit sum Q(n + 1) of the successor of n are divisible by 5.

In Maple 2025.1, why don't all the hotkeys show up when I press the Alt key?

Maple 2025.0 is OK,but update to 2025.1 all the hotkeys don't show up when I press the Alt key.Why?

During a birthday party, the birthday child realizes: In 1968, I was the same age as the sum of the digits of my birth year. How old will I be now at the end of 2025?

(Please no AI solution)

In the decimal system, specify the smallest natural number k that begins with the digit 7 and has the following additional property:
If you delete the first digit 7 and write it at the end, the newly created number z = (1/3)*k.

My name is William Cheung.  I am a Snapdragon X AIPC product manager.  I am reaching out to inquire about the current and future compatibility of Maple on Snapdragon X Windows PCs, which run on ARM64 architecture. With the growing adoption of ARM-based systems in enterprise and academic environments—especially those powered by Qualcomm’s Snapdragon X Elite platform—we are evaluating software solutions that can run efficiently on these devices.

Could you please confirm:

  1. Whether Maple, as a commercial product, currently supports Windows on ARM (either natively or via emulation)?
  2. If not, is there a roadmap or timeline for ARM-native support?
  3. Are there any known performance or compatibility issues when running Maple on Snapdragon X PCs using Microsoft’s Prism emulation?

Note that I have successfully installed and run Maple 2025 with sample worksheet files on my Snapdragon X Windows PC.

Thank you for your time and assistance. I look forward to your response.

Hello Maple Community,

I'm working on solving a PDE system for fluid flow in an L-shaped cavity (similar to previous work on H-shaped and square domains ) recently asked by some one here. I've implemented the governing equations in Maple but need help generating contour plots.

I'm attaching:
1. My Maple worksheet l_shape_cavity.mw
2. A PDF with the problem description and equations L_shape_cavity_work.pdf

The main issues I'm facing are:
- Setting up proper boundary conditions for the L-shaped domain
- Generating contour plots for velocity and pressure fields
- Ensuring the solution converges properly

Could you please help me with:
1. Correct implementation of the L-shaped domain geometry
2. Generating proper contour plots
3. Any suggestions for improving the numerical solution

Thank you for your assistance!

In Maple 2025 we are having problems with units.

In Maple 2024 we did not see same errors. Maple sheets working with Maple 2024 does not run in Maple 2025.

In Maple 2025, some units which are f.x. obviously 'J' is not recognices as 'J'.

In Maple 2025, part of units dissapears from expressions.

Just an example; the function v(t) below has correct units m/s, but when using simplify, the unit is changes to m

The Maple file: Error_with_units.mw

Is there a way to use alternative plotting packages in Maple such as open source options: Plotly or Matplotlib?

https://plotly.com/python/

https://matplotlib.org/

It's possible to carve a hole through a unit cube, without splitting it into pieces, so that another unit cube can pass through that hole.  This is know as the Prince Rupert problem and was first analyzed by John Wallis, a contemporary of Isaac Newton.  Here's what the result looks like:

If your computer can play audio, have a look at  Ruperts Cube with music!

Here is the worksheet that produced the cube and the animation: ruperts-cube.mw

Solving an ode, dsolve says it used exact method and gives two solutions as result (correct result).

But when asking dsolve to solve same ode but now specifying that it uses exact method (i.e. same one it used itself before), now it gives one solution only, not two as before.

Why is that? Should not both commands give same result? i.e. two solutions?

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

SupportTools:-Version()

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

Physics:-Version()

`The "Physics Updates" version in the MapleCloud is 1881 and is the same as the version installed in this computer, created 2025, October 7, 16:4 hours Pacific Time.`

restart;

ode:=3*y(x)^3*x^2+y(x)^4+(3*x^3*y(x)^2+y(x)^4+4*x*y(x)^3)*diff(y(x),x) = 0;
DEtools:-odeadvisor(ode);
infolevel[dsolve]:=5:

3*y(x)^3*x^2+y(x)^4+(3*x^3*y(x)^2+y(x)^4+4*x*y(x)^3)*(diff(y(x), x)) = 0

[_exact, _rational]

sol:=dsolve(ode); #gives two solutions

Classification methods on request

Methods to be used are: [exact]

----------------------------

* Tackling ODE using method: exact

--- Trying classification methods ---

trying exact

<- exact successful

y(x) = 0, x*y(x)^4+x^3*y(x)^3+(1/5)*y(x)^5+c__1 = 0

maple_sol:=dsolve(ode,[exact]);  #why y=0 solution do not show here??

Classification methods on request

Methods to be used are: [exact]

----------------------------

* Tackling ODE using method: exact

--- Trying classification methods ---

trying exact

<- exact successful

x*y(x)^4+x^3*y(x)^3+(1/5)*y(x)^5+c__1 = 0

sol:=dsolve(ode); #gives two solutions again

Classification methods on request

Methods to be used are: [exact]

----------------------------

* Tackling ODE using method: exact

--- Trying classification methods ---

trying exact

<- exact successful

y(x) = 0, x*y(x)^4+x^3*y(x)^3+(1/5)*y(x)^5+c__1 = 0

 

 

Download why_different_solutions_maple_2025_1_oct_22_2025.mw

I've been using the $include statement when using code written in Visual Basic Editor for packages quite extensively, and this works very nice.

Would it be possible to outsource code to external files for code not in packages / libraries as well? If yes, how would this be done?

In the below I was looking for an explanation why D is not printed in italic whereas I does.

I

I

(1)

whattype(I)

complex(extended_numeric)

(2)

_local(I)

I

(3)

whattype(I)

symbol

(4)

D

D

(5)

whattype(D)

procedure

(6)

_local(D)

D

(7)

whattype(D)

symbol

(8)

NULL

After the local call D and I are both symbols but not printed the same way. Why is that?

Download local_D_still_roman.mw

Just found strange bug in Maple. 

If I do solve(eq1...) then solve(eq2,...) then the second solve now gives internal error.

But If I do (from clean start)  just solve(eq2,...) then it times out ok, No internal error. So clearly the first call to solve changed something in internal memory/cache to do this.

Is there a way to correctly "clear" solve cache or its memory tables after each call, so that earlier calls to solve do not affect future calls behavior? Was going to add this to my collection of maple bugs post, but thought a new question will be better to make it easier to reply.

The problem is that these internal Maple errors can not be cought by try/catch. Which means the whole program crashes and there is no way to continue automatically.

restart;

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1881 and is the same as the version installed in this computer, created 2025, October 7, 16:4 hours Pacific Time.`

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

solve timesout OK if no call made before it

 

restart;

kernelopts('assertlevel'=2):

eq_2:=1 = -X*((1/2/(2*X*Y*a+2*X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*Y*(2*Y
*a+2*a*y0)+1/2/(2*X*Y*a+2*X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*y0*(
2*Y*a+2*a*y0))/(-2*X*Y*a-2*X*a*y0-2*Y*a*x0-2*a*x0*y0+Y^2+2*Y*y0+a^2+y0^2)-(-a*Y
-a*y0+(2*X*Y*a+2*X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*Y+(2*X*Y*a+2*
X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*y0)/(-2*X*Y*a-2*X*a*y0-2*Y*a*
x0-2*a*x0*y0+Y^2+2*Y*y0+a^2+y0^2)^2*(-2*Y*a-2*a*y0))/Y/((-a+1/2/(2*X*Y*a+2*X*a*
y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*Y*(2*X*a+2*a*x0-2*Y-2*y0)+(2*X*Y*a
+2*X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)+1/2/(2*X*Y*a+2*X*a*y0+2*Y*a
*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*y0*(2*X*a+2*a*x0-2*Y-2*y0))/(-2*X*Y*a-2*X*
a*y0-2*Y*a*x0-2*a*x0*y0+Y^2+2*Y*y0+a^2+y0^2)-(-a*Y-a*y0+(2*X*Y*a+2*X*a*y0+2*Y*a
*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*Y+(2*X*Y*a+2*X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2
-2*Y*y0-y0^2)^(1/2)*y0)/(-2*X*Y*a-2*X*a*y0-2*Y*a*x0-2*a*x0*y0+Y^2+2*Y*y0+a^2+y0
^2)^2*(-2*X*a-2*a*x0+2*Y+2*y0)):
try
    timelimit(30,solve(identity(eq_2,X),[x0,y0]));
catch:
    print("good. No crash");
end try;

"good. No crash"

 

 

Same solve gives internal error when calling another solve before it

 

restart;

kernelopts('assertlevel'=2):

 

 

eq_1:=1 = -X*(1/2*(2*X*Y^3-4*X*Y^2*a+6*X*Y^2*y0-8*X*Y*a*y0+6*X*Y*y0^2-4*X*a*y0^2+2*X*
y0^3+2*Y^3*x0-4*Y^2*a*x0+6*Y^2*x0*y0-8*Y*a*x0*y0+6*Y*x0*y0^2-4*a*x0*y0^2+2*x0*
y0^3+2*X*Y+2*X*y0+2*Y*x0+2*x0*y0)/(X^2*Y*a+X^2*a*y0+2*X*Y*a*x0+2*X*a*x0*y0+Y*a*
x0^2+a*x0^2*y0+X^3+3*X^2*x0+3*X*x0^2+x0^3)-1/2*(X^2*Y^3-2*X^2*Y^2*a+3*X^2*Y^2*
y0-4*X^2*Y*a*y0+3*X^2*Y*y0^2-2*X^2*a*y0^2+X^2*y0^3+2*X*Y^3*x0-4*X*Y^2*a*x0+6*X*
Y^2*x0*y0-8*X*Y*a*x0*y0+6*X*Y*x0*y0^2-4*X*a*x0*y0^2+2*X*x0*y0^3-Y^3*a^2+Y^3*x0^
2-3*Y^2*a^2*y0-2*Y^2*a*x0^2+3*Y^2*x0^2*y0-3*Y*a^2*y0^2-4*Y*a*x0^2*y0+3*Y*x0^2*
y0^2-a^2*y0^3-2*a*x0^2*y0^2+x0^2*y0^3+X^2*Y+X^2*y0+2*X*Y*x0+2*X*x0*y0+Y*x0^2+x0
^2*y0)/(X^2*Y*a+X^2*a*y0+2*X*Y*a*x0+2*X*a*x0*y0+Y*a*x0^2+a*x0^2*y0+X^3+3*X^2*x0
+3*X*x0^2+x0^3)^2*(2*X*Y*a+2*X*a*y0+2*Y*a*x0+2*a*x0*y0+3*X^2+6*X*x0+3*x0^2))/Y/
(1/2*(3*X^2*Y^2-4*X^2*Y*a+6*X^2*Y*y0-4*X^2*a*y0+3*X^2*y0^2+6*X*Y^2*x0-8*X*Y*a*
x0+12*X*Y*x0*y0-8*X*a*x0*y0+6*X*x0*y0^2-3*Y^2*a^2+3*Y^2*x0^2-6*Y*a^2*y0-4*Y*a*
x0^2+6*Y*x0^2*y0-3*a^2*y0^2-4*a*x0^2*y0+3*x0^2*y0^2+X^2+2*X*x0+x0^2)/(X^2*Y*a+X
^2*a*y0+2*X*Y*a*x0+2*X*a*x0*y0+Y*a*x0^2+a*x0^2*y0+X^3+3*X^2*x0+3*X*x0^2+x0^3)-1
/2*(X^2*Y^3-2*X^2*Y^2*a+3*X^2*Y^2*y0-4*X^2*Y*a*y0+3*X^2*Y*y0^2-2*X^2*a*y0^2+X^2
*y0^3+2*X*Y^3*x0-4*X*Y^2*a*x0+6*X*Y^2*x0*y0-8*X*Y*a*x0*y0+6*X*Y*x0*y0^2-4*X*a*
x0*y0^2+2*X*x0*y0^3-Y^3*a^2+Y^3*x0^2-3*Y^2*a^2*y0-2*Y^2*a*x0^2+3*Y^2*x0^2*y0-3*
Y*a^2*y0^2-4*Y*a*x0^2*y0+3*Y*x0^2*y0^2-a^2*y0^3-2*a*x0^2*y0^2+x0^2*y0^3+X^2*Y+X
^2*y0+2*X*Y*x0+2*X*x0*y0+Y*x0^2+x0^2*y0)/(X^2*Y*a+X^2*a*y0+2*X*Y*a*x0+2*X*a*x0*
y0+Y*a*x0^2+a*x0^2*y0+X^3+3*X^2*x0+3*X*x0^2+x0^3)^2*(X^2*a+2*X*a*x0+a*x0^2)):
timelimit(30,solve(identity(eq_1,X),[x0,y0]));

[]

 

 

#now try same solve as in first example

eq_2:=1 = -X*((1/2/(2*X*Y*a+2*X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*Y*(2*Y
*a+2*a*y0)+1/2/(2*X*Y*a+2*X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*y0*(
2*Y*a+2*a*y0))/(-2*X*Y*a-2*X*a*y0-2*Y*a*x0-2*a*x0*y0+Y^2+2*Y*y0+a^2+y0^2)-(-a*Y
-a*y0+(2*X*Y*a+2*X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*Y+(2*X*Y*a+2*
X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*y0)/(-2*X*Y*a-2*X*a*y0-2*Y*a*
x0-2*a*x0*y0+Y^2+2*Y*y0+a^2+y0^2)^2*(-2*Y*a-2*a*y0))/Y/((-a+1/2/(2*X*Y*a+2*X*a*
y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*Y*(2*X*a+2*a*x0-2*Y-2*y0)+(2*X*Y*a
+2*X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)+1/2/(2*X*Y*a+2*X*a*y0+2*Y*a
*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*y0*(2*X*a+2*a*x0-2*Y-2*y0))/(-2*X*Y*a-2*X*
a*y0-2*Y*a*x0-2*a*x0*y0+Y^2+2*Y*y0+a^2+y0^2)-(-a*Y-a*y0+(2*X*Y*a+2*X*a*y0+2*Y*a
*x0+2*a*x0*y0-Y^2-2*Y*y0-y0^2)^(1/2)*Y+(2*X*Y*a+2*X*a*y0+2*Y*a*x0+2*a*x0*y0-Y^2
-2*Y*y0-y0^2)^(1/2)*y0)/(-2*X*Y*a-2*X*a*y0-2*Y*a*x0-2*a*x0*y0+Y^2+2*Y*y0+a^2+y0
^2)^2*(-2*X*a-2*a*x0+2*Y+2*y0)):
try
    timelimit(30,solve(identity(eq_2,X),[x0,y0]));
catch:
    print("good. No crash");
end try;

Error, (in is/duplicates:-Normal) too many levels of recursion

 

 

Download second_sovle_fail_maple_2025_1_oct_21_2025.mw

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