Написать программу для вычисления СЛАУ методом простой итерации с использованием динамических массивов

//////////////////////////////////////////////////////////////////////////////
//  Solving nonlinear equations (iterations method)
//  (c) Johna Smith, 1996
//
//  Method description:
//    This method is used for solving equations like x=f(x)
//    We should set first approximation of the root and then
//    make some iterations x=f(x). When the difference between
//    two consequental iterations will be less than epsilon we
//    can say that we found the root.
//    To apply this method to equations like F(x)=0 we should
//    transform it into iterational form x=f(x), where f(x):
//      1) defined on [a;b]
//      2) for every point on [a;b] exists f'(x)
//      3) for every x from [a;b] f(x) is in [a;b]
//      4) exists number q : |f'(x)|<=q<1 for all x from [a;b]
//
//////////////////////////////////////////////////////////////////////////////
#include <stdio.h>
#include <math.h>
// this function returns value of f(x)
double f(double x)
{
  // x = x + 0.37(sin 2x - ln x)
  return x+0.37*(sin(2*x)-log(x));
}
double Solve(double x, double epsilon)
{
  double x1;
  do
  {
    x1=x;
    x=f(x);
  } while (fabs(x1-x)>epsilon);
  return x;
}
void main(void)
{
  printf("sin 2x - ln x = 0");
  printf("\nx0=%f",Solve(1.4,1e-6));
}

//////////////////////////////////////////////////////////////////////////////
//
//  Solving system of linear equations (iterations method)
//  (c) Johna Smith, 1996
//
//  Method description:
//  This program automatically transforms system of equations to
//  iterative form and solves it.
//  Iterative form:
//     x1=a11*x1+a12*x2+... + a1n*xn+b1
//     x2=a21*x1+a22*x2+... + a2n*xn+b2
//     ...
//     xn=an1*x1+an2*x2+... + ann*xn+bn
//  Setting first iteration of vector X we can get next iteration from
//  such form of system of equations. Iterations will be finished when
//  difference between two consequitive iterations will be less than epsilon
//
//////////////////////////////////////////////////////////////////////////////
#include <stdio.h>
#include <math.h>
#define N      3     // size of matrix
#define N1     N+1
float matrix[N][N1]=
    {{14.38,-2.41,1.39,5.86},
  {1.84,25.36,-3.31,-2.28},
  {2.46,-3.49,16.37,4.47}
    };
int maxiterations=10;  // maximum number of iterations
float epsilon=0.0001;  // required accuracy
void ShowMatrix(void)
{
  for (int i=0;i<N;i++)
  {
    for (int j=0;j<N;j++)
      printf("%+f*x%d",matrix[i][j],i+1);
    printf("=%f\n",matrix[i][N]);
  }
}
void main(void)
{
  // Variables declaration
  float x[N],y[N],t;
  register short int i,j,k;
  int iterations=0;
  // Printing given matrix
  ShowMatrix();
  // setting first iteration of vector X
  for (i=0;i<N;i++) x[i]=matrix[i][N];
  do
  {
    for (i=0;i<N;i++)
    {
      t=-matrix[i][N];
      for (j=0;j<N;j++)
        t+=matrix[i][j]*x[j];
        y[i]=(-t+matrix[i][i]*x[i])/matrix[i][i];
    }
    iterations++;
    k=0;
    // checking solution
    while (fabs(x[k]-y[k])<epsilon && k<N) k++;
    // new iteration becomes old
    for(i=0;i<N;i++) x[i]=y[i];
  } while (k!=N && iterations<maxiterations);
  if (iterations==maxiterations)
  {
    printf("Iterations are very slow...");
  } else
  {
    // Printing solution
    printf("\nSolution:\n");
    for (i=0;i<N;i++)
      printf("x%d=%f\n",i+1,x[i]);
    printf("%d iterations were made",iterations);
  }
}

#include <iostream>
#include <windows.h>
#include <math.h>
 
using namespace std;
 
