A program on matrix multiplication

#include<stdio.h>
#include<stdlib.h>
#define N 4

void multiply(int *p, int *q);

int a[N][N]={{0,1,1,1}, {0,0,0,1}, {0,0,0,1}, {0,0,0,0}};
int b[N][N];
int c[N][N];
int A2[N][N]={0};

int main()
{
   int i,j,k;

   printf("Please input matrix a：\n");  
   for(i=0;i<N;i++)
   {
       for(j=0;j<N;j++)
  {
  b[i][j]=a[i][j];
  printf("%d  ",a[i][j]);   
  }   
  printf("\n");
   }


/*
   printf("Please input a matrix：\n");   
   for(i=0;i<N;i++)
        for(j=0;j<N;j++)
  scanf("%d",&a[i][j]);
*/

   for(i=1;i<N-1;i++)
   {
  printf("The %d time multiplication:\n",i);
  multiply(a[0],b[0]);
  for(j=0;j<N;j++)
  for(k=0;k<N;k++)
  A2[j][k]=A2[j][k]+c[j][k];
     
   }

   printf("Please input matrix A2：\n");  
   for(i=0;i<N;i++)
   {
       for(j=0;j<N;j++)
  printf("%d  ",A2[i][j]);
  printf("\n");
   }


  system("pause");
  return 0;  
}


void multiply(int *p,int *q)
{
   int i,j,k;

   for(i=0;i<N;i++)
        for(j=0;j<N;j++)
  c[i][j]=0;

   for(i=0;i<N;i++)
        for(j=0;j<N;j++)
  {
    for(k=0;k<N;k++)
c[i][j]=c[i][j]+a[i][k]*b[k][j];  
  }

  for(i=0;i<N;i++)
   {
       for(j=0;j<N;j++)
  {
  b[i][j]=c[i][j];
  printf("%d  ",c[i][j]);   
  }   
  printf("\n");
   }
}