void rank(unsigned long n, unsigned long indx[], unsigned long irank[])
{
	unsigned long j;

	for (j=1;j<=n;j++) irank[indx[j]]=j;
}
/* (C) Copr. 1986-92 Numerical Recipes Software $!6)D$|2$3. */