function x=sup_solve(A, b) n = length(b); x=zeros(n, 1); for i=n:-1:1 %numeratore=b(i); %for j=i+1:n % numeratore=numeratore-A(i, j)*x(j); %end numeratore=b(i)-A(i, i+1:n)*x(i+1:n); x(i)=numeratore/A(i,i); end