Возвести матрицу в степень

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
 
 
 
class matrix
{
 public:
 float** m;
 int sz;
 
 matrix(int n)
 {
  int i;
  m = new float*[n];
  for(i=0;i<n;i++)
   m[i]=new float[n];
  sz = n;
 }
 
 void random()
 {
  int i,j;
  srand(time(0));
  for(i=0;i<sz;i++)
   for(j=0;j<sz;j++)
    m[i][j]=1.0/(rand()%9+1);
 }
 
 void print()
 {
  int i,j;
  for(i=0;i<sz;i++, printf("\n"))
   for(j=0;j<sz;j++)
    printf("%.4f ", m[i][j]);
 }
 
 
};
 
void multiplicate_matrix(matrix* a1, matrix* a2, matrix* result)
{
 float y;
 int i,j,k;
 int size;
 
 size = a1->sz;
 y = 0.0;
 for(i=0;i<size;i++)
 {
  for(j=0;j<size;j++)
  {
    y = 0.0;
    for(k=0;k<size;k++)
    {
     printf("m[%d][%d]=%f * m[%d][%d]=%f\n", i,k,a1->m[i][k], k,j,a2->m[k][j], a1->m[i][k]*a2->m[k][j]);
     y+=a1->m[i][k]*a2->m[k][j];
    }
    printf("%f\n", y);
    result->m[i][j]=y;
  }
 }
}
 
void stepen_matrix(int k, matrix* a1, matrix* a2)
{
 static int y =0;
 y++;
 matrix t4(3);
 multiplicate_matrix(a2, a1, &t4);
 t4.print();
 if (y+1==k)
  return;
 stepen_matrix(k, a1, &t4);
}
 
int main()
{
 matrix t1(3);
 
 t1.random();
 t1.print();
 printf("\n");
 
 stepen_matrix(4, &t1, &t1);
 return 0;
}