M =MatrixSpace(QQ,3,3)
A=M([0,6,13,6,0,-8,13,-8,0])
b=vector([61,-38,79])
B=matrix([b])
C=A.augment(B.transpose())
print C
C.rescale_row(2,1/13)
C.rescale_row(1,1/6)
C.add_multiple_of_row(1,2,-1)
C.add_multiple_of_row(2,1,1)
C.rescale_row(0,1/6)
C.rescale_row(1,-13/8)
C.add_multiple_of_row(0,1,1)
C.rescale_row(0,3/13)
z=C[0,3]
x=C[2,3]+C[2,2]*z
y=C[1,2]*z-C[1,3]
print x,y,z

Kreyszig-18.1-7-U (last edited 2010-12-18 10:46:55 by K.Swetha)