2维数组矩阵相乘

#include <stdio.h>
#define M 3
#define N 3
void main()
{
int A[M]
,B[M]
,C[M]
;
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维数组。。