next up previous
Next: 5 典型的なスキーム (2) 線形多段法 Up: 4.7.4 RKF45 自作プログラム開発 Previous: 4.7.4.4 testrkf1 の実行結果

4.7.4.5 testrkf1 の実行結果の分析

区間 $ [0,1]$ $ 1$ 等分した場合 ($ h=1$ ) でもそれなりに解けているのがすごい。 $ h$ が極端に大きいときと極端に小さいときを除くと、 $ h$ を半分にするごとに誤差は $ 32$ 分の $ 1$ 程度になっていて、 公式が $ 5$ 次の精度であることを裏付けている。 一方、推定される誤差は $ 4$ 次公式と $ 5$ 次公式の差を使っているので、 $ 4$ 次公式の誤差を推定したものなので、 $ h$ が半分になるごとに $ 16$ 分の $ 1$ 程度になっている。


next up previous
Next: 5 典型的なスキーム (2) 線形多段法 Up: 4.7.4 RKF45 自作プログラム開発 Previous: 4.7.4.4 testrkf1 の実行結果
桂田 祐史
2015-05-30