ルンゲ＝クッタ法（RK4法）を３次元へ適用 （２階微分方程式：３次元）
刻み幅dtの結果と刻み幅dt/2で2回進めた結果を比較し、エラーが許容値を
超える場合は刻み幅を半分にして計算しなおす。
又、エラーが許容値内ならそのままプロットし、次回の刻み幅は２倍にして
計算させる。
※もちろん、次回の刻み幅は２倍である必要はなく(整数である必要もなく)倍数を変えて計算させ、結果を比較することで
簡易な精度チェックとして使えると思われる。
（計算機内部が２進扱いであること、すなわち、丸め誤差にもこだわらないとして。）
（刻み幅自動調整なしの場合） ↓少しづつズレている
（刻み幅自動調整ありの場合） ↓ズレがなくなる 許容誤差：0.000001に設定
※きれいな楕円になりました。
参考図書：Ｃによる数値計算（朝倉書店）p183-184
数値計算プログラミング（岩波書店）p124-126
パソコンで見る天体の動き(地人書館) p47-52