《数字信号处理》教学参考资料（Numerical Recipes in C，The Art of Scientific Computing Second Edition）Chapter 21.1 References

Table of Prototype Declarations (nr.h) 931 unsigned long *ilob,*iupb,*ncumfq,jdif,nc,minint,nch,ncum,nrad; arithcode; #define_ARITHCODE_DECLARE_T_ #endif /*ARITHCODE_DECLARE_T_/ #ifndef_HUFFCODE_DECLARE_T_ typedef struct unsigned long *icod,*ncod,*left,*right,nch,nodemax; huffcode; #define_HUFFCODE_DECLARE_T_ #endif /*_HUFFCODE_DECLARE_T_*/ http://www.n read able files Permission is #include <stdio.h> #if defined(STDC)II defined(ANSI)I defined(NRANSI)/ANSI/ 83g void addint(double *uf,double **uc,double **res,int nf); void airy(float x,float *ai,float *bi,float *aip,float *bip); (including this one) void amebsa(float **p,float y],int ndim,float pb[], float *yb, float ftol,float (*funk)(float [])int *iter,float temptr); 11-800-872 void amoeba(float **p,float y[],int ndim,float ftol, float (*funk)(float [])int *iter); 7423 to any server computer, tusers to make one paper 1988-1992 by Cambridge University Press.Programs from NUMERICAL RECIPES IN float amotry(float **p,float y[],float psum[],int ndim, float (*funk)(float [])int ihi,float fac); float amotsa(float **p,float y[],float psum[],int ndim,float pb[], float *yb,float (*funk)(float [])int ihi,float *yhi,float fac); void anneal(float x[],float y[],int iorder[],int ncity); (North America THE double anorm2(double **a,int n); ART void arcmak(unsigned long nfreg[], unsigned long nchh,unsigned long nradd, arithcode *acode); only). void arcode(unsigned long *ich,unsigned char **codep,unsigned long *lcode, unsigned long *lcd,int isign,arithcode tacode); void arcsum(unsigned long iin[],unsigned long iout[],unsigned long ja, int nwk,unsigned long nrad,unsigned long nc); void asolve(unsigned long n,double b[],double x[],int itrnsp); void atimes(unsigned long n,double x],double r[],int itrnsp); to dir void avevar(float data[],unsigned long n,float *ave,float *var); void balanc(float *a,int n); void banbks(float **a, unsigned long n,int m1,int m2,float **al, unsigned long indx,float b[]); ectcustser void bandec(float *a,unsigned long n,int mi,int m2,float *al, OF SCIENTIFIC COMPUTING(ISBN 0-521 unsigned long indx☐，float*d); void banmul(float **a,unsigned long n,int mi,int m2,float x[],float b[]); void bcucof(f1oaty[],f1oaty1[☐，f1oaty2[☐，f1oaty120,f1oatd1, float d2,float **c); void bcuint(float y[],float yi[],float y2[],float y12[], 1988-1992 by Numerical Recipes -431085 float xil,float xiu,float x21,float x2u,float x1, float x2,float *ansy,float *ansy1,float *ansy2); void beschb(double x,double *gam1,double *gam2,double *gampl, double *gammi); float bessi(int n,float x); float bessio(float x); float bessi1(float x); (outside North America) Software. ying of void bessik(float x,float xnu,float *ri,float *rk,float *rip, float *rkp); float bessj(int n,float x); float bessjo(float x); float bessi1(float x); void bessjy(float x,float xnu,float *rj,float *ry,float *rjp, float *ryp); float bessk(int n,float x); float bessko(float x); float bessk1(float x); float bessy(int n,float x);

