R2=matrix([[10,1,1],[1,10,1],[1,1,10]])
print R2
p2=vector([6,6,6])
x2=R2.solve_right(p2)
print x2

Strang-18.3 3-U (last edited 2010-12-17 12:25:32 by 172)