Runge-kutta method를 이용한 상미분방정식의 풀이

소개글

Runge-kutta method 를 이용하여 포트란으로 코딩하고
상미분방정식 dy/dx=x^2+y,dy/dx=2y/(x+1) 풀이

목차

1. 변수 설명표
2. 순서도
3. 코딩
4. 변수 설명표
5. 순서도
6. 코딩
7. 실행결과
8. 변수 설명표
9. 순서도
10. 코딩
11. 실행 결과

본문내용

변수 설명표
변수명
내용
M
1 이면 2계 (b=0.5)로 계산
2 이면 2계 (b=1)로 계산
3 이면 3계로 계산
4 이면 4계로 계산
X0, XF
X의 하한 및 상한치 [x0,xf]
Y0
Y의 초기값(y0)
H0
초기 미소증분
HOUT
출력간격
TOL
허용오차
X, Y
출력시의 x 및 y값
H
허용오차 이내의 y값을 얻기 위한 미소증분
ITMAX
최대 반복횟수
EPS
허용오차
순서도
Runge-Kutta 법

참고 자료

C및 fortran에 의한 컴퓨터 수치해석