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