Подстроение алгоритма определителя

    
int  a[10][10];
int  I=0, J=0;
 
        cout <<"\nВведите кол-во строк в матрице :\n";scanf("%d",&I);
    cout <<"\nВведите кол-во столбцов в матрице :\n";scanf("%d",&J);
    
    cout << "Введите элементы матрицы:\n";
 
        for(int i=0;i<I;i++)
            for(int j=0;j<J;j++)
                scanf("%d",&a[i][j]);
 
    cout << "Введенная матрица имеет вид:\n";
    
 
        for (int i=0;i<I;i++)
        {
            for(int j=0;j<J;j++)
    
            printf("\t%6.1d",a[i][j]);
            printf("\n");
        }
 
//Определитель
int j,x,y;
double k;
for (x=0;x<(I-1);x++)
{
 
for (y=(x+1);y<I;y++)
{
k=-a[y][x]/a[x][x];
for (j = 0;j< I;j++)
{
        a[y][j]+=k*a[x][j];
        if (abs(a[y][j]) < EPS) a[y][j]=0;
}
}
}
int det=1;
for(int i=0;i<I;i++)
det*=a[i][i];
 
 cout <<"\nОпределитель:" << det << "\n";

#include "stdafx.h"
#include <iostream>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
#define ESP 0.000000001
 
using namespace std;
 
int n;
double d;
 
void imas(double** X,int n);
void omas(double** X,int n);
void det(double** a,int n);
//void max(double** a,int x);
 
void main()
{
cout << "\nrazmernoct kvadratnoi matrix A(nxn)";
 
cin >> n;
double **A;
A=(double**)calloc(n,sizeof(double)); //Amn
for(int i=0;i<n;i++)
A[i]=(double*)calloc(n,sizeof(double));
 
imas(A,n);
omas(A,n);
det(A,n);
delete A;
getch();
}
//-----------------------------------
void imas(double** X,int n)
{
cout << "\nVvod matrix "<<n <<" x "<<n <<"\n";
int i,j;
for(i=0;i<n;i++)
for(j=0;j<n;j++)
cin >> X[i][j];
}
 
void omas(double** X,int n)
{
int i,j;
for(i=0;i<n;i++)
{
cout <<"\n";
for( j=0;j<n;j++)
cout<< X[i][j]<<"\t";
}
cout << "\n";
}
 
void det(double** a,int n)
{
int j,x,y;
double k;
for (x=0;x<(n-1);x++)
{
 
for (y=(x+1);y<n;y++)
{
k=-a[y][x]/a[x][x];
for (j = 0;j< n;j++)
{
a[y][j]+=k*a[x][j];
if (fabs(a[y][j]) < ESP) a[y][j]=0;
}
}
}
double det=1;
for(int i=0;i<n;i++)
det*=a[i][i];
cout <<"\ndet= "<< det;
}

int det=1;

double det=1.0;

scanf("%lf",&a[i][j]);

printf("\t%lf",a[i][j]);

        int  I=0, J=0, check=0, ri=0;
    double a[10][10];
    double rate=0,ratio=0,b[x][x],t;
    bool F;
 
//Проверка на нулевую строку с ненулевым вектором свободных членов данной строки
for(int i=0;i<I;i++)
  {
  for(int j=0;j<J-1;j++)
    if ((a[i][j]==0)&&(a[i][J-1]!=0))
      check++;
  if (check==J-1)
    {
    cout <<"Матрица не совместна.\n";
    
    F=false;
    system("pause"); //Пауза системы
    } 
  check=0;
  }
 
//Проверка на нулевой столбец - система имеет бесконечно много решений на данном столбце
for(int j=0;j<J-1;j++)
  {
  for(int i=0;i<I;i++)
    if (a[i][j]==0)
      check++;
  if (check==I)
    {
    printf("Матрица имеет бесконечно много решений на %d столбце. \n",j+1);
    
    F=false;
    getch();
    char c;
    c = getch();
    if (c==13) exit(0);
 
    system("pause"); //Пауза системы, скушай твикс
    }
  check=0;
  }
 
if (I<J-1)
  {
    printf("Матрица имеет бесконечно много решений, т.к. строк меньше, чем столбцов.\n");
 
    F=false;
    getch();
    char c;
    c = getch();
    if (c==13) exit(0);
  system("pause");//Пауза системы, скушай снигерс
  }
 
//Дублирование проверки
for (int i=0;i<I;i++)
  for(int j=0;j<J;j++)
    b[i][j]=a[i][j]*1.0;
 
