A=matrix([[1,1,1],[1,2,1],[2,3,2]])
x=var('x');y=var('y');z=var('z')
D=matrix([[x],[y],[z]])
B=matrix([[2],[3],[5]])
x=A.solve_right(B);x
A*x