residuo.f90

subroutine residuo(aa, bb, u)

!Cálculo del residuo para comprobar la solución obtenida
!Entrada:	aa (matriz de coeficientes original)
!		bb (términos independientes originales)
!		b (solución del sistema)
!Salida:	bb (residuo)

implicit none
real(8), intent (in) :: aa(:,:), u(:)
real(8), intent (inout) :: bb(:)
integer :: i, n

n = size(bb)

do i= 1, n
	bb(1:n) = bb(1:n) - aa(1:n, i) * u(i)	
end do
bb=-bb

return

end subroutine