next up previous contents
Next: 4.3 計算機の精度 Up: 第4章 数値解法の基本: Euler method Previous: 4.1 Euler 法---一変数の場合

4.2 例:線形方程式の場合

とりあえず、答がわかっている場合にオイラー法がどう振舞うかということを みてみることにしよう。初期値問題

t=1 での数値解を考える。いうまでもないが、厳密解が で 与えられるわけである。

なので、

である。指数関数は定義により

と書けるので、 問題は、この「数値解」がどれくらい正しいかということであるが、

となり、の極限で数値解が厳密解に一致するというこ とが示せたことになる。上の証明は任意の係数、初期条件と積分区間に一般化 できる。

なお、多変数の場合にもまったく同様に厳密解に収束することが証明できるが、 これは計算練習ということにしておく。各自証明を試みること。



Jun Makino
Thu Aug 13 14:18:16 JST 1998