Как реализовать проверку решения системы линейных уравнений по методу Зейделя?

#include "stdafx.h"
#include <fstream>
#include <iostream>
#include <conio.h>
#include <stdlib.h>
#include <string>
#include <iostream>
#include <iomanip>
using namespace std;
void ArrayInit(float **A, int N, int M);// Заполнение массива с клавиатуры
void entermass(float *A, int N);
void ArrayShow(float **A, int N, int M);
void privodVid(float **A, float *B, int N, int M);
void uravnenieVid(float **A, float *B, int N, int M);
void matrixShow(float **A, float *B, int N, int M);
void Itershow(float **A, float *X0, float *B, int N, int M, int iter, int t);
void nullout(float *X0, int M);
int main()
{
    setlocale(LC_ALL, "RUS");
    char simvol;
    int matrica,N,M, iterac, tochnost,numbers=0;
    float **A;
    cout << "\n Введите размеры Уравнения: ";
    cin >> N >> M;
    A = new float*[N];
    for (int i = 0; i < N; i++) A[i] = new float[M];
    //создаем матрицу
    ArrayInit(A, N, M);
    cout << "Вы ввели \n";
    ArrayShow(A, N, M);
    float *B = new float[M];
    cout << "\n Введите точность знаков после запятой: "; cin >> tochnost;
    cout << "\n Введите количество итераций: "; cin >> iterac;
    //Заполняем матрицу
    cout << "____________________________\n";
    //Заполняем матрицу
    cout << "\n Введите В \n\n";
    entermass(B, M);
    cout << "____________________________\n\n";
    cout << " Уравнение имеет вид: \n\n";
    //вывод матрицы
    uravnenieVid(A, B, N, M);
    cout << "\n Приведем эту систему к виду, удобному для итерации\n\n";
    //-----------------------------------------------------—
    privodVid(A, B, N, M);
    //-----------------------------------------—
    //вывод матрицы
    matrixShow(A, B, N, M);
    //--------------------------------------------------
    float *X0;
    X0 = new float[N];
    X0[1 - N] = 0;
    for (int i = 0; i < N; i++) 
    {
        X0[i] = 0;
    }
    X0[0] = B[0];
    cout << "\n В качестве нулевых приближений корней возьмем :\n\n";
    //вывод нулевых
    nullout(X0, M);
    cout << "======================\n";
    Itershow(A,X0, B, N, M, iterac, tochnost);
    cout << " \n \n";
    system("Pause");
    return 0;
 
}
void ArrayInit(float **A, int N, int M)
{
    for (int i = 0; i < M; i++)
    {
        for (int j = 0; j < N; j++)
        {
            cout << "Введите элемент матрицы: " << "[" << i << "][" << j << "]";//Ввод чисел в двумерный массив А
            cin >> A[i][j];
        }
    }
}
void entermass(float *A, int N)
{
    for (int i = 0; i < N; i++)
    {
        cout << "Введите элемент строки: " << "[" << i << "]";//Заполнение строки одномерного массива
        cin >> A[i];
    }
}
void ArrayShow(float **A, int N, int M)
{
    for (int i = 0; i < M; i++)
    {
        for (int j = 0; j < N; j++)
        {
            cout << " " << A[i][j];
        }
        cout << "\n";
    }
}
void privodVid(float **A, float *B, int N, int M)
{
    float vrem = 0;
    int numbers = 0;
    //приводим к правильному виду
    for (int i = 0; i < M; i++) {
        vrem = A[i][numbers];
        numbers++;
        for (int j = 0; j < N; j++) {
            A[i][j] = A[i][j] / vrem;
            if (j == M - 1)
                B[i] = B[i] / vrem;
        }
    }
}
void uravnenieVid(float **A, float *B, int N, int M)
{
    for (int i = 0; i < M; i++)
    {
        for (int j = 0; j < N; j++)
        {
            if (j == 0)
            {
                cout << "\t" << A[i][j] << "x" << j + 1 << " ";
            }
            else
            {
                if (A[i][j] >= 0)
                    cout << "+\t" << A[i][j] << "x" << j + 1 << " ";
                else
                {
                    cout << "+\t" << A[i][j] << "x" << j + 1 << " ";
                }
            }
            if (j == M - 1)
            {
                cout << "\t= " << B[i] << " ";
            }
        }
        cout << "\n";
    }
}
void matrixShow(float **A, float *B, int N, int M)
{
    int numbers;
    for (int i = 0; i < M; i++)
    {
        cout << "\tx" << i + 1 << " = ";
        numbers = 0;
        for (int j = 0; j < N; j++)
        {
            if (i != j)
            {
                A[i][j] *= -1;
                if (numbers == 0)
                {
                    numbers++;
                    cout << "\t" << A[i][j] << "x" << j + 1 << " ";
                }
            }
            if (j == M - 1)
                if (B[i] >= 0)
                    cout << "+\t" << B[i] << " ";
                else
                    cout << "+\t(" << B[i] << ") ";
        }
        cout << "\n";
    }
}
void Itershow(float **A,float *X0, float *B, int N, int M, int iter,int t)
{
    float prov = 0;
    float prov2 = 0;
    float otvet = 0;
    float vivod = 0;
    float buf=0;
    int p = 0;
    float vrem1;
    for (int k = 0; k < iter; k++)
    {
        cout << " Итерация №" << k + 1 << ":\n";
        for (int i = 0; i < M; i++)
        {
            vivod = 0;
            for (int j = 0; j < N; j++)
            {
                if (i != j)
                {
                    vrem1 = A[i][j] * X0[j];
                    vivod += vrem1;
                }
            
            }
            vivod += B[i];
            X0[i] = vivod;
    
            std::cout << " X" << i + 1 << "|" << k << "|=" << std::fixed << std::setprecision(t) << vivod << std::endl;
        }
        
        
        
        cout << "======================\n\n";
    }
    
}
void nullout(float *X0, int M)
{
    for (int i = 0; i < M; i++)
    {
        cout << " X" << i << "|0|=" << X0[i] << "\n";
    }
}