Table of Prototype Declarations (nr.h) 931 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). unsigned long *ilob,*iupb,*ncumfq,jdif,nc,minint,nch,ncum,nrad; } arithcode; #define _ARITHCODE_DECLARE_T_ #endif /* _ARITHCODE_DECLARE_T_ */ #ifndef _HUFFCODE_DECLARE_T_ typedef struct { unsigned long *icod,*ncod,*left,*right,nch,nodemax; } huffcode; #define _HUFFCODE_DECLARE_T_ #endif /* _HUFFCODE_DECLARE_T_ */ #include <stdio.h> #if defined(__STDC__) || defined(ANSI) || defined(NRANSI) /* ANSI */ void addint(double **uf, double **uc, double **res, int nf); void airy(float x, float *ai, float *bi, float *aip, float *bip); void amebsa(float **p, float y[], int ndim, float pb[], float *yb, float ftol, float (*funk)(float []), int *iter, float temptr); void amoeba(float **p, float y[], int ndim, float ftol, float (*funk)(float []), int *iter); float amotry(float **p, float y[], float psum[], int ndim, float (*funk)(float []), int ihi, float fac); float amotsa(float **p, float y[], float psum[], int ndim, float pb[], float *yb, float (*funk)(float []), int ihi, float *yhi, float fac); void anneal(float x[], float y[], int iorder[], int ncity); double anorm2(double **a, int n); void arcmak(unsigned long nfreq[], unsigned long nchh, unsigned long nradd, arithcode *acode); void arcode(unsigned long *ich, unsigned char **codep, unsigned long *lcode, unsigned long *lcd, int isign, arithcode *acode); void arcsum(unsigned long iin[], unsigned long iout[], unsigned long ja, int nwk, unsigned long nrad, unsigned long nc); void asolve(unsigned long n, double b[], double x[], int itrnsp); void atimes(unsigned long n, double x[], double r[], int itrnsp); void avevar(float data[], unsigned long n, float *ave, float *var); void balanc(float **a, int n); void banbks(float **a, unsigned long n, int m1, int m2, float **al, unsigned long indx[], float b[]); void bandec(float **a, unsigned long n, int m1, int m2, float **al, unsigned long indx[], float *d); void banmul(float **a, unsigned long n, int m1, int m2, float x[], float b[]); void bcucof(float y[], float y1[], float y2[], float y12[], float d1, float d2, float **c); void bcuint(float y[], float y1[], float y2[], float y12[], float x1l, float x1u, float x2l, float x2u, float x1, float x2, float *ansy, float *ansy1, float *ansy2); void beschb(double x, double *gam1, double *gam2, double *gampl, double *gammi); float bessi(int n, float x); float bessi0(float x); float bessi1(float x); void bessik(float x, float xnu, float *ri, float *rk, float *rip, float *rkp); float bessj(int n, float x); float bessj0(float x); float bessj1(float x); void bessjy(float x, float xnu, float *rj, float *ry, float *rjp, float *ryp); float bessk(int n, float x); float bessk0(float x); float bessk1(float x); float bessy(int n, float x);

