Download this file
1 2 3 4 5 6 7 8 9 10 11
extern long idum; void ranpt(float pt[], float regn[], int n) { float ran1(long *idum); int j; for (j=1;j<=n;j++) pt[j]=regn[j]+(regn[n+j]-regn[j])*ran1(&idum); } /* (C) Copr. 1986-92 Numerical Recipes Software )&%. */