[SciPy-user] integrating a system of differential equations

[Joshua Stults]

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 at 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 at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>

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