Метод простых итераций

#include <iostream>
#include <cstdlib>
#include <iomanip>
#include <cmath>
 
#define EPS 1e-10
 
 
 
 
using namespace std;
using myType = double;
 
class SLE
{
 
private:
    myType** _arguments;
    myType* _freeMembers;
    int _rowsCount;
    int _columnsCount;
 
public:
    SLE(int rowsC, int columnsC);
    void setUserValus();
    void setRandomValues();
    void showSLE();
    bool checkConvergence();
    bool checkInfelicity(myType* decisions, myType* tempDecisions);
    myType* findDecision();
 
 
 
 
};
 
 
int main()
{
    SLE newSLE(4,4);
    newSLE.setUserValus();
    newSLE.showSLE();
    newSLE.checkConvergence();
    newSLE.findDecision();
    return 0;
}
 
 
 
 
SLE::SLE(int rowsC, int columnsC)
{
    _columnsCount = columnsC;
    _rowsCount = rowsC;
    _arguments = new myType* [_rowsCount];
    _freeMembers = new myType [_columnsCount];
 
    for(int i = 0; i < _rowsCount; i++)
        _arguments[i] = new myType [_columnsCount];
}
 
void SLE::setRandomValues()
{
    for(int i = 0; i <_rowsCount; i++)
    {
        for(int j = 0; j < _columnsCount; j++)
        {
            _arguments[i][j] = rand()%100 - 50;
        }
 
        _freeMembers[i] = rand()%100 - 50;
    }
 
}
 
void SLE::showSLE()
{
    for(int i = 0; i < _rowsCount; i++)
    {
        for(int j = 0; j < _columnsCount; j++)
        {
            cout << setw(10) << _arguments[i][j];
        }
        cout << "   " << "(" << setw(5)  <<_freeMembers[i] << setw(5) << ")" << endl;
    }
}
 
bool SLE::checkConvergence()
{
    int temp;
    for(int i = 0; i < _rowsCount; i++)
    {
        myType sum = 0;
        for(int j = 0; j < _columnsCount; j++)
        {
            if(i != j)
                sum += abs(_arguments[i][j]);
        }
 
        if(abs(_arguments[i][i]) < sum)
        {
            cout << "не определяется" << endl;
            return false;
        }
 
    }
    cout << "Определяется" << endl;
 
    return true;
}
 
void SLE::setUserValus()
{
    cout << "введите аргументы" << endl;
 
    for(int i = 0; i < _rowsCount; i++)
        for(int j = 0; j < _columnsCount; j++)
        {
            cout << i << "  " << j<< " : ";
            cin >> _arguments[i][j];
            cout << endl;
        }
    cout << "введите свободные члены" << endl;
    for(int i = 0; i < _columnsCount; i++)
    {
        cout << i << " : ";
        cin >> _freeMembers[i];
        cout << endl;
 
    }
 
}
 
 
 
myType* SLE::findDecision()
{
    cout << "FIND" << endl;
    myType* decision = new myType [_rowsCount];
    myType* tempDecision = new myType [_rowsCount];
    int temp;
    bool flag = false;
 
    for(int i = 0; i < _rowsCount; i++)
    {
        decision[i] = 0;
        tempDecision[i] = 0;
    }
 
    do
    {
        for(int i = 0; i < _rowsCount; i++)
        {
            decision[i] = _freeMembers[i];
 
            for(int j = 0; j < _columnsCount; j++)
            {
                if( i != j)
                    decision[i] -= _arguments[i][j] * tempDecision[j];
 
            }
 
            decision[i] = decision[i]/_arguments[i][i];
        }
 
 
        flag = checkInfelicity(decision, tempDecision);
 
        for(int i = 0; i < _rowsCount; i++)
            tempDecision[i] = decision[i];
 
 
    }while(!flag);
 
 
    for(int i = 0; i < _rowsCount; i++)
    {
        cout << "x" << i << "  = " << decision[i] << "  " ;
    }
    cout << endl;
 
    return decision;
 
}
 
bool SLE::checkInfelicity(myType* decisions, myType* tempDecisions)
{
 
 
    for(int i = 0; i < _rowsCount; i++)
    {
        if(abs(decisions[i] - tempDecisions[i]) > EPS)
            return false;
    }
 
    return true;
}