Решение системы уравнений

#include <iostream>
#include <fstream>
using namespace std;
 
int main()
{
    int n, m,i,j,i2,j2;
    double** mas;
    double k1, k2;
    string test = "test.txt";
    ifstream file(test);
    file >> n >> m;
    m++;
    mas = new double* [n];
    int *s = new int [m];
    for (i = 0; i < n; i++) {
        mas[i] = new double[m];
        for (j = 0; j < m; j++) {
            file >> mas[i][j];
            s[j] = 0;
        }
    }
    for (i = 0; i < n; i++)
        for (j = 0; j < m; j++) {
            if (mas[i][j] != 0 && s[j] == 0) {
                for (i2 = 0; i2 < n; i2++) mas[i2][j] = 0;
                for (j2 = 0; j2 < m; j2++) mas[i][j2] = mas[i][j2]/mas[i][j];
 
                for (i2 = 0; i2 < n; i2++)
                    for (j2 = 0; j2 < m; j2++) 
                        if (i2 != i && j2 != j) {
                            k1 = mas[i2][j];
                            k2 = mas[i][j2];
                            if ((i <= i2 && j <= j) || (i <= i && j <= j2) || (i2 <= i && j2 <= j2) || (i2 <= i2 && j2 <= j)) {
                                mas[i2][j2] = k1 * k2 - mas[i2][j2] * mas[i][j];
                            }
                            else mas[i2][j2] = mas[i2][j2] * mas[i][j] - k1 * k2;
                        }
                mas[i][j] = 1;
                s[j] = 1;
            }
        }
    for (i = 0; i < n; i++){
        for (j = 0; j < m; j++)
            cout << mas[i][j]<<" ";
        cout << "\n";
}