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

Strang-page18,22-U (last edited 2010-12-16 12:45:45 by tejaswi)