Il Linguaggio Fortran 90/95

242 Unità di Programma
– abbia un parametro formale che sia un array fittizio di forma presunta, un puntatore
oppure un target.
– faccia parte di una lista di argomenti (in questo caso l’interface block è raccomandato
ma non obbligatorio).
Il programma che segue, concepito per la risoluzione di una equazione non lineare con il classico
metodo di Newton-Raphson, rappresenta un esempio di utilizzo di funzioni esterne la cui
interfaccia sia resa esplicita attraverso l’uso di un interface block. Si ricorda che il metodo di
Newton-Raphson si basa sull’applicazione della formula iterativa:
xi+1 = xi + f(xi)
f ′ (xi)
la quale, a partire da un valore iniziale x0, produce una successione di valori x1, x2, . . . che,
in molti casi, si caratterizza per una rapida convergenza verso un’approssimazione della radice
x della funzione f(x). L’efficacia e la rapidità del metodo dipendono, comunque, anche dalla
scelta del valore iniziale e dalle caratteristiche della funzione nell’intorno della radice.
PROGRAM main
IMPLICIT NONE
INTEGER, PARAMETER :: dp=KIND(1D0)
REAL(KIND=dp) :: x, eps
INTEGER :: nmax
! LEGENDA:
! eps: approssimazione richiesta
! nmax: numero massimo di iterazioni
! x: (IN) valore di tentativo / (OUT) radice della funzione
INTERFACE
FUNCTION f(x)
IMPLICIT NONE
INTEGER, PARAMETER :: dp=KIND(1D0)
REAL(KIND=dp), INTENT(IN) :: x
REAL(KIND=dp) :: f
END FUNCTION f
!
FUNCTION df(x)
IMPLICIT NONE
INTEGER, PARAMETER :: dp=KIND(1D0)
REAL(KIND=dp), INTENT(IN) :: x
REAL(KIND=dp) :: df
END FUNCTION df
!
SUBROUTINE newton(f,df,x,nmax,eps)