sistuspf.f90

subroutine sistuspf(a, b, u, ip)

!Resolución del sistema triangular superior salvo permutación de filas ip
!Entrada:	a (matriz triangular superior salvo permutación de filas)
!		b (términos independientes)
!		ip (puntero de las permutaciones)
!Salida:	b (solución del sistema)


implicit none
real(8), intent (in) :: a(:,:), b(:)
integer, intent (in) :: ip(:)
real(8), intent (out) :: u(:)
real(8) :: sum
integer :: i, n, j

n = size(b)

do i = n, 1, -1
	sum = 0d0
	do j = i+1, n
		sum = sum+a(ip(i), j)*u(j)
	end do
	u(i) = (b(ip(i))-sum)/a(ip(i), i)
end do

return
end subroutine sistuspf