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]
