課題７

ルンゲクッタ法を用いて、微分方程式


dx/dt ={x3,x4,-x1/r3乗,-x2/r3乗}
x={x1,x2,x3,x4}
r=√(x1の2乗+x2の2乗)
x(0)={1,0,0,1}


の解を求めるプログラムを作成し、1(Eular法)、４次の場合の比較検討を行う。