Fortran을 이용한 여러 적분법 인토딩 예

소개글

본 보고서는 Fortran으로 사다리꼴 적분법, Romberg 적분법, Simpson 법칙, Gauss 적분법을 직접 인코딩 해보고 그 결과를 정리하였다.

목차

Question 6.1, 사다리꼴 적분공식
Question 6.2, Romberg 적분법
Question 6.3, Simpson 법칙
Question 6.4, Gauss 적분법

본문내용

Question 6.1, 사다리꼴 적분공식
TRAPEZOID RULE
F(X) = EXP(X)
OLD = 1.E35
WRITE(6,200)
200 FORMAT( /` ENTER ITMAX, A, B, EPS`/)
READ(5,*)ITMAX, A, B, EPS
IF(A .GT. B) THEN
TEMP = A
A = B
B = TEMP
ENDIF
H = B - A
ENDS = (F(A) + F(B)) / 2
SUM = 0.0
DO 10 ITN = 1, ITMAX
H = H / 2
M = 2**ITN - 1
DO 20 I = 1, M, 2
SUM = SUM + F(A + I * H)
20 CONTINUE
AREA = (ENDS + SUM) * H
IF( ABS((AREA - OLD) / AREA) .LT. EPS) GOTO 1000
OLD = AREA
10 CONTINUE
WRITE(6, 210) AREA
210 FORMAT(//` SORRY, ITMAX EXCEEDED.`//,
+ ` THE BEST APPROXIMATION TO THE AREA IS `, E15.7/)
STOP
1000 WRITE(6, 220) AREA
220 FORMAT(//` THE AREA IS `, E15.7/)
STOP
END

참고 자료

없음