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] 。 集成的結果是相同類型的狀態對列表。