932 Appendix A float bessyo(float x); float bessy1(float x); float beta(float z,float w); float betacf(float a,float b,float x); float betai(float a,float b,float x); float bico(int n,int k); void bksub(int ne,int nb,int jf,int ki,int k2,float ***c); float bnldev(float pp,int n,long *idum); float brent(float ax,float bx,float cx, float (*f)(float),float tol,float *xmin); void broydn(float x[],int n,int *check, void(*vecfunc)(int,f1oat0,f1oat[☐))； http://www.nr read able files 后Per void bsstep(float y[],float dydx[],int nv,float *xx,float htry, float eps,float yscal[],float *hdid,float *hnext, void (*derivs)(float,float []float [])) 83g void caldat(long julian,int *mm,int *id,int *iyyy); void chder(float a,float b,float c[],float cder[],int n); float chebev(float a,float b,float c[],int m,float x); (including this one) void chebft(float a,float b,float c,int n,float (+func)(float)); void chebpc(float c[],float d[],int n); 11-800-872 void chint(float a,float b,float c],float cint,int n); float chixy(float bang); 7423 from NUMERICAL RECIPES IN void choldc(float *a,int n,float p); void cholsl(float **a,int n,float p[],float b[],float x[]); void chsone(float bins[],float ebins],int nbins,int knstrn, float *df,float *chsq,float *prob); void chstwo(float binsi],float bins2[],int nbins,int knstrn, (North America 州Me:Is to tusers to make one paper 1988-1992 by Cambridge University Press.Programs THE float *df,float *chsq,float *prob); ART void cisi(float x,float *ci,float *si); void cntabi(int **nn,int ni,int nj,float *chisq only). float *df,float *prob,float *cramrv,float *ccc); void cntab2(int **nn,int ni,int nj,float *h,float *hx,float why, float *hygx,float *hxgy,float *uygx,float tuxgy,float *uxy); void convlv(float data[],unsigned long n,float respns[],unsigned long m, st st int isign,float ans[]); void copy(double **aout,double **ain,int n); to dir Copyright (C) void correl(float datal,float data2],unsigned long n,float ans[]); void cosft(f1oaty[☐，intn,int isign); void cosfti(float y[],int n); OF SCIENTIFIC COMPUTING(ISBN void cosft2(f1oaty[☐，intn,int1sigr); void covsrt(float **covar,int ma,int ia[],int mfit); void crank(unsigned long n,floatw],float *s); void cyclic(float a[],float b[],float c[],float alpha,float beta 12-.621 float r[],float x,unsigned long n); void daub4(float a[],unsigned long n,int isign); float dawson(float x); To order Numerical Recipes books or personal use.Further reproduction, 1988-1992 by Numerical Recipes 43108 float dbrent(float ax,float bx,float cx float (*f)(float),float (*df)(float),float tol,float *xmin); void ddpoly(float c[],int nc,float x,float pd[],int nd); int decchk(char string[],int n,char *ch); void derivs(float x,float y[],float dydx[]); Software. float dfidim(float x); void dfour1(double data],unsigned long nn,int isign); (outside North America) ying of void dfpmin(float p[],int n,float gtol,int *iter,float *fret, float(*func)(f1oat[门)，void(*dfunc)(f1oat☐，float[J)); float dfridr(float (*func)(float),float x,float h,float err); void dftcor(float w,float delta,float a,float b,float endpts[], float *corre,float *corim,float *corfac); void dftint(float (*func)(float),float a,float b,float w, float *cosint,float *sinint); void difeq(int k,int k1,int k2,int jsf,int is1,int isf, int indexv[],int ne,float **s,float **y); void dlinmin(float p,float xi[],int n,float *fret, float (*func)(float [])void (*dfunc)(float []float[]));

932 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). float bessy0(float x); float bessy1(float x); float beta(float z, float w); float betacf(float a, float b, float x); float betai(float a, float b, float x); float bico(int n, int k); void bksub(int ne, int nb, int jf, int k1, int k2, float ***c); float bnldev(float pp, int n, long *idum); float brent(float ax, float bx, float cx, float (*f)(float), float tol, float *xmin); void broydn(float x[], int n, int *check, void (*vecfunc)(int, float [], float [])); void bsstep(float y[], float dydx[], int nv, float *xx, float htry, float eps, float yscal[], float *hdid, float *hnext, void (*derivs)(float, float [], float [])); void caldat(long julian, int *mm, int *id, int *iyyy); void chder(float a, float b, float c[], float cder[], int n); float chebev(float a, float b, float c[], int m, float x); void chebft(float a, float b, float c[], int n, float (*func)(float)); void chebpc(float c[], float d[], int n); void chint(float a, float b, float c[], float cint[], int n); float chixy(float bang); void choldc(float **a, int n, float p[]); void cholsl(float **a, int n, float p[], float b[], float x[]); void chsone(float bins[], float ebins[], int nbins, int knstrn, float *df, float *chsq, float *prob); void chstwo(float bins1[], float bins2[], int nbins, int knstrn, float *df, float *chsq, float *prob); void cisi(float x, float *ci, float *si); void cntab1(int **nn, int ni, int nj, float *chisq, float *df, float *prob, float *cramrv, float *ccc); void cntab2(int **nn, int ni, int nj, float *h, float *hx, float *hy, float *hygx, float *hxgy, float *uygx, float *uxgy, float *uxy); void convlv(float data[], unsigned long n, float respns[], unsigned long m, int isign, float ans[]); void copy(double **aout, double **ain, int n); void correl(float data1[], float data2[], unsigned long n, float ans[]); void cosft(float y[], int n, int isign); void cosft1(float y[], int n); void cosft2(float y[], int n, int isign); void covsrt(float **covar, int ma, int ia[], int mfit); void crank(unsigned long n, float w[], float *s); void cyclic(float a[], float b[], float c[], float alpha, float beta, float r[], float x[], unsigned long n); void daub4(float a[], unsigned long n, int isign); float dawson(float x); float dbrent(float ax, float bx, float cx, float (*f)(float), float (*df)(float), float tol, float *xmin); void ddpoly(float c[], int nc, float x, float pd[], int nd); int decchk(char string[], int n, char *ch); void derivs(float x, float y[], float dydx[]); float df1dim(float x); void dfour1(double data[], unsigned long nn, int isign); void dfpmin(float p[], int n, float gtol, int *iter, float *fret, float (*func)(float []), void (*dfunc)(float [], float [])); float dfridr(float (*func)(float), float x, float h, float *err); void dftcor(float w, float delta, float a, float b, float endpts[], float *corre, float *corim, float *corfac); void dftint(float (*func)(float), float a, float b, float w, float *cosint, float *sinint); void difeq(int k, int k1, int k2, int jsf, int is1, int isf, int indexv[], int ne, float **s, float **y); void dlinmin(float p[], float xi[], int n, float *fret, float (*func)(float []), void (*dfunc)(float [], float[]));

