v1=vector([3,1,1,2]);v2=vector([1,6,2,5]);v3=vector([1,0,0,0]);v4=vector([0,9,3,6])
A=column_matrix([v1, v2, v3, v4])
w=vector([8, 2, 2, 5])
B=column_matrix([w])
aug=A.augment(B, subdivide=True)
Q=aug.rref()
print "RREF(A)= "
print Q
|
RREF(A)=
[ 1 0 0 0| 2]
[ 0 1 0 0| 1]
[ 0 0 1 0| 1]
[ 0 0 0 1|-2/3]
RREF(A)=
[ 1 0 0 0| 2]
[ 0 1 0 0| 1]
[ 0 0 1 0| 1]
[ 0 0 0 1|-2/3]
|