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

Strang-1.3-13-U (last edited 2010-12-18 12:27:56 by Jrachit)