# .rref() over the reals?

Hi everyone,

I just discovered sage and was trying to obtain the reduced echelon form of a real augmented matrix with the following code:

```
A = matrix(RR,[[27.6, 30.2], [3100, 6400], [250, 360]]);
B = column_matrix(RR, [162 , 23610, 1623]);
C = A.augment(B, subdivide = True);
C.rref();
```

The output is the identity matrix:

```
[1.0 0.0|0.0]
[0.0 1.0|0.0]
[0.0 0.0|1.0]
```

This is not right over the reals, as A\B yields

```
[3.90000000000000]
[1.80000000000000]
```

Is there any way I obtain the correct output with rref(), that is

```
[ 1.00000000000000 0.000000000000000| 3.90000000000000]
[0.000000000000000 1.00000000000000| 1.80000000000000]
[0.000000000000000 0.000000000000000|0.000000000000000]
```