Algoritmo Fortran -

Olá Pessoal, este algoritmo em fortran 90 para encontrar as raízes de uma função pelo método de Newton-Raphson

T(x)=x-f(x) / f'(x)

--------------------------------------------------
PROGRAM newton
IMPLICIT NONE
REAL(8):: precisao, x, y
INTEGER(2):: n
PRINT*, "digite valor inicial para x"
READ*, x
n=0
precisao=1D0

DO WHILE (precisao>0.000001D0)
y = x - ((x**2-3*x+2)/(2*x-3))
precisao = ABS(x-y)
x=y
n=n+1
ENDDO
PRINT* , "A RAIZ DA FUNCAO E",y
PRINT* , "com", n, "interacoes"

END PROGRAM newton

------------