Marco Caliari schrieb:
> Moreover, I noticed what I consider a bug in ode45.m (as well as in all
> the other solvers): at line 336, it should be
>
> vdelta = abs(y5 - y4);
By example:
If I solve the "Van der Pol" oscillator between [0, 2] having initial values
[2, 0] and using your suggestion with ode23.m then the stats information is
ans =
{
nsteps = 184
nfailed = 2
nfevals = 555
npds = 0
ndecomps = 0
nlinsols = 0
}
Without your suggestion I get the following results
ans =
{
nsteps = 58
nfailed = 2
nfevals = 177
npds = 0
ndecomps = 0
nlinsols = 0
}
So number of evaluations and steps increased dramatically. I can see the same
for ode45..ode78. Can you send an example where I can see results becoming
better if I use your modification?
Thanks in advance,
Thomas
-------------------------------------------------------------------------
This SF.net email is sponsored by the 2008 JavaOne(SM) Conference
Don't miss this year's exciting event. There's still time to save $100.
Use priority code J8TL2D2.
http://ad.doubleclick.net/clk;198757673;13503038;p?http://java.sun.com/javaone
_______________________________________________
Octave-dev mailing list
Octave-dev@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/octave-dev
