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

 int range;
    cin >> range;
    // Colum of free members
    double *free = new double[range];
    ifstream ifsf("B.txt");
    double q;
    int i = 0;
    int j = 0;
    while ( ifsf >> free[i] )
    {
        i++;
    }
    ifsf.close();
    ofstream ofs("B-0.txt");
    for ( int i = 0; i < range; i++ )
    {
        ofs << free[i] << endl;
    }
    ofs.close();

ifsf>>mass[i];

#include <iostream.h>
#include <fstream.h>
 
int main()
{
    int range;
    cin >> range;
    // Colum of free members
    double *mass = new double[range];
    ifstream ifsf("B.txt");
    int i = 0;
    int j = 0;
    while (i < range)
    {
        ifsf>>mass[i];
        cout<<mass[i]<<"\r\n";
        i++;
    }
    ifsf.close();
 
    ofstream ofs("B-0.txt");
    for (j = 0; j < range; j++ )
        ofs<<mass[i]<<"\r\n";
    ofs.close();
    return 0;
}

    int range;
    cout << "Please enter tne range os ypur system" << endl;
    cin >> range;
    // Colum of free members
    double *free = new double[range];
    ifstream ifsf("B.txt");
    double q;
    int i = 0;
    int j = 0;
    while ( i < range )
    {
        ifsf >> free[i];
        i++;
    }
    ifsf.close();
    ofstream ofs("B-0.txt");
    for ( int w = 0; w < range; w++ )
    {
        ofs << free[w] << endl;
    }
    ofs.close();

    int range;
    cout << "Please enter tne range os ypur system" << endl;
    cin >> range;
 
    // Allocate main matrix
    double **matrix = new double*[range];
    for (int i = 0; i < range; i++) 
    {
        matrix[i] = new double[range];
    }
 
    // adjacency matrix
    double **matrix2 = new double*[range];
    for (int i = 0; i < range; i++) 
    {
        matrix2[i] = new double[range];
    }
    // reading matrix from file
    ifstream ifs("A.txt");
 
    while ( ifs >> q )
    {
        matrix[i][j] = q;
        j++;
        if ( j % range == 0 ) { i++; j = 0;}
    }
    ifs.close();