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

Kreyszig-6.3-2 (last edited 2010-12-15 12:38:36 by vasanthi emani)