indir.f90

subroutine indir(a, b)

!Introducción de los coeficientes de la matriz y del término independiente
!Salida: 	a (matriz de coeficientes)
!		b (términos independientes)

implicit none
real(8), intent (out) :: a(:,:)
real(8), intent (out) :: b(:)
integer :: i, n

n = size(b)


do i = 1, n
	read*, a(i,1:n), b(i)
end do

return
end subroutine indir