오일러 방법

오일러 방법(Euler's Method)은 수치해석학에서 상미분방정식(ODE)의 초기값 문제를 근사적으로 푸는 가장 기본적인 1차 수치적 방법이다. 18세기 스위스의 수학자 레온하르트 오일러(Leonhard Euler)의 이름에서 유래하였다.

개념 및 원리

오일러 방법은 테일러 급수(Taylor series) 전개에서 유도된다. 주어진 미분방정식이 다음과 같은 형태일 때,

$$ \frac{dy}{dt} = f(t, y), \quad a \leq t \leq b, \quad y(a) = \alpha $$

구간 $[a, b]$를 $N$개의 등간격으로 나누어 $t_i = a + ih$ ($h = (b-a)/N$)라 하면, 테일러 급수 전개로부터 다음 관계를 얻는다.

$$ y(t_{i+1}) = y(t_i) + h f(t_i, y(t_i)) + \frac{h^2}{2} y''(\xi_i) $$

여기서 2차 이상의 항(나머지 항)을 제거하면 다음과 같은 근사식을 얻으며, 이것이 오일러 방법의 핵심이다.

$$ \omega_0 = \alpha, \quad \omega_{i+1} = \omega_i + h f(t_i, \omega_i), \quad i = 0, 1, \dots, N-1 $$

여기서 $\omega_i$는 실제 해 $y(t_i)$에 대한 수치적 근사값이다.

기하학적 의미

기하학적으로 오일러 방법은 초기점에서 출발하여 각 점에서의 접선(미분방정식에 의해 결정된 기울기)을 따라 일정 간격 $h$만큼 전진하는 과정을 반복함으로써, 미지의 곡선을 다각형(polygonal) 근사선으로 대체하는 방법이다. 이때 $h$를 충분히 작게 설정하면 근사 곡선이 실제 곡선에서 크게 벗어나지 않는다.

오차 특성

오일러 방법은 1차 방법(global error가 $O(h)$)에 속한다. 즉, 간격 $h$를 절반으로 줄이면 전체 오차가 약 절반으로 감소한다. 이는 동일한 계산 비용 대비 더 높은 정확도를 제공하는 룽게-쿠타 방법(Runge-Kutta Method) 등 고차 방법에 비해 상대적으로 오차가 큰 편이다. 국소 절단 오차(local truncation error)는 $O(h^2)$이다.

활용 및 한계

오일러 방법은 구현이 매우 간단하고 직관적이어서 수치해석 교육 및 기초적인 시뮬레이션에 널리 사용된다. 물리학, 공학, 생물학 등 다양한 분야에서 미분방정식의 수치적 해를 구하는 데 활용된다. 그러나 오차가 상대적으로 크고, 특히 강성(stiff) 미분방정식의 경우 수치적 불안정성을 보일 수 있어 실제 고정밀 계산에서는 개선된 방법(예: 역 오일러 방법, 룽게-쿠타 방법 등)이 선호된다.

관련 방법

오일러 방법을 변형한 여러 파생 방법이 존재한다. 역 오일러 방법(Implicit Euler Method)은 수치적 안정성이 향상된 음해법이며, 반-암시적 오일러 방법(Semi-implicit Euler Method), 지수 오일러 방법(Exponential Euler Method) 등이 있다. 또한 호인의 방법(Heun's Method)이나 중간점 방법(Midpoint Method)은 오일러 방법을 개선한 2차 방법에 해당한다.