Advanced Search

Maple 2025 Questions and Posts

These are Posts and Questions associated with the product, Maple 2025
View unanswered posts
Maple 2025 Questions and Posts Feed

why Maple 2025 hangs on this separable ode?...

Asked by: nm 10998
ode dsolve differential-equations
8 hours ago
2  0

This looks like regression in dsolve.

In Maple 2024.2, dsolve solves this with no problem and very quickly. 

In Maple 2025 it just hangs.

Any one could find why this is the case? infolevel does not show why. Below is Maple 2024.2 worksheet and Maple 2025 worksheet.  This is Maple 2024.2 NO HANG

> 

interface(version);

`Standard Worksheet Interface, Maple 2024.2, Windows 10, October 29 2024 Build ID 1872373`

> 

ode:=diff(y(x),x)+((y(x)+1)*(y(x)-1)*(y(x)-2))/(x+1)=0;
IC:=y(1)=0;

diff(y(x), x)+(y(x)+1)*(y(x)-1)*(y(x)-2)/(x+1) = 0

y(1) = 0

> 

DEtools:-odeadvisor(ode);

[_separable]

> 

infolevel[dsolve]:=5;
dsolve([ode,IC])

5

 -> Computing symmetries using: way = 3

Methods for first order ODEs:

--- Trying classification methods ---

trying a quadrature

trying 1st order linear

trying Bernoulli

trying separable

<- separable successful

y(x) = RootOf(-2048+(x^6+6*x^5+15*x^4+20*x^3+15*x^2+6*x+257)*_Z^18+(-6*x^6-36*x^5-90*x^4-120*x^3-90*x^2-36*x-1542)*_Z^12+(9*x^6+54*x^5+135*x^4+180*x^3+135*x^2+54*x+3081)*_Z^6)^6-1

 

 

Download dsolve_2024_no_hang_april_20_2025.mw

This is Maple 2025. HANGed. Had to terminate it after 15 minutes. It seems to hang on resolving initial conditions. 

> 

interface(version);

`Standard Worksheet Interface, Maple 2025.0, Linux, March 24 2025 Build ID 1909157`

> 

ode:=diff(y(x),x)+((y(x)+1)*(y(x)-1)*(y(x)-2))/(x+1)=0;
IC:=y(1)=0;

diff(y(x), x)+(y(x)+1)*(y(x)-1)*(y(x)-2)/(x+1) = 0

y(1) = 0

> 

DEtools:-odeadvisor(ode);

[_separable]

> 

infolevel[dsolve]:=5;
dsolve([ode,IC])

5

 -> Computing symmetries using: way = 3

Methods for first order ODEs:

--- Trying classification methods ---

trying a quadrature

trying 1st order linear

trying Bernoulli

trying separable

<- separable successful

 

 

Download dsolve_2025_on_linux_hangs_april_20_2025.mw

unable to evaluate the function...

Asked by: Mapleliquid 20
syntax plotting
15 hours ago
1  3

Hello everyone,

How I get a plot for this function ? What means numeric values ? 

Thanks !

 
> 

restart
> 

T[S] := 290
> 

eta := 17; lambda := 24
> 

h := .2; `&ohm;` := 2*Pi*10; R := 2
> 

T(x[2]) := eta*`&ohm;`^2*R^2*[x[2]/h-(1/2)*(x[2]/h)^2]/lambda+T[S]

eta*`&ohm;`^2*R^2*[x[2]/h-(1/2)*(x[2]/h)^2]/lambda+T[S]

 (1)
> 

plot(T(x[2]), x[2] = .1 .. .2)

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

  

 


Download TSL_bung_9.mwTSL_bung_9.mw

Bug of "MatrixExponential" Function in t...

Asked by: MathGuru 5
linearalgebra matlab matrixexponential
Yesterday at 11:23 AM
1  7

Dear Maple Support Team and Community,

    I would like to report what seems to be a bug in Maple's MatrixExponential function in the LinearAlgebra package.

    Consider the following 8×8 skew-Hermitian matrix B:
 

B := Matrix(8, 8, [[0, I, 0, -I, 0, I, 0, -I], 
                   [-I, 0, I, 0, -I, 0, I, 0], 
                   [0, -I, 0, I, 0, -I, 0, I], 
                   [I, 0, -I, 0, I, 0, -I, 0], 
                   [0, I, 0, -I, 0, I, 0, -I], 
                   [-I, 0, I, 0, -I, 0, I, 0], 
                   [0, -I, 0, I, 0, -I, 0, I], 
                   [I, 0, -I, 0, I, 0, -I, 0]]);

    Now compute the matrix exponential exp(I*t*B) using MatrixExponential

