Ask Your Question

Revision history [back]

.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]