Fortran: False position, Newton 법

소개글

본 보고서는 Fortran으로 False position, Newton 법을 직접 인코딩 해보고 그 결과를 정리하였다.

목차

False Position Method, Question 1.2
Newton’s Second Order, Question 1.4
Newton Raphson Method, Question 1.6

본문내용

False Position Method, Question 1.2
F(X)=X*ALOG(X)-1
READ(5,*)X1,X2,ITM
EPS=1.0E-4
IT=0
10 IT=IT+1
FX1=F(X1)
FX2=F(X2)
X3=(X1*FX2-X2*FX1)/(FX2-FX1)
FX3=F(X3)
IF(IT.GT.ITM)GOTO 1000
IF(FX1*FX3)11,13,12
11 IF(ABS(X2-X3).LE.EPS)GOTO 13
X2=X3
GOTO 10
12 IF(ABS(X3-X1).LE.EPS)GOTO 13
X1=X3
GOTO 10
13 WRITE(6,200)X3,IT
200 FORMAT(1X,`ROOT=`,F10.3/1X,`ITERATION=`,I3)
STOP
1000 WRITE(6,210) IT
210 FORMAT(1X,5X,`ITERATION=`,I2,5X,`Check Error`)
STOP
END

Newton’s Second Order, Question 1.4
F(X)=X**2-11-2*EXP(-X)
D1(X)=2*X+2*EXP(-X)
D2(X)=2-2*EXP(-X)
READ(5,*)XN,ITM
EPS=1.0E-4
IT=0
10 IT=IT+1
FX=F(XN)
D1X=D1(XN)
D2X=D2(XN)
XNI=XN-FX/(D1X-D2X*FX/(2.0*D1X))
FX=F(XNI)
IF(IT.GT.ITM) GOTO 1000
IF(ABS(FX).LE.EPS) GOTO 20
XN=XNI

참고 자료

없음