A = matrix([[3,1,-1],[2,2,-3],[-1,1,-2]])
print A

B = vector([3,1,-2])
x= A.solve_right(B)
print x



A * x

x, y, z = var('x, y, z')
solve([3*x + y - z     == 3
       2*x + 2*y - 3*z == 1
        -x +   y - 2*z == -2], x, y, z)

SageDays/Kreyszig-7.11-21 (last edited 2010-08-09 14:30:38 by shambulingayya.N.D)