Question: ODESteps gives division by zero for basic first order quadrature ode

Any workaround or why ODESteps gives this internal Maple error?  Is this known limitation of ODESteps?

ps. reported also to Maplesoft just in case.

17128

interface(version);

`Standard Worksheet Interface, Maple 2024.0, Windows 10, March 01 2024 Build ID 1794891`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1750 and is the same as the version installed in this computer, created 2024, May 31, 10:47 hours Pacific Time.`

restart;

20456

ode:=diff(y(x), x) - 2*y(x) = 2*sqrt(y(x));
DEtools:-odeadvisor(ode);
 

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

[_quadrature]

ic:=y(0) = 1;
Student:-ODEs:-ODESteps([ode,ic]);

y(0) = 1

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

restart;

20456

ode:=diff(y(x), x) - 2*y(x) = 2*sqrt(y(x));
ic:=y(0) = 1;
dsolve([ode,ic]);

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

y(0) = 1

y(x) = 4*exp(2*x)-4*exp(x)+1

 

 

Download bug_in_odesteps_divide_by_zero_may_31_2024.mw

Please Wait...