function gaussa(A)
% GAUSSA(A) produces the RREF
% of the matrix A.
X = rref(A);
disp(['RREF of A:'])
disp(X)
end
