ode86 对以下形式的常微分方程组进行积分dy/dx=f(x,y), y(x0)=y0, 使用12阶,8阶和6阶龙格-库塔公式对。 该方法使用高阶公式(使用局部外推法)进行改进。 对于比 1e-6 严格的公差,结果预计将优于 ODE45。 另见 ODE23 ODE45 和 ODEDEMO.M。
基于代码 ODE45 CB Moler,25-3-1987,MathWorks, Inc.
误差控制方法和系数取自通道Tsitouras 和 SN Papakostas,“Runge-Kutta 方法的廉价误差估计”,SIAM J. Sci。 计算。 20(1999) 2067-2088。
已测试 MATLAB 版本:6.1
2022-12-20 16:13:40
3KB
matlab
1