A = matrix([[1, -1, 1], [-1, 1, -1], [0, 10, 25], [20, 10, 0]])
b = vector([0, 0, 90, 80])
x = A.solve_right(b)
print x