Solution

C=MATMUL(A,B) performs `proper' mathematical matrix multiplication. Element C(i,j) equals DOT_PRODUCT(A(i,:),B(:,j))

C=A*B, however, performs element-by-element multiplication. Element C(i,j) equals A(i,j)*B(i,j).