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

#include <windows.h>
#include <stdio.h>
#include <conio.h>
 
int ** a = (int **)malloc(sizeof(int));
int * b = (int *)malloc(sizeof(int));
 
void out_massiv(int n, int *mass);
int compress(int m, int **a, int *b);
int * compress(int m, int i, int * vec);
 
void main()
{
    srand(0);
    int n,m,i,j,k;
    do
    {
        printf("Enter num of rows in massiv: ");
        scanf("%d",&m);
        a = (int **)realloc((void *)a,m*sizeof(int));
        b = (int *)realloc((void *)b,m*sizeof(int));
        printf("\tMatrix A\r\n");
        for(i = 0,k = 0; i < m; i++)
        {
            b[i] = /*rand()*/i;
            a[i] = (int *)malloc(m*sizeof(int));
            for(j = 0; j < m; j++)
            {
/*              a[i][j] = rand();
                while(1000 < a[i][j])
                    a[i][j] = rand();*/
                a[i][j] = k++;
            }
            if(i % 2 == 0)
                a[i][i] = 0;
            out_massiv(m, a[i]);
        }
        printf("\tVector B\r\n");
        out_massiv(m, b);
        printf("Compressing\r\n");
        printf("Moved %d rows\r\n",n = compress(m, a, b));
        m -= n;
        printf("\tMatrix A\r\n");
        for(i = 0; i < m; i++)
            out_massiv(m, a[i]);
        printf("\tVector B\r\n");
        out_massiv(m, b);
        printf("[Y/N] Y - enter new input data\r\n");scanf("%d",m);
    }
    while(toupper(getch()) == 'Y');
}
 
void out_massiv(int n, int *mass)
{
    for(int i = 0;i < n; i++)
        printf("%02d ",mass[i]);
    printf("\r\n");
}
 
int compress(int m, int **a, int *b)
{
    int n = 0;
    for(int i = m - 1,j; 0 <= i; i--)
    {
        if(a[i][i] == 0)
        {
            printf("TEMP A\r\n");
            for(j = 0; j < m; j++)
            {
                a[j] = compress(m, i, a[j]);
                out_massiv(m - 1, a[j]);
            }
            for(j = i; j < m; j++)
                memmove((void *)&a[j],(void *)&a[j + 1],m);
            b = compress(m, i, b);
            printf("TEMP B\r\n");
            out_massiv(m, b);
            n++;
            m = m - n;
        }
    }
    return n;
}
 
int * compress(int m, int i, int * vec)
{
    memmove((void *)&vec[i],(void *)&vec[i + 1],m - (i + 1));
    return vec;
}

