Il Linguaggio Fortran 90/95

5.18 Visibilità 307
! Parametri formali
REAL,INTENT(IN) :: a, b
REAL,OPTIONAL,INTENT(IN) :: tolerance
! Interface block per le funzioni esterne argomento
INTERFACE
FUNCTION f(x)
IMPLICIT NONE
REAL :: f
REAL,INTENT(IN) :: x
END FUNCTION f
FUNCTION df(x)
IMPLICIT NONE
REAL :: df
REAL,INTENT(IN) :: x
END FUNCTION df
END INTERFACE
OPTIONAL df
! Variabili locali
REAL :: tol
! *** Sezione esecutiva ***
! Inizializzazione di tol
IF (PRESENT(tolerance)) THEN
tol = tolerance
ELSE
tol = default_tolerance
END IF
! Scelta del metodo da impiegare
IF (PRESENT(df)) THEN ! Usa il metodo di Newton
findRoot = newton((a+b)/2.,f,df,tol)
ELSE ! Usa il metodo delle secanti
findRoot = secant(a,b,f,tol)
END IF
END FUNCTION findRoot
RECURSIVE FUNCTION secant(a,b,f,tol) RESULT(root)
! *** Sezione dichiarativa ***
IMPLICIT NONE
! Tipo della funzione
REAL :: root
! Parametri formali
REAL, INTENT(IN) :: a, b, tol
! Interface block per la funzione esterna argomento
INTERFACE