python中开普勒轨道数值解的微分方程？

所以，我正试图写一个代码来解开普勒势V（r）=-1/r中轨道的微分方程

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

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

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

现在我试着用数值方法求解，首先我说du/dfi=y 然后定义一个函数（我取任意的M和M）

并从scipy.integrate导入odeint导入 然后放进去

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

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

ys=[fi,u(fi)] 

---

还是别的什么？

订单一系统的状态是

[价值，衍生]

。集成的结果是一个相同类型的状态对列表。

您是否尝试读取

help（odeint）

？我尝试了，但不幸的是仍然没有完全理解：/“返回…数组，其中包含t中每个所需时间的y值，初始值

y0

，在第一行。”因此结果是（t，y（t））而我想得到你（t）？

ys=[fi,u(fi)]