int main()
{
    SetConsoleOutputCP(1251);
    char con;
    int n;
    double delta=0.0, eps=0.0;
    do
    {
        cout << " Введите количество неизвестных: " << endl;
        cin >> n;
        cout << " Введите точность: " << endl;
        cin >> eps;
 
        float **a = new float*[n];//выдиление динам. памяти для массива
        for(int i=0; i<n; i++)
        {
            a[i] = new float[n];
        }
        float *b = new float[n];// выдиление памяти для свободных членов
        float *x = new float[n];
        float *xp = new float[n];
 
        for(int i=0; i<n; i++)// ввод данных
        {
            for(int j=0; j<n; j++)
            {
                cout << " Введите " << (j+1) << "й коэффициент " << (i+1) << "го уравнения: ";
                cin >> a[i][j];
            }
            cout << " Введите " << (i+1) << "й свободный член " << (i+1) << "го уравнения: ";
            cin >> b[i];
        }
 
        for(int i=0; i<n; i++)// вывод СЛАУ
        {
            cout << " | ";
            for(int j=0; j<n; j++)
            {
                cout << a[i][j] << "*x(" << (i+1) << ") + " ;
            }
            cout << " = " << b[i];
            cout << endl;
        }
 
        do
        {
            for(int i=0; i<n; i++) // НЕ РАБОТАЕТ(((
            {
                for(int j=1; j<n; j++)
                {
                    if(j!=i)
                    {
                        xp[i]=b[i]/a[i][i]; // начальное значение Х
                        x[i]=(b[i]-(a[i][j]*xp[j]))/a[i][i];
                        if(fabs(x[i]-xp[i])>delta)
                            delta=fabs(x[i]-xp[i]);
                    }
                }
            }
        }
        while(delta>eps);
 
        for(int i=0; i<n; i++)// вывод результата
        {
            cout << "  x(" << i << ") = " << x[i] << endl;
        }
 
 
 
 
        for(int i=0; i<n; i++)// освобождение выдиленой памяти
        {
               delete [] a[i];
        }
        delete [] a;
        delete [] b;
        delete [] x;
        delete [] xp;
        cout << " Для продолжения программы, введите - 'y'" << endl;
        cout << " для завершения, введите любой символ..." << endl;
        cin >> con;
    }
    while(con=='y');
    return 0;
}

#include <iostream>
#include <windows.h>
#include <math.h>
 
using namespace std;
 
int main()
{
    SetConsoleOutputCP(1251);
    char con;
    int n;
    double delta=0.0, eps=0.0;
    do
    {
        cout << " Введите количество неизвестных: " << endl;
        cin >> n;
        cout << " Введите точность: " << endl;
        cin >> eps;
 
        double **a = new double*[n];//выдиление динам. памяти для массива
        for(int i=0; i<n; i++)
        {
            a[i] = new double[n];
        }
        double *b = new double[n];// выдиление памяти для свободных членов
        double *x = new double[n];
        double *xp = new double[n];
 
        for(int i=0; i<n; i++)// ввод данных
        {
            for(int j=0; j<n; j++)
            {
                cout << " Введите " << (j+1) << "й коэффициент " << (i+1) << "го уравнения: ";
                cin >> a[i][j];
            }
            cout << " Введите " << (i+1) << "й свободный член " << (i+1) << "го уравнения: ";
            cin >> b[i];
        }
 
        for(int i=0; i<n; i++)// вывод СЛАУ
        {
            cout << " | ";
            for(int j=0; j<n; j++)
            {
                cout << a[i][j] << "*x(" << (i+1) << ") + " ;
            }
            cout << " = " << b[i];
            cout << endl;
        }
 
        do
        {
            for(int i=0; i<n; i++) // НЕ РАБОТАЕТ(((
            {
                xp[i]=b[i]/a[i][i]; // начальное значение Х
                for(int j=0; j<n; j++)
                {
                    if(j!=i)
                    {
                        x[i]=(b[i]-(a[i][j]*xp[j]))/a[i][i];
                    }
                }
                delta=fabs(xp[0]-x[0]);
                if(fabs(xp[i]-x[i])>delta)
                    delta=fabs(xp[i]-x[i]);
                xp[i]=x[i];
            }
        }
        while(delta>eps);
 
        for(int i=0; i<n; i++)// вывод результата
        {
            cout << "  x(" << i << ") = " << x[i] << endl;
        }
 
 
 
 
        for(int i=0; i<n; i++)// освобождение выдиленой памяти
        {
               delete [] a[i];
        }
        delete [] a;
        delete [] b;
        delete [] x;
        delete [] xp;
        cout << " Для продолжения программы, введите - 'y'" << endl;
        cout << " для завершения, введите любой символ..." << endl;
        cin >> con;
    }
    while(con=='y');
    return 0;
}

