format short
A = floor(10*rand(6));
b = floor(20*rand(6,1))-10;
x = A\b;
error = A*x-b;
Aug = [A b]
U=rref(Aug)
error2 = x-U(:,7)