输入两个4*4矩阵，求矩阵乘积

#include <stdio.h>
int main()
{
int i,j,k;
int a1[4][4]={0},a2[4][4]={0};
int b[4][4]={0};

printf("第一个矩阵:\n");
for(i=0;i<4;i++)
{
    printf("第%d行>:",i+1);
    for(j=0;j<4;j++)
    {
    scanf("%d",&a1[i][j]);
    }
}

printf("\n");
printf("第二个矩阵:\n");
for(i=0;i<4;i++)
{
    printf("第%d行>:",i+1);
    for(j=0;j<4;j++)
    {
    scanf("%d",&a2[i][j]);
    }
}
printf("\n");

for(i=0;i<4;i++)
{
   for(j=0;j<4;j++)
   {
      for(k=0;k<4;k++)
      {
         b[i][j]+=a1[i][k]*a2[k][j];
      }
      printf("%d  ",b[i][j]);
      if((j+1)%4==0)
      printf("\n");
   }
}
printf("\n");
return 0;

}