我用拉格朗日方程求解出行星轨道方程，用 Scipy 求数值解并绘制出行星轨道的模拟。

图形的绘制使用的是 Pygame，Pygame 万岁。

代码。