Il Linguaggio Fortran 90/95

320 Array Processing
matrice quadrata un valore scalare λ è detto autovalore di A se esiste un vettore non nullo u
tale che
Au = λu
Tale vettore u è chiamato autovettore di A associato a λ.
PROGRAM eigen
IMPLICIT NONE
! Programma per testare gli autovalori e gli autovettori
! di una matrice A di tipo intero
REAL, ALLOCATABLE, DIMENSION (:,:) :: A
REAL, ALLOCATABLE, DIMENSION (:) :: u
REAL, ALLOCATABLE, DIMENSION (:) :: zero
REAL :: lamda
REAL, PARAMETER :: eps=0.001
INTEGER :: i, n, stat
LOGICAL :: correct=.TRUE.,answer=.TRUE.
!
PRINT*, "Dimensione di A: "
READ(*,*) n
ALLOCATE(A(1:n,1:n),u(1:n),zero(1:n),STATUS=stat)
IF (stat/=0)
PRINT*, "Allocazione degli array fallita!"
PRINT*, "Codice dell’errore: ", stat
STOP
END IF
PRINT*, "Introdurre A per righe: "
DO i=1,n
READ(*,*) A(i,1:n)
END DO
PRINT*, "Introdurre un autovalore di A: "
READ(*,*) lamda
PRINT*, "Introdurre il corrispondente autovettore: "
READ(*,*) u(1:n)
PRINT*, "A: "
DO i=1,n
PRINT*, A(i,1:n)
END DO
FORALL (i=1:n)
A(i,i) = A(i,i) - lamda
END FORALL
zero = MATMUL(A,u)
i = 1
DO WHILE (correct.AND.i