python中kepler轨道数值解的微分方程？

Question

所以，我正在尝试编写一个代码，用于解决kepler势V（r）= - 1 / r

中轨道的（我们称之为）轨道微分方程

当你进行数学运算时，你会得到一个如下所示的微分方程：

d ^ 2u / d（fi）^ 2 + u - m / M ^ 2 = 0

其中u = 1 / r 我们最终都在寻找r（fi）

现在我尝试使用数值方法解决它，首先我说du / dfi = y 然后定义一个函数（我采取了一些任意的M和m）

 def func(y,fi): 
    m=4
    M=5
    return [y[1],m/M^2-y[0]]$

并从scipy.integrate import odeint导入 然后放入

 ts = np.linspace(0,15,150)
 ys = odeint(func, y0, ts)

现在这让我得到150个两个数组的数组 我真的不明白第一个数字意味着什么，第二个数字是什么意思

ys=[fi,u(fi)] 

或其他什么？

Answer 1

订单状态为[value, derivative]。集成的结果是相同类型的状态对列表。