function [x] = my_triangular_upper_01_04(a,b)
%calcola ax=b con a triangolare superiore
n=length(b);
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