Il Linguaggio Fortran 90/95

6.5 Array allocabili 327
!
INTEGER,INTENT(IN) :: n
REAL,INTENT(IN),DIMENSION(n,n) :: A
END FUNCTION det
END INTERFACE
OPEN(UNIT=11,FILE=’matrix.dat’,STATUS=’OLD’,ACTION=’READ’)
READ(11,*) n
ALLOCATE(A(n,n),IOSTAT=stato)
IF (stato==0) THEN
READ(11,*) A
CLOSE(11)
! Calcolo del determinante e stampa del risultato
WRITE(*,’(1X,A22,1X,G10.5)’) "Il determinante vale: ",det(n,A)
ELSE
WRITE(*,*) "L’allocazione della matrice e’ fallita"
END IF
END PROGRAM determinante
!
!----------------------------------------------------------
!
RECURSIVE FUNCTION det(n,A) RESULT(detA)
! N.B.: Questa subroutine valuta il determinante della
! matrice quadrata A(n,n) in maniera ricorsiva secondo
! la regola di Laplace applicata alla prima riga
!
! Sezione dichiarativa
REAL :: detA ! determinante di A (valore restituito dalla funzione)
! Parametri formali
IMPLICIT NONE
INTEGER,INTENT(IN) :: n
REAL,INTENT(IN),DIMENSION(n,n) :: A ! N.B.: array fittizio
! Variabili locali
INTEGER :: j
REAL,DIMENSION(n-1,n-1) :: m ! N.B.:array automatico
REAL :: detm ! determinante del generico minore
! A e’ la matrice di cui valutare il determinante
! m e’ il minore associato al generico elemento A(1,j)
! Sezione esecutiva
IF (n==1) THEN
detA = A(1,1)
ELSE
detA = 0.0
DO j=1,n