[SciPy-user] integrating a system of differential equations

Joshua Stults joshua.stults@gmail....
Wed May 27 09:16:53 CDT 2009


Massimo,

If you're having stability problems, usually going to an implicit
integration scheme will help.  Looks like the scipy.integrate.ode
class lets you choose an integration scheme based on backward
difference formulas, which should be unconditionally stable.  You'll
just need to supply a system function and a Jacobian function.

Hope that helps.

On Wed, May 27, 2009 at 8:24 AM, ms <devicerandom@gmail.com> wrote:
> Hello,
>
> I have to integrate a huge system of differential equations. The system
> is such that for each time step, the solution of equation j-1 is a
> parameter for equation j, so they have to be all integrated together at
> the same time.
>
> I tried to do it myself but it seems whatever I do is by no means stable.
>
> Can one do that using odeint? I tried to look for odeint documentation
> but it's not very clear to me if it is possible (especially for someone
> like me not being exactly accustomed to numerical resolution of ODEs).
>
> Thanks!
> Massimo
>
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>

-- 
Joshua Stults
Website: http://j-stults.blogspot.com


More information about the SciPy-user mailing list