太陽を周回する地球の運動について見てみる。 解くべき運動方程式は時刻tでの軌道面における地球の直交座標を(x(t),y(t))として、 である。太陽は原点から動かないものとする。GMは万有引力定数と太陽質量の積である。 program ode use runge_kutta_4 implic…