934 Appendix A void gaucof(int n,float a[],float b[],float amuo,float x[],float w[]); void gauher(float x[],float w[],int n); void gaujac(float x[],float w[],int n,float alf,float bet); void gaulag(float x[],float w[],int n,float alf); void gauleg(float x1,float x2,float x[],float w[],int n); void gaussj(float *a,int n,float **b,int m); void gcf(float *gammcf,float a,float x,float *gln); float golden(float ax,float bx,float cx,float (*f)(float),float tol, float *xmin); void gser(float *gamser,float a,float x,float *gln); void hpsel(unsigned long m,unsigned long n,float arr[],float heap); void hpsort(unsigned long n,f1 oat ra☐)； 常 void hgr(f1oat**a,intn,f1 oat wr☐，f1 oat wi☐)； void hufapp(unsigned long index[],unsigned long nprob],unsigned long n, unsigned long i); void hufdec(unsigned long *ich,unsigned char *code,unsigned long lcode, 鱼 18881892 unsigned long *nb,huffcode *hcode); void hufenc(unsigned long ich,unsigned char **codep,unsigned long *lcode, 11800 nsig ed long *nb,huffcode *hcode) void hufmak(unsigned long nfreg[],unsigned long nchin,unsigned long *ilong, 3 unsigned long *nlong,huffcode *hcode); from NUMERICAL RECIPES I void hunt(float xx[],unsigned long n,float x,unsigned long *jlo); void hypdrv(f1oats,f1 oat yy☐，f1 oat dyyds▣)； fcomplex hypgeo(fcomplex a,fcomplex b,fcomplex c,fcomplex z); (Nor to make void hypser(fcomplex a,fcomplex b,fcomplex c,fcomplex z, fcomplex *series,fcomplex *deriv); America server computer, one paper Cambridge University Press. THE unsigned short icrc(unsigned short crc,unsigned char *bufptr, unsigned long len,short jinit,int jrev); ART unsigned short icrcl(unsigned short crc, unsigned char onech); unsigned long igray(unsigned long n,int is); void iindexx(unsigned long n,long arr[],unsigned long indx[]); 9 void indexx(unsigned long n,float arr[],unsigned long indx[]); strictly prohibited. for their Programs void interp(double **uf,double **uc,int nf); int irbitl(unsigned long *iseed); int irbit2(unsigned long *iseed); void jacobi(float **a,int n,float d],float **v,int *nrot); to dir void jacobn(float x,float y[],float dfdx[],float **dfdy,int n); long julday(int mm,int id,int iyyy); void kendl1(float datal,float data2[],unsigned long n,float *tau,float *z, 18881920 OF SCIENTIFIC COMPUTING(ISBN float *prob); void kendl2(float **tab,int i,int j,float *tau,float *z,float *prob); void kermom(double w☐，double y,intm); void ks2dis(float xi,float y1[],unsigned long ni, void (*quadvl)(float,float,float *float *float *float *) float *d1,float *prob); void ks2d2s(f1oatx10,f1oaty1[],nsigned1ongn1,f1oatx2☐，f1oaty2[☐， idge.org unsigned long n2,float *d,float *prob); Numerical Recipes -43106 void ksone(float data[],unsigned long n,float (*func)(float),float *d, float *prob); (outside void kstwo(float datal[],unsigned long n1,float data2[],unsigned long n2, float *d,float *prob); North Software. void laguer(fcomplex a],int m,fcomplex *x,int *its); void lfit(float x[],float y[],float sig[],int ndat,float a],int ia], int ma,float **covar,float *chisq,void (*funcs)(float,float []int)); void linbcg(unsigned long n,double b[],double x[],int itol,double tol, int itmax,int *iter,double err); void linmin(float p[],float xi[],int n,float *fret, f1oat(*func)(f1oat[☐))； void Insrch(int n,float xold[],float fold,float g],float p[],float x[], float *f,float stpmax,int *check,float (*func)(float [])) void load(f1oatx1,f1oatv[☐，f1oaty0); void1oad1(f1oatx1,f1oatv1▣，f1oaty[]); void load2(float x2,float v2,float y[]); void locate(float xx[],unsigned long n,float x,unsigned long *j);934 Appendix A Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copyin Copyright (C) 1988-1992 by Cambridge University Press. Programs Copyright (C) 1988-1992 by Numerical Recipes Software. Sample page from NUMERICAL RECIPES IN C: THE ART OF SCIENTIFIC COMPUTING (ISBN 0-521-43108-5) g of machine￾readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website http://www.nr.com or call 1-800-872-7423 (North America only), or send email to directcustserv@cambridge.org (outside North America). void gaucof(int n, float a[], float b[], float amu0, float x[], float w[]); void gauher(float x[], float w[], int n); void gaujac(float x[], float w[], int n, float alf, float bet); void gaulag(float x[], float w[], int n, float alf); void gauleg(float x1, float x2, float x[], float w[], int n); void gaussj(float **a, int n, float **b, int m); void gcf(float *gammcf, float a, float x, float *gln); float golden(float ax, float bx, float cx, float (*f)(float), float tol, float *xmin); void gser(float *gamser, float a, float x, float *gln); void hpsel(unsigned long m, unsigned long n, float arr[], float heap[]); void hpsort(unsigned long n, float ra[]); void hqr(float **a, int n, float wr[], float wi[]); void hufapp(unsigned long index[], unsigned long nprob[], unsigned long n, unsigned long i); void hufdec(unsigned long *ich, unsigned char *code, unsigned long lcode, unsigned long *nb, huffcode *hcode); void hufenc(unsigned long ich, unsigned char **codep, unsigned long *lcode, unsigned long *nb, huffcode *hcode); void hufmak(unsigned long nfreq[], unsigned long nchin, unsigned long *ilong, unsigned long *nlong, huffcode *hcode); void hunt(float xx[], unsigned long n, float x, unsigned long *jlo); void hypdrv(float s, float yy[], float dyyds[]); fcomplex hypgeo(fcomplex a, fcomplex b, fcomplex c, fcomplex z); void hypser(fcomplex a, fcomplex b, fcomplex c, fcomplex z, fcomplex *series, fcomplex *deriv); unsigned short icrc(unsigned short crc, unsigned char *bufptr, unsigned long len, short jinit, int jrev); unsigned short icrc1(unsigned short crc, unsigned char onech); unsigned long igray(unsigned long n, int is); void iindexx(unsigned long n, long arr[], unsigned long indx[]); void indexx(unsigned long n, float arr[], unsigned long indx[]); void interp(double **uf, double **uc, int nf); int irbit1(unsigned long *iseed); int irbit2(unsigned long *iseed); void jacobi(float **a, int n, float d[], float **v, int *nrot); void jacobn(float x, float y[], float dfdx[], float **dfdy, int n); long julday(int mm, int id, int iyyy); void kendl1(float data1[], float data2[], unsigned long n, float *tau, float *z, float *prob); void kendl2(float **tab, int i, int j, float *tau, float *z, float *prob); void kermom(double w[], double y, int m); void ks2d1s(float x1[], float y1[], unsigned long n1, void (*quadvl)(float, float, float *, float *, float *, float *), float *d1, float *prob); void ks2d2s(float x1[], float y1[], unsigned long n1, float x2[], float y2[], unsigned long n2, float *d, float *prob); void ksone(float data[], unsigned long n, float (*func)(float), float *d, float *prob); void kstwo(float data1[], unsigned long n1, float data2[], unsigned long n2, float *d, float *prob); void laguer(fcomplex a[], int m, fcomplex *x, int *its); void lfit(float x[], float y[], float sig[], int ndat, float a[], int ia[], int ma, float **covar, float *chisq, void (*funcs)(float, float [], int)); void linbcg(unsigned long n, double b[], double x[], int itol, double tol, int itmax, int *iter, double *err); void linmin(float p[], float xi[], int n, float *fret, float (*func)(float [])); void lnsrch(int n, float xold[], float fold, float g[], float p[], float x[], float *f, float stpmax, int *check, float (*func)(float [])); void load(float x1, float v[], float y[]); void load1(float x1, float v1[], float y[]); void load2(float x2, float v2[], float y[]); void locate(float xx[], unsigned long n, float x, unsigned long *j);