Euler/RK4 pythonのOrbiteシミュレーション

Question

私は8つの惑星（太陽が（0,0）にある）のorbiteシミュレーションのコードを開始しました 私はXとYで加速度を得ました。 XとYのスピード、まずオイラーを使い、次にRK4を使い、私は立ち往生した。私はどこから始めるべきか分からない。 numpyによって提供されるように

Answer 1

あなたが実際にベクトル演算を使用する必要があります。ここに私のコードです。

パラメータ配列を構成する方法を統一する必要があります。位置と速度を区別する明確な理由はありません。おそらく、リストを拡張する手順add_body(m,x,y,vx,vy)を専用にすることさえできます。

あなたの加速計算が間違っていたため、あなたは1体の加速を計算するコードを書き始めましたが、パラメータの不一致を生成するすべての本体の加速リストを計算するように切り替えました。あなたはゼロ除算のエラーを得ているはずです。

オイラー法は定性的誤差を非常に速く生成しますが、それによって落胆することはありません。