function [x] = my_triangular_lower_31_03(t, b)
%Implementa la sostituzione in avanti 
% per la risoluzione di sistemi triangolari inferiori
%tx=b
n=length(t);
%%%%%%%test di singolarita'%%%%%%%%%%%%%%
for k=1:n
    if(t(k,k)==0)
        disp('matrice singolare');
        return
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x=zeros(n,1);
x(1)=b(1)/t(1,1);
for k=2:n
    s=0;
    for j=1:k-1
        s=s+t(k,j)*x(j);
    end
    x(k)=(b(k)-s)/t(k,k);
end
end