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)