#include <iostream>
#include <conio.h>
#include <locale>
 
using namespace std;
 
int main()
{
    setlocale(LC_ALL, "rus");
    int n;
    double delta, eps;
 
    do
    {
        cout << " Введите количество неизвестных: " << endl;
        cin >> n;
        cout << " Введите точность: " << endl;
        cin >> eps;
 
        double **a = new double *[n];//выдeление динам. памяти для массива
        for(int i=0; i<n; i++)
        {
            a[i] = new double[n];
        }
        double *b = new double[n];// выдeление памяти для свободных членов
        double *x = new double[n];
        double *xp = new double[n];
 
        for(int i=0; i<n; i++)// ввод данных
        {
            for(int j=0; j<n; j++)
            {
                cout << " Введите " << (j+1) << "й коэффициент " << (i+1) << "го уравнения: ";
                cin >> a[i][j];
            }
            cout << " Введите свободный член " << (i+1) << "го уравнения: ";
            cin >> b[i];
        }
 
        for(int i=0; i<n; i++)// вывод СЛАУ
        {
            cout << " | ";
            for(int j=0; j<n; j++)
            {
                cout << a[i][j] << "*x(" << (i+1) << ") + " ;
            }
            cout << " = " << b[i];
            cout << endl;
        }
 
        for(int i=0; i<n; i++)
        {
            double temp = 0.0;
            for(int j=0; j<n; j++)
            {
                if(j != i) temp += a[i][j];
            }
            if(a[i][i] < temp)
            {
                cout << " Необходимое условие сходимости не выполняется" << endl;
                break;
            }
        }
 
        for(int i=0; i<n; i++)
        {
            xp[i] = b[i]/a[i][i]; // начальное значение Х
        }
 
        do
        {
            for(int i=0; i<n; i++)
            {
                double temp = 0.0;
                for(int j=0; j<n; j++)
                {
                    if(j != i) temp += a[i][j]*xp[j];
                }
                x[i] = (b[i] - temp)/a[i][i];
                delta = fabs(x[0] - xp[0]);
                for(int j=1; j<n; j++)
                {
                    if(fabs(x[j] - xp[j]) > delta)
                        delta = fabs(x[j] - xp[j]);
                }
                xp[i] = x[i];
            }
        }
        while(delta > eps);
 
        for(int i=0; i<n; i++)// вывод результата
        {
            cout << "  x(" << i << ") = " << x[i] << endl;
        }
 
        cout << " Для продолжения программы, введите 'y'" << endl;
        cout << " Для завершения введите любой другой символ..." << endl;
        int ch = getche();
        if((char)ch != 'y') break;
 
        for(int i=0; i<n; i++)// освобождение выдeленной памяти
        {
            delete [] a[i];
        }
        delete [] a;
        delete [] b;
        delete [] x;
        delete [] xp;
 
    } while(1);
 
 
    return 0;
}

  for(int i=0; i<n; i++)
        {
            double temp = 0.0;
            for(int j=0; j<n; j++)
            {
                if(j != i) temp += a[i][j];
            }
            if(a[i][i] < temp)
            {
                cout << " Необходимое условие сходимости не выполняется" << endl;
                break;
            }
        }
 
        for(int i=0; i<n; i++)
        {
            xp[i] = b[i]/a[i][i]; // начальное значение Х
        }
 
        do
        {
            for(int i=0; i<n; i++)
            {
                double temp = 0.0;
                for(int j=0; j<n; j++)
                {
                    if(j != i) temp += a[i][j]*xp[j];
                }
                x[i] = (b[i] - temp)/a[i][i];
                delta = fabs(x[0] - xp[0]);
                for(int j=1; j<n; j++)
                {
                    if(fabs(x[j] - xp[j]) > delta)
                        delta = fabs(x[j] - xp[j]);
                }
                xp[i] = x[i];
            }
        }
        while(delta > eps);

if(fabs(a[i][i]) < fabs(temp))
            {
                cout << " Необходимое условие сходимости не выполняется" << endl;
                break;
            }