"puissance niéme d'un matrise", j'ai besoin de votre aide...

Bonsoir,
Je veux que vous m'aidiez à trouver une solution "puissance d'un matrice" en utilise la fonction produit.
A,A^2,A^3...........A^n
le fichier "s1":
4
1 2 3 4
1 1 1 2
1 4 3 2

2 3 6 3
1 2 5 0
2 3 1 2
1 0 4 4

J'ai essayé, mais je ne pouvais pas trouver une solution.

#include<stdio.h>
#include<stdlib.h>

int lecture(FILE *f,int m1[][20] )
{
int i ,j , n;

fscanf(f , "%d" , &n) ;

for( i = 0 ; i < n ; i ++ )
for( j = 0 ; j <n ; j++ )
fscanf(f , "%d" , &m1[i][j]) ;

return n ;
}
//***********************************************************
int affiche( FILE *f , int m[20][20] , int n )
{
int i ,j ;

for( i = 0 ; i < n ; i ++ )
{
fprintf(f,"\n") ;
for( j = 0 ; j <n ; j++ )
fprintf(f, "\t%d" , m[i][j]) ;
}
return 0;
}
//*****************************************************************//
void produit (int m1[][20],int m2[][20],int m[][20],int n)
{
int i,j,k;
for( i = 0 ; i < n ; i ++ )
{ for(j=0;j<n;j++ )
{
m[i][j]=0;
for(k=0;k<n;k++)
{
m[i][j]+=m1[i][k]*m2[k][j];
}
}
}

}

//***************************************************************
void puiss_mat(m[][20],int n ,int p)
{
int i,j,l;
for(l=2;l<k-1;l++)
{
if(k<3) return produit(m1,m1,m,n);
?????????
}
}
//****************************************************************
int main()
{
int m1[20][20] ;
int m2[20][20] ;
int m[20][20];
int n;
char s1[20] ;
printf(" Fichier de données : " );
scanf("%s" , s1 ) ;
FILE *graphe = fopen(s1, "r+" ) , *resultat;
n= lecture(graphe,m1 );
lecture(graphe,m2 );
produit(m1,m2,m,n);
resultat = fopen("res.txt" ,"w+") ;
affiche(resultat, m, n ) ;
system("pause") ;
return 0 ;
}