A=matrix(-1,1,2],[3,-1,1],[-1,3,4) B=vector([2,6,4]) x=A.solve_right(B) print x