Метод Гаусса-Зейделя

#include "stdafx.h"
#include "math.h"
#include <iostream>
#include <iomanip>
#include <time.h>
using namespace std;
 
bool converge(double *xk, double *xkp, double eps, int n)
{
    for (int i = 0; i < n; i++) 
    {
        if (fabs(xk[i] - xkp[i]) >= eps) 
            return false;
    }
    return true;
}
 
int _tmain(int argc, _TCHAR* argv[])
{
    srand((unsigned)time(NULL));
    double *p, *x, **a, *b, eps;
    int n,i,j,k=0;
    cin >> n;
    cin >> eps,
    p = new double[n];
    x = new double[n];
    a = new double*[n];
    b = new double[n];
    for (i = 0; i < n; i++)
    {
        a[i] = new double[n];
        b[i] = ((double)rand()/RAND_MAX)*100;
        x[i] = 0.0;
        for (j = 0; j < n; j++)
        a[i][j] = ((double)rand()/RAND_MAX)*100;
    
    }
    cout << setw (13) << "A" << setw (16) << "B";
    cout<< endl;
    for (i = 0; i < n; i++)
    {
        cout<< endl;
        for (j = 0; j < n; j++)
            cout << a[i][j] << " ";
        cout << " " << b[i];
    }
    do
    {
        k++;
        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 < i; j++)
                var += (a[i][j] * x[j]);
            for (int j = i + 1; j < n; j++)
                var += (a[i][j] * p[j]);
            x[i] = (b[i] - var) / a[i][i];
        }
    }
    while (!converge(x, p, eps, n));
    cout<< endl;
    for (i = 0; i < n; i++)
    {
        cout << x[i] << endl;
    }
    cout << k << endl;
    system ("pause");
 
    return 0;
}