Python可以调用C。行业标准是ODEPACK中的LSODE。它是公共领域。您可以下载C版本。这些求解器非常棘手，因此最好使用一些经过良好测试的代码。
补充：确保您确实有一个僵硬的系统，即比率（特征值）相差2或3个数量级以上。同样，如果系统是刚性的，但您仅在寻找稳态解，则这些求解器为您提供了代数求解某些方程的选项。否则，像DVERK这样的Runge- Kutta求解器将是一个很好的解决方案，并且更为简单。
在此处添加内容是因为它不适合注释：这来自DLSODE标头文档：
C T :INOUT Value of the independent variable. On return it
C will be the current value of t (normally TOUT).
C
C TOUT :IN Next point where output is desired (.NE. T).
同样，是的，Michaelis- Menten动力学是非线性的。不过，Aitken加速器可以使用它。（如果需要简短的解释，请首先考虑Y是标量的简单情况。运行系统可获得3个Y（T）点。通过它们拟合指数曲线（简单代数）。然后将Y设置为渐近线，重复，现在将Y推广为向量，假设这3个点在一个平面上- 否则就可以。）此外，除非您具有强制功能（例如恒定的IV滴注），否则MM消除将衰减离开，系统将接近线性。希望能有所帮助。
602392714
清零编程群