R1=matrix([[2,1,1],[1,2,1],[1,1,2]])
print R1
p1=vector([4,4,4])
x1=R1.solve_right(p1)
print x1

Kreyszig-18.3 2 (last edited 2010-12-18 07:30:23 by pratyureepu)