A=matrix([[1,2,3],[7,3,1],[1,6,2]])
b=vector([20,13,0])
x = A.solve_right(b)
print x
A*x

Kreyszig-6.6-18-U (last edited 2010-12-18 12:09:47 by 10)