この解は楕円軌道だと仮定しを代入して計算してやればすぐに半径1、中心0の円であることがわかります。ソースは変数が増えるのでややこしくなりますがこうです。0;}このソースでの計算はを繰り返すだけです。メッシュの幅は0.2、t=20まで計算しました。赤い線が理論値で緑が計算結果です。これははっきりいって悲惨ですね。外に膨らんでいくのは変わりません。これはオイラー法の欠陥なのです。年後を計算すれば太陽系から離脱してしまいます。