C语言-两n阶矩阵相乘

#include<stdio.h>
#define M 20
#define N 20

int multip(int a[M][N],int b[M][N],int c[M][N],int n)
{
	int i,j,k;
	for(i=0;i<n;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];
			}
		}
	}
	return c[M][N];
}

int main()
{
	int a[M][N],b[M][N],x[M][N];
	int i,j,n;
	printf("请输入矩阵阶数："); 
	scanf("%d",&n);
	printf("请输入矩阵A："); 
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
			scanf("%d",&a[i][j]);
	printf("请输入矩阵B："); 
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
			scanf("%d",&b[i][j]);
	multip(a,b,x,n);
	printf("所得矩阵C为：\n"); 
	for(i=0;i<n;i++)
	{
		for(j=0;j<n;j++)
		{
			printf("%d ",x[i][j]);
		}
		printf("\n");
	}
	return 0;	

}