Алгоритм вычисления определителя квадратной матрицы

#include <stdio.h>
 
double x[50][50];
int rang;
double com=1;
 
//Деление всех элементов строки на указанный элемент этой же строки 
void delst(int stroka, int stolb)
{
int j;
double t=x[stroka][stolb];
 
        for (j=0;j<rang;j++)
            {
 
            x[stroka][j]=x[stroka][j]/t;
 
            }
 
        com=com*t;
}
//--------------------------------------------------------------------
 
//Умножение всех элементов указанной строки на элемент указанной строки , указанного столбца 
void multst(int stroka_out, int stroka_in, int stolb)
{
int j;
double t=x[stroka_in][stolb];
 
        for (j=0;j<rang;j++)
            {
 
            x[stroka_out][j]=x[stroka_out][j]*t;
 
            }
 
        com=com/t;
}
//--------------------------------------------------------------------
 
 
//Вычитание элементов строк ---------------
void minusst(int stroka_in, int stroka_out)
{
int j;
        for (j=0;j<rang;j++)
            {
 
            x[stroka_out][j]=x[stroka_out][j]-x[stroka_in][j];
 
            }
}
//--------------------------------------------------------------------
 
//Перестановка двух уравнений 
void perestanovka(int stroka_1, int stroka_2)
{
int j;
double str1[50], str2[50];
 
        for (j=0;j<rang;j++)
            {
 
            str1[j]=x[stroka_1][j];
            str2[j]=x[stroka_2][j];
 
            }
 
        for (j=0;j<=rang;j++)
            {
 
            x[stroka_1][j]=str2[j];
            x[stroka_2][j]=str1[j];
 
            }
}
//--------------------------------------------------------------------
 
//Выбор ведущего элемента ----------------------------------------------------------------------
int scan_nenol(int stroka, int stolb)
{
int i;
        for (i=stroka+1;i<rang;i++)
            {
 
            if (x[i][stolb]!=0)
               {
               return (i);
               break;
               }
            else
               {
               if (i==rang-1) return 0;
               }
 
            }
}
//--------------------------------------------------------------------
 
//Элементарные преобразования ----------------------------------------------------------------------
void nul()
{
int i,j,ss;
 
        for (j=0;j<rang-1;j++)
            {
 
            for (i=j+1;i<rang;i++)
                {
 
                if (x[i][j]==0)
                   {
                   continue;
                   }
                else
                   {//*******************************************************************
 
                   //----------------//
                   if (x[j][j]==0)
                      {
                      ss = scan_nenol(j,j);
                      if (ss==0)
                         {
                         continue;                   
                         }
                      else
                         {
                         perestanovka(j,ss);
                         }
                      }
                   //----------------//
                   
                   else
                       {
                       delst(j,j);
                       multst(j,i,j);
                       minusst(j,i);
                       }
 
                   }//********************************************************************
 
                }