Метод Гаусса - Зейделя, программа работает неправильно

#include <iostream>
#include <cmath>
#include <vector>
 
using namespace std;
 
int main()
{
    setlocale(LC_ALL, "Russian");
    int n, p = 0;
    cout << "Размер квадратной матрицы = ";
    cin >> n;
    double e, k, C;
    std::vector<double> b(n), F(n, 0), x(n, 0), T(n, 0);
    std::vector<std::vector<double>> a(n, std::vector<double>(n));
    cout << "[a].[x]=[b]" << "\n";
    cout << "Заполните элемент a:" << "\n";
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; j++) {
            cout << "a[" << i << "," << j << "] = ";
            cin >> a[i][j];
        }
    }
    cout << "Заполните элемент дополнительной матрицы b:" << "\n";
    for (int j = 0; j < n; j++)
    {
        cout << "b[" << j << "] = ";
        cin >> b[j];
    }
    cout << "Введите точность = ";
    cin >> e;
 
    k = 1e10;
    while (k >= e) {
        p = p + 1;
        for (int i = 0; i < n; i++) {
            double B = 0;
            T[i] = b[i] / a[i][i];
            C = a[i][i];
            for (int j = 0; j < n; j++) {
                if (j != i)
                    B += a[i][j] * x[j];
            }
            x[i] = T[i] - B / C;
        }
        cout << x[0] << "\n";
        cout << x[1] << "\n";
 
        for (int i = 0; i < n; i++)
        {
            F[i] = 0.0;
            for (int j = 0; j < n; j++) {
                F[i] = F[i] + a[i][j] * x[j];
            }
        }
        for (int i = 0; i < n; i++) {
            F[i] = F[i] - b[i];
        }
        k = 0.0;
        for (int i = 0; i < n; i++) {
            k += F[i] * F[i];
        }
        k = std::sqrt(k);
    }
 
    cout << "Решение :\n";
    for (int i = 0; i < n; i++) {
        cout << "x" << i << "=" << x[i] << "\n";
    }
 
    cout << "Количество итераций: " << p << "\n";
    for (int i = 0; i < n; i++) {
        cout << "F" << i << "=" << F[i] << "\n";
    }
    cout << "Итоговая погрешность = " << k << "\n";
    return 0;
}