编号类型领域子领域问题问题的数学分析算法算法逐步推理思考的数学方程式算法实现代码参数列表及数值范围及边界条件关联知识1应用研究空间物理天体力学/轨道动力学给定初始位置和速度，计算卫星在三维空间中的轨道运动基于牛顿万有引力定律，建立二阶常微分方程：$$\frac{d^2 \mathbf{r}}{dt^2} = -\frac{GM}{\mathbf{r}^3} \mathbf{r}$br其中\mathbf{r}为位置矢量，G为引力常数，M$ 为中心天体质量四阶Runge-Kutta（RK4）数值积分方法1. 定义状态向量 y=