function gaussab(A,b)
% GAUSS(A,b) produces the RREF
% of the augmented matrix (A|b).
X = rref(horzcat(A,b));
disp(['RREF of (A|b):'])
disp(X)
end
