MaplePrimes Questions

Currently, when I have solution to an ode, say y(x)=sin(x)+_C1 and have some initial condition, then to solve for _C1,  I manually substitute the solution into the IC and replace each x by x0 and replace each derivative manually and so on.

This is because I could not find automatic way to do this. Using another software, it is possible to automate this by writing the solution using the  y -> Function[{x}, ...] syntax. But in Maple, I was not sure how to do the same.

Here is a simple made up example. 

sol:= y(x) = sin(x)+_C1;
IC := a*D(y)(x0)+c*y(x0)= b*y0+exp((D@@2)(y)(x0));

The goal is to replace the solution (which is function y(x)) into the IC, and automatically replace all its derivatives and replace x by x0 then solve for _C1 from the equation that results.

Now I do this manually like this

eval(IC,[ y(x0)=eval(rhs(sol),x=x0), 
          D(y)(x0)=eval(diff(rhs(sol),x),x=x0), 
         (D@@2)(y)(x0)= eval(diff(rhs(sol),x$2),x=x0) ])

which gives

But this is too much work.

Using the other software, I can do the above much more easily like this

sol = y -> Function[{x}, Sin[x] + C[1]]
ic = a*y'[x0] + c*y[x0] == b*y0 + Exp[y''[x0]];
ic /. sol

I looked at algsubs, dchange, or making the solution as function instead, and so on but could not emulate the y -> Function[{x}, Sin[x] + C[1]] method in Maple.

What would be similar method in Maple to do the above automatically?  May be there is already builtin function in Maple?

Variation of (a) Skin friction ∂W/∂Z​, (b) Heat Transfer ∂θ/∂Z​, and (c) Mass Transfer ∂ϕ/∂Z​ for γ=10.0, Pr=7.0, ε=1.0, Nt=0.4, and Nb=0.2.

10.0   0.03301 1.90406 0.21772
20.0   0.01212 1.90403 0.20269
30.0   0.00727 1.90402 0.19325
40.0   0.00522 1.90400 0.18645


how to get this values by solving the PDE by using the pdsolve method 

 

 # Equation 37 (continuity): Simplified

    pde1 := diff(U(X,Z), X) + diff(W(X,Z), Z) = 0;

    

    # Equation 38 (momentum): Simplified

    pde2 := U(X,Z)*diff(U(X,Z), X) + W(X,Z)*diff(U(X,Z), Z) = 

           diff(U(X,Z), Z, Z) - theta(X,Z) - phi(X,Z);

    

    # Equation 39 (energy): Simplified

    pde3 := U(X,Z)*diff(theta(X,Z), X) + W(X,Z)*diff(theta(X,Z), Z) = 

           (1/Pr)*diff(theta(X,Z), Z, Z) + 

           Nb*diff(phi(X,Z), Z)*diff(theta(X,Z), Z) + 

           Nt*diff(theta(X,Z), Z)^2;

    

    # Equation 40 (concentration): Simplified

    pde4 := U(X,Z)*diff(phi(X,Z), X) + W(X,Z)*diff(phi(X,Z), Z) = 

           (1/Sc_val)*diff(phi(X,Z), Z, Z) + 

           (Nt/Nb)*diff(theta(X,Z), Z, Z);

    

    # Boundary conditions

    bcs := {

        # At Z = 0

        U(X,0) = 0,

        W(X,0) = 0,

        theta(X,0) = 1,

        phi(X,0) = 1,

        

        # At Z = Z_max

        U(X,Z_max) = 0,

        theta(X,Z_max) = 0,

        phi(X,Z_max) = 0,

      

    };


Variation of (a)W(b) θ and (c) ϕ for different value of γ when ε = 10.0, Nb = 0.4, ε = 10.0, Sc = 0.5 and Pr = 7.0.
consider X as 0.1

I frequently receive this message after installing the most recent Maple update:

Can anyone help me with this? Thank you so much for your valuable support.

is it possible to collect using pattern? For example, given 

How to tell Maple to collect on  r^power terms to produce

This came up in another forum here  and using that other software, it is possible to ask collect to collect on pattern r^_

Is there a way in Maple to collect on all powers of r in the above? Here is worksheet

A:=r^(2*a)+r^2*(1+a+r^(2*a)) + r + a*r;
B:=(1+a)*r+(1+a)*r^2+r^(2*a)+r^(2+2*a);

