nm

11908 Reputation

20 Badges

13 years, 177 days

MaplePrimes Activity


These are replies submitted by nm

@delvin 

I want to define the expression diff(H(xi), xi, xi)/H(xi)^2 as a new variable

Then why you said in your question you wanted the coefficient of the above to be zero?

Ok,, try this. if you want to eliminate H(xi), then try solving the ode  H'/H^2=Z for H and then use dchange to chnage the dependent variable H to Z.

The new ode will be in term of Z(xi) only and no H(xi) in it. 

Again, do not know if this is what you meant. If not, then may be someone else can better guess than I am able to.

restart;

df := diff(H(xi), xi, xi)= 2*diff(H(xi), xi)^2/H(xi) + A *H(xi)^2 + B*(diff(H(xi), xi)/H(xi))^2 + C*diff(H(xi), xi):
eq1:= -216*a1*(2*diff(H(xi), xi)^2/H(xi))+2*a3*(2*diff(H(xi), xi)^2/H(xi) + A + B*diff(H(xi), xi)^2/H(xi)^2 + C*diff(H(xi), xi))*c^2/H(xi)^2+4*diff(H(xi), xi)^3/H(xi)^2+5*diff(H(xi), xi)^2/H(xi)^3+7*diff(H(xi), xi)^3/H(xi)^5:         
eq2 := expand(subs(df, eq1));
 

-432*a1*(diff(H(xi), xi))^2/H(xi)+4*a3*c^2*(diff(H(xi), xi))^2/H(xi)^3+2*a3*c^2*A/H(xi)^2+2*a3*c^2*B*(diff(H(xi), xi))^2/H(xi)^4+2*a3*c^2*C*(diff(H(xi), xi))/H(xi)^2+4*(diff(H(xi), xi))^3/H(xi)^2+5*(diff(H(xi), xi))^2/H(xi)^3+7*(diff(H(xi), xi))^3/H(xi)^5

ode:=diff(H(xi), xi)/H(xi)^2=Z(xi);

(diff(H(xi), xi))/H(xi)^2 = Z(xi)

sol:=dsolve(ode);

H(xi) = 1/(Int(-Z(xi), xi)+c__1)

PDEtools:-dchange({sol},eq2,[Z(xi)])

-432*a1*Z(xi)^2/(Int(-Z(xi), xi)+c__1)^3+4*a3*c^2*Z(xi)^2/(Int(-Z(xi), xi)+c__1)+2*a3*c^2*(Int(-Z(xi), xi)+c__1)^2*A+2*a3*c^2*B*Z(xi)^2+2*a3*c^2*C*Z(xi)+4*Z(xi)^3/(Int(-Z(xi), xi)+c__1)^4+5*Z(xi)^2/(Int(-Z(xi), xi)+c__1)+7*Z(xi)^3/(Int(-Z(xi), xi)+c__1)

 

 

Download repl_may_19_2025_V2.mw

@delvin 

 I want to remove the derivatives in eq 2 and the new function will not have any derivative terms at all.

If this is what you want, then why not just set diff(H(xi), xi) to zero in each term that has it? Looking at your eq2, this will leave only 2*a3*c^2*A/H(xi)^2 term.

No problem for me on Linux, Maple 2025 with ST 17.

I also was able to stop the kernel by clicking on lower left red botton. Timelimit also works on this.

worksheet below

