a=matrix([[1,1,1],[1,3,3],[1,3,5]])
b=vector([2,0,2])
x=a.solve_right(b)
print x