СЛАУ: метод простой итерации и метод Зейделя

#include <iostream>
#include <cmath>
using namespace std;
 
// Условие окончания
bool converge(double xk[10], double xkp[10], int n, double eps)
{
    double norm = 0;
    for (int i = 0; i < n; i++)
        norm += (xk[i] - xkp[i]) * (xk[i] - xkp[i]);
    return (sqrt(norm) < eps);
}
 
double okr(double x, double eps)
{
    int i = 0;
    double neweps = eps;
    while (neweps < 1)
    {
        i++;
        neweps *= 10;
    }
    int okr = pow(double(10), i);
    x = int(x * okr + 0.5) / double(okr);
 
    return x;
}
 
bool diagonal(double** a, int n)
{
    int i, j, k = 1;
    double sum;
    for (i = 0; i < n; i++) {
        sum = 0;
        for (j = 0; j < n; j++) sum += abs(a[i][j]);
        sum -= abs(a[i][i]);
        if (sum > a[i][i])
        {
            k = 0;
            std::cout << a[i][i] << " < " << sum << std::endl;
        }
        else
        {
            std::cout << a[i][i] << " > " << sum << std::endl;
        }
 
 
    }
 
    return (k == 1);
 
}
 
 
 
void Jacobi(int N, double** A, double* F, double* X, double eps)
{
    double* TempX = new double[N];
    double norm; // норма, определяемая как наибольшая разность компонент столбца иксов соседних итераций.
 
    do {
        for (int i = 0; i < N; i++) {
            TempX[i] = F[i];
            for (int g = 0; g < N; g++) {
                if (i != g)
                    TempX[i] -= A[i][g] * X[g];
            }
            TempX[i] /= A[i][i];
        }
        norm = fabs(X[0] - TempX[0]);
        for (int h = 0; h < N; h++) {
            if (fabs(X[h] - TempX[h]) > norm)
                norm = fabs(X[h] - TempX[h]);
            X[h] = TempX[h];
        }
    } while (norm > eps);
    delete[] TempX;
}
 
 
int main()
{
    setlocale(LC_ALL, "");
 
    double eps;
    double** a, * b, * x, * p;
    int n, i, j, m = 0;
    int method;
    std::cout << "Введите размер квадратной матрицы: ";
    std::cin >> n;
    a = new double* [n];
    for (int i = 0; i < n; i++) {
        a[i] = new double[n];
    }
    b = new double[n];
    x = new double[n];
    p = new double[n];
 
    std::cout << "Введите точность вычислений: ";
    std::cin >> eps;
    std::cout << "Заполните матрицу А: " << std::endl << std::endl;
    for (i = 0; i < n; i++)
        for (j = 0; j < n; j++)
        {
            cout << "A[" << i+1 << "][" << j+1 << "] = ";
            cin >> a[i][j];
        }
    std::cout << std::endl << std::endl;
    std::cout << "Ваша матрица А: " << std::endl << std::endl;
    for (i = 0; i < n; i++)
    {
        for (j = 0; j < n; j++)
            cout << a[i][j] << " ";
        std::cout << std::endl;
    }
 
    std::cout << std::endl;
 
    std::cout << "Заполните столбец свободных членов: " << std::endl << std::endl;
    for (i = 0; i < n; i++)
    {
        std::cout << "В[" << i + 1 << "] = ";
        std::cin >> b[i];
    }
 
    std::cout << "Система уравнений" << std::endl;
 
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            std::cout << a[i][j] << "*x" << j + 1;
            if (j + 1 < n)
                std::cout << " + ";
 
        }
        std::cout << " = " << b[i] << std::endl;
    }
 
    std::cout << std::endl << std::endl;
 
 
    /*
    Ход метода, где:
    a[n][n] - Матрица коэффициентов
    x[n], p[n] - Текущее и предыдущее решения
    b[n] - Столбец правых частей
    Все перечисленные массивы вещественные и
    должны быть определены в основной программе,
    также в массив x[n] следует поместить начальное
    приближение столбца решений (например, все нули)
    */
 
    for (int i = 0; i < n; i++)
        x[i] = 1;
 
    std::cout << "Диагональное преобладание: "  << std::endl;
    if (diagonal(a, n)) {
        do
        {
            for (int i = 0; i < n; i++)
                p[i] = x[i];
            for (int i = 0; i < n; i++)
            {
                double var = 0;
                for (int j = 0; j < n; j++)
                    if (j != i) var += (a[i][j] * x[j]);
 
                x[i] = (b[i] - var) / a[i][i];
                
            }
            std::cout << x[i] << std::endl;
            m++;
        } while (!converge(x, p, n, eps));
        
        std::cout << "Условие сходимости выпоняется" << std::endl << std::endl;
 
 
        
        std::cout << "Решение системы:" << std::endl << std::endl;
        std::cout << "Метод Зейделя" << std::endl;
        for (i = 0; i < n; i++) std::cout << "x" << i << " = " << okr(x[i], eps) << "" << std::endl;
        std::cout << "Итераций: " << m << std::endl << std::endl;
        std::cout << "Метод простых итераций: " << std::endl;
        double* xb = new double[n];
        Jacobi(n, a, b, xb, eps);
        for (int i = 0; i < n; i++)cout << "x" << i << " = " << xb[i] << "" << std::endl;
    }
    else {
        std::cout << "Условие сходимости не выполняется" << std::endl;
    }
 
    //  #перенаправление [[Название целевой страницы]]
    system("pause");
    return 0;
}