Add an augment method to Matrix.
Add functionality for finding reduced row echelon form.
Only sort the rows once.