R3=matrix([[5,1,2],[1,4,-2],[2,3,8]])
print R3
p3=vector([19,-2,39])
x3=R3.solve_right(p3)
print x3