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

Strang-ex-1.2-17-U (last edited 2010-12-17 12:22:36 by prasad)