function [x]=sostit_avanti(L,b) 
% Risoluzione di un sistema lineare Lx=b con le sostituzioni in 
% avanti: versione per righe.

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