小行星运动轨迹的RungeKutta法模拟

小行星运动的Runge-Kutta法模拟背景介绍由于两个恒星作用下行星运动问题没有解析解，只能用数值方法求解微分方程。但是在用一阶近似求解微分方程的时候存在严重的误差累积。当只考虑一个恒星引力影响时的