Fortran program Birge vieta

subrutine birge_vieta (a, m, x0, n, tol, Raiz, clave)! -------------------------------------------------- ----------------- IMPLICIT NONEINTEGER, hensigt (IN) :: m! Grado del polinomioREAL, DIMENSION (0: m), hensigt (IN) :: en! Vector DE M + 1 Elementos! de los coef. del polinomio REAL, hensigt (IN) :: x0! Aproximación inicial a la raízREAL, hensigt (IN) :: tol! Tolerancia para el fejl absolutINTEGER, hensigt (INOUT) :: n! Número máximo de iteraciones /! Iteraciones realizadas REAL, hensigt (OUT) :: Raiz! Aproximación a la raizINTEGER, hensigt (OUT) :: clave! Clave de Éxito:! 0: Éxito! > 0: iteraciones excedidas! -------------------------------------------------- ----------------- INTEGER :: i, jREAL :: xx0, b, c! -------------------------------------------------- ----------------- xx0 = x0DO i = 1, n! -------------------------! ESQUEMA de Horner! ------------------------- B = a (m) c = a (m) DO j = m-1,1, -1 b = b * xx0 + a (j) C = C * xx0 + b ENDDO b = b * xx0 + a (0)! -------------------------! Método de Newton! ------------------------- Raiz = xx0 - b /c IF (ABS ((Raiz-xx0) /Raiz) & lt; tol) SÅ clave = 0 n = i RETURN SLUT HVIS xx0 = raizEND DOclave = 1RETURNEND subrutine birge_vieta