обратная матрица методом гауса

//Пусть матрица А имеет размер NxN.
//Заведем матрицу а размера Nx2N для хранения (А|E).
//Последовательность действий должна быть такая:
    // Спросить у пользователя А.
     //Дописать справа Е.
     //Применить метод Гаусса.
     //Напечатать ответ.
#include<stdio.h>
#include<conio.h>
#include<math.h>
 
#define N 4
double a[N][2*N];
 
void gauss(void);
void sprositMatr(void);
void dopisatE(void);
void otvet(void);
void delimStr(int s,double lam);
void zanulimStlb(int ns);
void vichestStr(int s1,int s2,double lam);
void nomerStrSNenulElem(void);
void meniaemStr(int s1,int s2);
 
void main(void){
  //clrscr();
  sprositMatr();
  dopisatE();
  gauss();
  otvet();
  getch();
  //delimStr();
  //zanulimStlb();
  //vichestStr();
  nomerStrSNenulElem();
  //meniaemStr();
}
void sprositMatr(void){   //Спросить у пользователя А
  int i,j;
  for(i=0;i<N;i++){
      for(j=0;j<N;j++){
         printf("A%d%d=",i,j);
         scanf("%lf",&a[j]);
      }
  }
}
 
void otvet(void){         //Напечатать ответ
  int i,j;
  for(i=0;i<N;i++){
      for(j=N;j<2*N;j++)
         printf("%6.2f",a[j]);
       printf("\n");
  }
}
//Так как мы объявили матрицу а глобально, она автоматически заполнилась нулями.
//Поэтому, чтобы дописать единичную матрицу к нашей, достаточно поставить единицы на диагонали:
void dopisatE(void){
  int i;
  for (i=0;i<N;i++)
    a[N+i]=1;
}
//Реализуем функцию gauss.
//Сделаем слева диагональную матрицу.
//Сначала, пользуясь нулевой строкой матрицы,
//добьемся того, чтобы все элементы нулевого столбика,
//кроме а00, стали равны 0.
//Далее, с помощью первой строки зануляем все элементы первого столбика, кроме а11.
//Продолжая эту процедуру дальше, мы сделаем слева диагональную матрицу.
//После этого остается только разделить нулевую строку на  а00, первую на а11 и т.д. Программируем:
void gauss(void){
  int i;
  for(i=0;i<N;i++)
    zanulimStlb(i);
  for(i=0;i<N;i++)
    delimStr(i,a);
}
 
//Функция delimStr(s,lam) делит строку номер s на число lam.
//Реализуем ее:
void delimStr(int s,double lam){
 
  int j;
 
  for(j=0;j<2*N;j++)
    a[j]/=lam;
}
 
 
//Функция zanulimStlb(ns), пользуясь строкой номер ns, делает нулями все элементы столбца ns, кроме ans,ns . Чтобы этого добиться, нужно проделать следующие действия:
//(0)-(ns)*a0,ns/ ans,ns
//(1)-(ns)*a1,ns/ ans,ns
//…
//(N-1)-(ns)*aN-1,ns/ ans,ns
//Необходимо только при вычитании не забыть, что из самой себя строку ns вычитать не надо. Программируем:
//Если аns,ns=0, то нужно найти ненулевой элемент из чисел аns+1,ns , аns+2,ns, …, аN-1,ns, а затем строку с этим элементом поменять местами со строкой ns.
//Если такого элемента не окажется, то A-1 не существует.
void zanulimStlb(int ns){
  double lam;
  int nst;
  if( a[ns][ns]==0 ){
     //nst=nomerStrSNenulElem(ns);
      if(nst==N){
           printf("Обратной матрицы нет");
           getch();
           exit(1);
       }
        meniaemStr(ns,nst);
  }
  int i;
  for(i=0;i<N;i++){
    if(i!=ns){
      lam=a[ns]/a[ns][ns];
      vichestStr(i,ns,lam);
    }
  }
}
 
//Функция nomerStrSNenulElem(ns) ищет ненулевой элемент из чисел аns+1,ns , аns+2,ns, …, аN-1,ns и возвращает номер строки с этим элементом.
//Если элемента не оказалось, она возвращает N:
void nomerStrSNenulElem(void){
   int ns;
   int i;
   for(i=ns+1;i<N && fabs(a[ns])<1E-12;i++)
   return(i);
}
 
//Функция meniaemStr(s1,s2) меняет местами строки с номерами s1 и s2:
void meniaemStr(int s1,int s2){
 
 
    double temp;
    int j;
    for(j=0;j<2*N;j++){
        temp=a[s1][j];
        a[s1][j]=a[s2][j];
        a[s2][j]=temp;
    }
}
 
//функция vichestStr(s1,s2,lam), которая из строки s1 вычитает строку s2, умноженную на число lam:
void vichestStr(int s1,int s2,double lam){
  int j;
  for(j=0;j<2*N;j++)
    a[s1][j]-=a[s2][j]*lam;
}