Reduced Row Echelon Form?



Is there a function to change a matrix to reduced row echelon form? (eg. 
Gauss-Jordan algorithm).  Similar to ECHELON(M), but not quite.  Or something 
similar to octave's \ or, I think its Mathematica's, lsolve().

Thanks,

Nate