4个耦合一阶微分方程的Runge-Kutta四阶

计算科学 龙格库塔
2021-12-16 05:09:15

我必须解决 4 个耦合的一阶微分方程f(t),g(t),h(t)w(t)女巫只是t,但对于每个参考链接,假设有 3 个变量的函数(即http://www.phy.davidson.edu/FacHome/dmb/py200/RungeKuttaMethod.htm

我想知道是否有人可以在这个问题上启发我,方程式是:

dfdt=αfβf+θg(f+h)f

dgdt=ψfϕg

dhdt=ξf+μhτh+ϵw(f+h)h

dwdt=νhχh

1个回答

是一个使用scipy.integrate.odeint. 它不使用 Runge-Kutta 方法,但通过更改scipy.integrate.odeinttoscipy.integrate.ode并使用set_integrator方法将积分器设置为“dopri5”,您可以使用 Runge-Kutta 4(5) 方法来求解方程。