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