Table of Prototype Declarations (nr.h) 933 double dpythag(double a,double b); void drealft(double data],unsigned long n,int isign); void dsprsax(double sa[],unsigned long ija[],double x[],double b[], unsigned long n); void dsprstx(double sa[],unsigned long ija[],double x[],double b[] nsig知ed long n)； void dsvbksb(double **u,double w[],double **v,int m,int n,double b[], double x[]); void dsvdcmp(double **a,int m,int n,double w[],double **v); void eclass(int nf[],int n,int lista[],int listb[],int m); void eclazz(int nf[],int n,int (*equiv)(int,int)); float ei(float x); http://www.nr read able files Permission is void eigsrt(float d[],float **v,int n); float elle(float phi,float ak); float ellf(float phi,float ak); 839 float ellpi(float phi,float en,float ak); granted fori void elmhes(float 来来a.1ntn)t float erfcc(float x); (including this one) float erff(float x); 11-800-872 internet float erffc(float x); void eulsum(float *sum, float term,int jterm,float wksp[]); float evlmem(float fdt,float d[],int m,float xms); to any server computer,is tusers to make one paper 1988-1992 by Cambridge University Press.Programs from NUMERICAL RECIPES IN float expdev(long *idum); float expint(int n,float x); float f1(float x); float fidim(float x); float f2(float y); (North America THE float f3(float z); ART float factln(int n); float factrl(int n); s. void fasper(float x[],float y[],unsigned long n,float ofac,float hifac, float wk1[],float wk2[],unsigned long nwk,unsigned long *nout, unsigned long *jmax,float *prob); void fdjac(int n,float x[],float fvec[],float **df st st void(*vecfunc)(int,f1oat☐，f1oat[U)); void fgauss(float x,float a[],float *y,float dyda,int na); to dir Copyright (C) void fillo(double *u.int n): void fit(float x,float y[],int ndata,float sig[],int mwt, float *a,float *b,float *siga,float *sigb,float *chi2,float *q) rectcustser umencal fecpes o00 1988-1992 by Numerical Recipes OF SCIENTIFIC COMPUTING(ISBN void fitexy(float x[],float y,int ndat,float sigx[],float sigy[, float *a,float *b,float *siga,float *sigb,float *chi2,float q); void fixrts(f1oatd☐，intm); void f1eg(f1oatx,f1 oat pl[☐，int nl); 10-621 void flmoon(int n,int nph,long *jd,float *frac); f1 oat fmin(f1oatx[□)； void four1(float data[],unsigned long nn,int isign); -43108 void fourew(FILE *file[5],int *na,int *nb,int *nc,int *nd); void fourfs(FILE *file[5],unsigned long nn[],int ndim,int isign); void fourn(float data[],unsigned long nn[],int ndim,int isign); void fpoly(float x,float p[],int np); void fred2(int n,float a,float b,float t[],float f[],floatw (outside North Software. float (g)(float),float (*ak)(float,float)); float fredin(float x,int n,float a,float b,float t[],float f[],floatw], float (g)(float),float (*ak)(float,float)); void frenel(float x,float *s,float *c); void frprmn(float p[],int n,float ftol,int *iter,float *fret, float (*func)(float [])void (*dfunc)(float []float [])) void ftest(float datal[],unsigned long n1,float data2[],unsigned long n2, float *f,float *prob); float gamdev(int ia,long *idum); float gammln(float xx); float gammp(float a,float x); float gammq(float a,float x); float gasdev(long *idum);

