Il Linguaggio Fortran 90/95

360 Array Processing
ELSE
m(k,k:n+1) = m(k,k:n+1)/m(k,k)
m(k+1:n,k+1:n+1) = m(k+1:n,k+1:n+1) - &
SPREAD(m(k,k+1:n+1),1,n-k)*SPREAD(m(k+1:n,k),2,n-k+1)
END IF
END DO triang_loop
! Fase di sostituzione all’indietro
IF (error) THEN
x = 0.0
ELSE
DO k=n,1,-1
x(k) = m(k,n+1)-SUM(m(k,k+1:n)*x(k+1:n))
END DO
END IF
DEALLOCATE(m,temp_row)
END IF
END SUBROUTINE solve_linear_equations
END MODULE solve_module
PROGRAM p_solve
USE solve_module
REAL, DIMENSION(3,3) :: A
REAL, DIMENSION(3) :: b, x
INTEGER :: i, j
LOGICAL :: error
! costruzione di una matrice di coefficienti di prova
DO i=1,3
DO j=1,3
a(i,j) = i+j
END DO
END DO
A(3,3) = -A(3,3)
! costruzione della colonna dei termini noti
b = (/20,26,-4/)
CALL solve_linear_equations(A,x,b,error)
IF (error==.FALSE.) THEN
PRINT*, "La soluzione e’: ", x
ELSE
PRINT*, "Errore nella matrice dei coefficienti"
END IF
END PROGRAM p_solve