Написать программу которая находит решение СЛАУ методом Гаусса - C (СИ) - Ответ 3499172

#include<stdio.h>
#include<conio.h>
#include<math.h>
#define N 5
#define Q -2.23
#define D 2
 
int main()
{
  double a[5][5]={Q, 1, 0, 0, 0,     //исходная матрица
                   1, -2, 1, 0, 0,
                   0, 1, -2, 1, 0,
                   0, 0, 1, -2, 1,
                   0, 0, 0, 1, Q};
 
   double b[5]={0, D, D, D, 0};    //вектор свободных коофициэнтов
   double x[5]={0};
   double k, m, p, j, i;  //счётчики
   double r, c, s;
 
   for(k=0;k<N-1;k++)
      {
       p=k;
       //printf("%f\n\n", p);
          //printf("k=%f\n\n", k);
       for(m=k+1;m<N;m++)
          {
           //printf("%f and %f\n\n", a[p][k], a[m][k]);
         // printf("m=%f\n\n",m );
           if(fabs(a[p][k])<fabs(a[m][k]))
              {
               p=m;
               //printf("%f and %f\n\n", a[p][k], a[m][k]);
              }
          }
 
       for(j=k;j<N;j++)
          {
           //printf("j=%f\n\n", j);
           r=a[k][j];
           //printf("%f\n\n", a[k][j]);
           a[k][j]=a[p][j];
           //printf("%f\n\n", a[p][j]);
           a[p][j]=r;
          }
 
      /*printf("\n\n\n\n");
      for(i=0;i<n;i++)
         {
         printf("\n");
         for(j=0;j<n;j++)
            {
             printf("%10f", a[i][j]);
            }
          }*/
 
 
       r=b[k];
       b[k]=b[p];
       b[p]=r;
      /* printf("\n\n\n");
       for(i=0;i<n;i++)
          printf("%10f", b[i]); */
 
       for(m=k+1;m<N;m++)
          {
           c=a[m][k]/a[k][k];
          // printf("c=%10f\n", c);
           b[m]=b[m]-c*b[k];
 
           for(i=k;i<N;i++)
              {
               a[m][i]-=c*a[k][i];
              }
          }
        }
      x[N-1]=b[N-1]/a[N-1][m-1];
 
      for(k=N-2;k!=0;k--)
         {
          s=0;
 
          for(i=k-1;i<N-1;i++)
             {
              s=s+a[k][i]*x[i];
             }
 
          x[k]=(b[k]-s)/a[k][k];
         }
    //printf("s=%f", s);
    printf("\n\n\n");
   for(i=0;i<N;i++)
      printf("%f  ", x[i]);
 
   getch();
   return 0;
}