/* CAUTION: This is the ANSI C (only) version of the Numerical Recipes   utility file nr.h.  Do not confuse this file with the same-named   file nr.h that is supplied in the 'misc' subdirectory.   *That* file is the one from the book, and contains both ANSI and   traditional K&R versions, along with #ifdef macros to select the   correct version.  *This* file contains only ANSI C.               */#ifndef _NR_H_#define _NR_H_#ifndef _FCOMPLEX_DECLARE_T_typedef struct FCOMPLEX {float r,i;} fcomplex;#define _FCOMPLEX_DECLARE_T_#endif /* _FCOMPLEX_DECLARE_T_ */#ifndef _ARITHCODE_DECLARE_T_typedef struct {	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>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);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);float brent_arg(float ax, float bx, float cx,	float (*f)(float, void*), float tol, float *xmin, void *arg);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[]));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 f1dim_arg(float x, void *arg);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 frprmn_arg(float p[], int n, float ftol, int *iter, float *fret,	float (*func)(float [], void*), void (*dfunc)(float [], float [], void*), void* arg);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);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,