A=matrix([[1,0],[2,1]])
B=matrix([[1,3],[0,1]])
Q=A*B
print A*B
Q.transpose()== B.transpose()*A.transpose()