text: /* Driver for routine eigsrt */
  
  #include <stdio.h>
  #define NRANSI
  #include "nr.h"
  #include "nrutil.h"
  
  #define NP 10
  
  int main(void)
  {
  	int i,j,nrot;
  	static float c[NP][NP]=
  		{5.0,4.3,3.0,2.0,1.0,0.0,-1.0,-2.0,-3.0,-4.0,
  		4.3,5.1,4.0,3.0,2.0,1.0,0.0,-1.0,-2.0,-3.0,
  		3.0,4.0,5.0,4.0,3.0,2.0,1.0,0.0,-1.0,-2.0,
  		2.0,3.0,4.0,5.0,4.0,3.0,2.0,1.0,0.0,-1.0,
  		1.0,2.0,3.0,4.0,5.0,4.0,3.0,2.0,1.0,0.0,
  		0.0,1.0,2.0,3.0,4.0,5.0,4.0,3.0,2.0,1.0,
  		-1.0,0.0,1.0,2.0,3.0,4.0,5.0,4.0,3.0,2.0,
  		-2.0,-1.0,0.0,1.0,2.0,3.0,4.0,5.0,4.0,3.0,
  		-3.0,-2.0,-1.0,0.0,1.0,2.0,3.0,4.0,5.0,4.0,
  		-4.0,-3.0,-2.0,-1.0,0.0,1.0,2.0,3.0,4.0,5.0};
  	float *d,**v,**e;
  
  	d=vector(1,NP);
  	v=matrix(1,NP,1,NP);
  	e=convert_matrix(&c[0][0],1,NP,1,NP);
  	printf("****** Finding Eigenvectors ******\n");
  	jacobi(e,NP,d,v,&nrot);
  	printf("unsorted eigenvectors:\n");
  	for (i=1;i<=NP;i++) {
  		printf("eigenvalue %3d = %12.6f\n",i,d[i]);
  		printf("eigenvector:\n");
  		for (j=1;j<=NP;j++) {
  			printf("%12.6f",v[j][i]);
  			if ((j % 5) == 0) printf("\n");
  		}
  		printf("\n");
  	}
  	printf("\n****** Sorting Eigenvectors ******\n\n");
  	eigsrt(d,v,NP);
  	printf("sorted eigenvectors:\n");
  	for (i=1;i<=NP;i++) {
  		printf("eigenvalue %3d = %12.6f\n",i,d[i]);
  		printf("eigenvector:\n");
  		for (j=1;j<=NP;j++) {
  			printf("%12.6f",v[j][i]);
  			if ((j % 5) == 0) printf("\n");
  		}
  		printf("\n");
  	}
  	free_convert_matrix(e,1,NP,1,NP);
  	free_matrix(v,1,NP,1,NP);
  	free_vector(d,1,NP);
  	return 0;
  }
  #undef NRANSI
  /* (C) Copr. 1986-92 Numerical Recipes Software 7&X*. */
  