M = MatrixSpace(QQ, 3,4)
A = M([1,3,1,2,2,6,4,8,0,0,2,4])
A
b=vector([1,3,1])
x = A.solve_right(b)
print x
A * x

Kreyszig-2.2 - 36(b)-U (last edited 2010-12-18 08:10:21 by Archana)