Система линейных алгебраических уравнений

#include <vector>
#include <iomanip>
#include <iostream>
#include <cmath>
#include <locale.h> 
#include <stdio.h>
#include <fstream>
#include <Windows.h>
 
using namespace std;
 
int main()
{
    int k=1;
    char *str = new char [256];
    int n=0;
    ifstream base("file.txt");
    while (!base.eof())
    {
        base.getline(str, 256, '\n');
        n++;
    }
    base.close();
    delete str;
    int m=n+1;
    fstream fi;
    fi.open("file.txt", ios::in);
    float ** matr=new float *[n];
    for (int i=0; i<n; i++) 
    {
        matr[i]=new float [m];
    }
    int i, j;
    if (!fi)
    { 
        cerr << "error file doesn't exist\n"; 
        exit(1);
    }
    for(i=0;i<n;i++)
    {
        for(j=0;j<m;j++)
        {
            fi>>matr[i][j];
        }
    }
    for (i=0; i<n; i++)
    {
        for(j=0; j<m; j++)
        {
            cout<<setw(10)<<matr[i][j]<<"    ";
        }
        cout<<"\n";
    }
    fi.close();
    long double eps=0.001;
    // Введем вектор значений неизвестных на предыдущей итерации,
    // размер которого равен числу строк в матрице, т.е. size,
    // причем согласно методу изначально заполняем его нулями
    vector <long double> previousVariableValues (n, 0);
    // Будем выполнять итерационный процесс до тех пор, 
    // пока не будет достигнута необходимая точность    
    while (true)
    {
        
        // Введем вектор значений неизвестных на текущем шаге       
        vector <long double> currentVariableValues (n);
        // Посчитаем значения неизвестных на текущей итерации
        // в соответствии с теоретическими формулами
        cout<<"Reshenie na"<<' '<<k<<' '<< "iteracii:"<<endl;
        for (int i = 0; i < n; i++)
        {
            // Инициализируем i-ую неизвестную значением 
            // свободного члена i-ой строки матрицы
            currentVariableValues[i] = matr[i][n];
            // Вычитаем сумму по всем отличным от i-ой неизвестным
            for (int j = 0; j < n; j++)
            {       
                    currentVariableValues[i] = currentVariableValues[i] + matr[i][j] * previousVariableValues[j];
            }
        
        cout<<"x"<< i+1<< "="<< currentVariableValues[i]<<endl;
        }
        
        // Посчитаем текущую погрешность относительно предыдущей итерации
        long double max = 0;
        for (int i = 0; i < n; i++)
        {
            max=fabs (currentVariableValues[0] - previousVariableValues[0]);
            if(max<fabs (currentVariableValues[i] - previousVariableValues[i]))
            {   
                max=fabs (currentVariableValues[i] - previousVariableValues[i]);
            }
        }
        // Если необходимая точность достигнута, то завершаем процесс
        if (max < eps)
        {
            break;
        }
        // Переходим к следующей итерации, так 
        // что текущие значения неизвестных 
        // становятся значениями на предыдущей итерации
        previousVariableValues = currentVariableValues;
        k++;
    }
    // Выводим найденные значения неизвестных с 8 знаками точности
    cout<<endl;
    cout<< "Reshenie s neobxodimoi tochnost'u"<<endl;
    for ( i = 0; i < n; i++)
    {
        cout<<"x"<< i+1<< "=" << previousVariableValues[i]<<endl;
    }
    system("pause");
    return 0;
}