How to solve a stiff ode with Python?

Question

I'm a Python beginner. I'm trying to switch some programs that I have in matlab. I need solve a stiff ode equation, whose inputs are all matrices. In matlab I use

[ttT,uT] = ode23s('SST',t,fT);

Answer 1

For most things you do in Matlab, you can do them with the NumPy module in Python. It can be found here.

You might also find the related module SciPy useful as well.

PyDSTool might also be of relevance to you. It's a wrapper around the Radau solver.

Then you might like to try matplotlib for plotting. It works quite like Matlab's plotting thing.

The following links might help, too:

• http://www.ews.uiuc.edu/~mrgates2/ode/

• http://wiki.python.org/moin/NumericAndScientific?action=show&redirect=SciPy

• Integrate stiff ODEs with Python