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