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-1.4-28-U (last edited 2010-12-18 12:07:53 by 10)