for(int col=0;col<J-1;col++)
  {
  for(int i=col;i<I;i++)
    if (b[i][col]!=0)
      {
      rate=b[i][col];
      ri=i;
      break;
      }
  for(int j=0;j<J;j++)
    b[ri][j]/=rate;
  for(int j=0;j<J;j++) //Поднимаем строку с ненулевым элементом наверх на позиции кол
    {
    t=b[col][j];
    b[col][j]=b[ri][j];
    b[ri][j]=t;
    }
  for(int i=0;i<I;i++)
    if (i!=col)
      {
      ratio=b[i][col]/b[col][col];
      for(int j=0;j<J;j++)
        b[i][j]-=ratio*b[col][j];
      }
  }
cout << "Матрица преобразований имеет вид:\n";
gr;
for (int i=0;i<I;i++)
  {
  for(int j=0;j<J;j++)
    
    printf("\t%6.1f",b[i][j]);//Вывод матрицы преобразований
    printf("\n");
    
}
gr;
for(int i=0;i<I;i++)
  {
  for(int j=0;j<J-1;j++)
    if ((b[i][j]==0)&&(b[i][J-1]!=0))
      check++;
  if (check==J-1)
    {
    printf("Матрица не совместна.\n");
    
    F=false;
    getch();
    char c;
    c = getch();
    if (c==13) exit(0);
    system("pause"); //Пауза системы, скушай баунти
    }
  check=0;
  }
 
cout << "\nПолученные решения:\n";
for(int i=0;i<J-1;i++)
  printf("\nx%d = %.3lf\n",i+1,b[i][J-1]);

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#define EPS 0.000000001
 
int main()
{
    double **A, **Am, *X, *b, mu, t;  // AX = b
    int i, j, n, k, x, y;
 
    printf ("Введите размерность матрицы:\n");
    scanf("%d", &n);
 
    A = (double**)malloc(n * sizeof(double*));
    for (i = 0; i < n; i ++)
        A[i] = (double*)malloc(n * sizeof(double));
    Am = (double**)malloc(n * sizeof(double*));
    for (i = 0; i < n; i ++)
        Am[i] = (double*)malloc(n * sizeof(double));
    X = (double*)malloc(n * sizeof(double));
    b = (double*)malloc(n * sizeof(double));
 
    printf ("Заполните матрицы A i b:\n");
 
    for (i = 0; i < n; i ++) {
        for (j = 0; j < n; j ++) {
            printf ("Введите A[%d][%d]=", i, j);
            scanf ("%lf", &A[i][j]);
            Am[i][j] = A[i][j];
        }
        printf ("Введите b[%d]=", i);
        scanf("%lf", &b[i]);
    }
 
    printf ("Исходная матрица:\n");
    for (i = 0; i < n; i ++) {
        for (j = 0; j < n; j ++)
            printf ("%9.2lf", A[i][j]);
        printf ("\n");
    }
 
    double kk;
 
    for (x = 0; x < (n - 1); x ++) {
        for (y = (x + 1); y < n; y ++) {
            kk =- Am[y][x] / Am[x][x];
            for (j = 0; j < n; j ++) {
                Am[y][j] += kk * Am[x][j];
                if ( fabs(A[y][j] ) < EPS) Am[y][j] = 0;
            }
        }
    }
 
    double det = 1.0;
    for(int i = 0; i < n; i ++)
        det *= Am[i][i];
 
    printf("\nОпределитель = %lf\n", det);
 
    printf ("Метод Гаусса:\n");
    for (k = 0; k < n; k ++)
        for (i = k + 1 ; i < n; i ++) {
            mu = A[i][k] / A[k][k];
            for (j = 0; j < n; j ++)
                A[i][j] -= A[k][j] * mu;
            b[i] -= b[k] * mu;
        }
 
    printf ("Матрица А, приведённая к треугольному виду; матрица b :\n");
    for (i = 0; i < n; i ++) {
        for (j = 0; j < n; j ++)
            printf ("%9.2lf", A[i][j]);
        printf ("%25.2lf", b[i]);
        printf ("\n");
    }
 
    for (i = n - 1; i >= 0; i --) {
        t = b[i];
        for(j = i + 1; j < n; j++)
            t -= A[i][j] * X[j];
        X[i] = t / A[i][i];
    }
 
    printf ("Ответ для метода Гаусса\n");
    for (k = 0; k < n; k ++) printf ("%lf\n", X[k]);
 
    for (i = 0; i < n; i ++) free(A[i]);
    free (A);
 
    for (i = 0; i < n; i ++) free(Am[i]);
    free (Am);
 
    free (X);
    free (b);
    return 0;
}