A=matrix([[1,1,1],[1,0,-1]])
B=matrix([[1],[2]])
A.solve_right(B)