3
#define SWAP(g,h) {y=(g);(g)=(h);(h)=y;}
5
hess_elim(double **a, int n)
10
for(m=1; m < n-1 ; m++) {
13
for(j=m; j < n ; j++) {
14
if(fabs(a[j][m-1]) > fabs(x)) {
20
for(j=m-1; j < n ; j++) SWAP(a[i][j],a[m][j])
21
for(j=0; j < n ; j++) SWAP(a[j][i],a[j][m])
24
for(i=m+1; i < n ; i++) {
28
for(j=m; j < n ; j++) a[i][j] -= y*a[m][j];
29
for(j=0; j < n; j++) a[j][m] += y*a[j][i];