NyanCoder

2.c

Apr 4th, 2022 (edited)
413
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C 1.29 KB | None | 0 0
  1. #include "stdlib.h"
  2. #include "stdio.h"
  3.  
  4. int main(int, char**)
  5. {
  6.     system("clear");
  7.     int n = 0;
  8.     printf("Введите размер матриц\n");
  9.     scanf("%i", &n);
  10.     float** a = (float**)malloc(sizeof(void*) * n);
  11.     printf("Введите матрицу A\n");
  12.     for (int i = 0; i < n; i++)
  13.     {
  14.         a[i] = (float*)malloc(sizeof(float) * n);
  15.         for (int j = 0; j < n; j++)
  16.             scanf("%f", &a[i][j]);
  17.     }
  18.     printf("Введите матрицу B\n");
  19.     float* b = (float*)malloc(sizeof(float) * n);
  20.     for (int i = 0; i < n; i++)
  21.         scanf("%f", &b[i]);
  22.     for (int k = 0; k < n; k++) /* Прямой ход Гаусса */
  23.     {
  24.         for (int j = k + 1; j < n; j++)
  25.         {
  26.             const float r = a[j][k] / a[k][k];
  27.             for (int i = k; i < n; i++)
  28.                 a[j][i] = a[j][i] - r * a[k][i];
  29.             b[j] = b[j] - r * b[k];
  30.         }
  31.     }
  32.     for (int k = n - 1; k > -1; k--) /* Обратный ход Гаусса */
  33.     {
  34.         float r = 0;
  35.         for (int j = k + 1; j < n; j++)
  36.             r += a[k][j] * b[j];
  37.         b[k] = (b[k] - r) / a[k][k];
  38.     }
  39.     printf("Корни системы:\n");
  40.     for (int i = 0; i < n; i++)
  41.         printf("x[%i]=%f   ", i + 1, b[i]);
  42.     return 0;
  43. }
Add Comment
Please, Sign In to add comment