Python 科技研究之 04 使用 Python 简化二维三体问题模拟

n体问题是行星轨道研究中的一个显着问题。从本质上讲，使用牛顿力学分析单个物体的动力学很容易。两个物体在引力作用下相互作用的系统会变得有点困难，但我们仍然可以找到解析解。“混乱”始于向系统中添加第三个物体，并且随着我们向系统中添加更多物体，它变得更加复杂。但是废话少说，让我们动手吧。

根据牛顿运动方程，每个物体之间的力由以下二阶非线性微分方程定义：

这个等式捕捉到一个质量（表达式的左侧，F=m · a）周围的其他质量（右侧）感受到的引力。

对于等式 1，r是质量为m的物体的二维位置矢量 ( r = (x, y)) 。为了计算简单，引力常数G以及物体的质量都将设置为 1。

现在，让我们深入研究方程式 1的数值解。我们将把二阶 ODE 分解为两个一阶微分方程，其中对r相对于时间进行微分得到速度矢量v（Eq 2），对v进行类似的操作得到加速度（Eq 3），即自然是等式 1 的右侧。

如果我们在这里所做的导致的问题多于实际帮助