6 double det(
double **aa,
int n)
9 double d, big, dum, sum, tem;
17 if ((tem=fabs(aa[i][j])) > big) big=tem;
18 if (big == 0.0)
nrerror (
"Singular matrix in routine dim");
24 for (k=1;k<i;k++) sum -= aa[i][k]*aa[k][j];
31 sum -= aa[i][k]*aa[k][j];
33 if ( (dum=vv[i]*fabs(sum)) >= big) {
47 if (aa[j][j] == 0.0) aa[j][j]=
TINY;
50 for (i=j+1;i<=n;i++) aa[i][j] *= dum;
double det(double **aa, int n)
void free_vector(double *v, int nl, int nh)
void nrerror(char *error_text)
double * vector(int nl, int nh)