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