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

Strang-page 17_22-U (last edited 2010-12-16 12:27:16 by sarada)