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