function x=Triug(A,b);
n=length(b);
x=b;
for k=n:-1:1
    x(k)=(x(k)-A(k,k+1:n)*x(k+1:n))/A(k,k);
end