我是Python初学者.我正在尝试切换我在matlab中拥有的某些程序.
我需要解决一个刚性ode方程,其输入都是矩阵.在matlab中我使用

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

解决方法:

对于您在Matlab中所做的大多数事情,都可以使用Python中的NumPy模块来完成.可以找到here.

您可能还会发现相关的模块SciPy也很有用.

PyDSTool也可能与您相关.它是Radau求解器的包装.

然后,您可能要尝试matplotlib进行绘图.它的工作原理与Matlab的绘图类似.

以下链接也可能有帮助：

> http://www.ews.uiuc.edu/~mrgates2/ode/
> http://wiki.python.org/moin/NumericAndScientific?action=show&redirect=SciPy
> Integrate stiff ODEs with Python