Table of Prototype Declarations (nr.h) 933 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). double dpythag(double a, double b); void drealft(double data[], unsigned long n, int isign); void dsprsax(double sa[], unsigned long ija[], double x[], double b[], unsigned long n); void dsprstx(double sa[], unsigned long ija[], double x[], double b[], unsigned long n); void dsvbksb(double **u, double w[], double **v, int m, int n, double b[], double x[]); void dsvdcmp(double **a, int m, int n, double w[], double **v); void eclass(int nf[], int n, int lista[], int listb[], int m); void eclazz(int nf[], int n, int (*equiv)(int, int)); float ei(float x); void eigsrt(float d[], float **v, int n); float elle(float phi, float ak); float ellf(float phi, float ak); float ellpi(float phi, float en, float ak); void elmhes(float **a, int n); float erfcc(float x); float erff(float x); float erffc(float x); void eulsum(float *sum, float term, int jterm, float wksp[]); float evlmem(float fdt, float d[], int m, float xms); float expdev(long *idum); float expint(int n, float x); float f1(float x); float f1dim(float x); float f2(float y); float f3(float z); float factln(int n); float factrl(int n); void fasper(float x[], float y[], unsigned long n, float ofac, float hifac, float wk1[], float wk2[], unsigned long nwk, unsigned long *nout, unsigned long *jmax, float *prob); void fdjac(int n, float x[], float fvec[], float **df, void (*vecfunc)(int, float [], float [])); void fgauss(float x, float a[], float *y, float dyda[], int na); void fill0(double **u, int n); void fit(float x[], float y[], int ndata, float sig[], int mwt, float *a, float *b, float *siga, float *sigb, float *chi2, float *q); void fitexy(float x[], float y[], int ndat, float sigx[], float sigy[], float *a, float *b, float *siga, float *sigb, float *chi2, float *q); void fixrts(float d[], int m); void fleg(float x, float pl[], int nl); void flmoon(int n, int nph, long *jd, float *frac); float fmin(float x[]); void four1(float data[], unsigned long nn, int isign); void fourew(FILE *file[5], int *na, int *nb, int *nc, int *nd); void fourfs(FILE *file[5], unsigned long nn[], int ndim, int isign); void fourn(float data[], unsigned long nn[], int ndim, int isign); void fpoly(float x, float p[], int np); void fred2(int n, float a, float b, float t[], float f[], float w[], float (*g)(float), float (*ak)(float, float)); float fredin(float x, int n, float a, float b, float t[], float f[], float w[], float (*g)(float), float (*ak)(float, float)); void frenel(float x, float *s, float *c); void frprmn(float p[], int n, float ftol, int *iter, float *fret, float (*func)(float []), void (*dfunc)(float [], float [])); void ftest(float data1[], unsigned long n1, float data2[], unsigned long n2, float *f, float *prob); float gamdev(int ia, long *idum); float gammln(float xx); float gammp(float a, float x); float gammq(float a, float x); float gasdev(long *idum);

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);

