2维数组矩阵相乘

#include <stdio.h>
#define M 3
#define N 3
void main()
{
	int A[M][N],B[M][N],C[M][N];
    int i,j,k;
		for(i=0;i<M;i++)
			for(j=0;j<N;j++)
				scanf("%d",&A[i][j]);
			for(i=0;i<M;i++)
			for(j=0;j<N;j++)
				scanf("%d",&B[i][j]);
			for(i=0;i<M;i++)
			for(j=0;j<N;j++)
			{	
				C[i][j]=0;
				for(k=0;k<N;k++)
					C[i][j]+=A[i][k]*B[k][j];
			}
			for(i=0; i<M; i++)
    {
        for(j=0; j<N; j++)
            printf("%d\t", C[i][j]);
        printf("\n");
    }
}

写得极丑无比，且只能2个矩阵相乘，行与列必须在程序中改，对于很久没用c语言的我，很久没用2维数组。。