Умножить одну матрицу на вторую - оптимизация циклов

#include<stdio.h>
#include<conio.h>
main()
{
clrscr();
FILE *A, *B, *C;
int i,j;
int a[5][5], b[5][5], c[5][5], aT[5][5], aTb[5][5];
 
A=fopen("A.TXT", "r");
for(i=1; i<=4; i++)
  {
   for(j=1; j<=4; j++)
   fscanf(A, "%d ", &a[i][j]);
  }
fclose(A);
 
B=fopen("B.TXT", "r");
for(i=1; i<=4; i++)
  {
   for(j=1; j<=4; j++)
   fscanf(B, "%d ", &b[i][j]);
  }
fclose(B);
 
printf("Matrycja A ");
for(i=1; i<=4; i++)
  {
   printf("\n");
   for(j=1; j<=4; j++)
     printf("%d ", a[i][j]);
  }
 
printf("\n\nMatrycja B ");
for(i=1; i<=4; i++)
  {
   printf("\n");
   for(j=1; j<=4; j++)
     printf("%d ", b[i][j]);
  }
 
for(i=1; i<=4; i++)
  {
   for(j=1; j<=4; j++)
   aT[j][i]=a[i][j];
  }
 
printf("\n\nMatrycja AT ");
for(i=1; i<=4; i++)
  {
   printf("\n");
   for(j=1; j<=4; j++)
     printf("%d ", aT[i][j]);
  }
 
aTb[1][1]=aT[1][1]*b[1][1]+aT[1][2]*b[2][1]+aT[1][3]*b[3][1]+aT[1][4]*b[4][1];
aTb[1][2]=aT[1][1]*b[1][2]+aT[1][2]*b[2][2]+aT[1][3]*b[3][2]+aT[1][4]*b[4][2];
aTb[1][3]=aT[1][1]*b[1][3]+aT[1][2]*b[2][3]+aT[1][3]*b[3][3]+aT[1][4]*b[4][3];
aTb[1][4]=aT[1][1]*b[1][4]+aT[1][2]*b[2][4]+aT[1][3]*b[3][4]+aT[1][4]*b[4][4];
 
aTb[2][1]=aT[2][1]*b[1][1]+aT[2][2]*b[2][1]+aT[2][3]*b[3][1]+aT[2][4]*b[4][1];
aTb[2][2]=aT[2][1]*b[1][2]+aT[2][2]*b[2][2]+aT[2][3]*b[3][2]+aT[2][4]*b[4][2];
aTb[2][3]=aT[2][1]*b[1][3]+aT[2][2]*b[2][3]+aT[2][3]*b[3][3]+aT[2][4]*b[4][3];
aTb[2][4]=aT[2][1]*b[1][3]+aT[2][2]*b[2][3]+aT[2][3]*b[3][3]+aT[2][4]*b[4][3];
 
aTb[3][1]=aT[3][1]*b[1][1]+aT[3][2]*b[2][1]+aT[3][3]*b[3][1]+aT[3][4]*b[4][1];
aTb[3][2]=aT[3][1]*b[1][2]+aT[3][2]*b[2][2]+aT[3][3]*b[3][2]+aT[3][4]*b[4][2];
aTb[3][3]=aT[3][1]*b[1][3]+aT[3][2]*b[2][3]+aT[3][3]*b[3][3]+aT[3][4]*b[4][3];
aTb[3][4]=aT[3][1]*b[1][3]+aT[3][2]*b[2][3]+aT[3][3]*b[3][3]+aT[3][4]*b[4][3];
 
aTb[4][1]=aT[4][1]*b[1][1]+aT[4][2]*b[2][1]+aT[4][3]*b[3][1]+aT[4][4]*b[4][1];
aTb[4][2]=aT[4][1]*b[1][2]+aT[4][2]*b[2][2]+aT[4][3]*b[3][2]+aT[4][4]*b[4][2];
aTb[4][3]=aT[4][1]*b[1][3]+aT[4][2]*b[2][3]+aT[4][3]*b[3][3]+aT[4][4]*b[4][3];
aTb[4][4]=aT[4][1]*b[1][3]+aT[4][2]*b[2][3]+aT[4][3]*b[3][3]+aT[4][4]*b[4][3];
 
printf("\n\nMatrycja AT*B");
for(i=1; i<=4; i++)
  {
   printf("\n");
   for(j=1; j<=4; j++)
     printf("%d ", aTb[i][j]);
  }
 
for(i=1; i<=4; i++)
  {
   for(j=1; j<=4; j++)
   c[i][j]=3*(aTb[i][j]-b[i][j])-b[i][j];
  }
 
printf("\n\nObchysljujemo matrycju C=3(ATB-B)-B");
printf("\nMatrycja C dorivnjuje:\n");
for(i=1; i<=4; i++)
  {
   printf("\n");
   for(j=1; j<=4; j++)
     printf("%d ", c[i][j]);
  }
printf("\n\nVykonala studentka grupy KNz11 Dyka Oksana");
 
C=fopen("c.txt","w");
for(i=1; i<=4; i++)
  {
   for(j=1; j<=4; j++)
     fprintf(C, "%d ", c[i][j]);
   fprintf(C, "\n");
  }
return 0;
}