簡単な常微分方程式の数値積分

簡単な常微分方程式の例として、ここでは

をとりあげる。 初期値として、x = 0 dx/dt = 1 を t=0で与えると 解は sin(t) となることは自明であるが、それを もっとも簡単な矩形公式の陽解法で積分する。

ソースコードはここである

これを、コンパイルし、 出力結果を gnuplot で図化したのち、convert (ImageMagick) で gif に変換したものが 以下の図である

戻る

1999.Oct 玉川