function x = indietro(U,y)
% x = indietro(U,y)  calcola la soluzione  x  di un sistema
% triangolare superiore  U * x = y
n = length(y);
x = y(n)/U(n,n);
for k = n-1 : -1 : 1
x = [(y(k)-U(k,k+1 : n)*x)/U(k,k); x];
end
end
