Re: Translating ODE from Matlab to Octave

[Bård Skaflestad]

On Tue, 2013-11-26 at 15:19 +0100, RICHARD Dominique wrote:

> f = @(t,y) daesystem(t,y,params);  
> [t,y] = ode15s(f,t,y0,options);
> 
> which I transformed to :

[...]

> f = @(t,y) daesystem(t,y,params);
> [y, t]= lsode(@(y,t) f, y0, t)

Try f = @(y, t) daesystem(t, y, params)

The call

        lsode(@(y,t) f, ...)

makes no sense because 'f' is a function handle so

        @(y,t) f

is a function that, for all parameter pairs (y,t) returns the original
function handle ('f') rather than the value of the function (daesystem)
at (y,t).

Alternatively, you *could* use

        [y, t] = lsode(@(y, t) f(t, y), y0, t)

if you don't want to change the definition of 'f'.


Sincerely,
-- 
Bård Skaflestad <address@hidden>
SINTEF ICT, Applied mathematics