Нахождения детерминанта

double thre(double a,double  x,double  c,double a1,double x1,double c1,double a2,double x2,double c2)
{
    SetConsoleCP(1251);
    SetConsoleOutputCP(1251);
    double delt = (a * x1 * c2) + (a1 * x2 * c) + (x * c1 * a2) - (c * x1 * a2) - (x * a1 * c2) - (c1 * x2 * a);
    return delt;
}

    #include <iostream>
    #include <cstdlib>
    #include <windows.h>
    #include "Header3.h"
    using namespace std;
    void forxfor()
    {
        SetConsoleCP(1251);
        SetConsoleOutputCP(1251);
        double a, x, c, b, f;
        double a1,x1,c1,b1,f1;
        double a2,x2,c2,b2,f2;
        double a3,x3,c3,b3,f3;
        cout << "Введіть коефіцієнт x y z b і це, що після знака дорівнює : " <<endl;
        cin >> a >> x >> c >> b >> f; 
        cout << "Введіть коефіцієнт x1 y1 z1 b1 і це, що після знака дорівнює : " <<endl;
        cin >> a1 >> x1 >> c1 >> b1 >> f1;
        cout << "Введіть коефіцієнт x2 y2 z2 b2 і це, що після знака дорівнює : " <<endl;
        cin >> a2 >> x2 >> c2 >> b2 >> f2;
        cout << "Введіть коефіцієнт x3 y3 z3 b3 і це, що після знака дорівнює : " <<endl;
        cin >> a3 >> x3 >> c3 >> b3 >> f3;
        double delta = a * thre(x1,c1,b1,x2,c2,b2,x3,c3,b3) - a1 * thre(x, c, b,x2,c2,b2,x3,c3,b3) + a2 * thre(x, c, b,x1,c1,b1,x3,c3,b3) - a3 * thre(x, c, b,x1,c1,b1,x2,c2,b2);
        double delta1 = f * thre(x1,c1,b1,x2,c2,b2,x3,c3,b3) - f1 * thre(x, c, b,x2,c2,b2,x3,c3,b3) + f2 * thre(x, c, b,x1,c1,b1,x3,c3,b3) - f3 * thre(x, c, b,x1,c1,b1,x2,c2,b2);
        double rozv_x = delta1 / delta;
        double delta2 = a * thre(f1,c1,b1,f2,c2,b2,f3,c3,b3) - a1 * thre(f, c, b,f2,c2,b2,f3,c3,b3) + a2 * thre(f, c, b,f1,c1,b1,f3,c3,b3) - a3 * thre(f, c, b,f1,c1,b1,f2,c2,b2);
        double rozv_y = delta2 / delta;
        double delta3 = a * thre(x1,f1,b1,x2,f2,b2,x3,f3,b3) - a1 * thre(x, f, b,x2,f2,b2,x3,f3,b3) + a2 * thre(x, f, b,x1,f1,b1,x3,f3,b3) - a3 * thre(x, f, b,x1,f1,b1,x2,f2,b2);
        double rozv_z = delta3 / delta;
        double delta4 = a * thre(x1,c1,f1,x2,c2,f2,x3,c3,f3) - a1 * thre(x, c, f,x2,c2,f2,x3,c3,f3) + a2 * thre(x, c, f,x1,c1,f1,x3,c3,f3) - a3 * thre(x, c, f,x1,c1,f1,x2,c2,f2);
        double rozv_b = delta4 / delta;
        if ( delta == 0 )
        {
            cout << "Визначник = 0, помилка при коефіцієнтах"<<endl;
        }
        else
        {
        cout<< "x = " << delta1 << " / " << delta <<" = "<< rozv_x <<endl;
        cout << "y = " << delta2 << " / " << delta  <<" = "<<rozv_y <<endl;
        cout << "z = " << delta3 << " / " << delta  <<" = "<<rozv_z <<endl;
        cout << "b = " << delta4 << " / " << delta  <<" = "<<rozv_b <<endl;
        }
    }

#include <iostream>
#include <cstdlib>
#include <cmath>
#include <windows.h>
#include "det2.h"
#include "header1.h"
#include "header2.h"
#pragma once;
#include <string>
using namespace std;
int main()
{
    SetConsoleCP(1251);
    SetConsoleOutputCP(1251);
    int otvet;
    cout << "Введіть порядок матриці 2 | 3 | 4 : ";
    cin >> otvet;
    if (otvet == 2)
    {
        twoxtwo();
    }
    if (otvet == 3)
    {
        threxthre();
    }
    if (otvet == 4)
    {
        forxfor();
    }
    return 0;
}