python中的僵硬颂歌

Question

你能就下列问题给出一些一般性的建议或建议吗？你知道吗

在Python中，我试图求解一个由大约100个非线性复杂的刚性ODEs组成的系统。我及时地发展我的方程，只找到一个变量的稳态解。结果表明，对于某些参数，解是非常光滑的，但对于另一些参数，解会强烈振荡，在某个时间点变得非物理（降到零下，而物理上是不可能的）。你知道吗

我使用积分器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，但没有什么真正的帮助，例如，有时我可以解决这些数值上的困难，但对于其他一些参数-不能。你知道吗