Il Linguaggio Fortran 90/95

310 Unità di Programma
x = findRoot(a,b,g,tolerance=1.0E-10)
serve ancora ad invocare la funzione secant, ma con un valore della tolleranza ben preciso.
Si noti che, in quest’ultimo caso, se non si fosse usata la parola chiave tolerance si sarebbe
presentata una condizione ambigua di errore in quanto il valore1.0E-10 sarebbe stato associato
al termine df (che nella lista dei parametri formali di findRoot occupa la quarta posizione)
mentre per la tolleranza si sarebbe erroneamente utilizzato il valore di default.
Di seguito si riporta un programma per testare il modulo RootFinders, una funzione di
prova di cui valutare la radice e la sua derivata:
PROGRAM Test_Root
! Scopo: verificare l’efficacia delle funzioni definite
! nel modulo RootFinders
USE RootFinders
IMPLICIT NONE
REAL a, b
REAL, PARAMETER :: tol = 1.0E-6
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
!
PRINT*, "Inserisci gli estremi sinistro e destro &
&dell’intervallo di tentativo: "
READ(*,*) a, b
PRINT*, "Metodo di Newton - La radice di f vale: ", findRoot(a,b,f,df)
PRINT*, "Metodo delle secanti - La radice di f vale: ", findRoot(a,b,f)
END PROGRAM Test_Root
FUNCTION f(x)
! Funzione esterna di cui valutare la radice
IMPLICIT NONE
REAL :: f
REAL,INTENT(IN) :: x
f = x+EXP(x)
END FUNCTION f