function [x]=sostit_indietro(U,b) 
% Risoluzione di un sistema lineare Ux=b con le sostituzioni 
% all'indietro: versione per righe.

n = size(U,1);
x = zeros(n,1);
x(n,1) = b(n)/U(n,n); 
for k=(n-1):-1:1 
  x(k,1)=(b(k) - sum(U(k,k+1:n).*x(k+1:n,1)'))/U(k,k); 
end; 