オイラー法

オイラー法（ーほう）(Euler's Method) とは、1階常微分方程式の数値解法の一つ。

1階常微分方程式 $\frac{dy}{dx} = f(x,y)$ と初期値x₀,y₀(=Y₀)が与えられた場合に、x_n=x_n-1+h対するy_nの近似値Y_nは、以下のようにあらわされる。（ここでhは、0<hで小さい数である。）

Y n = Y n - 1 + h f (x n - 1, Y n - 1)

（n=1,2,3,...）

この方法は、数学的に理解しやすく、プログラム的にも簡単なので、数値解析の初歩的な学習問題としてよく取りあげられる。しかし、1階常微分方程式の数値解法としては精度が悪く、元の微分方程式によってはいかなるhをとっても元の方程式の解に収束しないこともある方法なので、学習目的以外であまり使われない。

方程式によってはhを小さく取れば、計算時間はかかるもののある程度は精度が出る、しかし、極端にhを小さく取ると丸め誤差が累積してきて逆に精度が落ちることもあり得る。

[編集] 説明

微分の定義から

$\frac{dy}{dx} = f(x,y)= \lim_{h \rightarrow 0} \frac{y(x+h)-y(x)}{h}$

hが小さいときは、

$f(x,y) \simeq \frac{y(x+h)-y(x)}{h}$

すなわち

$y(x+h) \simeq y(x)+hf(x,y)$

x,yの初期値をx₀,y₀とし、 x₁=x₀+h,x₂=x₁+h,...y₁=y(x₁),y₂=y(x₂),...などとすると、

$y_{n} \simeq y_{n-1}+hf(x_{n-1},y_{n-1})$

(n=1,2,3,...)

ここで

Y n = Y n - 1 + h f (x n - 1, Y n - 1)

とおくと、Y_nは、y_nの近似値になる。