オイラー法・RK4 比較シミュレーター

構造動的解析（時刻歴応答解析）：地震や衝撃による建物や橋の振動をシミュレーションします。運動方程式（2階の常微分方程式）を解くために、オイラー法やより高精度なニューマークβ法、ルンゲ＝クッタ法が用いられます。計算の安定性が重要です。

宇宙機・ミサイルの軌道計算：惑星や衛星の軌道は重力の影響を受ける常微分方程式で記述されます。長期間にわたる高精度な計算が必要なため、RK4やそれ以上の高次ルンゲ＝クッタ法、あるいは可変刻み幅の手法が標準的に使われています。

電気回路の過渡現象解析：スイッチを入れた瞬間からの電流や電圧の時間変化を求めます。回路方程式は常微分方程式系になるため、SPICEなどの回路シミュレータの内部でこれらの数値積分法が駆使されています。

CFD（数値流体力学）における粒子追跡：流体中を漂う微粒子や汚染物質の動きを追跡します。粒子の速度場が与えられた時、その位置は速度を積分することで求まり、これにオイラー法やルンゲ＝クッタ法が適用されます。特に複雑な流れでは精度が重要です。

まず、「RK4は常にオイラー法より優れている」と思いがちですが、実はそうとも限りません。確かに精度は圧倒的に高いですが、1ステップあたり4回も関数$f(x, y)$を計算するため、計算コストが約4倍になります。例えば、単純な振動現象を短時間だけシミュレーションする場合、オイラー法で刻み幅を十分小さくすれば、RK4と同等の精度をより短い計算時間で得られることもあります。常に「精度 vs 計算時間」のバランスを考えましょう。

次に、「刻み幅$h$を小さくすればするほど良い」というのも落とし穴です。確かに誤差は理論通り小さくなりますが、計算ステップ数が爆発的に増え、丸め誤差が蓄積し始めます。特にオイラー法では、$h$を極端に小さくしても、あるところから誤差が逆に増大する現象が見られることがあります。このシミュレーターで$h$を0.001や0.0001と極端に小さくして試してみると、理論と実測が完全には一致しない部分を観察できるはずです。

最後に、「この比較は1次元の常微分方程式だけの話」と限定して考えないでください。実務で扱うCAE問題、例えば構造物の振動（連立微分方程式）や流体の熱対流（偏微分方程式の時間発展）では、これらの手法が時間積分の基本部品として組み込まれています。RK4の考え方は、より高次元で複雑な系にも拡張されているのです。