# echelon_form not giving reduced echelon form

According to the document(emphasis mine):

`echelon_form(algorithm='default', cutoff=0, **kwds)`

Return the echelon form of self.

OUTPUT:The

reduced row echelon formof self, as an immutable matrix.

Here is what I was doing:

```
sage: A = Matrix([[1,0,3,1,2],[-1,3,0,-1,1],[2,1,7,2,5],[4,2,14,0,6]])
sage: A.echelon_form()
[1 0 3 1 2]
[0 1 1 0 1]
[0 0 0 4 4]
[0 0 0 0 0]
```

While the result is in row echelon form, it's not in **reduced** row echelon form. What I expect is something like this:

```
[1 0 3 0 1]
[0 1 1 0 1]
[0 0 0 1 1]
[0 0 0 0 0]
```

Also, if this is merely a documentation bug, how can I get the reduced row echelon form of a matrix?

Cross posted at http://stackoverflow.com/questions/40...