Решть систему уравнений - C++ - Ответ 1578472

#include <stdio.h>
#include <stdlib.h>
#include <iostream.h>
#include <conio.h>
 
// Функция обращения матрицы
// A - матрица
// N - размерность матрицы
void inversion(double **A, int N)
{
    double **E;  // Единичная матрица
    double temp; // Временный элемент
    int i, j, k; // Счётчики
 
    // Выделяем память под массив указателей на строки матрицы
    E = (double **)malloc(N * sizeof(double *));
    
    // Выделяем память под каждую строку матрицы
    for (i = 0; i < N; i++)
        E[i] = (double *)malloc(N * sizeof(double));
    
    // Заполняем единичную матрицу
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < N; j++)
        {
            E[i][j] = 0.0;
 
            if (i == j)
                E[i][j] = 1.0;
        }
    }
    
    // Проход по матрице вниз (полчаем нижнюю треугольную матрицу)
    for (k = 0; k < N; k++)
    {
        temp = A[k][k]; // Запоминаем очередной диагональный элемент
 
        // Делим каждый элемент очередной строки на запомненный
        // ранее диагональный элемент (как в исходной матрице,
        // так и в единичной)
        for (j = 0; j < N; j++)
        {
            A[k][j] /= temp;
            E[k][j] /= temp;
        }
        
        // Проходим по всем строкам ниже данной
        for (i = k + 1; i < N; i++)
        {
            // Запоминаем элемент соответствующей строки,
            // находящийся под диагональным
            temp = A[i][k];
            
            // Производим над каждым элементом соответствующей
            // строки необходимые операции (подробнее см. 
            // метод Гаусса-Жордана)
            for (j = 0; j < N; j++)
            {
                A[i][j] -= A[k][j] * temp;
                E[i][j] -= E[k][j] * temp;
            }
        }
    }
 
    // Проход по матрице вверх (получаем единичную матрицу на месте
    // исходной и обратную на месте единичной)
    for (k = N - 1; k > 0; k--)
    {
        // Проходим по всем строкам выше данной
        for (i = k - 1; i >= 0; i--)
        {
            // Запоминаем элемент соответствующей строки,
            // находящийся над диагональным
            temp = A[i][k];
 
            // Производим над каждым элементом соответствующей
            // строки необходимые операции
            for (j = 0; j < N; j++)
            {
                A[i][j] -= A[k][j] * temp;
                E[i][j] -= E[k][j] * temp;
            }
        }
    }
    
    // Переписываем значения из единичной матрицы в исходную
    // (поскольку в данный момент именно в единичной матрице
    // хранится обратная)
    for (i = 0; i < N; i++)
        for (j = 0; j < N; j++)
            A[i][j] = E[i][j];
 
    // Очищаем память из-под единичной матрицы
    for (i = 0; i < N; i++)
        free(E[i]);
 
    free(E);
}
 
int main()
{
    double **matrix,s; // Матрица
    double **c;
    double *x,*b;
    int N;           // Размерность
    int i, j;        // Счётчики
 
    // Вводим размерность
    printf("Enter N: ");
    scanf("%d", &N);
 
    // Выделяем память
    matrix = (double **)malloc(N * sizeof(double *));
    с== (double **)malloc(N * sizeof(double *));
    for (i = 0; i < N; i++)
        matrix[i] = (double *)malloc(N * sizeof(double));
    
    // Вводим матрицу
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < N; j++)
        {
            printf("matrix[%d][%d] = ", i, j);
            scanf("%lf", &matrix[i][j]);
        }
    }
    cout<<"vvedite vektor b";
     for (i=0;i<N;i++)
     cin>>b[i];
    
    // Обращаем матрицу
    inversion(matrix, N);
 
 
   //ищем х через умножение обратной матр на вектор  
 { for(i=0;i<N;i++)
  {s=0;
   for(int j=0;j<N;j++)
    {s+=matrix[i][j]*b[j];}
   x[i]=s;}}
 for(i=0;i<N;i++) cout<<"\nx="<<x[i];
 
   // умножаем х на х транспон.
 for(i=0;i<N;i++)
 {cout<<"\n";
  for(int j=0;j<N;j++)
 
   {
    c[i][j]=x[i]*x[j];
   }
   }
 //H=E-(X*Xt).цикл для вычитания матрицы Е. 
  for(i=0;i<N;i++) 
   for(j=0;j<N;j++)
    {c[i][j]*=-1;
      if (i==j)
       {c[i][j]+=1;
       }
     }
 
    // Выводим  матрицу
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < N; j++)
            printf("%lf  ", matrix[i][j]);
 
        printf("\n");
    }
 
    // Освобождаем память
    for (i = 0; i < N; i++)
        free(matrix[i]);
 
    free(matrix);
 
    return 0;
}