A = [1 -2 1 3; 2 -3 4 8; 3 -8 0 6];

% rref(A)
A(2,:) = A(2,:)-2*A(1,:);
A(3,:) = A(3,:)-3*A(1,:);

A(3,:) = A(3,:) + 2*A(2,:);

A(1,:) = A(1,:)-A(3,:);
A(2,:) = A(2,:)-2*A(3,:); 

A(1,:) = A(1,:) + 2*A(2,:);

A