Echelon Form

[Ariel Pacetti <apacetti@math.utexas.edu>]

Given a matrix (say with integer coefficients), is there a way to find its 
reduced Echelon form? (I looked through the manual but couldn't find this 
specific routine).
To fixed conventions, the reduced Echelon form is the usual Gaussian 
elimination with the condition that given a pivot position all entries 
above it are zero. For example:

[1 1]
[0 1]

is not Echelon form while

[1 0]
[0 1]

is. Cheers,

Ariel