function [x] = my_triangular_upper_31_03(t, b)
%Implementa la sostituzione all'indietro 
% per la risoluzione di sistemi triangolari superiori
%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(n)=b(n)/t(n,n);
for k=n-1:-1:1
    s=0;
    for j=k+1:n
        s=s+t(k,j)*x(j);
    end
    x(k)=(b(k)-s)/t(k,k);
end
end