Table of Prototype Declarations (nr.h) 935 void lop(double **out,double **u,int n); void lubksb(float *a,int n,int *indx,float b[]); void ludcmp(float **a,int n,int *indx,float *d); void machar(int *ibeta,int it,int *irnd,int *ngrd, int *machep,int *negep,int *iexp,int *minexp,int *maxexp float *eps,float *epsneg,float *xmin,float *xmax); void matadd(double *a,double **b,double **c,int n); void matsub(double +a,double **b,double **c,int n); void medfit(float x[],float y],int ndata,float *a,float *b,float *abdev); void memcof(float data[],int n,int m,float *xms,float d[]); int metrop(float de,float t); void mgfas(double **u,int n,int maxcyc); void mglin(double **u,int n,int ncycle); float midexp(float (*funk)(float),float aa,float bb,int n); float midinf(float (*funk)(float) float aa,float bb,int n); float midpnt(float (*func)(float),float a,float b,int n); granted for 18881892 float midsql(float (*funk)(float),float aa,float bb,int n); float midsqu(float (*funk)(float) float aa,float bb,int n); 1800 void miser(float (*func)(float [])float regn],int ndim,unsigned long npts, float dith,float ave,float *var); void mmid(float y[],float dydx[],int nvar,float xs,float htot, Cambridge from NUMERICAL RECIPES I int nstep,float yout []void (*derivs)(float,float []float[])); void mnbrak(float *ax,float *bx,float *cx,float *fa,float *fb, float *fc,float (*func)(float)); void mnewt(int ntrial,float x[],int n,float tolx,float tolf); (Nor void moment(float data[],int n,float ave,float *adev,float *sdev America server computer, e University Press. THE float *var,float *skew,float *curt); void mp2dfr(unsigned char a[],unsigned char s[],int n,int *m); one paper ART void mpadd(unsigned char w[],unsigned char u[],unsigned char v[],int n); void mpdiv(unsigned char g],unsigned char r[],unsigned char u[], unsigned char v[],int n,int m); Programs void mpinv(unsigned char u[],unsigned char v[],int n,int m); void mplsh(unsigned char un,int n); void mpmov(unsigned char u[],unsigned char v[],int n); void mpmul(unsigned char w[],unsigned char u[],unsigned char v[],int n, int m); to dir void mpneg(unsigned char un],int n); void mppi(int n); void mprove(float **a,float **alud,int n,int indx[],float b[] SCIENTIFIC COMPUTING(ISBN float x[); 188819920 void mpsad(unsigned char w[],unsigned char u[],int n,int iv); void mpsdv(unsigned char w[],unsigned char u[],int n,int iv,int *ir); void mpsmu(unsigned char w[], unsigned char u[],int n,int iv); 10-621 void mpsqrt(unsigned char w[],unsigned char u[],unsigned char v[],int n, int m); void mpsub(int*is,unsigned char w[],unsigned char u[☐，unsigned char v☐， Numerical Recipes -43106 int n); void mrqcof(float x[],float y[],float sig[],int ndata,float a, int ia[□，int ma,f1oat**alpha,f1 oat beta☐，f1oat*chisq, (outside void (*funcs)(float,float [float *float []int)); void mrqmin(float x[],float y[],float sig[],int ndata,float a[] Software. int ia[],int ma,float **covar,float *alpha,float *chisq, void (*funcs)(float,float []float *float []int),float *alamda); void newt(float x[],int n,int *check, void (*vecfunc)(int,float ]float [])) void odeint(float ystart[],int nvar,float x1,float x2, float eps,float hi,float hmin,int *nok,int *nbad, void (*derivs)(float,float []float []) void (*rkgs)(float []float []int,float *float,float, float [float *float *void (*)(float,float []float []))) void orthog(int n,float anu[,float alpha[],float beta[],float a] f1oatb[▣)； void pade(double cof☐，intn,float*resid); void pccheb(float d[],float c[],int n);

Table of Prototype Declarations (nr.h) 935 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 lop(double **out, double **u, int n); void lubksb(float **a, int n, int *indx, float b[]); void ludcmp(float **a, int n, int *indx, float *d); void machar(int *ibeta, int *it, int *irnd, int *ngrd, int *machep, int *negep, int *iexp, int *minexp, int *maxexp, float *eps, float *epsneg, float *xmin, float *xmax); void matadd(double **a, double **b, double **c, int n); void matsub(double **a, double **b, double **c, int n); void medfit(float x[], float y[], int ndata, float *a, float *b, float *abdev); void memcof(float data[], int n, int m, float *xms, float d[]); int metrop(float de, float t); void mgfas(double **u, int n, int maxcyc); void mglin(double **u, int n, int ncycle); float midexp(float (*funk)(float), float aa, float bb, int n); float midinf(float (*funk)(float), float aa, float bb, int n); float midpnt(float (*func)(float), float a, float b, int n); float midsql(float (*funk)(float), float aa, float bb, int n); float midsqu(float (*funk)(float), float aa, float bb, int n); void miser(float (*func)(float []), float regn[], int ndim, unsigned long npts, float dith, float *ave, float *var); void mmid(float y[], float dydx[], int nvar, float xs, float htot, int nstep, float yout[], void (*derivs)(float, float[], float[])); void mnbrak(float *ax, float *bx, float *cx, float *fa, float *fb, float *fc, float (*func)(float)); void mnewt(int ntrial, float x[], int n, float tolx, float tolf); void moment(float data[], int n, float *ave, float *adev, float *sdev, float *var, float *skew, float *curt); void mp2dfr(unsigned char a[], unsigned char s[], int n, int *m); void mpadd(unsigned char w[], unsigned char u[], unsigned char v[], int n); void mpdiv(unsigned char q[], unsigned char r[], unsigned char u[], unsigned char v[], int n, int m); void mpinv(unsigned char u[], unsigned char v[], int n, int m); void mplsh(unsigned char u[], int n); void mpmov(unsigned char u[], unsigned char v[], int n); void mpmul(unsigned char w[], unsigned char u[], unsigned char v[], int n, int m); void mpneg(unsigned char u[], int n); void mppi(int n); void mprove(float **a, float **alud, int n, int indx[], float b[], float x[]); void mpsad(unsigned char w[], unsigned char u[], int n, int iv); void mpsdv(unsigned char w[], unsigned char u[], int n, int iv, int *ir); void mpsmu(unsigned char w[], unsigned char u[], int n, int iv); void mpsqrt(unsigned char w[], unsigned char u[], unsigned char v[], int n, int m); void mpsub(int *is, unsigned char w[], unsigned char u[], unsigned char v[], int n); void mrqcof(float x[], float y[], float sig[], int ndata, float a[], int ia[], int ma, float **alpha, float beta[], float *chisq, void (*funcs)(float, float [], float *, float [], int)); void mrqmin(float x[], float y[], float sig[], int ndata, float a[], int ia[], int ma, float **covar, float **alpha, float *chisq, void (*funcs)(float, float [], float *, float [], int), float *alamda); void newt(float x[], int n, int *check, void (*vecfunc)(int, float [], float [])); void odeint(float ystart[], int nvar, float x1, float x2, float eps, float h1, float hmin, int *nok, int *nbad, void (*derivs)(float, float [], float []), void (*rkqs)(float [], float [], int, float *, float, float, float [], float *, float *, void (*)(float, float [], float []))); void orthog(int n, float anu[], float alpha[], float beta[], float a[], float b[]); void pade(double cof[], int n, float *resid); void pccheb(float d[], float c[], int n);