pythonのケプラー軌道数値解の微分方程式は？

Question

私はケプラーポテンシャルV（r）= - 1/rの軌道の微分方程式を解くコードを書こうとしています。V（r）= - 1/r

M/M^2 = 0

ここで、u = 1/R と我々は -

D^2U/D（FI）^ 2 + U：このようになり微分方程式最終的にr（fi）を探します

と今は

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

とscipy.integrateインポートodeintからインポート（私はいくつかの任意のMおよびMを取った）最初のiは、関数をdefinigデュ/ DFI = Y 前記、数値的方法を用いてそれを解決しようとしましたその後、 と

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

に入れ、今、これは私に二つの数字 の150の配列の配列を取得し、私は本当に最初の番号を意味dodesかを理解し、どのような二番目の数字は意味があることは

ありません

ys=[fi,u(fi)] 

などですか？

Answer 1

ご注文のシステムの状態は[value, derivative]です。統合の結果は、同じタイプの状態ペアのリストです。