A = matrix(4,4,[1,-1,0,0,0,1, -1, -1,0, 1, 3, 0,0, 1, 0, 6])
b = vector([0,0,9,9])
print "RREF([A: b]) ="
print A.augment(b).rref()
print "i = "
print A.solve_right(b)
|
RREF([A: b]) =
[1 0 0 0 3]
[0 1 0 0 3]
[0 0 1 0 2]
[0 0 0 1 1]
i =
(3, 3, 2, 1)
RREF([A: b]) =
[1 0 0 0 3]
[0 1 0 0 3]
[0 0 1 0 2]
[0 0 0 1 1]
i =
(3, 3, 2, 1)
|