HP Forums
odesolve - Printable Version

+- HP Forums (https://www.hpmuseum.org/forum)
+-- Forum: HP Calculators (and very old HP Computers) (/forum-3.html)
+--- Forum: HP Prime (/forum-5.html)
+--- Thread: odesolve (/thread-7395.html)



odesolve - roadrunner - 12-09-2016 03:54 PM

odesolve(sin(x),[x,y],[0,-1],0.01) returns [−0.999950000417]
odesolve(sin(x),[x,y],[0,-1],0.0001) returns [−0.999999995]
odesolve(sin(x),[x,y],[0,-1],0.00000001) returns [-1.]

but

odesolve(sin(x),[x,y],[0,-1],0.) returns [0]

perhaps a bug?

[attachment=4269]

-road


RE: odesolve - parisse - 12-10-2016 07:17 AM

Indeed, fixed in source.


RE: odesolve - GeJo - 11-13-2021 10:20 AM

Hi,
I have another issue with odesolve:

Entering a simple differential equation does not seem to lead to proper results, e.g
odesolve(t,[t,y],[0,0],3) should give 4.5[attachment=10025]

The following expression gives a correct result:
odesolve(y,[t,y],[0,0],3)
Gives the correct result 20.0855...

It holds true for a hardware Prime with latest software as well as for the Prime Virtual Calculator:
Prime Software Version 2.1.14596 (2021 10 1)
The same ho[0.,[1.005]]lds true for the virtual calculator Build: 2.1.14593 (2021 10 1)

I found an older thread in The HP Support forum where it is said tha it got fixed in version 11586, but for the latest version, it still or again seems to be the same.
https://h30434.www3.hp.com/t5/Calculators/limitations-of-the-odesolve-command/td-p/6515319

https://h30434.www3.hp.com/t5/Calculators/limitations-of-the-odesolve-command/td-p/6515319/page/2

[attachment=10025]


RE: odesolve - parisse - 11-14-2021 07:08 PM

The answer form [0.0,[0.005]] shows that the numerical algorithm fails. The reason is that the algorithm is unable to estimate the step by comparing errors with 2 different methods, because it's not really a differential equation and error is 0 for both methods.