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

Strang-page18-27-U (last edited 2010-12-16 12:45:31 by sarada)