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