%ODEOPTIONS itask=2 don't work with solver type="rk" or "rkf"
Reported by German Bresciano
BUG DESCRIPTION:
----------------
when ode is called with solver type="rk" or "rkf" and %ODEOPTIONS itask parameter=2, ode returns just the last calculated point instead of all mesh points.
ERROR LOG:
----------
HOW TO REPRODUCE THE BUG:
-------------------------
function ydot=f(t, y)
ydot=y^2-y*sin(t)+cos(t)
endfunction
%ODEOPTIONS=[2,0,0,%inf,0,2,500,12,5,0,-1,-1];
y=ode('adams',0,0,%pi,f)
y=ode('rk',0,0,%pi,f)
y=ode('rkf',0,0,%pi,f)
OTHER INFORMATION:
------------------