principal.f90

program principal
use indir_mod
use gauss_mod
use sistu_mod
use residuo_mod

implicit none
real(8), allocatable :: a(:,:), aa(:,:), b(:), bb(:)
integer :: n, i

read*, n

allocate (a(n,n), aa(n,n), b(n), bb(n))

call indir(a, b)
do i = 1, n
	aa(i,1:n) = a(i,1:n)
	bb(i) = b(i)
end do

call gauss(a, b)

call sistu(a, b)
print*, b(1:n)

call residuo(aa, bb, b)
print*, bb(1:n)


deallocate (a, aa, b, bb)
stop
end program principal