restart:
with(linalg): with(plots):
A:=matrix(2,3,[1,2,-1,0,5,2]);
B:=matrix(3,3,[1,2,0,0,3,-7,5,2,1]);
evalm(A&*B);
X:=matrix(3,1,[x,y,z]);
evalm(A&*X);
evalm(A*B);
evalm(B&*A);
print(A);
print(X);
evalm(A&*X);
B:=matrix(3,2,[6,3,-1,2,8,1]);
evalm(A&*B);
A1:=matrix(2,2,[1,2,-3,1]);
A2:=matrix(2,2,[3,-2,5,4]);
evalm(A1&*A2);
evalm(A2&*A1);
A:=matrix(2,2,[1,-2,2,-4]);
evalm(inverse(A));