interface(version);

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

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 17 and is the same as the version installed in this computer, created May 5, 2025, 12:37 hours Eastern Time.`

restart;

try
   timelimit(30,simplify(tan(Pi/5) + sin(Pi/15) - sqrt(15)));
catch:
   print(lastexception);
end try;

proc (x) local x1; option `Copyright (c) 1999 Waterloo Maple Inc. All rights reserved.`; `...` end proc, "time expired"

Download simplify_may_17_2025.mw

Always export to vector based format and not raster based format. So try to export to postscript or pdf and not to png or jpeg. 

this looks like a bug as help clearly says

"The convert(A, Vector) function converts the parameter A into a Vector.
"The orientation of the result can be specified via the parameter o. The default orientation is column."

Also below it says

"If A is a 1-dimensional array or Array, a list, a vector or a Vector, it is simply passed directly to the Vector constructor"

but Vector([1, 2, 3]) gives column vector.

@C_R 

I see. But if the reason it timedout for you when you run in B.M. (battery mode) is because Maple was slower and the 30 seconds was reached before entering the DLL, then this will not help me. It is then just by chance. It is like me setting timelimit of 2 seconds instead of 30 second. If you change your timelimit from 30 seconds to say 2 minutes, does it still not hang for you in B.M. mode?

What I do not understand is why timelimit does not work when code is inside DLL. Assuming for now that the hang happens inside this GMP DLL. does timelimit have no effect in this case because the code is inside dynamic linked library?

I can'f find definite answer or information on this now using google.

But there has to be a way to terminate a call, even if it is inside DLL which is faulty and hangs.

I am sure someone inside Maplesoft knows the answer for this and could comment and explain what is going on?

Also, I still go not understand why GMP was even used here for this simplify call, since everything is exact and symbolic.

@C_R 

         This time the system ran in optimized battery mode

Could you give instructions how to make Maple run in this mode? Never heard of this and googling does not show how to do it. Is this some option somewhere? I'd like to try it.

I run Maple now 2025 inside Linux which is installed inside virtual box.  

Or is this setting something done in the PC itself and not in Maple?  

I still do not understand what is going on. Why would running in "optimized battery" mode not make Maple hang on this call. Why would this make any difference in how Maple itself work?

@C_R 

great detective work. I assume this GMP is same as this https://gmplib.org/  The GNU
Multiple Precision Arithmetic Library

But why would simplify and Maple be using GMP here, since everything is symbolic and exact? 

look at parfrac in help pages

@janhardo 

Linear or non-linear odes:  is there difference for trying to get parameteric solution for both types of ode?

For linear ode's, parametric solution is not needed, since almost all such ode's can be solved directly (either explicit or implicit).

Only time to resort to parametric form of solution method is for non-linear ode's. Lets conside only first order for now. Parametric solution will have the form   

          y(x)=f(t,...)
          x(t)=g(t...)

In the above, "t" is the parameter. So we have two equations, and one parameter that holds them.

If we can eliminate the parameter "t", say using the eliminate command in Maple or most likely by solving for t as function of x from the second equation using RootOf and replace that in first equation, then we get explicit solution.

But most of the times this is not possible to do, or it can make the solution much more complicated. So in this case, some books leave the ode solution as above in parametric form. 

Also sometimes solving the ode parametric form can actually be simpler than using direct approach for nonlinear odes'. 

But the bottom line,  parametric solution is a way to make it possible to express solution of  nonlinear ode's where other direct method are not possible or produce much more complicated answers. 

Here is one example from the handbook of exact solutions

ode:=y(x)*diff(y(x),x)-y(x)=-6/25*x+2/25*A*(2*sqrt(x)+19*A+6*A^2/sqrt(x));

Book solution is

The parameter above is tau.  You see, the solution is easy to read (one line!) , but this is because it is parametric. Most such odes can be hard or not possible to rewrite as explicit or even implicit form. But possible as parametric. 

As far as I know, Maple dsolve have no option to obtain/ask for parametric solution, in the sense the above book gives. 

Maple solution for the above is below. You see it is more complicated (specially the explicit one) but on the other hand, it is not parametric, so no need to try to solve for the parameter

 

restart;

ode:=y(x)*diff(y(x),x)-y(x)=-6/25*x+2/25*A*(2*sqrt(x)+19*A+6*A^2/sqrt(x));

y(x)*(diff(y(x), x))-y(x) = -(6/25)*x+(2/25)*A*(2*x^(1/2)+19*A+6*A^2/x^(1/2))

dsolve(ode)

2*(-x^(1/2)+3*A)*((3*A^2+8*x^(1/2)*A-3*x+5*y(x))/(6*A^2-2*x^(1/2)*A+y(x)))^(3/2)*y(x)/((-(-x^(1/2)+3*A)^2/(6*A^2-2*x^(1/2)*A+y(x)))^(1/2)*A*(12*A^2+2*x^(1/2)*A-2*x+5*y(x)))-(1/2)*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+y(x)))^(1/2)*2^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*y(x))/(6*A^2-2*x^(1/2)*A+y(x)))^(1/2)*(6*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*y(x)+26*A*x^(3/2)-30*A*x^(1/2)*y(x)-3*x^2+5*x*y(x))/(6*A^2-2*x^(1/2)*A+y(x))^2)^(1/2)-5*arctanh((1/2)*(6*A^2-14*x^(1/2)*A+4*x-5*y(x))/((6*A^2-2*x^(1/2)*A+y(x))*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*y(x)+26*A*x^(3/2)-30*A*x^(1/2)*y(x)-3*x^2+5*x*y(x))/(6*A^2-2*x^(1/2)*A+y(x))^2)^(1/2))))/(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*y(x)+26*A*x^(3/2)-30*A*x^(1/2)*y(x)-3*x^2+5*x*y(x))/(6*A^2-2*x^(1/2)*A+y(x))^2)^(1/2)+c__1 = 0

dsolve(ode,useInt)

c__1+(2*(6*A^3-2*A^2*x^(1/2)+A*y(x))*Intat(exp(-25*(Int(1/(12*_a^3-4*_a^2-19*_a-3), _a))), _a = (6*x^(1/2)*A-2*x+3*y(x))/(12*A^2-4*x^(1/2)*A+2*y(x)))-3*(-(1/3)*x^(1/2)+A)*exp(-25*Intat(1/(12*_a^3-4*_a^2-19*_a-3), _a = (6*x^(1/2)*A-2*x+3*y(x))/(12*A^2-4*x^(1/2)*A+2*y(x))))*y(x))/(6*A^2-2*x^(1/2)*A+y(x)) = 0

dsolve(ode,explicit)

y(x) = RootOf(48*A^2*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(3/2)*_Z*x^(1/2)-8*A*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(3/2)*x*_Z-72*A^3*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(3/2)*_Z-12*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(3/2)*_Z^2*A+4*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(3/2)*_Z^2*x^(1/2)+10*A*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*arctanh((1/2)*(6*A^2-14*x^(1/2)*A+4*x-5*_Z)/((6*A^2-2*x^(1/2)*A+_Z)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)))*x*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)*_Z-48*A^2*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*x^(1/2)*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)*_Z-12*A*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*x*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)*_Z+40*A^2*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*arctanh((1/2)*(6*A^2-14*x^(1/2)*A+4*x-5*_Z)/((6*A^2-2*x^(1/2)*A+_Z)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)))*x^(1/2)*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)*_Z-8*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*x^(1/2)*A*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*c__1*_Z+432*A^5*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)-360*A^5*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*arctanh((1/2)*(6*A^2-14*x^(1/2)*A+4*x-5*_Z)/((6*A^2-2*x^(1/2)*A+_Z)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)))*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)+24*A^2*x^(3/2)*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)-20*A^2*x^(3/2)*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*arctanh((1/2)*(6*A^2-14*x^(1/2)*A+4*x-5*_Z)/((6*A^2-2*x^(1/2)*A+_Z)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)))*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)-210*A^3*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*arctanh((1/2)*(6*A^2-14*x^(1/2)*A+4*x-5*_Z)/((6*A^2-2*x^(1/2)*A+_Z)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)))*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)*_Z-25*A*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*arctanh((1/2)*(6*A^2-14*x^(1/2)*A+4*x-5*_Z)/((6*A^2-2*x^(1/2)*A+_Z)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)))*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)*_Z^2+60*A^4*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*arctanh((1/2)*(6*A^2-14*x^(1/2)*A+4*x-5*_Z)/((6*A^2-2*x^(1/2)*A+_Z)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)))*x^(1/2)*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)+80*A^3*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*arctanh((1/2)*(6*A^2-14*x^(1/2)*A+4*x-5*_Z)/((6*A^2-2*x^(1/2)*A+_Z)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)))*x*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)+252*A^3*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)*_Z+30*A*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)*_Z^2-72*A^4*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*x^(1/2)*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)-96*A^3*(-2*(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*((3*A^2+8*x^(1/2)*A-3*x+5*_Z)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*x*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*2^(1/2)+72*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*A^4*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*c__1+5*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*c__1*_Z^2-12*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*A^3*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*x^(1/2)*c__1+42*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*A^2*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*c__1*_Z-16*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*x*A^2*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*c__1+4*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*x^(3/2)*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*A*c__1-2*(-(27*A^4+54*A^3*x^(1/2)-72*x*A^2+45*A^2*_Z+26*A*x^(3/2)-30*A*x^(1/2)*_Z-3*x^2+5*x*_Z)/(6*A^2-2*x^(1/2)*A+_Z)^2)^(1/2)*x*(-(x-6*x^(1/2)*A+9*A^2)/(6*A^2-2*x^(1/2)*A+_Z))^(1/2)*c__1*_Z)

 


 

Download parametric_sol.mw

 

The above book has many such ode's I found which Maple 2025 dsolve can't solve, but book has solutions for in parametric form.

 

   

@janhardo 

Not sure I understand the question. If you mean how to ask Maple for implicit solution to an ode, then the option 'implicit' is used for this.

And to prevent Maple dsolve from evaluating integrals internally while it is solving an ode, the option 'useInt' is used.

@Kitonum 

I understand all that about expand.

The issue is this is done in code. Not looking at screen.

So code has to keep trying different things and then pick the smallest (in leaf count) obtained. That is why I wrote the above full_simplify() command.

My point is that simplify() should just have worked here.  One does not have to use tricks or try 100 different combinations of commands to obtain simpler result.

Simplify is the most important core command in any CAS system. And Maple is supposed to be the most powerful symbolic CAS there is. Yet, I find it is always a struggle to obtain simpler result using just the simplify command.

Maple needs a new full_simplify() command, where it tries much much harder so users do not have to resort to tricks and different options hoping to get simpler result.

@Kitonum 

Ok, will add expand to the list I have to try. 

But any idea why simplify did not do it? Given that B is "simpler" in the sense it is smaller in size. (and also more clear to read).

@Scot Gould 

There seems to be something not setup right on your Maple? I just now did the command and worked fine

@vv 

But if there is no algorithm, how would you suggest the other software did it? (hard to trace what the other software does internally).

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