Решение системы линейных уравнений третьего порядка по правилу Крамера

#include <iostream>
#include <locale>
 
using std::cout;
using std::cin;
using std::endl;
 
const int n = 3, m = 4;
 
double determinant(double (*mat)[m]) {
    double det;
    double h1 = 1, h2 = 0, h3 = 1, h4 = 0;
    for (int i = 0; i < n ;i++) {
        h1 = 1;
        h3 = 1;
        int l = (2 * n) - 1 - i;
        for (int j = 0; j < n; j++) {
            h3 *= mat[j][l % n];
            l--;
            h1 *= mat[j][(j + i) % (n)];
        }
        h4 += h3;
        h2 += h1;
    }
    det = h2-h4;
    return det;
}
 
int main() {
    setlocale(LC_ALL, "Russian");
    double a[n][m];
    double b[n];
    double x[n];
    cout << "Введите матрицу 3х4 (первые три элемента в строке - a1, a2, a3; четвертый - b)= " << endl;
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < m; j++) {
            cin >> a[i][j];
        }
    }
    cout << endl << "Ваша система линейных уравнений:" << endl;
    for (int i = 0; i < n; i++) {
        b[i] = a[i][m-1];
        for (int j = 0; j < n; j++) {
            cout << "(" << a[i][j] << ")" << "*x" << j << " + " ;
        }
        cout << "\b\b = " << b[i] << endl;
    }
 
    cout << endl << "Определитель матрицы (без вектора b)" << endl;
    double det = determinant(a);
    cout << det;
    double h1 = 1, h2 = 0, h3 = 1, h4 = 0;
    if (det != 0 ) {
        cout << endl << "Корни уравнений:" << endl;
        for (int k = 0; k < n; k++) {
            h1 = 1;
            h2 = 0;
            h3 = 1;
            h4 = 0;
            for (int i = 0; i < n; i++) {
                h1 = 1;
                h3 = 1;
                int l = (2 * n) - 1 - i;
                for (int j = 0; j < n; j++) {
                    if ((j+i) % (n) == k)
                        h1 *= b[j];
                      else
                        h1 *= a[j][(i + j) % n];
                    if (l % n == k)
                        h3 *= b[j];
                      else
                        h3 *= a[j][l % n];
                l--;
                }
                h4 += h3;
                h2 += h1;
            }
            x[k+1] = h2 - h4;
            cout << "x(" << k+1 << ") = " << x[k+1]/det << endl;
            }
    } 
      else cout << "Система не имеет решений!";
    return 0;
}

determinant(double mat[][m], int n)

double det = determinant(a, n);