r^(2*a)+r^2*(1+a+r^(2*a))+r+a*r

(1+a)*r+(1+a)*r^2+r^(2*a)+r^(2+2*a)

simplify(A-B)

0

collect(A,r)

r^2*(1+a+r^(2*a))+(1+a)*r+r^(2*a)

collect(A,r,'distributed')

r^2*(1+a+r^(2*a))+(1+a)*r+r^(2*a)

collect(A,r,'recursive')

r^2*(1+a+r^(2*a))+(1+a)*r+r^(2*a)

collect(A,r,expand)

(1+a+(r^a)^2)*r^2+(1+a)*r+(r^a)^2

collect(A,r^(n::anything))

Error, (in collect) cannot collect r^n::anything

 


 

Download collect_using_pattern.mw

Using that other software:

 

In the below plot switches between to solutions of a RootOf expression when the plot range starts at zero.

plot3d on the other hand sticks to one root.

Why is that and how to get a plot starting at zero showing only one root?

restart

a := RootOf(JacobiCN(sqrt(2)*sqrt(alpha), (1/2)*sqrt(2)*_Z)^2*_Z^2+_Z^2-2)

RootOf(JacobiCN(2^(1/2)*alpha^(1/2), (1/2)*2^(1/2)*_Z)^2*_Z^2+_Z^2-2)

(1)

allvalues(a)

RootOf(JacobiCN(2^(1/2)*alpha^(1/2), (1/2)*2^(1/2)*_Z)^2*_Z^2+_Z^2-2)

(2)

plot(a, alpha = 0 .. .5)

 

eval(a, [alpha = 1/20])

RootOf(JacobiCN((1/20)*2^(1/2)*20^(1/2), (1/2)*2^(1/2)*_Z)^2*_Z^2+_Z^2-2)

(3)

evalf(allvalues(RootOf(JacobiCN((1/20)*2^(1/2)*20^(1/2), (1/2)*2^(1/2)*_Z)^2*_Z^2+_Z^2-2)))

1.024662619, -1.024662619

(4)

_ValuesMayBeLost

true

(5)

plot3d(a)

 

NULL

Download plot_of_RootOf.mw

Hello,

I have upgraded to maple 2025, but the ui fonts are too small and very thin. I went to Files -> Options -> Interface -> Default Zoom, thet sat it up to 150%. It only changed the document area not the UI options. This solution used to work with Maple 2024. I am on ubuntu 22.04.

i did my try to sketch the best shape of graph by existing code but the 3D shape in matlab is not what i am looking and is  not intresting for this kind of plot so i want use and design a better shape of 3D plot for thus contour  i need help for that 

plot-help.mw

I do not remember if this came up before. And this is all done in code, without looking at the screen.

Given sqrt(1-cos(x)^2), Maple's simplify does not return sqrt(sin(x)^2), instead it returns ugly result csgn(sin(x))*sin(x) which is correct ofcourse, but why not just return sqrt(sin(x)^2)?  As sqrt(sin(x)^2) is much easier to read than csgn(sin(x))*sin(x).

length of sqrt(1-cos(x)^2) is 29 and length of sqrt(sin(x)^2) is 21.

What Maple simplify seems to do is simplify sqrt(1-cos(x)^2) to sqrt(sin(x)^2) internally, but instead of stopping there, it keeps going and "simplifies" sqrt(sin(x)^2)  to csgn(sin(x))*sin(x).

How to make it stop at sqrt(sin(x)^2)?

Again, this is done in code. Not interactive. code uses simplify() command on most things. 

interface(version)

`Standard Worksheet Interface, Maple 2025.2, Windows 10, November 11 2025 Build ID 1971053`

A:=sqrt(1-cos(x)^2)

(1-cos(x)^2)^(1/2)

B:=sqrt(sin(x)^2)

(sin(x)^2)^(1/2)

length(A)

29

length(B)

21

simplify(A)

csgn(sin(x))*sin(x)

simplify(A,trig)

csgn(sin(x))*sin(x)

simplify(A,size)

(1-cos(x)^2)^(1/2)

simplify(A) assuming x>0 and x<Pi

sin(x)

 

 

Download simplify_dec_3_2025.mw

That other system does it better

