【C++】四阶龙格库塔算法实现递推轨道飞行器位置速度

四阶龙格库塔法递推轨道飞行器位置速度

算法
本篇介绍求解精度更高的四阶龙格库塔法(Runge-Kutta)，其截断误差为O(h⁵)。
四阶龙格库塔法(Runge-Kutta)求解算法为：

上式是关于 y(k) 向 y(k+1) 的递推形式，可以根据初始条件按照递推关系依次求出 y(1)、y(2)、y(3)、… y(n) ，此离散序列即为微分方程的数值解。

原理
微分方程的数值解法，本质是使用数值积分来实现 y ˙ \dot{y}