H := LinearAlgebra:-MatrixExponential(I * t * B);

However, this result is inconsistent with the matrix exponential computed using MATLAB and SageMath, both of which give the same result, different from Maple's.

I would appreciate it if someone from the development team could look into this.

Error, (in PatternMatching:-AlgStruct:-Match) stri...

Asked by: nm 10998
patmatch
Yesterday at 2:25 AM
0  1

Encountered this error using patmatch with condition. I have changed my code since then  to avoid such cases.

But do you think this is valid error? It only happens when adding conditional. 

> 

interface(version);
> 

restart;
> 

RHS:=1/2/lambda(y)*f(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*a+1/2/lambda(y)*f(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*b-1/lambda(y)*f(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*y-1/2/lambda(y)^2*D(f)(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*a*b+1/2/lambda(y)^2*D(f)(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*a*y+1/2/lambda(y)^2*D(f)(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*b*y-1/2/lambda(y)^2*D(f)(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*y^2;

(1/2)*f(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*a/lambda(y)+(1/2)*f(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*b/lambda(y)-f(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*y/lambda(y)-(1/2)*(D(f))(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*a*b/lambda(y)^2+(1/2)*(D(f))(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*a*y/lambda(y)^2+(1/2)*(D(f))(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*b*y/lambda(y)^2-(1/2)*(D(f))(RootOf(f(_Z)*a*b-f(_Z)*a*y-f(_Z)*b*y+f(_Z)*y^2+lambda(y)^2))*y^2/lambda(y)^2

> 

patmatch(RHS,F::anything*lambda(y)^(n::anything)+H::anything,'la')

true

> 

patmatch(RHS,conditional(F::anything*lambda(y)^(n::anything)+H::anything, not (_has(H,lambda(y)) or _has(n,y))),'la')

Error, (in PatternMatching:-AlgStruct:-Match) string or symbol expected for substring

 

 

Download error_patmatch_april_18_2025.mw

design question: should Maple dsolve hang on ode w...

Asked by: nm 10998
ode dsolve differential-equations
April 18 2025
1  3

I wonder what is the general view on this.

Maple tries hard to find analytical solutions by trying different algorithms. Which is very good. But the question is, should it also hang doing this? Should not there be a circuit breaker to prevent the hang?

I mean there must be a limited number of algorithms it tries. So at one point one would expect it will finish and return either no solution or the solution it found.

For this Abel ode   y'=x+y^3, which is known not to be solvable, Maple hangs on 

           > Step 2: calculating resultants to eliminate F and get candidates for 

I waited for almost one hour. Clearly this indicates a problem internally. Right?

There should be some internal checks to prevent this hang I would think.  I do not know where it actually hangs, since trace only shows the last step above.

It will good to find out the cause of the hang and add code to prevent this in a future version of Maple dsolve to make it more robust.

btw, using that another software, it returns instantly on this ode with no solution. May be the other software did not try as hard, but at least it did not hang :)

> 

restart;
> 

interface(version);

`Standard Worksheet Interface, Maple 2025.0, Linux, March 24 2025 Build ID 1909157`

> 

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1861 and is the same as the version installed in this computer, created 2025, April 10, 15:58 hours Pacific Time.`

> 

restart;
> 

infolevel[dsolve]:=5;
ode:=diff(y(x),x)=x+y(x)^3;
sol:=dsolve(ode,y(x))

5

diff(y(x), x) = x+y(x)^3

Methods for first order ODEs:

--- Trying classification methods ---

trying a quadrature

trying 1st order linear

trying Bernoulli

trying separable

trying inverse linear

trying homogeneous types:

trying Chini

Chini's absolute invariant is: (1/27)/x^5

differential order: 1; looking for linear symmetries

trying exact

trying Abel

The relative invariant s3 is: x

The first absolute invariant s5^3/s3^5 is: 1/x^5

The second absolute invariant s3*s7/s5^2 is: 0

...checking Abel class AIL (45)

...checking Abel class AIL (310)

...checking Abel class AIR (36)

...checking Abel class AIL (301)

...checking Abel class AIL (1000)

...checking Abel class AIL (42)

...checking Abel class AIL (185)

...checking Abel class AIA (by Halphen)

...checking Abel class AIL (205)

...checking Abel class AIA (147)

...checking Abel class AIL (581)

...checking Abel class AIL (200)

...checking Abel class AIL (257)

...checking Abel class AIL (400)

...checking Abel class AIA (515)

...checking Abel class AIR (1001)

...checking Abel class AIA (201)

...checking Abel class AIA (815)

Looking for potential symmetries

... changing x -> 1/x, trying again

Looking for potential symmetries

The third absolute invariant s5*s7/s3^4 is: 0

 ->         ======================================

 ->             ...checking Abel class D (by Appell)

 -> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 1

*** No disqualifying factor on F was found ***

 -> Step 2: calculating resultants to eliminate F and get candidates for C

 

 

Download why_hangs_dsolve_april_18_2025.mw

division by zero calling dsolve with Abel ode...

Asked by: nm 10998
ode dsolve differential-equations
April 16 2025
3  6

Do you think could be a bug in dsolve? 

> 

restart;
> 

interface(version);

`Standard Worksheet Interface, Maple 2025.0, Linux, March 24 2025 Build ID 1909157`

> 

Physics:-Version()

`The "Physics Updates" version in the MapleCloud is 1861 and is the same as the version installed in this computer, created 2025, April 10, 15:58 hours Pacific Time.`

> 

restart;
> 

ode:=x^2+3*x*diff(y(x),x)=y(x)^3+2*y(x);
DEtools:-odeadvisor(ode);
dsolve([ode,y(1)=1],[Abel])

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

[_rational, _Abel]

Error, (in dsolve) numeric exception: division by zero

 

 

Download dsolve_division_by_zero_abel_april_16_2025.mw

If I do not tell it to use Abel, then dsolve does not give divison by zero.

Maple 2025 on Linux

where is Maple's help page for Abel second kin...

Asked by: nm 10998
ode differential-equations documentation
April 15 2025
3  1

I can't find the help page for Abel second kind, class B. 

Maple has help page for Abel second kind, class A and Abel second kind, class C. But not for class B. 

Here is an example of Abel second kind class B

ode:=(3*t*y(t)+y(t)^2)+(t^2+t*y(t))*diff(y(t),t)=0;
DEtools:-odeadvisor(ode)

I wanted to know the difference and the transformation used for class B to make it Abel first kind.

I googled and can't find it. Also local help skips over class B.

Is this documented somewhere else?

btw, find error on the help page for class A. Transformation used is wrong. Will leave this for another question.

Strange rendering of assignement operator: Has som...

Asked by: C_R 3207
assignment
April 14 2025
3  6

I have repeatedly seen this on two Windows PCs:

The assignement operator := is rendered as a roman d

This happens after using Maple for some time.
 Exiting Maple and restart of Maple is required. Has someone noticed the same?
All on Windows 10 and for sure in screen reader mode (my default, cannot report on the new GUI).

Ribbon interface ......

Asked by: vv 13560
interface gui ribbon
April 14 2025
2  3

It seems that the new Ribbon interface has several bugs (probably an update will come soon). So, not only the Export As is not working, but I see that (at least in Windows), opening a worksheet with a large output will display the output using the Maple input font.
Just save a .mw with the content:

expand((x+1)^200);

and then open it.

Magnification in Maple worksheets and help pages...

Asked by: Jean-Michel 35
worksheet initialization zoom
April 13 2025
2  4

hello everyone,

I am a bit shortsighted even with my glasses.

So my question is easy :: is there a way that when I open a (new) worksheet the magnification is 150%.

I can do it easily by hand but I am looking for 'something' in my maple.ini file to do this?

As well with the help pages.

Because even if I have a 150% magnification and I call a help page the display is still 100 % not 150%.

A part that I would like to set for every session a Digits:16 .

In the .ini file where I must input this? A ";" or not?

Thank you veru much.

Jean-Michel

workaround for Maple unable to verify its own solu...

Asked by: nm 10998
ode dsolve series
April 13 2025
2  8

I've reported this problem to Maple many years ago but still not fixed in Maple 2025.

Maple not able to verify its own solution on some series solutions using odetest.

Any one can suggest a workaround one can use other than odetest to verify ode solution when using series? Or have an idea why this happens?

The solution is correct but for some reason odetest returns FAIL 

Below is one example of many I have

> 

restart;
> 

interface(version);

`Standard Worksheet Interface, Maple 2025.0, Linux, March 24 2025 Build ID 1909157`

> 

 

> 

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1861 and is the same as the version installed in this computer, created 2025, April 10, 15:58 hours Pacific Time.`

> 

ode:=x*diff(y(x),x$2)+y(x)=0;
IC:=y(0)=1;
maple_sol:=dsolve([ode,IC],y(x),'series');
odetest(maple_sol,[ode,IC],'series','point'=0);

x*(diff(diff(y(x), x), x))+y(x) = 0

y(0) = 1

y(x) = c__1*x*(series(1-(1/2)*x+(1/12)*x^2-(1/144)*x^3+(1/2880)*x^4-(1/86400)*x^5+O(x^6),x,6))+ln(x)*(series(-x+(1/2)*x^2-(1/12)*x^3+(1/144)*x^4-(1/2880)*x^5+O(x^6),x,6))+(series(1-(3/4)*x^2+(7/36)*x^3-(35/1728)*x^4+(101/86400)*x^5+O(x^6),x,6))

Warning, unable to compute series necessary to test the given solution

FAIL

 

 

Download series_unable_to_odetest_april_12_2025.mw

Export As, Save As...

Asked by: vv 13560
export gui save
April 10 2025
3  4

In the File menu,  the Export As dialog is missing. Only the default (.pdf) can be used; no .mpl or .tex)
The Save As dialog is incomplete (e.g. save to .mws).

Or, is the installation corrupted?

Exploring the Mapping of Vertical Lines under 𝑤 ...

Asked by: Paras31 215
complex plotting conformal
April 10 2025
0  2

In the transformation w = e^z, vertical lines in the complex z-plane of the form Re(z)=a are mapped to circles of radius e^a in the w-plane.

Can I create a Maple visualization that does the following:

  1. Displays the result side by side, showing each vertical line in the z-plane and its corresponding mapped circle in the w-plane.

  2. Does it use a different color and style for each line and its image (e.g., dashed, dotted, solid, or different colors)?

> 

restart; interface(imaginaryunit = 'i')
> 

z := 2+2*I

2+2*I

 (1)
> 

with(plots); point1 := pointplot([[Re(z), Im(z)]], symbol = solidcircle, color = blue, axes = normal, labels = ["Re(z)", "Im(z)"], title = "Complex Number Plot")

 
> 

mod_z := abs(z)

2*2^(1/2)

 (2)
> 

vector1 := arrow([0, 0], [Re(z), Im(z)], color = red, shape = double_arrow, width = 0.5e-1, border = false, head_width = .1, head_length = .1); vector2 := arrow([0, 0], [0, Im(z)], color = green, shape = double_arrow, width = 0.5e-1, border = false, head_width = .1, head_length = .1); vector3 := arrow([0, 0], [Re(z), 0], color = yellow, shape = double_arrow, width = 0.5e-1, border = false, head_width = .1, head_length = .1)
> 

display([vector1, point1, vector2, vector3], axes = normal, labels = ["Re(z)", "Im(z)"], view = [0 .. 3, -1 .. 4], scaling = constrained, title = "Modulus of complex number")

 
> 

w := 2+3*I

2+3*I

 (3)
> 

z+w

4+5*I

 (4)
> 

z-w

-I

 (5)
> 

z.w

-2+10*I

 (6)
> 

z/w

10/13-(2/13)*I

 (7)
> 

restart
> 

with(plots); a := 2; complexplot(a+I*t, t = -Pi .. Pi, title = "Vertical Line Re(z)=2")

 
> 

A := [-3, -1, 0, 1, 2]; P := [seq(complexplot(x+I*y, y = -10 .. 10, color = red), `in`(x, A))]; display(P, title = "Multiple Vertical Lines in Complex Plane")

 
> 

with(plots); B := [-3, -1, 0, 1, 2]; H := [seq(complexplot(x+I*y, x = -10 .. 10, color = blue), `in`(y, B))]; display(H, title = "Multiple Horizontal Lines in Complex Plane")

 
 

NULL

Download operations.mw

Display quality with rotated screens...

Asked by: Anthrazit 880
resolution java orientation
April 09 2025
2  12

I need to take up an issue which is quite old, but still unsolved.

https://www.mapleprimes.com/questions/234249-Text-Quality-On-Screen

The problem is that in specific display configurations with one (rotated) vertical screen and multiple horizontal screens, the quality of fonts especially is significantly lower in Maple as with a configuration with just horizontal screens.

The affect is visible in Maple when it is on the horizontal (unrotated) screens.

Here are screenshots for comparision, and the layout of my desktop.

Screen layout:

Screenshot rotated

Screenshot unrotated

It seems that the problem is present when starting Maple. Changing the rotation of the screen when Maple is running does not affect the font quality in Maple.

Other (Windows) programs are not affected at all of this changes, so this seems to be a Java issue.

Open file dialog: How to list only Maple files...

Asked by: C_R 3207
gui
April 07 2025
2  3

With the new GUI:
I get all files (*.*) listed by default.
Other file types to filter are not listed.

Is this the same on other machines?

Can I do something about it?

1 2 3 Page 1 of 3