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