a=matrix([[1,1,1],[3,3,-1],[1,-1,1]])
print a
b=vector([-2,6,-1])
x=a.solve_right(b)
print x

Strang-page18,28-U (last edited 2010-12-16 12:49:28 by tejaswi)