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

Kreyszig-18.3 3-U (last edited 2010-12-18 05:21:31 by 10)