矩阵计算模板

for (i=0,j=0;i<N&&j<N;i++,j++)
{
    if (a[i][j]==0)
    {
        for (m=i;a[m][j]==0;m++);
        /* if (m==N)
        {
        sn=0;
        printf("detA=%lf\n",sn);
        exit(0);
        }
        else*/
        for (n=j;n<N;n++)
        {
             c=a[i][n];
             a[i][n]=a[m][n];
             a[m][n]=c;
        }
        k*=(-1);
    }
    for (s=N-1;s>i;s--)
    {
        x=a[s][j];
        for (t=j;t<N;t++)
        a[s][t]-=a[i][t]*(x/a[i][j]);
    }
}
for (i=0;i<N;i++)
f*=a[i][i];
sn=k*f;
printf ("%.lf\n",sn);