使用python建立物理模型并解常微分方程，普通方程组

Question

最近学校项目需要使用python重现并扩展一个物理模型，其主要方程包含5个常微分方程和7个普通方程组。
方程组如下：

其中S8、S4、S2、S、Sp、EH、EL、iH、iL、nH、nL、V为变量，共12个，且互相影响。

其他详细参数可以参考原文：https://doi.org/10.1039/C5CP05755H

之前尝试使用了scipy.intergrate.odeint函数求解，但是报错...

在文章中作者使用MATLAB对该方程组求解并且使用了jacobian矩阵（和dae函数？）结果如下：

所以想请教各位大佬有没有什么思路或者可用的python函数推荐？最好能给出一个大致的框架。谢谢！

Answer 1

from scipy.integrate import odeint  这个可以  但是方程太多了不知道行不行  之前解过一个三个等式的可以用