Метод Гаусса

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
float A[100][100];
float B [100];
int n, m;
 
void readmatr(){
    int i, j;
    for (i=1; i<=n; i++){
        for (j=1; j<=m; j++){
            scanf("%f", &A[i][j]);
        }
    }
}//функция считывает матрицу неизвестных
void readsol(){
    int i;
    for (i=1; i<=n; i++){
        scanf("%f", &B[i]);
    }
}//функция считывает столбец свободных членов
 
void writesys(){
    int i, j;
    for (i=1; i<=n; i++){
        printf("%.3f*X[%d]", A[i][1], 1);
        for (j=2; j<=m; j++){
             printf("%+.3f*X[%d]", A[i][j], j); 
        } 
        printf(" = %.3f\n", B[i]);
    } 
}//функция выводит систему на экран
 
void write(){
     int i, j;
     for (i=1; i<=n; i++){
          for (j=1; j<=m; j++){
              printf("%.3f ", A[i][j]);
          } 
          printf("| %.3f\n", B[i]);
      }
}
 
void liniesperm(int a, int b){
    int i;
    float t;
    for(i=1; i<=m; i++){
              t=A[a][i];
              A[a][i]=A[b][i];
              A[b][i]=t;
    }    
}//функция для перестановки двух строк местами
 
void permutation (int l, int c){
    int strok, stolb;
    float t;
    strok=l+1;
    stolb=c;
    while (A[l][c] == 0){
        if(strok > n){
            stolb++;
            strok=l;
        }
        if(stolb > n){
            break;
        }
        if (A[strok][stolb] != 0){
           liniesperm(strok, l);
           t = B[l];
           B[l]=B[strok];
           B[strok]=t;
       }
       if (A[l][c] == 0){
           strok++;
       }
      if (strok > n){break;}
    }
} //функция сортировки, обеспечивает что бы на диагонали стояли не нулевые элементы
 
 
void nulling (int k){
       int i, j;
       if (A[k][k] != 0){
          B[k] /= A[k][k];
          for ( j = n; j >= k; j-- ){
              A[k] [j] /= A[k] [k];
          }
          for ( i = k + 1; i <= n; i++ ){
               B[i]-=B[k]*A[i][k];
               for ( j = n; j >= k; j-- ){
                    A[i] [j] -= A[k] [j] * A[i] [k];
               }
          }
          for ( i = k - 1; i >= 1; i-- ){
               B[i]-=B[k]*A[i][k];
               for ( j = n; j >= k; j-- ){
                   A[i] [j] -= A[k] [j] * A[i] [k];
               }
          }
       }
}//функция для получения нулей над и под главной диагональю
 
int compatibility (){
    int i, j, flag=0;
    float t;
    for (i=1; i<=n; i++){
        t=0;
        for (j=1; j<=m; j++){
            if (A[i][j] !=0){
                t=A[i][j];
            }
        }
        if ((t==0) && (A[i] != 0)){
            flag = 1;
        }  
    }
    return(flag);
}//функция определяет совместность системы
 
 
int main() {
    int i, j, flag;
    begin: ;
    printf("\nVvedite kolichestvo uravneniy: ");
    scanf("%d", &n);
    printf("\nVvidite kolichestvo neizvestnih: ");
    scanf("%d", &m);
    printf("\nVvedite metricu neizvestnih:\n");
    readmatr();                                    // вводим систему в матричном виде
    printf("\nVvidine stolbec svobodnih chlenov:\n");
    readsol();
    printf("\nPoluchenaya sistema:\n");//выводим систему в стандартном виде
    writesys();
    printf("\nEsli sistema verna vvelite 0, esli net vvedite 1:");
    scanf("%d", &flag);
    if (flag == 1){
        system ("cls");
        goto begin;
    }//проверяем верно ли введена система, в случае ошибки в системе переходим к началу
    printf("\nProizvodimie preobrazovania:");
    for (i=1; i<=n; i++){
        permutation(i, i);
        nulling (i);
        printf("\nShag #%d:\n", i);//выводим преобразования на каждом шаге(для удобства проверки)
        write();  
    }//приводим к треугольному виду
    printf("\n\nReshenia:\n\n");
    if (compatibility() == 1){
        printf("\nSistena nesovmestna\n");
    }// проверяем систему на совместность
    else{
      if (n >= m){//проверяем наличие свободных переменных
          for (i=1; i<=n; i++){
              printf("x[%d] = %.3f\n", i, B[i]);
          }//в случае отсутствия свободных выводим главные
      } 
      else {
          for (i=1; i<=n; i++){
              printf("x[%d] = %.3f", i, B[i]);
              for (j=n+1; j<=m; j++){
                  printf(" %+.3fx[%d]", -A[i][j], j);
              }
              printf("\n");
          }//выводим главные переменные выраженные через свободные
      }
    }
    return 0;
}