python開普勒軌道數值解的微分方程？

Question

因此，我試圖編寫一個代碼來解開開普勒勢阱中的（我們稱之爲）軌道的微分方程，當你進行數學計算時，你將得到一個數學公式微分方程，看起來像這樣：

d^2U/d（FI）^ 2 + U - 米/ M^2 = 0

其中u = 1/R ，我們最終尋找r（fi）

和現在我嘗試使用數值方法來解決它，第一I所述的du/DFI = Y 然後definig的函數（I採取了一些任意的M和M）

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

和從scipy.integrate進口odeint進口把 然後在

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

現在這讓我兩個數字 150個數組的數組，我真的不明白dodes第一個數字是說又是什麼第二個數字的意思是是

ys=[fi,u(fi)] 

還是別的？

Answer 1

您訂購一個系統的狀態是[value, derivative]。整合的結果是一個相同類型的狀態對列表。