附加的脚本使用MATLAB支持的子程序ode45, ode23或ode23,也使用一个简单的不支持的子程序ode4(可在附件和文件交换的其他地方)解决两体轨道矢量常微分方程。金宝app所有这些子程序都使用龙格-库塔算法求解ODE。脚本比较结果与精确解析解是一个适当的椭圆:因此我们可以比较子程序的准确性。我们发现,支持的子程序显示的误差范金宝app围在0.1到0.5%;而ode4的峰值误差小于10^-7%,即小于10^6倍。这些结果表明,ode4可能是第一个尝试求解ODE的龙格-库塔子程序。所附的信息包含这些结果的图和复制结果所需的所有脚本。
引用作为
约翰Keevil(2021)。ODE4给出比ODE45, ODE23, ODE23更准确的结果(//www.tatmou.com/matlabcentral/fileexchange/59044-ode4-gives-more-accurate-results-than-ode45-ode23-ode23s), MATLAB中央文件交换。检索.