int main()
{
    int range = 5;
 
    // Allocate 2 matrix
    double **matrix = new double*[range]();
    for (int i = 0; i < range; i++) 
    {
        matrix[i] = new double[range]();
    }
 
    // reading from file
    ifstream ifs("aa.txt");
    double q;
    int i = 0;
    int j = 0;
    while ( ifs >> q )
    {
        matrix[i][j] = q;
        j++;
        if ( j % range == 0 ) { i++; j = 0;}
    }
    ifs.close();
    // delete nessesary colums and rows
 
    // счетчик количества нулевых строк/столбцов
    int k = 0;
    ofstream inf("inf.txt");
 
    // Тут я просто хочу удалить из массива нулевую строчку\столбец.
    for ( int i = 0; i < range; i++ )
    {
        if ( matrix[i][i] == 0 )
        {
            inf << i << endl;
            k++;
            delete matrix[i];
        }
    }
    inf.close();
 
    // output matrix 
    ofstream ofs("bb.txt");
    for (int i = 0; i < range-k; i ++)
    {
        for (int j = 0; j < range-k; j++)
        {
            ofs << matrix[i][j] << " ";
        } ofs << endl;
    }
    for ( int i =0; i < range-k; i++)
    delete[] matrix[i];
    delete[] matrix;

        for (int i = 0; i < range-k; i ++)
        {
                if ( matrix[i][i]!=NULL)
                {
                        for (int j = 0; j < range-k; j++)
                        {
                                ofs << matrix[i][j] << " ";
                        } ofs << endl;
                }
        }

#include <iostream>
#include <vector>
#include <limits>
#include <algorithm>
using namespace std;
const int INF = 123456789;
const double EPS = numeric_limits<double>::epsilon();
int gauss (vector < vector<double> > a, vector<double> & ans) {
    int n = (int) a.size();
    int m = (int) a[0].size() - 1;
    vector<int> where (m, -1);
    for (int col=0, row=0; col<m && row<n; ++col) {
        int sel = row;
        for (int i=row; i<n; ++i)
            if (abs (a[i][col]) > abs (a[sel][col]))
                sel = i;
        if (abs (a[sel][col]) < EPS)
            continue;
        for (int i=col; i<= m; ++i)
            swap (a[sel][i], a[row][i]);
        where[col] = row;
 
        for (int i=0; i<n; ++i)
            if (i != row) {
                double c = a[i][col] / a[row][col];
                for (int j=col; j<= m; ++j)
                    a[i][j] -= a[row][j] * c;
            }
        ++row;
    }
 
    ans.assign (m, 0);
    for (int i=0; i<m; ++i)
        if (where[i] != -1)
            ans[i] = a[where[i]][m] / a[where[i]][i];
    for (int i=0; i<n; ++i) {
        double sum = 0;
        for (int j=0; j<m; ++j)
            sum += ans[j] * a[i][j];
        if (abs (sum - a[i][m]) > EPS)
            return 0;
    }
 
    for (int i=0; i<m; ++i)
        if (where[i] == -1)
            return INF;
    return 1;
}
int main(){
    freopen("input.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
    int n;
    std::cin >> n;
    std::vector < std::vector<double> > vec(n);
    for (int i = 0; i < n; ++i)
    {
        vec[i].resize(n + 1);
        for (int j = 0; j <= n; ++j)
            cin >> vec[i][j];
    }
    vector<double> res(n);
    gauss(vec, res);
    for (vector<double>::iterator it = res.begin(); it != res.end(); ++it)
        std::cout << static_cast<int> (*it + (*it < 0 ? -0.5 : 0.5)) << ' ';
    return 0;
}

void PryamoiHod(int n, double **a, double *b)
{
        double v;
        for(int k = 0,i,j,im; k < n - 1; k++)
        {
                printf("UPDATING %04d ROW\r\n",k + 1);
                im = k;
                for(i = k + 1; i < n; i++)
                {
                        if(fabs(a[im][k]) < fabs(a[i][k]))
                        {
                                im = i;
                        }
                }
                if(im != k)
                {
                        for(j = 0; j < n; j++)
                        {
                            
                                v                = a[im][j];
                                a[im][j] = a[k][j];
                                a[k][j]  = v;
                        }
                        v     = b[im];
                        b[im] = b[k];
                        b[k]  = v;
                }
                for(i = k + 1; i < n; i++)
                {
                        v               = 1.0*a[i][k]/a[k][k];
                        a[i][k] = 0;
                        b[i]    = b[i] - v*b[k];
                        for(j = k + 1; j < n; j++)
                                a[i][j] -= v*a[k][j];
                }
        }
}
 
void ObratniHod(int n, double **a, double *b, double *x)
{
        double s = 0;
        x[n - 1] = 1.0*b[n - 1]/a[n - 1][n - 1];
        for(int i = n - 2, j; 0 <= i; i--)
        {
                printf("UPDATING %04d ROW\r\n",n - i - 1);
                s = 0;
                for(j = i + 1; j < n; j++)
                {
                        s = s+a[i][j]*x[j];
                }
                x[i] = 1.0*(b[i] - s)/a[i][i];
        }
}