오일러 방법 계산기

오일러 방법을 사용하여 모든 1계 상미분 방정식(ODE) y' = f(x, y)를 수치적으로 풉니다. 반복 테이블, 방향장 위에 겹쳐진 오일러 다각형, 그리고 h, h/2, h/4에서의 실시간 수렴 비교를 확인하세요. 엄밀해(closed-form solution)에 대한 선택적 오차 분석도 제공합니다.

오일러 방법 계산기

빠른 예제:

y' = x + y y' = y − x² y' = −2xy y' = sin x − y 로지스틱 y' = y(1−y)

미분 방정식

y' =

변수로 x와 y를 사용하세요. + − × ÷ ^ 및 함수 sin, cos, tan, exp, ln, log, sqrt, abs를 지원합니다.

초기 x₀

초기 y₀

단계 크기 h

단계 수 n

정확한 해 y(x) (선택 사항 — 오차 분석 활성화)

입력 시, 모든 오일러 지점이 y(x)와 비교되며 오차가 점선 형태의 주황색 스터브로 표시됩니다.

방향장 표시 h/2 및 h/4에서 비교

Embed 오일러 방법 계산기 Widget

오일러 방법 계산기 정보

오일러 방법 계산기는 고전적인 (전진) 오일러 방법을 사용하여 \( y' = f(x, y), \; y(x_0) = y_0 \) 형식의 모든 1계 초기값 문제를 수치적으로 해결합니다. 이 도구는 전체 반복 테이블을 제공하고, 실시간 방향장 위에 오일러 다각형을 그리며, 세 가지 다른 단계 크기에서 해를 비교하여 시각적으로 수렴 과정을 확인할 수 있게 해줍니다. 또한, 정확한 폐형식 해를 입력하면 단계별 오차 분석 결과도 생성합니다.

오일러 방법이란 무엇인가요?

오일러 방법은 초기값 문제의 해를 근사하기 위한 가장 단순한 알고리즘입니다. 해 곡선 위의 알려진 점 \( (x_0, y_0) \)에서 시작하여, 국소 기울기 \( f(x, y) \)를 따라 작은 단계 크기 h만큼 반복적으로 전진합니다.

y_n+1 = y_n + h · f(x_n, y_n), x_n+1 = x_n + h

기하학적으로 각 단계는 현재 지점에서의 미분 방정식 값과 동일한 기울기를 가진 짧은 직선 선분입니다. 그 결과로 나타나는 꺾은선인 오일러 다각형은 실제 (보통 곡선인) 해의 근사치가 됩니다.

얼마나 정확한가요?

오일러 방법은 1계 방법입니다. 각 단계에서의 국소 절단 오차는 \( O(h^2) \)이며, 고정된 구간에 대해 적분한 후의 전역 오차는 \( O(h) \)입니다. 실제적으로는 다음과 같습니다.

단계 크기를 절반으로 줄이면 전역 오차도 대략 절반으로 줄어듭니다.
오차는 적분 구간의 길이에 따라 선형적으로 증가합니다.
해의 곡률이 높은 곳에서 오차가 가장 심하게 발생합니다.

내장된 단계 크기 비교(h, h/2, h/4) 기능을 통해 이러한 선형 수렴을 직접 확인할 수 있습니다. 옵션을 활성화하고 세 최종 값이 공통 한계값에 접근하는지, 각 값이 이전 값보다 한계값에서 약 절반 정도 가까워지는지 확인해 보세요.

차트 읽는 법

시각화 도구는 단일 좌표 평면에 네 가지 정보를 계층화하여 표시합니다.

회색 방향장 (Slope Field) — 해당 지점에서 \( f(x, y) \)와 동일한 경사를 가진 짧은 선분들입니다. 이는 "미분 방정식이 지시하는 흐름"이라고 생각할 수 있습니다. 모든 해 곡선은 모든 지점에서 이 장에 접해야 합니다.
인디고 오일러 다각형 — 단계별 수치 해입니다. 각 선분은 이전 그리드 지점에서 시작하여 \( f(x_n, y_n) \) 방향으로 거리 h만큼 뻗어 나갑니다.
녹색 점선 정확한 해 곡선 — 폐형식 해를 입력했을 때만 나타납니다. 수직 방향의 주황색 점선 스터브는 부호가 있는 국소 오차 \( y_n - y_{\text{exact}}(x_n) \)를 나타냅니다.
주황색 및 녹색 비교 곡선 — 단계 크기 비교가 활성화되었을 때 표시되는 h/2 및 h/4에서의 동일한 문제 실행 결과입니다.

이 계산기 사용 방법

우변 입력: y' = 표시가 된 필드에 미분 방정식을 입력하세요. x와 y를 변수로 사용합니다. 지원되는 연산자는 + − × ÷ ^이며, 지원되는 함수에는 sin, cos, tan, asin, acos, atan, exp, ln, log, log10, log2, sqrt, abs 등이 포함됩니다.
초기 조건 설정: 시작 값 x₀, 해당 지점에서의 초기 y₀, 단계 크기 h(순방향 적분은 양수, 역방향 적분은 음수), 단계 수 n을 설정합니다.
(선택 사항) 정확한 해 제공: y(x)를 알고 있다면 입력하세요. 계산기가 매 단계에서 \( |y_n - y(x_n)| \)을 계산하고 최대 오차 및 최종 오차를 보고합니다.
시각화 옵션 전환: 방향장은 기본적으로 켜져 있습니다. 단계 크기 비교는 h/2 및 h/4에서의 두 곡선을 추가로 겹쳐서 보여줍니다.
실행 클릭: 결과 섹션에 요약 통계, 차트, 수렴 비교 패널 및 전체 반복 테이블이 표시됩니다. 테이블 행에 마우스를 올리면 차트의 해당 지점이 강조 표시됩니다(반대의 경우도 마찬가지).

계산 예시

\( y' = x + y, \; y(0) = 1 \)이고 h = 0.1로 10단계를 수행하는 경우를 가정해 봅시다. 정확한 해는 \( y(x) = -x - 1 + 2e^x \)입니다. 오일러 공식을 적용하면 다음과 같습니다.

n = 1: y₁ = 1 + 0.1 · (0 + 1) = 1.1 n = 2: y₂ = 1.1 + 0.1 · (0.1 + 1.1) = 1.22 n = 3: y₃ = 1.22 + 0.1 · (0.2 + 1.22) = 1.362 ⋮ n = 10: y₁₀ ≈ 3.1875, 정확한 y(1) = 2e − 2 ≈ 3.4366

최종 오차는 약 0.249입니다. h를 0.05로 절반 줄이면 최종 오차는 약 0.13으로 떨어지고, 다시 0.025로 줄이면 약 0.067로 떨어집니다. 이는 이론이 예측하는 대로 깔끔한 선형 수렴을 보여줍니다.

오일러 방법 vs 기타 수치 해석 방법

방법	차수	단계당 평가 횟수	전역 오차	참고 사항
오일러 (전진)	1	1	O(h)	가장 단순한 방법. 교육 및 프로토타이핑에 적합.
개선된 오일러 (호인)	2	2	O(h²)	단계의 시작과 끝에서의 기울기를 평균함.
중점 방법 (RK2)	2	2	O(h²)	각 단계의 중간 지점에서 기울기를 평가함.
룬게-쿠타 4 (RK4)	4	4	O(h⁴)	범용 솔버의 주력 모델. 단계당 매우 높은 정확도.
후진 (암시적) 오일러	1	1 (근 구하기 포함)	O(h)	무조건적 안정성. 뻣뻣한(stiff) ODE에 필수적.

오일러 방법이 잘못될 수 있는 경우

전진 오일러 방법은 다음 세 가지 상황에서 오작동할 수 있습니다.

단계 크기가 너무 큼 — 다각형이 진동하거나 발산합니다. 해결책은 h를 줄이는 것이며, h, h/2, h/4 비교 기능을 통해 이를 즉시 확인할 수 있습니다.
뻣뻣한(stiff) 미분 방정식 — 빠르게 감쇠하는 모드와 천천히 감쇠하는 모드가 공존하는 방정식은 안정성을 위해 h를 매우 작게 설정해야 합니다. 이 경우 암시적(후진 오일러) 또는 BDF 방법으로 전환하십시오.
f(x, y)의 특이점 — 0으로 나누기, 음수의 sqrt, 또는 0 이하의 ln 등은 적분을 중단시킵니다. 계산기는 문제가 된 단계를 명확하게 보고합니다.

일반적인 응용 분야

물리학 — 1계 연립 방정식으로 표현된 뉴턴의 제2법칙, 방사성 붕괴 \( \dot{N} = -\lambda N \), 뉴턴의 냉각 법칙.
생물학 및 역학 — 로지스틱 성장 \( \dot{y} = r\,y(1 - y/K) \), 구획 SIR 모델.
경제학 — 연속 복리 이자, 단순 솔로우 성장 모델.
화학 — 1차 반응 속도론 \( \dot{c} = -k c \).
교육 — RK4나 적응형 솔버로 넘어가기 전 수치 적분의 개념 도입.

자주 묻는 질문

오일러 방법이란 무엇인가요?

오일러 방법은 초기값 문제 y' = f(x, y), y(x0) = y0를 해결하기 위한 가장 간단한 수치적 절차입니다. 각 단계에서 y_{n+1} = y_n + h · f(x_n, y_n)을 통해 해를 전진시키며, 현재 지점의 기울기를 거리 h만큼 따라갑니다. 이는 1차 정확도를 가지며, 전역 오차는 O(h)입니다.

오일러 방법은 얼마나 정확한가요?

오일러 방법은 국소 절단 오차 O(h²)와 전역 오차 O(h)를 갖습니다. 단계 크기를 절반으로 줄이면 전역 오차도 대략 절반으로 줄어듭니다. 이 계산기에서 제공하는 h, h/2, h/4 수렴 비교를 통해 오차가 h에 따라 선형적으로 감소하는 것을 확인할 수 있습니다.

오일러 방법이 실패하는 경우는 언제인가요?

오일러 방법은 뻣뻣한 문제이거나 단계 크기가 해의 국소 곡률에 비해 너무 클 때 불안정해질 수 있습니다. 수치 해가 진동하거나 무한대로 발산하거나 실제 해에서 눈에 띄게 벗어날 수 있습니다. 대개 h를 줄이면 도움이 되며, 뻣뻣한 방정식의 경우 후진 오일러와 같은 암시적 방법이 선호됩니다.

단계 크기는 어떻게 선택해야 하나요?

관심 구간에 대해 약 10~50단계가 되는 h로 시작하세요. 오일러 다각형이 방향장이나 정확한 해에서 눈에 띄게 벗어나면 h를 절반으로 줄이고 다시 실행하세요. 내장된 h, h/2, h/4 비교 기능을 사용하여 세 곡선이 서로 수렴하는지 확인하세요.

오일러 방법과 룬게-쿠타(RK4)의 차이점은 무엇인가요?

4차 룬게-쿠타 방법은 단계당 4개 지점에서 기울기를 평가하고 (1, 2, 2, 1)/6의 가중치로 결합하여 전역 오차 O(h⁴)를 제공합니다. 이는 동일한 단계 수에서 오일러의 O(h)보다 수 차례 더 정밀합니다. 오일러 방법은 수치 적분의 개념을 교육하거나 매우 단순한 응용 분야에서 여전히 가치가 있습니다.

연립 미분 방정식에도 사용할 수 있나요?

이 계산기는 단일 스칼라 1계 상미분 방정식 y' = f(x, y)을 처리합니다. 연립 방정식이나 고계 미분 방정식의 경우, 방정식을 1계 연립 방정식으로 다시 작성하여 전용 솔버를 사용하거나, 2계 방정식을 두 개의 1계 방정식으로 변환하여 성분별로 해결해야 합니다.