Re: about lsode

[Geordie McBain]

On Fri, 2005-06-10 at 09:33 -0500, Pascal A. Dupuis wrote:
> Hello,
> a few questions about lsode:
> 1) is lsode_options() persistant ? I mean, if I set some param, will
> it stay forever of should the call be repeated each time before
> lsode() itself ?
> 2) I'm trying to integrate the function
>    xdot= a2 Il^2 + a1 Il + b3 x^3 + ...
>    where Il is a function of time, a2, a1, b3, ... are constant
>    about the jacobian: is this correct: 
>    jac = 3 b3 x^2 + 2 b2 x + ...
>    what about t ? what about terms like 2 a2 Il dIl/dt ?
> 
> Thanks in advance
> 
> Pascal    
> 

Hello.  I don't know the answer off hand, but you can find some
information that might be relevant in the upstream documentation for the
Fortran function LSODE which is in ODEPACK at netlib: 

        http://www.netlib.org/odepack/opkd-sum

        http://www.netlib.org/odepack/opkdmain.f

Hope this helps.




-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------