六个微分方程求状态方程,用递推的方式求卫星运动轨迹