function [x] = my_lower_sparse_03_04(a,b)
%risoluzione di sistemi triangolari inferiori bidiagonali
n=length(b);
% for k=1:n
%     if (a(k,k)==0)
%         disp('matrice singolare');
%         return
%     end
% end
x=zeros(n,1);
x(1)=b(1);
for k=2:n
    s=a(k-1)*x(k-1);
    x(k)=b(k)-s;
end
end