E=matrix(3,2,1],[2,1,1],[6,2,4) F=vector([3,0,6]) z = E.solve_right(F) print z