数値解析 第8回 (3) ホイン法
ホイン法のアイデア
オイラー法では、
区間左端の $x$ 座標を用いて $f(x,y(x))$ ≒ $f(x_n,y(x_n))$ と考えて、$(9.9)$ の積分を
$\dps{\int_{x_n}^{x_{n+1}}f(x, y(x))dx}$ ≒ $\dps{h \times f(x_n,y_n)}$
と近似していました。
$n=0$ の場合 
ホイン法ではまず、 $(9.9)$ の積分にある意味「台形公式」を用いて
$\dps{\int_{x_n}^{x_{n+1}}f(x, y(x))dx}$ ≒ $\dps{\frac{h}{2}\Big\{f(x_n,y(x_n)) + f(x_{n+1},y(x_{n+1})\Big\}}$
と考えます。
さらに $y(x_n)$ は $y_n$ で近似し、まだ近似値の求まっていない $y(x_{n+1})$ にはオイラー法の近似式
$\dps{\hat y = y_n + h \times f(x_n,y_n)}$
を代入します。( この $y_n$ はホイン法の近似値です。)
$n=0$ の場合 
( $\hat y$ は $y(x_1)$ の近似値なので、ふつう $f(x_1,\hat y)$ と $f(x_1,y(x_1))$ は違っています。)
ホイン法の公式
公式 $(9.12)$
$\dps{\hat y = y_n + h \times f(x_n,y_n)}$ として
$\dps{y_{n+1}=y_n + \frac{h}{2} \times \left\{f(x_n,y_n)+f(x_{n+1},\hat y)\right\}}$
$(\,n=0,1,\cdots\,)$
※ ホイン法は、オイラー法を改良した「修正オイラー法」のひとつで、この他にもホイン法と呼ばれる公式がいくつもあります。
実行例
教科書 p.193 図 9.3 では刻み幅 1/1024 で誤差は 0.00052% ほどになり、オイラー法より格段に精度が上がっています。