Vectorized 2-dimensional multiplication of multidimensional matrices

If you dont want to use slow for loops to apply a matrix multiplication between two n-dimensional matices, this submission could be helpful. Unfortunately, singelton expansion is not supported. Thus it is required that boths inputs are the same size.
Input 1: n-by-m-by-... Matrix
Input 2: m-by-p-by-... Matrix.

Output: n-by-p-by-... Matrix.