x,y,z=var('x,y,z')
A=matrix(QQ,[[1,1,1],[1,3,3],[1,3,5]])
b = vector([2,0,2])
C=A.augment(b.transpose())
print C
C.add_multiple_of_row(2,1,-1)
C.add_multiple_of_row(1,0,-1)
print C
z=C[2,3]/C[2,2]
y=(C[1,3]-C[1,2]*z)/C[1,1]
x=C[0,3]-C[0,2]*z-C[0,1]*y
print x,y,z

Strang-1.3-27-U (last edited 2010-12-18 12:12:26 by Shubham)