function y = avanti(L,b)
% y = avanti(L,b) calcola la soluzione y del sistema triangolare
% inferiore  L * y = b,  dove  L  e' una matrice triangolare
% inferiore la cui diagonale e' costituita da elementi tutti uguali a 1
y = b(1);
for k = 2 : length(b)
y = [y; b(k) - L(k, 1: k-1) * y];
end
end
