Function Reference
— Function File: [r, k] = rref (a, tol)

Returns the reduced row echelon form of a. tol defaults to eps * max (size (a)) * norm (a, inf).

Called with two return arguments, k returns the vector of "bound variables", which are those columns on which elimination has been performed.