A=matrix([[1,0],[2,4],[2,1]])
B=matrix([[3,3,0],[1,2,1]])
print A*B
p=A[::,:1]*B[:1]
q=A[::,1:]* B[1:]
print p+q
print A*B==(p+q)

#the given matrix AB using coloumn times rows