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