Детерминант матрицы

#include <iostream>
#include <cstdlib>
#include <fstream>
 
using namespace std;
 
void print(double **,const int);
double determinant(double** ,int );
int sign(const int);
double** createMatrix(double **,const int ,const int);
 
int main()
{
   int size=9;
   double **a=new double*[size];
   ofstream array("array.txt",ios::out);
   
   for(int k=0;k<size;k++)
      a[k]=new double[size];
   
   for(int i=0;i<size;i++)
      for(int j=0;j<size;j++)
      {   
         a[i][j]=rand()%10;
         array<<a[i][j]<<endl;
      }  
   print(a,size);
         
   cout<<"\nDeterm is "<<determinant(a,size)<<endl;
   
   system("pause");
   return 0;   
}
 
double determinant(double **b,int rank)
{
   if(rank==2)
      return b[0][0]*b[1][1]-b[0][1]*b[1][0];
   else
   {
      double temp_sum=0; 
       
      for(int i=0;i<rank;i++)  
         temp_sum=temp_sum+sign(i)*b[i][0]*determinant(createMatrix(b,i,rank),rank-1);   
         
      return temp_sum;
   }
}
 
inline int sign(const int k)
{
   if(k%2==0)
      return 1;
   else
      return -1;       
}
 
double** createMatrix(double **matrix,const int p,const int r)
{
   double **newMatrix=new double*[r-1];
   
   for(int k=0;k<r;k++)
      newMatrix[k]=new double[r-1];
      
   for(int i=0;i<r-1;i++)
   {        
        if(i<p) 
           for(int j=0;j<r-1;j++)
              newMatrix[i][j]=matrix[i][j+1];    
        else
            for(int j=0;j<r-1;j++)
              newMatrix[i][j]=matrix[i+1][j+1];   
   }  
   return newMatrix;
}
 
void print(double **matrix,const int size)
{
   for(int i=0;i<size;i++)
      for(int j=0;j<size;j++)
         cout<<matrix[i][j]<<(j==size-1 ? "\n": "  ");
         
   cout<<endl;
}

    for(int k=0;k<r;k++)
        newMatrix[k]=new double[r-1];