結果 コード 太陽、地球、月と書いてありますが上の結果は太陽でも地球でも月でもありません。 初期条件と書かれたstateの部分には太陽、地球、月のx, y座標およびx, y方向の初期速度が入ります。 import numpy as np import matplotlib.pyplot as plt from matplotlib.animation import FuncAnimation import math # ルンゲクッタ法 def runge_kutta(h, f, y): k1 = h * f(y) k2 = h * f(y + k1 / 2) k3 = h * f(y + k2…