标签: python jupyter-notebook ipython runge-kutta
求解微分方程的耦合系统时,我必须求助于Dormand-Prince方法,因为我需要一种自适应阶跃方法(使用四阶Runge-Kutta方法,积分在积分间隔即将结束时会爆炸) 。
但是,我还必须为步长选择一个最大值和最小值,我不知道选择哪个值。我的初始步长为h = 0.002133404446814894,公差为0.001。
我认为没有必要编写方程组,但如果您认为有帮助,我会做的。