#include <iostream>
#include <cmath>
using namespace std;
const int n = 4;
 // функция оценки прогрешности
double Delta(double *X, double *x, double a);
void PrintVector(double *X);
int main()
{
  // исходные данные ===================================================
  // Матрица А
    double A[n][n] = {  
                        { 0.53,       0.011,  0.035,  0.09 },
                        { -0.73,       0.58,  0,  0 },
                        { 0.154,  0.012,       0.42,  0.03 },
                        { 0,  -0.061,  0.02,       0.32 },
                     }  ;              
    // вектор б
    double b[n] = { 1.39,  1.74,  -2.05,  -1.73 };  
    //  норма
    double a =0.299;
    // эпсилон
    double eps = 0.001;
  // вычисления ========================================================
    // счетчик итераций
    int k = 0;
    // вектор приблиения
    double x[n];   
    // вектор с решением системы
    double X[n];
    // итерации
    while (true)
    {     
        // если первая итерация-в качествеприближения берем вектор б
        if (!k) 
            memcpy(x, b, n*sizeof(double));
        else
            memcpy(x, X, n*sizeof(double));
        // обнуляем вектор с решением системы
        memset(X, 0, n*sizeof(double));    
        // подсчитываем значения
        for (int i=0; i<n; i++)
        {
            for (int j=0; j<n; j++)
            {
                if (i==j) continue;
                X[i] += A[i][j]*x[j];
            }
            X[i] += b[i];
        }
        // выводим решение системы
        PrintVector(X); 
        k++;  
        // оцениваем погрешность
        // если удовлетворяет заданной - прекращаем вычисления
        if (Delta(X, x, a) <= eps) break;
    }
    // выводи количество итераций
    std::cout << "\n Count of iteration is " << k << std::endl;
    system("pause");
    return 0;
} 
double Delta(double *X, double *x, double a)
{
    double max = 0;
    for (int i=0; i<n; i++)
    {
       double f = fabs(X[i]-x[i]);
       if (f > max) max = f;
    }
    return (a/(1-a))*max;
} 
void PrintVector(double *X)
{
    for (int i=0; i<n; i++)
        std::cout << " X["<< i+1 << "] = "  << X[i] << std::endl;
    std::cout << "------" << std::endl;        
}