function x = backward_subst(U,y)
%
% Input:
%      U (matrice triangolare superiore di dimensione NxN)
%      y (vettore dei termini noti di dimensione N)
% Output:
%      x (soluzione del sistema Ux = y)

N = length(y);
x = zeros(N,1);

x(N) = y(N)/U(N,N);
for i=N-1:-1:1
      x(i) = y(i);
      for j=i+1:N
            x(i) = x(i)-U(i,j)*x(j);
      end
      x(i)=x(i)/U(i,i);
end