936 Appendix A void pcshft(float a,float b,float d[],int n); void pearsn(float x[],float y[],unsigned long n,float *r,float *prob, float *z); void period(float x[],float y[],int n,float ofac,float hifac, float px[],float py[],int np,int *nout,int +jmax,float *prob); void piksr2(intn,f1 oat arr☐，f1 oat brr☐)； void piksrt(int n,float arr[]); void pinvs(int iel,int ie2,int jel,int jsf,int jc1,int k, float ***c,float **s); float plgndr(int 1,int m,float x); float poidev(float xm,long *idum); void polcoe(float x[],float y[],int n,float cof[); http://www.nr. void polcof(float xa☐，f1 oat ya],intn,float cof[☐)； void poldiv(float u[],int n,float v[],int nv,float q[],float r[]); void polin2(float xia[],float x2a[],float **ya,int m,int n, 83g float x1,float x2,float *y,float *dy); granted for 19881992 void polint(float xa],float ya[],int n,float x,float *y,float *dy); void powell(float p[],float **xi,int n,float ftol,int *iter,float *fret, 1-600 f1oat(*func)(f1oat☐))； void predic(f1 oat data☐，int ndata,f1oatd[☐，intm,float future[☐，int nfut); float probks(float alam); void psdes(unsigned long *lword,unsigned long *irword); /Cambridge University Press. from NUMERICAL RECIPES IN void pwt(float a[],unsigned long n,int isign); void pwtset(int n); float pythag(float a,float b); void pzextr(int iest,float xest,float yest[],float yz[],float dy], int nv); (North America server computer, to make float qgaus(float (*func)(float),float a,float b); one paper THE ART void grdcmp(float **a,int n,float *c,float *d,int *sing); float gromb(float (*func)(float),float a,float b); float gromo(float (*func)(float),float a,float b, strictly prohibited. Programs float (*choose)(float (*)(float),float,float,int)); void groot(float p[],int n,float *b,float *c,float eps); void grsolv(float **a,int n,float c[],float d[],float b[]); void grupdt(float **r,f1oat**gt,intn,f1oatu[门，f1oatv门)； float qsimp(float (*func)(float),float a,float b); float qtrap(float (*func)(float),float a,float b); to dir float quad3d(float (*func)(float,float,float),float x1,float x2); void quadct(float x, float y,float xx[],float yy],unsigned long nn OF SCIENTIFIC COMPUTING(ISBN float *fa,float *fb,float *fc,float *fd); rectcustser void quadmx(float **a,int n); void quadvl(float x,float y,float *fa,float *fb,float *fc,float *fd); float ran0(long *idum); 10-621 float ran1(long *idum); float ran2(long *idum); float ran3(long *idum); idge.org 1988-1992 by Numerical Recipes -43108 float ran4(long *idum); void rank(unsigned long n,unsigned long indx[],unsigned long irank[]); void ranpt(float pt[],float regn[],int n); (outside void ratint(float xa[],float ya[],int n,float x,float *y,float *dy); void ratlsq(double (*fn)(double),double a,double b,int mm,int kk, North Software. double cof],double *dev); double ratval(double x,double cof[],int mm,int kk); Amer float rc(float x,float y); float rd(float x,float y,float z); void realft(float data[],unsigned long n,int isign); void rebin(float rc,int nd,float r[],float xin[],float xi[]); void red(int iz1,int iz2,int jz1,int jz2,int jm1,int jm2,int jmf int ic1,int jci,int jcf,int kc,float ***c,float **s); void relax(double *u,double **rhs,int n); void relax2(double *u,double **rhs,int n); void resid(double **res,double **u,double **rhs,int n); float revcst(float x[,float y[],int iorder[],int ncity,int n[]); void reverse(int iorder[],int ncity,int n[]);936 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 pcshft(float a, float b, float d[], int n); void pearsn(float x[], float y[], unsigned long n, float *r, float *prob, float *z); void period(float x[], float y[], int n, float ofac, float hifac, float px[], float py[], int np, int *nout, int *jmax, float *prob); void piksr2(int n, float arr[], float brr[]); void piksrt(int n, float arr[]); void pinvs(int ie1, int ie2, int je1, int jsf, int jc1, int k, float ***c, float **s); float plgndr(int l, int m, float x); float poidev(float xm, long *idum); void polcoe(float x[], float y[], int n, float cof[]); void polcof(float xa[], float ya[], int n, float cof[]); void poldiv(float u[], int n, float v[], int nv, float q[], float r[]); void polin2(float x1a[], float x2a[], float **ya, int m, int n, float x1, float x2, float *y, float *dy); void polint(float xa[], float ya[], int n, float x, float *y, float *dy); void powell(float p[], float **xi, int n, float ftol, int *iter, float *fret, float (*func)(float [])); void predic(float data[], int ndata, float d[], int m, float future[], int nfut); float probks(float alam); void psdes(unsigned long *lword, unsigned long *irword); void pwt(float a[], unsigned long n, int isign); void pwtset(int n); float pythag(float a, float b); void pzextr(int iest, float xest, float yest[], float yz[], float dy[], int nv); float qgaus(float (*func)(float), float a, float b); void qrdcmp(float **a, int n, float *c, float *d, int *sing); float qromb(float (*func)(float), float a, float b); float qromo(float (*func)(float), float a, float b, float (*choose)(float (*)(float), float, float, int)); void qroot(float p[], int n, float *b, float *c, float eps); void qrsolv(float **a, int n, float c[], float d[], float b[]); void qrupdt(float **r, float **qt, int n, float u[], float v[]); float qsimp(float (*func)(float), float a, float b); float qtrap(float (*func)(float), float a, float b); float quad3d(float (*func)(float, float, float), float x1, float x2); void quadct(float x, float y, float xx[], float yy[], unsigned long nn, float *fa, float *fb, float *fc, float *fd); void quadmx(float **a, int n); void quadvl(float x, float y, float *fa, float *fb, float *fc, float *fd); float ran0(long *idum); float ran1(long *idum); float ran2(long *idum); float ran3(long *idum); float ran4(long *idum); void rank(unsigned long n, unsigned long indx[], unsigned long irank[]); void ranpt(float pt[], float regn[], int n); void ratint(float xa[], float ya[], int n, float x, float *y, float *dy); void ratlsq(double (*fn)(double), double a, double b, int mm, int kk, double cof[], double *dev); double ratval(double x, double cof[], int mm, int kk); float rc(float x, float y); float rd(float x, float y, float z); void realft(float data[], unsigned long n, int isign); void rebin(float rc, int nd, float r[], float xin[], float xi[]); void red(int iz1, int iz2, int jz1, int jz2, int jm1, int jm2, int jmf, int ic1, int jc1, int jcf, int kc, float ***c, float **s); void relax(double **u, double **rhs, int n); void relax2(double **u, double **rhs, int n); void resid(double **res, double **u, double **rhs, int n); float revcst(float x[], float y[], int iorder[], int ncity, int n[]); void reverse(int iorder[], int ncity, int n[]);