常微分方程式のプログラム例

最終更新日: 2000-11-17


できれば常微分方程式の解き方のページも見てください。 それと、私はC言語ユーザーですが、 諸般の事情により、主にFortran77版を公開しています。

振動子
減衰も外力もない、ただの振動子の数値計算プログラム。 2次のRunge-Kutta法(修正Euler法)を用いているので、 dt=0.05程度の粗い刻みで10周期にわたって計算を続けても、 生じる誤差は0.1%にも満たない。
強制振動子 (応用数学演習I 夏休み課題 1997)
2階Runge-Kutta法(修正Euler法)を用いて、 強制振動の運動方程式を解く。 oscil.fforce1.fに分かれている。 外力の種類を変更するには、 たとえば force1.fforce2.fに差し替えて コンパイルしなおす。

[「常微分方程式」に戻る] [「数値計算プログラム例」に戻る]