Как выводить степень матрицы

#include <stdio.h>
#include <conio.h>
#include <math.h>
#define Max 100
#define N 5
int main ()
{int min;
int i,j,C[N][N],A[N][N],B[N][N],B_new[N][N],k,fl,temp,stepen=1;
printf("Enter the elements of the matrix C of dimension N, with the infinity take the number 100\n");
for (int i = 0; i<N; i++){
    for (int j = 0; j<N; j++){
         scanf("%i", &C[i][j]);
}
}
for (i=0;i<N;++i)
    for (j=0;j<N;++j)
        A[i][j]=B[i][j]=C[i][j];
for ( i=0; i<N; i++)
   C[i][i]= A[i][i] =B[i][i]=0;//замена диагональных элементов матрицы на 0
//B_new[i][j]
do{fl=0;
++stepen;
for (i=0;i<N;++i)
    for (j=0;j<N;++j){
     min = A[i][0]+B[0][j];
     for (k=1;k<N;++k){
          temp=A[i][k]+B[k][j];
          if( min >temp) min = temp;
     }
     B_new[i][j]= min;
        }
    for(i=0;i<N;i++){
        for(j=0;j<N;j++){
              if(B[i][j]!=B_new[i][j])
              {fl=1;
                for (i=0;i<N;++i)
                    for (j=0;j<N;++j)
                        B[i][j]=B_new[i][j];
              break;
          }}
          if(fl) break;}
          //B^stepen
          }while(fl);
 
 
//exit
printf("Matrix B\n");
for ( i = 0; i < N; i ++ ) {
for ( j = 0; j < N; j ++ )
    printf ( "\t%d" ,B[i][j] );
printf("\n");
 
}
return 0;
}