Вычислить определитель матрицы

#include <iostream>
#include <math.h>
 
#define MAX_BUF 1000
 
double readValue()
{
  while (true)
  {
    double a;
    std::cin >> a;
    if (std::cin.fail())
    {
      std::cin.clear();
      std::cin.ignore(MAX_BUF, '\n');
      std::cout << "invalid input, try again" <<  std::endl;
    }
    else
    {
      std::cin.ignore(MAX_BUF, '\n');
      return a;
    }
  }
}
 
double determinant(double matrix[4][4], int n)
{   
  double det = 0;
  double submatrix[4][4];
  if (n == 2)
    return ((matrix[0][0] * matrix[1][1]) - (matrix[1][0] * matrix[0][1]));
  else
  {
    for (int x = 0; x < n; x++)
    {
      int subi = 0; 
      for (int i = 1; i < n; i++)
      {
        int subj = 0;
        for (int j = 0; j < n; j++)
        {
          if (j == x)
            continue;
          submatrix[subi][subj] = matrix[i][j];
          subj++;
        }
        subi++;
      }
      det = det + (pow(-1, x) * matrix[0][x] * determinant( submatrix, n - 1 ));
    }
 }
 return det;
}
 
int main()
{
  int i, j;
  double matrix[4][4];
  std::cout << "Enter the elements of the matrix:" << std::endl;
  for (i = 0; i < 4; i++)
  {
    std::cout << i+1 << " row:" << std::endl;
    for (j = 0; j < 4; j++)
    {
      std::cout << "element " << i+1 << ", " << j+1 << "= ";
      matrix[i][j] = readValue();
    }
  }
  std::cout << "The entered matrix is:" << std::endl;
  for (i = 0; i < 4; ++i)
  {
    for (j = 0; j < 4; ++j)
      std::cout << matrix[i][j] << " ";
    std::cout << std::endl;
  }
  std::cout << "Determinant of the matrix is " << determinant(matrix, 4) << std::endl;
  return 0;
}

#include <iostream>
#include <math.h>
 
#define MAX_BUF 1000
 
bool repeat()
{
  while (true)
  {
    std::cout << "do you want to countinue (y/n)? ";
    char sm;
    std::cin >> sm;
    std::cin.ignore(MAX_BUF,'\n');
    if (sm == 'y')   
      return true;
    else if (sm == 'n')
      return false;
  } 
}
 
double readValue()
{
  while (true)
  {
    double a;
    std::cin >> a;
    if (std::cin.fail())
    {
      std::cin.clear();
      std::cin.ignore(MAX_BUF, '\n');
      std::cout << "invalid input, try again" <<  std::endl;
    }
    else
    {
      std::cin.ignore(MAX_BUF, '\n');
      return a;
    }
  }
}
 
double determinant(double matrix[4][4], int n)
{   
  double det = 0;
  double submatrix[4][4];
  if (n == 2)
    return ((matrix[0][0] * matrix[1][1]) - (matrix[1][0] * matrix[0][1]));
  else
  {
    for (int x = 0; x < n; x++)
    {
      int subi = 0; 
      for (int i = 1; i < n; i++)
      {
        int subj = 0;
        for (int j = 0; j < n; j++)
        {
          if (j == x)
            continue;
          submatrix[subi][subj] = matrix[i][j];
          subj++;
        }
        subi++;
      }
      det = det + (pow(-1, x) * matrix[0][x] * determinant( submatrix, n - 1 ));
    }
 }
 return det;
}
 
int main()
{
  int i, j;
  double matrix[4][4];
  while (true)
  {
    std::cout << "Enter the elements of the matrix:" << std::endl;
    for (i = 0; i < 4; i++)
    {
      std::cout << i+1 << " row:" << std::endl;
      for (j = 0; j < 4; j++)
      {
        std::cout << "element " << i+1 << ", " << j+1 << "= ";
        matrix[i][j] = readValue();
      }
    }
    std::cout << "The entered matrix is:" << std::endl;
    for (i = 0; i < 4; ++i)
    {
      for (j = 0; j < 4; ++j)
        std::cout << matrix[i][j] << " ";
      std::cout << std::endl;
    }
    std::cout << "Determinant of the matrix is " << determinant(matrix, 4) << std::endl;
    if (!repeat())
      return 0;
  }
}