Код вычисления определителя матрицы до 10-го порядка

#include <iostream>
#include <cmath>
#include <ctime>
 
using namespace std;
 
void Gauss(double *a, double *b, double *x,int n) // а - матрица коэффициентов, строки подряд, b - свободные члены, x - решение.
{
 int i,j,k,t;
 double kof,s;
 double temp;
 for (i=n-1; i>0; --i) // Цикл по уравнениям - вычитаемым
 {
  for (t=i, j=i-1; j>=0; --j) //Ищем строку с максимальным в i-том столбце коэффициентом.
  {
   if (fabs(a[i*n+t])<fabs(a[i*n+j]))
   {
    t=j;
   }
  }
  if (a[i*n+t]==0.0)
  {
   return;
  }
  if (t!=i) // Если она не i-тая
  {
   for (k=n-1; k>=0; --k) // Меняем её с i-той
   {
    temp=a[k*n+t];
    a[k*n+t]=a[k*n+i];
    a[k*n+i]=temp;
   }
  }
  for (j=i-1; j>=0; --j) // Цикл пo уравнениям-уменьшаемым
  {
   kof=a[i+j*n]/a[i*n+i]; // kof=a[j][i]/a[i][i];
   for (a[i+j*n]=0.0, b[j]-=b[i]*kof, k=i-1; k>=0; --k) // Цикл по столбцам
   {
    a[k+j*n]-=a[k+i*n]*kof; // a[k][j]-=a[k][i]*kof
   }
  }
 }
 x=new double [n];
 if (x)
 {
  for (i=0; i<n; ++i)
  {
   for (s=0.0, j=i-1; j>=0; --j)
   {
    s+=a[j+i*n]*x[j]; //s+=a[j][i]*x[j]
   }
   x[i]=(b[i]-s)/a[i*n+i]; // x[i]=(b[i]-s)/a[i][i];
   //cout << x[i] << "\n";
  }
 }
 cout << "\n";
 for (i=0; i<n; ++i)
cout << "x " << x[i] << endl;
}
 
 
main()
{
srand((unsigned)time(NULL));
      int n = 6;
     double gauss2[6] ={1,5,3,7,6,3};
      double M[36]={0, 0, 1,0, 0, 0,0, 0, 0, -1, 0, 0,0, 8, 0, 0, 0, 4,0, 0, 0, 0, 6,0, 0, 6, 0, 0, 0, 0,1, 0, 0, 0, 0, 0 };
        
      double *y = new double[n];
 
    /*  for(int i=0; i<n; i++)
        {
            for(int j=0; j<n; j++)
                {
                    //cin >> M[i+j*n];
                     M[i+j*n] = rand()/100;
//cout << "gauss1["<< i << "][" << j << "]=";
//cin >> gauss1[i][j];
                }
        }
*/
for(int i=0; i<n; i++)
        {
            for(int j=0; j<n; j++)
                {
                     cout << M[i*n+j] << "  ";
                }
                cout << gauss2[i] << endl;
        }   
      Gauss(M,gauss2,y,n);
       cin.get();
}

#include <stdarg.h>
#include <iostream>
 
using namespace std;
 
int main()
{
    setlocale (LC_ALL,"RUS");
    int i,j,N,M;double *a;
 
    a = new double [N*M];
    for (i=0;i<N;i++)
    for (j=0;j<M;j++)
 
        cout<<"N="; cin>>N;
    
    if (N<=10)
    {
        cout<<"M="; cin>>M;
    }
    else
    {
        cout<<"число строк не должно превышать 10, введите другое число"; cin>>N;
    }
    
    cout<<"M="; cin>>M;
    if (M=N)
    {
        cout<<"введите матрицу А"<<endl;
    }
    else
    {
        cout<<"вы ввели не правильное число столбцов, оно должно равняться числу строк"; cin>>M;
    }
    {
        cin>>*(a+i*M+j);
    }
 
    delete []a;
    system("pause");
    return 0;
}

#include <iostream>
#include <cmath>
using namespace std;
 
int Gauss(double *aa, int n, int m, double *det, double eps)
#define a(i,j) aa[n*(j)+(i)]
{
    double e,p;
    int i,j,k,nm;
 
    *det=1.0;
    nm=n+m;
    for (i=0; i<n; i++) {
        e=0.0;
        for (j=i; j<n; j++)
        if ((p=fabs(a(j,i)))>e) {e=p; k=j;}
        if (e < eps) return i+1;
        if (k != i) {
            for (j=i; j<nm; j++) {
                p=a(i,j);
                a(i,j)=a(k,j);
                a(k,j)=p;
            }
            *det=-*det;
        }
        *det *= a(i,i);
        for (k=i+1; k<nm; k++) a(i,k) /= a(i,i);
        for (j = m!=0 ? 0 : i+1; j<n; j++) if (j != i) {
            p=a(j,i);
            for (k=i+1; k<nm; k++) a(j,k) -=a(i,k)*p;
        }
    }
    return 0;
}
 
double a[16]={
       1, -2,  3,  4,
       2,  1, -4,  3,
       3, -4, -1, -2,
       4,  3,  2, -1};
 
int main() {
    double det;
    int ret;
    ret=Gauss(a, 4, 0, &det, 0.000001);
    if(ret)
        det=0.0;
    cout << det << endl;
    return 0;
}