Решение СЛАУ методом простых итераций

#include "stdafx.h"
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
#define n 4
 
void del(double **A, double *X, double *Xk,   double *b);
 
int main()
 
{
 
    int i,j;
    double **A;     //matrix equations
    double *X;      //communications solution and iteration (k+1)
    double *Xk;     //iteration (k)
    double *b;      //free to members of
    double sumd = 0;//additional variables for computation
    double eps;
 
//--------------dynamic memory allocation------------------------------//
    A  = new double* [n];
    X  =  new double [n];
    Xk =  new double [n];
    b  =  new double [n];
    for(i = 0; i < n; i++)
        A[i] = new double [n];
//-------------input system of equations-------------------------------//
    printf("Type elements in the variables in the system: \n");
    for(i = 0; i < n; i++)
       for(j = 0; j < n; j++)
        scanf("%lf",&A[i][j]);
 
 
    printf("Members enter free system: \n");
    for(i = 0; i < n; i++)
        scanf("%lf",&b[i]);
    printf("enter the calculation accuracy: \n");
        scanf(" %lf",&eps);
//-------------test system on the diagonal elements--------------------//
    for(i = 0; i < n; i++)
       if (A[i][i] == 0) {del( A, X, Xk, b); printf("!On the main diagonal of the system is zero"); exit(0);}
//-------------checking for convergence of-----------------------------//
    for(i = 0; i < n; i++)
      {
       for(j = 0; j < n; j++)
          if (i != j)       sumd += fabs(A[i][j]);
       if (fabs(A[i][i]) < sumd)  {del( A, X, Xk, b);  printf("!Diverging methods for this system"); exit(0);}
       sumd = 0;
      }
//----- approximate solutions of the system-//
    for(i = 0; i < n; i++)
         X[i]=0.0;
    
//-------------calculation system of linear algebraic equations---------//
    int  count = 0, flag = 1;
    double x, zh = 0;
    for(i = 0; i < n; i++)
    Xk[i] = X[i];
    do
      {
       count++;
       for(i = 0; i < n; i++)
     {
       x=0;
       for(j = 0; j < n; j++)
          {
         if (i != j)  x += Xk[j] * A[i][j];
         if (i == j) zh = A[i][j];
          }
       x = (b[i] - x) / zh;
       X[i] = x;
       if ((fabs(X[i]-Xk[i]))<=eps) flag=0;
     }
     for(i = 0; i < n; i++)
        Xk[i] = X[i];
      }while(flag);
    printf("solution call method of iterations: \n");
    for(i = 0; i < n; i++)
       printf("x[%d]=%0.5lf\n",i+1,X[i]);
       printf("count=%d\n",count);
    
       del( A, X, Xk, b);
       return 0;
}
 
 
void del(double **A, double *X, double *Xk,   double *b)
  {
    for(int i = 0; i < n; i++)
      delete[]A[i];
     delete[]A;
     delete[]X;
     delete[]Xk;
     delete[]b;
  }