Метод Гаусса-Зейделя для решения СЛАУ

#include <iostream>
#include <iomanip>
#include <cmath>
 
using namespace std;
class Matr
{
private:
    double eps = 0.001;
    static const int n = 4;
    double masA[n][n] = { { 0.05, -0.06, -0.12, 0.14 },
                  { 0.04, -0.12, 0.68, 0.11 },
                  { 0.34, 0.08, -0.06, 0.44 },
                  { 0.11, 0.12, 0, -0.03 } };
 
    double masB[n] = { -2.17, 1.4, -2.1, -0.8 };
    void Zeidel()
    {
        double x[n];       // решения n-ой итерации
        double x0[n];      // решения прошлой итерации
        double E[n];       // разница между соответствующими эл-ми
        double max = 0;    // для сравнения с eps
 
        for (int i = 0; i < n; i++)
       {
          x0[i] = masB[i]; // для нахождения первого приближения
        }
 
        int counter = 0;   // переменная для подсчёта итераций
 
        do
        {
            for (int i = 0; i < n; i++)   // цикл для нахождения неизвестных
            {
                x[i] = 0;                      
 
                for (int j = 0; j < n; j++)    
                {
                    x[i] += masA[i][j] * x0[j];     
                }
                x[i] += masB[i];    // прибавляем согласно формуле каждому x значение свободного члена
                E[i] = fabs(x[i] - x0[i]);  // разница между найденным и прошлым значением итерации
                x0[i] = x[i];    // присваивание значения переменной x[0]
            }
 
            max = 0;
            for (int i = 0; i < n; i++)
            {
                if (max < E[i])  // сравнение max и разницы
                    max = E[i];         // приравнение 
            }
            counter++;
        } while (max > eps);
 
        cout << endl << "Количество итераций: " << counter << endl << endl;
 
        cout << "Найденые иксы:" << endl;
        for (int i = 0; i < n; i++)
        cout << "x" << i + 1 << " = " << x[i] << " " << endl;