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

Question

所以，我正在尝试编写一个代码来解决开普勒势 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] 。 集成的结果是相同类型的状态对列表。