c语言实现求行列式值

double hanglieshi(double *h,int i,int j)
{
double *A = (double *)malloc(sizeof(double)*i*j);
int k,l,s;
double t = 1.0;
int count=0;//因为可能出现除数为0的情况，所以需要换行
for(k = 0;k<i;k++)
{
for(l = 0;l<j;l++)
{
*(A + k*j +l) = *(h + k*j +l);
}
}
if(i!=j)
{
printf("Error");
return 0;
}
if(i==1&&j==1)
{
return *(A);
}
    else
{
   for(k = 1;k<i;k++)
{
for(s = k;s < i;s++)
{
int tag =0;
if(*(A+(k-1)*j+k-1) == 0)
{
for(l=k;l<i;l++)
{
if(*(A+l*j)!=0)
break;
}
if(l == i)
{
return 0;
}
tag =l;
for(l = k-1;l<j;l++)
{
double B = 0.0;
B = *(A+(k-1)*j+l);
*(A+(k-1)*j+l) = *(A+tag*j+l);
*(A+tag*j+l) = B;
count++;
}
}
double B = *(A + s*j+k-1);
for(l = k-1;l <j;l++)
       {
      *(A + s*j+l) = *(A + s*j+l)-B*(*(A+(k-1)*j+l))/(*(A+(k-1)*j+k-1));
       }
}
   
}
}
/*for(k = 0;k < i;k++)
{
for(l = 0;l < j; l++)
{
printf("%lf   ",*(A+k*j+l));
}
printf("\n");
}*/
for(k = 0,l = 0;k<i&&l<j;k++,l++)
{
t = t*(*(A+k*j+l));
}
free(A);
if(count%2 == 0)
   return t;
else
return -t;
}