Умножить матрицу на транспонированную - C (СИ)

#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <time.h>
# define n 3
int main()
{
 srand(time(NULL));
 int i,j,k, a[n][n],b[n][n],c[n][n];
 printf("Random matrix A\n");
 for (i=0; i<n; i++)
  {
   for(j=0; j<n; j++)
    {
      a[i][j]=rand()%10;
      printf("%3d",a[i][j]);
    }
   printf("\n");
  }
 printf("\nTransposition matrix A*\n");
 for (j=0; j<n; j++)
  for(i=0; i<n; i++)
   b[i][j]=a[j][i];
 for (i=0; i<n; i++)
  {
   for(j=0; j<n; j++)
    {
     printf("%3d",b[i][j]);
    } 
   printf("\n");
  }
 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]);
    }
  }
 printf("\nProduct A & A*\n");
 for (i=0; i<n; i++)
  {
   for(j=0; j<n; j++)
    {
      printf("%5d",c[i][j]);
    }
   printf("\n");
  }
 getch();
 return 0;
}