Is it possible to make simplify do the same in Maple? 

DO not know if this is new in Maple 2025.2 or not as I have not looked yet if it is possible to install Maple 2025.1 on my new PC given Maple 2025.2 is already installed.

Found That Maple 2025.2 gives "Error, (in dsolve) numeric exception: division by zero" when asked to solve this ode using Lie symmetry.  The error comes when it tried to 

                        Computing canonical coordinates for the symmetry 

Even if this problem was in earlier Maple versions, this ofcourse should not happen.

interface(version);

`Standard Worksheet Interface, Maple 2025.2, Windows 10, November 11 2025 Build ID 1971053`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1882 and is the same as the version installed in this computer, created 2025, December 1, 6:58 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.`

libname;

"C:\Users\me\maple\toolbox\2025\Physics Updates\lib", "C:\Users\me\maple\toolbox\2025\Maple Customer Support Updates\lib", "C:\Program Files\Maple 2025\lib"

restart;

ode:=sin(x)*diff(y(x), x, x) + (2*sin(x) - cos(x))*diff(y(x), x) + (sin(x) - cos(x))*y(x) = exp(-x);

sin(x)*(diff(diff(y(x), x), x))+(2*sin(x)-cos(x))*(diff(y(x), x))+(sin(x)-cos(x))*y(x) = exp(-x)

dsolve(ode)

y(x) = exp(arcsin(cos(x)))*c__2+exp(arcsin(cos(x)))*cos(x)*c__1+exp(arcsin(cos(x)))*(-(Int(csc(x)^2*exp(-arcsin(cos(x))-x), x))*cos(x)+Int(cot(x)*csc(x)*exp(-arcsin(cos(x))-x), x))

dsolve(ode,Lie)

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

 

 

Download division_by_zero_dsolve_maple_2025_2_on_dec_2_2025.mw

I would like to change the document and worksheet colour in Maple 2025. I would like to darken or change the color of the icons and text in the tool bar. I have searched everywhere and can't seem to locate it. As currently configured, it is hard on the eyes; my eyesight is not the best. I beleive there is a softer yellow colour that can be picked, which I believe is called classic. Thank you

I am using Maple 2025.  I have the following in my ini file:

with(LinearAlgebra):
    with(VectorCalculus):
    with(plottools):
    with(plots):
    with(DocumentTools):
    with(Units[Simple]):
    with(StringTools):
    stoperror(all):
   interface(displayprecision = 4):
    _EnvUseHeavisideAsUnitStep := true:  
    Digits := 10:

I want to set my default number formatting to "Engineering."  I cannot find where or how to do this, either in the menues or the an ini file command.  I can do it in the context panel line by line, but want it to be global and automatic.

Thanks in advance for your assistance.

I have an expression with y,y',y''. I found that I had to use collect first to make Maple simplifies it more.

Why is that? Is this expected?

interface(version);

`Standard Worksheet Interface, Maple 2025.2, Windows 10, November 11 2025 Build ID 1971053`

B:=-(-t^2+1)*(-diff(y(t),t$2)*(-t^2+1)^(1/2)+diff(y(t),t)/(-t^2+1)^(1/2)*t)-(2*(-t^2+1)^(1/2)-t)*diff(y(t),t)*(-t^2+1)^(1/2)+((-t^2+1)^(1/2)-t)*y(t);

-(-t^2+1)*(-(diff(diff(y(t), t), t))*(-t^2+1)^(1/2)+(diff(y(t), t))*t/(-t^2+1)^(1/2))-(2*(-t^2+1)^(1/2)-t)*(diff(y(t), t))*(-t^2+1)^(1/2)+((-t^2+1)^(1/2)-t)*y(t)

simplify(B);

((diff(diff(y(t), t), t))*t^4+2*(diff(y(t), t))*t^2*(-t^2+1)^(1/2)-2*(diff(diff(y(t), t), t))*t^2-y(t)*t^2-y(t)*(-t^2+1)^(1/2)*t-2*(diff(y(t), t))*(-t^2+1)^(1/2)+diff(diff(y(t), t), t)+y(t))/(-t^2+1)^(1/2)

B:=collect(B,{diff(y(t),t$2),diff(y(t),t),y(t)},'distributed'):
simplify(B);

2*(t^2-1)*(diff(y(t), t))+((-t^2+1)^(1/2)-t)*y(t)+(-t^2+1)^(3/2)*(diff(diff(y(t), t), t))


Download why_collect_is_needed_to_simplify_more_maple_2025_2_nov_30_2025.mw

Notice how much simpler the result when doing collect first.

I've seen many here recommend using the screen reader for Maple 2025. So this is first time I tried it.

Clicked on this second icon in my windows 10 start menu

Maple came up using the old UI, which is good:

But here comes the big problem. The file->Open menu does not work. It automatically opened in location which is not what I want. I wanted to navigate to different place on my PC. But anything I click on just generates loud beep and does nothing. Can't type anything to change location either.

So can't open any file. Not only that, It is stuck and can't even close the windows. Can't get out. Everything is stuck. 

It could be because I have 2 monitors. I had this open on the second monitor on the right. Everything on the second monitor now do not respond like Maple. Luckily, on the left monitor I was able to start task manager and kill Maple in order to get out of this hang.    

Basically screen reader does not work at all for me. If I can't open a file.

It could be due to using 2 monitors. I do not know.

Does file open work for others?

Here is small movie.

I could not get the max() function to work with units.  I had to strip out the units, do the max() and then add them back in.  It seems overly complicated and wrong.  I am sure there is a better way.

Could someone please check if Maple 2025.1 gives this error?  I only have Maple 2025.2 (is there a way to install Maple 2025.1 on same PC, after installing Maple 2025.2? Will Maple complain or remove Maple 2025.2 if I do this? Can I use same activation key to do this? I only have one.

Will call support and ask them on monday. This way I can check myself from now on.

Calling dsolve on an ode gives Error, (in dsolve) improper op or subscript selector which comes from 

             Computing canonical coordinates for the symmetry [x*(x^2+y^2+a), -y*(-x^2-y^2+a)]

Worksheet below. 

restart;

interface(version);

`Standard Worksheet Interface, Maple 2025.2, Windows 10, November 11 2025 Build ID 1971053`

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

libname;

"C:\Users\me\maple\toolbox\2025\Physics Updates\lib", "C:\Users\me\maple\toolbox\2025\Maple Customer Support Updates\lib", "C:\Program Files\Maple 2025\lib"

ode:=x*y(x)*diff(y(x),x)^2+(a+x^2-y(x)^2)*diff(y(x),x)-y(x)*x = 0;

x*y(x)*(diff(y(x), x))^2+(a+x^2-y(x)^2)*(diff(y(x), x))-y(x)*x = 0

dsolve(ode);

Error, (in dsolve) improper op or subscript selector

infolevel[dsolve]:=5;

5

dsolve(ode);

Methods for first order ODEs:

   *** Sublevel 2 ***

   Methods for first order ODEs:

   -> Solving 1st order ODE of high degree, 1st attempt

   trying 1st order WeierstrassP solution for high degree ODE

   trying 1st order WeierstrassPPrime solution for high degree ODE

   trying 1st order JacobiSN solution for high degree ODE

   trying 1st order ODE linearizable_by_differentiation

   trying differential order: 1; missing variables

   trying simple symmetries for implicit equations

   Successful isolation of dy/dx: 2 solutions were found. Trying to solve each resulting ODE.

      *** Sublevel 3 ***

      Methods for first order ODEs:

      --- Trying classification methods ---

      trying homogeneous types:

      trying exact

      Looking for potential symmetries

      trying an equivalence to an Abel ODE

      trying 1st order ODE linearizable_by_differentiation

   -> Solving 1st order ODE of high degree, Lie methods, 1st trial

    -> Computing symmetries using: way = 3

[x^3+x*y^2+a*x, x^2*y+y^3-a*y]

    <- successful computation of symmetries.

 1st order, trying reduction of order with given symmetries:

[x*(x^2+y^2+a), -y*(-x^2-y^2+a)]

   1st order, trying the canonical coordinates of the invariance group

   -> Computing canonical coordinates for the symmetry [x*(x^2+y^2+a), -y*(-x^2-y^2+a)]

Error, (in dsolve) improper op or subscript selector

 

 

Download dsolve_gives_parsing_error_maple_2025_2_nov_20_2025.mw

This internal error happens even if Physics or Support tools are not in libname.

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