938 Appendix A void spread(float y,float yy],unsigned long n,float x,int m); void sprsax(f1 oat sa0,unsigned long ija☐，f1oatx[l,f1oatb[☐ unsigned long n); void sprsin(float **a,int n,float thresh, unsigned long nmax,float sa], nsigned1ong1ja[☐)； void sprspm(float sa],unsigned long ija[],float sb[],unsigned long ijb[], float sc[],unsigned long ijc[]); void sprstm(float sa[],unsigned long ija],float sb[],unsigned long ijb], float thresh,unsigned long nmax,float sc[],unsigned long ijc[]); void sprstp(float sa[],unsigned long ija[],float sb[],unsigned long ijb[]); void sprstx(float sa],unsigned long ija[],float x[],float b[], unsigned long n); void stifbs(float y[],float dydx[],int nv,float *xx, float htry,float eps,float yscal[,float *hdid,float *hnext, void (*derivs)(float,float []float [])) void stiff(float y[],float dydx[],int n,float *x, granted for from NUMERICAL RECIPES I 18881892 float htry,float eps,float yscal[],float *hdid,float *hnext, void(*derivs)(f1oat,f1oat[，f1oat)); 1.800 void stoerm(float y[],float d2y],int nv,float xs, float htot,int nstep,float yout[], void (*derivs)(float,float []float [])) void svbksb(float *u,float w[],float **v,int m,int n,float b[], to any Cambridge float x); void svdcmp(float **a,int m,int n,float w[],float **v); (North America to make 州bMe se e University Press. 2 void svdfit(float x[],float y[],float sig[],int ndata,float a], int ma,float **u,float **v,float w[],float *chisq, one paper THE void(*funcs.)(f1oat,f1oat☐，int)); void svdvar(float v,int ma,float w[,float **cvm); 是 ART void toeplz(float r[],float x[],float y[],int n); void tptest(float datal,float data2,unsigned long n,float *t,float *prob); void tqli(float d[],float e[],int n,float **z); Programs float trapzd(float (*func)(float),float a,float b,int n); void tred2(float **a,int n,float d[],float e[]); void tridag(float a☐，f1oatb[],f1oatc☐，f1oatr0,f1oatu▣ unsigned long n); float trncst(float x[],float y[],int iorder[],int ncity,int n[]); to dir void trnspt(int iorder[],int ncity,int n[]); void ttest(float datai],unsigned long n1,float data2[],unsigned long n2, float *t,float *prob); 19881982 OF SCIENTIFIC COMPUTING(ISBN void tutest(float datal] unsigned long n1,float data2[],unsigned long n2, float *t,float *prob); void twofft(float data1☐，f1 oat data2[☐，f1 oat fft1[0,f1 oat fft:2[], unsigned long n); 10-621 void vander(double x[],double w[],double g[],int n); void vegas(f1 oat regn[☐，int ndim,f1oat(wfxn)(f1oat[☐，f1oat),int init, unsigned long ncall,int itmx,int nprn,float *tgral,float *sd, -43106 float *chi2a); Numerical Recipes void voltra(int n,int m,float to,float h,float *t,float **f float (g)(int,float),float (*ak)(int,int,float,float)); (outside void wt1(float a[],unsigned long n,int isign, void (*wtstep)(float ]unsigned long,int)); Software. void wtn(f1oata☐，unsigned long nn☐，int ndim,int1sign, void (*wtstep)(float []unsigned long,int)); Amer ying of void wwghts(float wghts[],int n,float h, void(*kermom)(double☐，double,int)); int zbrac(float (*func)(float),float *x1,float *x2); void zbrak(float (*fx)(float),float x1,float x2,int n,float xbi[], float xb2[],int *nb); float zbrent(float (*func)(float),float x1,float x2,float tol); void zrhgr(float a[],int m,float rtr[],float rti[]); float zriddr(float (*func)(float),float x1,float x2,float xacc); void zroots(fcomplex a[□，1ntm,fcomplex roots[□，int polish); #else/米AWSI*/938 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 spread(float y, float yy[], unsigned long n, float x, int m); void sprsax(float sa[], unsigned long ija[], float x[], float b[], unsigned long n); void sprsin(float **a, int n, float thresh, unsigned long nmax, float sa[], unsigned long ija[]); void sprspm(float sa[], unsigned long ija[], float sb[], unsigned long ijb[], float sc[], unsigned long ijc[]); void sprstm(float sa[], unsigned long ija[], float sb[], unsigned long ijb[], float thresh, unsigned long nmax, float sc[], unsigned long ijc[]); void sprstp(float sa[], unsigned long ija[], float sb[], unsigned long ijb[]); void sprstx(float sa[], unsigned long ija[], float x[], float b[], unsigned long n); void stifbs(float y[], float dydx[], int nv, float *xx, float htry, float eps, float yscal[], float *hdid, float *hnext, void (*derivs)(float, float [], float [])); void stiff(float y[], float dydx[], int n, float *x, float htry, float eps, float yscal[], float *hdid, float *hnext, void (*derivs)(float, float [], float [])); void stoerm(float y[], float d2y[], int nv, float xs, float htot, int nstep, float yout[], void (*derivs)(float, float [], float [])); void svbksb(float **u, float w[], float **v, int m, int n, float b[], float x[]); void svdcmp(float **a, int m, int n, float w[], float **v); void svdfit(float x[], float y[], float sig[], int ndata, float a[], int ma, float **u, float **v, float w[], float *chisq, void (*funcs)(float, float [], int)); void svdvar(float **v, int ma, float w[], float **cvm); void toeplz(float r[], float x[], float y[], int n); void tptest(float data1[], float data2[], unsigned long n, float *t, float *prob); void tqli(float d[], float e[], int n, float **z); float trapzd(float (*func)(float), float a, float b, int n); void tred2(float **a, int n, float d[], float e[]); void tridag(float a[], float b[], float c[], float r[], float u[], unsigned long n); float trncst(float x[], float y[], int iorder[], int ncity, int n[]); void trnspt(int iorder[], int ncity, int n[]); void ttest(float data1[], unsigned long n1, float data2[], unsigned long n2, float *t, float *prob); void tutest(float data1[], unsigned long n1, float data2[], unsigned long n2, float *t, float *prob); void twofft(float data1[], float data2[], float fft1[], float fft2[], unsigned long n); void vander(double x[], double w[], double q[], int n); void vegas(float regn[], int ndim, float (*fxn)(float [], float), int init, unsigned long ncall, int itmx, int nprn, float *tgral, float *sd, float *chi2a); void voltra(int n, int m, float t0, float h, float *t, float **f, float (*g)(int, float), float (*ak)(int, int, float, float)); void wt1(float a[], unsigned long n, int isign, void (*wtstep)(float [], unsigned long, int)); void wtn(float a[], unsigned long nn[], int ndim, int isign, void (*wtstep)(float [], unsigned long, int)); void wwghts(float wghts[], int n, float h, void (*kermom)(double [], double ,int)); int zbrac(float (*func)(float), float *x1, float *x2); void zbrak(float (*fx)(float), float x1, float x2, int n, float xb1[], float xb2[], int *nb); float zbrent(float (*func)(float), float x1, float x2, float tol); void zrhqr(float a[], int m, float rtr[], float rti[]); float zriddr(float (*func)(float), float x1, float x2, float xacc); void zroots(fcomplex a[], int m, fcomplex roots[], int polish); #else /* ANSI */