Integrate stiff ODEs with Python

Question

I\'m looking for a good library that will integrate stiff ODEs in Python. The issue is, scipy\'s odeint gives me good solutions sometimes, but the slightest change

Answer 1

PyDSTool wraps the Radau solver, which is an excellent implicit stiff integrator. This has more setup than odeint, but a lot less than PyGSL. The greatest benefit is that your RHS function is specified as a string (typically, although you can build a system using symbolic manipulations) and is converted into C, so there are no slow python callbacks and the whole thing will be very fast.