Решение СЛАУ методом Гауса-Зейделя. - C (СИ) - Ответ 2977806

//Метод Гауса-Зейделя для решенcия СЛАУ
//------------------------------------------------------------------------------
#include <stdio.h>
#include <alloc.h>
#include <conio.h>
#include <math.h>
//------------------------------------------------------------------------------
 
int main ()
{
  int i,j;
  int n;
  double **A, *B, *U, *X, *Xp, *Xn, max_s, max_n, eps;
 
  printf("Enter the size of the matrix:\nn = ");
  scanf("%d",&n);
 
  printf("\nEnter the values of the marix A:\n");
 
  A=(double**)malloc(n*sizeof(double*));
 
  for(i=0;i<n;i++)
  {
        A[i]=(double*)malloc(n*sizeof(double));
 
        for(j=0;j<n;j++)
        {
                printf("A[%d][%d] = ", i,j);
                scanf("%lf", &(A[i][j]));
        }
  }
 
  printf("\nEnter the values of the marix B:\n");
 
  B=(double*)malloc(n*sizeof(double));
 
  for(i=0;i<n;i++)
  {
        printf("B[%d] = ", i);
        scanf("%lf", &(B[i]));
  }
 
  printf("\nEnter pogreshnost:\neps = ");
  scanf("%lf", &eps);
 
  clrscr();
 
  X=(double*)malloc(n*sizeof(double));
  Xn=(double*)malloc(n*sizeof(double));
  Xp=(double*)malloc(n*sizeof(double));
 
  for(i=0;i<n;i++)
  {
        X[i] = 0;
        Xn[i] = 0;
        Xp[i] = 0;
  }
 
  U=(double*)malloc(n*sizeof(double));
 
  for(i=0;i<n;i++)
  {
        for(j=0;j<n;j++)
        {
                if(i==j)
                {
                        U[i]=A[i][j];
                }
        }
  }
 
  for(i=0;i<n;i++)
  {
        for(j=0;j<n;j++)
        {
                if(i==j)
                {
                        A[i][j] = 0;
                }
                else
                {
                        A[i][j] = -(A[i][j])/U[i];
                }
        }
        B[i] = (B[i])/U[i];
  }
 
  free(U);
 
  do
  {
        for(i=0;i<n;i++)
        {
                Xn[i] = X[i];
                X[i] = 0;
        }
 
        for (i=0;i<n;i++)
        {
                for (j=0;j<n;j++)
                {
                        X[i] = A[i][j]*Xp[j] + X[i];
                }
                X[i] = X[i] + B[i];
                Xp[i] = X[i];
        }
 
        max_s = 0;
        max_n = 0;
 
        for(i=0;i<n;i++)
        {
                max_n = fabs(X[i] - Xn[i]);
                if(max_n > max_s)
                {
                        max_s = max_n;
                }
        }
 
        for(i=0;i<n;i++)
        {
        printf("%.3lf ", (X[i]));
        }
        printf("  %lf\n", max_s);
 
  }while(max_s > eps);
 
  free(X);
  free(A);
  free(B);
 
  getch();
  return 0;
}
//------------------------------------------------------------------------------

double proverka (double **A, int n)
//проверка преобладания диагонали
if (...) return p = 1;
else reurn p = 0.
 
double zeydel (double **A, double *B, double *X, double eps, int n)
//решение СЛАУ
return X[i];
 
int main()
//введение с клавиатуры **A, *В, eps, n
//вывод результатов
1) СЛАУ не может быть решена
2) Решение