Метод Зейделя СЛАУ

#include <iostream>
#include <math.h>
using namespace std;
 
double A[5][5], T[5][5], X[5], Y[5], C[5][5], D[5];
int n;
double E;
 
void vvodsort()
{
    int i, j, k, max;
    cout << "Похибка: ";
    cin >> E;
    cout << "Розмiрнiсть матрицi: ";
    cin >> n;
    cout << "Коефiцiєнти матрицi A:\n";
    for (i = 0; i<n; i++)
        for (j = 0; j<n; j++){
            cin >> A[i][j];
            T[i][j] = A[i][j];
        }
    cout << "Матриця вiльних членiв B:\n";
    for (i = 0; i<n; i++) cin >> A[i][j];
 
    for (i = 0; i<n; i++){
        for (j = i; j<n; j++){
            if (fabs(A[j][i])>fabs(A[i][i]))max = j;
        }
 
        for (k = 0; k<n; k++)
            Y[k] = A[i][k];
        int y;
        y = A[i][n];
        A[i][n] = A[max][n];
        A[max][n] = y;
        for (k = 0; k<n; k++)
        {
            A[i][k] = A[max][k];
            A[max][k] = Y[k];
        }
    }
}
 
int proverka()
{
    int i, j, k, z = 0; double sum;
 
    for (i = 0; i<n; i++){
        sum = 0;
        for (j = 0; j<n; j++)
            sum = sum + fabs(A[i][j]);
        if (fabs(A[i][i])<fabs(sum - fabs(A[i][i]))){ z = 1; break; }
    }
    if (z == 0)return 1; else return 0;
}
 
void zegeli()
{
    int i, j;
    for (i = 0; i<n; i++)
    {
        Y[i] = X[i] = A[i][n] / A[i][i];
        D[i] = A[i][n] / A[i][i];
    }
 
    for (i = 0; i<n; i++)
        for (j = 0; j<n; j++)
        {
            if (i == j)C[i][j] = 0;
            else C[i][j] = A[i][j] / (-A[i][i]);
        }
 
    int z = 0, c = 0;
    while (z != n)
    {
        for (i = 0; i<n; i++)
            Y[i] = X[i];
        double t = 0;
        for (i = 0; i<n; i++)
        {
            for (j = 0; j<n; j++)
                t = t + C[i][j] * X[j];
            X[i] = t + D[i];
            t = 0;
        }
        z = 0;
        for (i = 0; i<n; i++)
            if (X[i] - Y[i]<E)z++;
        c++;
    }
 
    cout << "\nKorni \n ";
    for (i = 0; i<n; i++)
        cout <<"x["<<i<<"]="<< X[i] << " ";
    cout << "\nКiлькiсть iтерацiй = " << c;
    cout << "\nПеревiрка ";
    z = n;
    double t = 0;
    for (i = 0; i<n; i++){
        cout << "\n";
        for (j = 0; j<n; j++){
            cout << A[i][j] << "X" << j + 1;
            if (A[i][j + 1]<0);
            else if (j != 2)cout << "+";
            t = t + A[i][j] * X[j];
        }
        cout << "=" << t;
        t = 0;
    }
}
 
void main()
{
    setlocale(0, "");
    system("CLS");
    vvodsort();
    if (proverka()) {
        cout << "\nUslovie vipolniaetsea ! ";
        zegeli();
    }
    else cout << "\nUslovie ne vipolniaetsea ! ";
}
/*
1.10 0.12 -0.13
0.12 1.71 0.15
-0.13 0.15 1.63 */