A=[2 -5 7; 0 6 -8; 0 0 5];
b=[1 2 3]';
for kolumn= size(A,2):-1:1
    b(kolumn)=b(kolumn)/A(kolumn,kolumn);
    A(kolumn,kolumn)=1;
    for rad=kolumn-1:-1:1
        b(rad)=b(rad)-A(rad,kolumn)*b(kolumn)
        A(rad,kolumn)=0;
    end
end
b
A=[2 -5 7; 0 6 -8; 0 0 5];
b=[1 2 3]';
A\b