function [x] = my_upper_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(n)=b(n)/a(n,n);
for k=n-1:-1:1
    s=0;
    for j=k+1:n
        s=s+a(k,j)*x(j);
    end
    x(k)=(b(k)-s)/a(k,k);
end
end