常微分方程式のプログラム例
最終更新日: 2000-11-17
できれば常微分方程式の解き方のページも見てください。
それと、私はC言語ユーザーですが、
諸般の事情により、主にFortran77版を公開しています。
- 振動子
- 減衰も外力もない、ただの振動子の数値計算プログラム。
2次のRunge-Kutta法(修正Euler法)を用いているので、
dt=0.05程度の粗い刻みで10周期にわたって計算を続けても、
生じる誤差は0.1%にも満たない。
- 強制振動子 (応用数学演習I 夏休み課題 1997)
- 2階Runge-Kutta法(修正Euler法)を用いて、
強制振動の運動方程式を解く。
oscil.fと
force1.fに分かれている。
外力の種類を変更するには、
たとえば force1.fを
force2.fに差し替えて
コンパイルしなおす。
[「常微分方程式」に戻る]
[「数値計算プログラム例」に戻る]