python中的僵硬ODE

Question

您能否提供一些一般性建议或建议在哪里查看以下问题？

在 Python 中我试图解决一个大约100个非线性复杂刚性ODE的系统。我及时发展了我的方程，只找到一个变量的稳态解。事实证明，对于某些参数，解决方案非常平滑，但对于其他解决方案 - 解决方案在某个时间点振荡强烈并变得非物理（降至零以下，而物理上则不可能）。

我使用积分器 zvode 和方法 bdf ：

    r = ode(diff).set_integrator('zvode', method = 'bdf', atol = ..., retol 
       = ...)
    r.set_initial_value(initial, t0)
    while r.successful() and r.t < tend:
        output = r.integrate(r.t+dt)

我试图让我的时间变小，我试图改变 atol 和 retol ，但没有什么真正有用，即有时候我可以解决这些数字难题，但对于其他一些参数 - 不能。