Newton Raphson Method | FORTRAN 95

PROGRAM nr
IMPLICIT NONE

REAL::f,g,ea=6,es=1,x1,x0
INTEGER::c=0

PRINT *, "============================================"
PRINT *, "PROGRAM TO FIND ROOTS USING - NEWTON RAPHSON METHOD [BY - www.BottomScience.com]"
PRINT *, "============================================"

PRINT *,'INITIAL APPROXIMATION?'
READ(*,*)x0

DO WHILE(ea>es)
x1=x0-((f(x0))/(g(x0)))
ea=abs(((x1-x0)/x1)*100)
x0=x1
c=c+1
IF(c>50) EXIT
PRINT *,'CURRENT ROOT IS',X1
END DO

PRINT *,'FINAL ROOT IS',X1

END PROGRAM

REAL function f(x1)
REAL::x1

f=(x1**2)-3

RETURN
END FUNCTION

REAL function g(x2)
REAL::x2

g=2*x2

RETURN
END FUNCTION

OUTPUT

============================================
PROGRAM TO FIND ROOTS USING - NEWTON RAPHSON METHOD [BY - www.BottomScience.com]
============================================
INITIAL APPROXIMATION?
1
CURRENT ROOT IS 2.00000000
CURRENT ROOT IS 1.75000000
CURRENT ROOT IS 1.73214281
CURRENT ROOT IS 1.73205078
FINAL ROOT IS 1.73205078


------------------
(program exited with code: 0)

Press any key to continue . . .