A=matrix(QQ, [[0,8,2],[3,5,2],[6,2,8]])
print A
b=vector([-7,8,26])
B=matrix([b])
C=A.augment(B.transpose())
print C
C.add_multiple_of_row(2,1,-2)
print C
C.add_multiple_of_row(2,0,1)
print C
C.rescale_row(2, 1/6)
print C
C.rescale_row(0,1/8)
print C
C.rescale_row(1,1/3)
print C
z=C[2,3]
y=C[0,3]-C[0,2]*z
x=C[1,3]-C[1,1]*y-C[1,2]*z
print x,y,z

Kreyszig-18.1-1 (last edited 2010-12-17 12:26:18 by 172)