		float [], float *, float *, void (*)(float, float [], float [])));
		float [], float *, float *, void (*)(float, float [], float [])));
		float [], float *, float *, void (*)(float, float [], float [])));
		float yout[], float yerr[], void (*derivs)(float, float [], float []));
		int nstep, float yout[], void (*derivs)(float, float[], float[]));
		void (*derivs)(float, float [], float []));
		void (*derivs)(float, float [], float []));
		void (*derivs)(float, float [], float []));
		void (*derivs)(float, float [], float []));
		void (*derivs)(float, float [], float []));
		void (*derivs)(float, float [], float []),
		void (*derivs)(float, float [], float []),
		void (*derivs)(float, float [], float []),
		void (*funcs)(float, float [], float *, float [], int));
		void (*kermom)(double [], double ,int));
		void (*rkqs)(float [], float [], int, float *, float, float,
		void (*rkqs)(float [], float [], int, float *, float, float,
		void (*rkqs)(float [], float [], int, float *, float, float,
		void (*vecfunc)(int, float [], float []));
		void (*vecfunc)(int, float [], float []));
		void (*vecfunc)(int, float [], float []));
		void (*vecfunc)(int, float [], float []));
	float (*func)(float []), void (*dfunc)(float [], float []))
	float *, float *, void (*)(float, float [], float [])))
	float yerr[], void (*derivs)(float, float [], float []))
	float yout[], void (*derivs)(float, float [], float []))
	float yout[], void (*derivs)(float, float[], float[]))
	float(*func)(float []), void (*dfunc)(float [], float []))
	int ma, float **covar, float *chisq, void (*funcs)(float, float [], int))
	void (*derivs)(float, float [], float []))
	void (*derivs)(float, float [], float []))
	void (*derivs)(float, float [], float []))
	void (*derivs)(float, float [], float []))
	void (*derivs)(float, float [], float []))
	void (*derivs)(float, float [], float []))
	void (*derivs)(float, float [], float []))
	void (*derivs)(float, float [], float []),
	void (*dfunc)(float [], float []))
	void (*funcs)(float, float [], float *, float [], int))
	void (*funcs)(float, float [], float *, float [], int), float *alamda)
	void (*funcs)(float, float [], int))
	void (*kermom)(double [], double ,int))
	void (*quadvl)(float, float, float *, float *, float *, float *),
	void (*rkqs)(float [], float [], int, float *, float, float, float [],
	void (*vecfunc)(int, float [], float []))
	void (*vecfunc)(int, float [], float []))
	void (*vecfunc)(int, float [], float []))
	void (*wtstep)(float [], unsigned long, int))
	void (*wtstep)(float [], unsigned long, int))
	void addint(double **uf, double **uc, double **res, int nf);
	void arcsum(unsigned long iin[], unsigned long iout[], unsigned long ja,
	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 avevar(float data[], unsigned long n, float *ave, float *var);
	void avevar(float data[], unsigned long n, float *ave, float *var);
	void avevar(float data[], unsigned long n, float *ave, float *var);
	void avevar(float data[], unsigned long n, float *ave, float *var);
	void avevar(float data[], unsigned long n, float *ave, float *var);
	void avevar(float data[], unsigned long n, float *ave, float *var);
	void balanc(float **a, int n);
	void bcucof(float y[], float y1[], float y2[], float y12[], float d1,
	void beschb(double x, double *gam1, double *gam2, double *gampl,
	void beschb(double x, double *gam1, double *gam2, double *gampl,
	void bessik(float x, float xnu, float *ri, float *rk, float *rip,
	void bessjy(float x, float xnu, float *rj, float *ry, float *rjp,
	void bessjy(float x, float xnu, float *rj, float *ry, float *rjp,
	void bksub(int ne, int nb, int jf, int k1, int k2, float ***c);
	void bsstep(float y[], float dydx[], int nv, float *xx, float htry,
	void copy(double **aout, double **ain, int n);
	void copy(double **aout, double **ain, int n);
	void covsrt(float **covar, int ma, int ia[], int mfit);
	void covsrt(float **covar, int ma, int ia[], int mfit);
	void crank(unsigned long n, float w[], float *s);
	void derivs(float x, float y[], float dydx[]);
	void derivs(float x, float y[], float dydx[]);
	void dfour1(double data[], unsigned long nn, int isign);
	void dftcor(float w, float delta, float a, float b, float endpts[],
	void difeq(int k, int k1, int k2, int jsf, int is1, int isf,
	void drealft(double data[], unsigned long n, int isign);
	void dsprsax(double sa[], unsigned long ija[], double x[], double b[],
	void dsprstx(double sa[], unsigned long ija[], double x[], double b[],
	void dsvbksb(double **u, double w[], double **v, int m, int n, double b[],
	void dsvdcmp(double **a, int m, int n, double w[], double **v);
	void eigsrt(float d[], float **v, int n);
	void exit();
	void fdjac(int n, float x[], float fvec[], float **df,
	void fdjac(int n, float x[], float fvec[], float **df,
	void fill0(double **u, int n);
	void fill0(double **u, int n);
	void fill0(double **u, int n);
	void fit(float x[], float y[], int ndata, float sig[], int mwt,
	void flmoon(int n, int nph, long *jd, float *frac);
	void four1(float data[], unsigned long nn, int isign);
	void four1(float data[], unsigned long nn, int isign);
	void four1(float data[], unsigned long nn, int isign);
	void fourew(FILE *file[5], int *na, int *nb, int *nc, int *nd);
	void fourn(float data[], unsigned long nn[], int ndim, int isign);
	void gauleg(float x1, float x2, float x[], float w[], int n);
	void gaussj(float **a, int n, float **b, int m);
	void gaussj(float **a, int n, float **b, int m);
	void gcf(float *gammcf, float a, float x, float *gln);
	void gcf(float *gammcf, float a, float x, float *gln);
	void gser(float *gamser, float a, float x, float *gln);
	void gser(float *gamser, float a, float x, float *gln);
	void hqr(float **a, int n, float wr[], float wi[]);
	void hufapp(unsigned long index[], unsigned long nprob[], unsigned long n,
	void hypdrv(float s, float yy[], float dyyds[]);
	void hypser(fcomplex a, fcomplex b, fcomplex c, fcomplex z,
	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);
	void interp(double **uf, double **uc, int nf);
	void interp(double **uf, double **uc, int nf);
	void jacobn(float x, float y[], float dfdx[], float **dfdy, int n);
	void jacobn(float x, float y[], float dfdx[], float **dfdy, int n);
	void kermom(double w[], double y, int m);
	void laguer(fcomplex a[], int m, fcomplex *x, int *its);
	void linmin(float p[], float xi[], int n, float *fret,
	void linmin(float p[], float xi[], int n, float *fret,
	void lnsrch(int n, float xold[], float fold, float g[], float p[], float x[],
	void lnsrch(int n, float xold[], float fold, float g[], float p[], float x[],
	void lnsrch(int n, float xold[], float fold, float g[], float p[], float x[],
	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 lop(double **out, double **u, int n);
	void lubksb(float **a, int n, int *indx, float b[]);
	void lubksb(float **a, int n, int *indx, float b[]);
	void lubksb(float **a, int n, int *indx, float b[]);
	void lubksb(float **a, int n, int *indx, float b[]);
	void lubksb(float **a, int n, int *indx, float b[]);
	void lubksb(float **a, int n, int *indx, float b[]);
	void lubksb(float **a, int n, int *indx, float b[]);
	void lubksb(float **a, int n, int *indx, float b[]);
	void lubksb(float **a, int n, int *indx, float b[]);
	void lubksb(float **a, int n, int *indx, float b[]);
	void ludcmp(float **a, int n, int *indx, float *d);
	void ludcmp(float **a, int n, int *indx, float *d);
	void ludcmp(float **a, int n, int *indx, float *d);
	void ludcmp(float **a, int n, int *indx, float *d);
	void ludcmp(float **a, int n, int *indx, float *d);
	void ludcmp(float **a, int n, int *indx, float *d);
	void ludcmp(float **a, int n, int *indx, float *d);
	void ludcmp(float **a, int n, int *indx, float *d);
	void ludcmp(float **a, int n, int *indx, float *d);
	void matadd(double **a, double **b, double **c, int n);
	void matsub(double **a, double **b, double **c, int n);
	void mmid(float y[], float dydx[], int nvar, float xs, float htot,
	void mnbrak(float *ax, float *bx, float *cx, float *fa, float *fb,
	void mnbrak(float *ax, float *bx, float *cx, float *fa, float *fb,
	void mnbrak(float *ax, float *bx, float *cx, float *fa, float *fb,
	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 mpinv(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 mplsh(unsigned char u[], int n);
	void mplsh(unsigned char u[], int n);
	void mpmov(unsigned char u[], unsigned char v[], int n);
	void mpmov(unsigned char u[], unsigned char v[], int n);
	void mpmov(unsigned char u[], unsigned char v[], 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,
	void mpmul(unsigned char w[], unsigned char u[], unsigned char v[], int n,
	void mpmul(unsigned char w[], unsigned char u[], unsigned char v[], int n,
	void mpmul(unsigned char w[], unsigned char u[], unsigned char v[], int n,
	void mpneg(unsigned char u[], int n);
	void mpneg(unsigned char u[], int n);
	void mprove(float **a, float **alud, int n, int indx[], float b[],
	void mpsdv(unsigned char w[], unsigned char u[], int n, int iv, int *ir);
	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,
	void mpsub(int *is, unsigned char w[], unsigned char u[], unsigned char v[],
	void mrqcof(float x[], float y[], float sig[], int ndata, float a[],
	void newt(float x[], int n, int *check,
	void newt(float x[], int n, int *check,
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void nrerror(char error_text[]);
	void odeint(float ystart[], int nvar, float x1, float x2,
	void odeint(float ystart[], int nvar, float x1, float x2,
	void odeint(float ystart[], int nvar, float x1, float x2,
	void pearsn(float x[], float y[], unsigned long n, float *r, float *prob,
	void pearsn(float x[], float y[], unsigned long n, float *r, float *prob,
	void pinvs(int ie1, int ie2, int je1, int jsf, int jc1, int k,
	void poldiv(float u[], int n, float v[], int nv, float q[], float r[]);
	void polint(float xa[], float ya[], int n, float x, float *y, float *dy);
	void polint(float xa[], float ya[], int n, float x, float *y, float *dy);
	void polint(float xa[], float ya[], int n, float x, float *y, float *dy);
	void polint(float xa[], float ya[], int n, float x, float *y, float *dy);
	void polint(float xa[], float ya[], int n, float x, float *y, float *dy);
	void psdes(unsigned long *lword, unsigned long *irword);
	void pzextr(int iest, float xest, float yest[], float yz[], float dy[],
	void pzextr(int iest, float xest, float yest[], float yz[], float dy[],
	void qrdcmp(float **a, int n, float *c, float *d, int *sing);
	void qrupdt(float **r, float **qt, int n, float u[], float v[]);
	void quadct(float x, float y, float xx[], float yy[], unsigned long nn,
	void quadct(float x, float y, float xx[], float yy[], unsigned long nn,
	void quadmx(float **a, int n);
	void ranpt(float pt[], float regn[], int n);
	void realft(float data[], unsigned long n, int isign);
	void realft(float data[], unsigned long n, int isign);
	void realft(float data[], unsigned long n, int isign);
	void realft(float data[], unsigned long n, int isign);
	void realft(float data[], unsigned long n, int isign);
	void realft(float data[], unsigned long n, int isign);
	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,
	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);
	void reverse(int iorder[], int ncity, int n[]);
	void rk4(float y[], float dydx[], int n, float x, float h, float yout[],
	void rkck(float y[], float dydx[], int n, float x, float h,
	void rkqs(float y[], float dydx[], int n, float *x,
	void rkqs(float y[], float dydx[], int n, float *x,
	void rlft3(float ***data, float **speq, unsigned long nn1,
	void rlft3(float ***data, float **speq, unsigned long nn1,
	void rlft3(float ***data, float **speq, unsigned long nn1,
	void rotate(float **r, float **qt, int n, int i, float a, float b);
	void rsolv(float **a, int n, float d[], float b[]);
	void rsolv(float **a, int n, float d[], float b[]);
	void rstrct(double **uc, double **uf, int nc);
	void rstrct(double **uc, double **uf, int nc);
	void score(float xf, float y[], float f[]);
	void score(float xf, float y[], float f[]);
	void shell(unsigned long n, float a[]);
	void shoot(int n, float v[], float f[]);
	void shootf(int n, float v[], float f[]);
	void simp1(float **a, int mm, int ll[], int nll, int iabf, int *kp,
	void simp2(float **a, int n, int l2[], int nl2, int *ip, int kp, float *q1);
	void simp3(float **a, int i1, int k1, int ip, int kp);
	void simpr(float y[], float dydx[], float dfdx[], float **dfdy,
	void slvsm2(double **u, double **rhs);
	void slvsml(double **u, double **rhs);
	void solvde(int itmax, float conv, float slowc, float scalv[],
	void sort(unsigned long n, float arr[]);
	void sort(unsigned long n, float arr[]);
	void sort(unsigned long n, float arr[]);
	void sort2(unsigned long n, float arr[], float brr[]);
	void spline(float x[], float y[], int n, float yp1, float ypn, float y2[]);
	void spline(float x[], float y[], int n, float yp1, float ypn, float y2[]);
	void splint(float xa[], float ya[], float y2a[], int n, float x, float *y);
	void spread(float y, float yy[], unsigned long n, float x, int m);
	void svbksb(float **u, float w[], float **v, int m, int n, float b[],
	void svdcmp(float **a, int m, int n, float w[], float **v);
	void tqli(float d[], float e[], int n, float **z);
	void tridag(float a[], float b[], float c[], float r[], float u[],
	void trnspt(int iorder[], int ncity, int n[]);
	void twofft(float data1[], float data2[], float fft1[], float fft2[],
	void twofft(float data1[], float data2[], float fft1[], float fft2[],
	void wwghts(float wghts[], int n, float h,
	void zroots(fcomplex a[], int m, fcomplex roots[], int polish);
extern void (*nrdfun)(float [], float []);
extern void (*nrfuncv)(int n, float v[], float f[]);
float rtnewt(void (*funcd)(float, float *, float *), float x1, float x2,
float rtsafe(void (*funcd)(float, float *, float *), float x1, float x2,
int main(void)		/* example1 */
int main(void)		/* example2 */
int main(void)		/* example3 */
int main(void)  /* Program badluk */
int main(void)  /* Program fredex */
int main(void)  /* Program sfroid */
int main(void)  /* Program sphfpt */
int main(void)  /* Program sphoot */
void (*nrdfun)(float [], float []);
void (*nrfuncv)(int n, float v[], float f[]);
void (*nrfuncv)(int n, float v[], float f[]);
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,
void amoeba(float **p, float y[], int ndim, float ftol,
void anneal(float x[], float y[], int iorder[], int ncity)
void arcmak(unsigned long nfreq[], unsigned long nchh, unsigned long nradd,
void arcode(unsigned long *ich, unsigned char **codep, unsigned long *lcode,
void arcsum(unsigned long iin[], unsigned long iout[], unsigned long ja,
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,
void bandec(float **a, unsigned long n, int m1, int m2, float **al,
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,
void bcuint(float y[], float y1[], float y2[], float y12[], float x1l,
void beschb(double x, double *gam1, double *gam2, double *gampl, double *gammi)
void bessik(float x, float xnu, float *ri, float *rk, float *rip, float *rkp)
void bessjy(float x, float xnu, float *rj, float *ry, float *rjp, float *ryp)
void bksub(int ne, int nb, int jf, int k1, int k2, float ***c)
void broydn(float x[], int n, int *check,
void bsstep(float y[], float dydx[], int nv, float *xx, float htry, float eps,
void caldat(long julian, int *mm, int *id, int *iyyy)
void chder(float a, float b, float c[], float cder[], int n)
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)
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,
void chstwo(float bins1[], float bins2[], int nbins, int knstrn, float *df,
void cisi(float x, float *ci, float *si)
void cntab1(int **nn, int ni, int nj, float *chisq, float *df, float *prob,
void cntab2(int **nn, int ni, int nj, float *h, float *hx, float *hy,
void convlv(float data[], unsigned long n, float respns[], unsigned long m,
void copy(double **aout, double **ain, int n)
void correl(float data1[], float data2[], unsigned long n, float ans[])
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,
void daub4(float a[], unsigned long n, int isign)
void ddpoly(float c[], int nc, float x, float pd[], int nd)
void derivs(float x, float y[], float dydx[])
void derivs(float x, float y[], float dydx[])
void dfour1(double data[], unsigned long nn, int isign)
void dfpmin(float p[], int n, float gtol, int *iter, float *fret,
void dftcor(float w, float delta, float a, float b, float endpts[],
void dftint(float (*func)(float), float a, float b, float w, float *cosint,
void difeq(int k, int k1, int k2, int jsf, int is1, int isf, int indexv[],
void dlinmin(float p[], float xi[], int n, float *fret, float (*func)(float []),
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))
void eigsrt(float d[], float **v, int n)
void elmhes(float **a, int n)
void eulsum(float *sum, float term, int jterm, float wksp[])
void fasper(float x[], float y[], unsigned long n, float ofac, float hifac,
void fdjac(int n, float x[], float fvec[], float **df,
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,
void fitexy(float x[], float y[], int ndat, float sigx[], float sigy[],
void fixrts(float d[], int m)
void fleg(float x, float pl[], int nl)
void flmoon(int n, int nph, long *jd, float *frac)
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[],
void free_convert_matrix(b,nrl,nrh,ncl,nch)
void free_convert_matrix(float **b, long nrl, long nrh, long ncl, long nch)
void free_cvector(unsigned char *v, long nl, long nh)
void free_cvector(v,nl,nh)
void free_dmatrix(double **m, long nrl, long nrh, long ncl, long nch)
void free_dmatrix(m,nrl,nrh,ncl,nch)
void free_dvector(double *v, long nl, long nh)
void free_dvector(v,nl,nh)
void free_f3tensor(float ***t, long nrl, long nrh, long ncl, long nch,
void free_f3tensor(t,nrl,nrh,ncl,nch,ndl,ndh)
void free_imatrix(int **m, long nrl, long nrh, long ncl, long nch)
void free_imatrix(m,nrl,nrh,ncl,nch)
void free_ivector(int *v, long nl, long nh)
void free_ivector(v,nl,nh)
void free_lvector(unsigned long *v, long nl, long nh)
void free_lvector(v,nl,nh)
void free_matrix(float **m, long nrl, long nrh, long ncl, long nch)
void free_matrix(m,nrl,nrh,ncl,nch)
void free_submatrix(b,nrl,nrh,ncl,nch)
void free_submatrix(float **b, long nrl, long nrh, long ncl, long nch)
void free_vector(float *v, long nl, long nh)
void free_vector(v,nl,nh)
void frenel(float x, float *s, float *c)
void frprmn(float p[], int n, float ftol, int *iter, float *fret,
void ftest(float data1[], unsigned long n1, float data2[], unsigned long n2,
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)
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,
void hufdec(unsigned long *ich, unsigned char *code, unsigned long lcode,
void hufenc(unsigned long ich, unsigned char **codep, unsigned long *lcode,
void hufmak(unsigned long nfreq[], unsigned long nchin, unsigned long *ilong,
void hunt(float xx[], unsigned long n, float x, unsigned long *jlo)
void hypdrv(float s, float yy[], float dyyds[])
void hypser(fcomplex a, fcomplex b, fcomplex c, fcomplex z, fcomplex *series,
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)
void jacobi(float **a, int n, float d[], float **v, int *nrot)
void jacobn(float x, float y[], float dfdx[], float **dfdy, int n)
void kendl1(float data1[], float data2[], unsigned long n, float *tau,
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 ks2d2s(float x1[], float y1[], unsigned long n1, float x2[], float y2[],
void ksone(float data[], unsigned long n, float (*func)(float), float *d,
void kstwo(float data1[], unsigned long n1, float data2[], unsigned long n2,
void laguer(fcomplex a[], int m, fcomplex *x, int *its)
void lfit(float x[], float y[], float sig[], int ndat, float a[], int ia[],
void linbcg(unsigned long n, double b[], double x[], int itol, double tol,
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[],
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)
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,
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[])
void mgfas(double **u, int n, int maxcyc)
void mglin(double **u, int n, int ncycle)
void miser(float (*func)(float []), float regn[], int ndim, unsigned long npts,
void mmid(float y[], float dydx[], int nvar, float xs, float htot, int nstep,
void mnbrak(float *ax, float *bx, float *cx, float *fa, float *fb, float *fc,
void mnewt(int ntrial, float x[], int n, float tolx, float tolf)
void moment(float data[], int n, float *ave, float *adev, float *sdev,
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[],
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,
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,
void mpsub(int *is, unsigned char w[], unsigned char u[], unsigned char v[],
void mrqcof(float x[], float y[], float sig[], int ndata, float a[], int ia[],
void mrqmin(float x[], float y[], float sig[], int ndata, float a[], int ia[],
void newt(float x[], int n, int *check,
void nrerror(char error_text[])
void nrerror(error_text)
void odeint(float ystart[], int nvar, float x1, float x2, float eps, float h1,
void orthog(int n, float anu[], float alpha[], float beta[], float a[],
void pade(double cof[], int n, float *resid)
void pccheb(float d[], float c[], int n)
void pcshft(float a, float b, float d[], int n)
void pearsn(float x[], float y[], unsigned long n, float *r, float *prob,
void period(float x[], float y[], int n, float ofac, float hifac, float px[],
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,
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,
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,
void predic(float data[], int ndata, float d[], int m, float future[],
void psdes(unsigned long *lword, unsigned long *irword)
void pwt(float a[], unsigned long n, int isign)
void pwtset(int n)
void pzextr(int iest, float xest, float yest[], float yz[], float dy[], int nv)
void qrdcmp(float **a, int n, float *c, float *d, int *sing)
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[])
void quadct(float x, float y, float xx[], float yy[], unsigned long nn,
void quadmx(float **a, int n)
void quadvl(float x, float y, float *fa, float *fb, float *fc, float *fd)
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,
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,
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)
void reverse(int iorder[], int ncity, int n[])
void rk4(float y[], float dydx[], int n, float x, float h, float yout[],
void rkck(float y[], float dydx[], int n, float x, float h, float yout[],
void rkdumb(float vstart[], int nvar, float x1, float x2, int nstep,
void rkqs(float y[], float dydx[], int n, float *x, float htry, float eps,
void rlft3(float ***data, float **speq, unsigned long nn1, unsigned long nn2,
void rotate(float **r, float **qt, int n, int i, float a, float b)
void rsolv(float **a, int n, float d[], float b[])
void rstrct(double **uc, double **uf, int nc)
void rzextr(int iest, float xest, float yest[], float yz[], float dy[], int nv)
void savgol(float c[], int np, int nl, int nr, int ld, int m)
void score(float xf, float y[], float f[])
void score(float xf, float y[], float f[])
void scrsho(float (*fx)(float))
void shell(unsigned long n, float a[])
void shoot(int n, float v[], float f[])
void shootf(int n, float v[], float f[])
void simp1(float **a, int mm, int ll[], int nll, int iabf, int *kp,
void simp2(float **a, int n, int l2[], int nl2, int *ip, int kp, float *q1)
void simp3(float **a, int i1, int k1, int ip, int kp)
void simplx(float **a, int m, int n, int m1, int m2, int m3, int *icase,
void simpr(float y[], float dydx[], float dfdx[], float **dfdy, int n,
void sinft(float y[], int n)
void slvsm2(double **u, double **rhs)
void slvsml(double **u, double **rhs)
void sncndn(float uu, float emmc, float *sn, float *cn, float *dn)
void sobseq(int *n, float x[])
void solvde(int itmax, float conv, float slowc, float scalv[], int indexv[],
void sor(double **a, double **b, double **c, double **d, double **e,
void sort(unsigned long n, float arr[])
void sort2(unsigned long n, float arr[], float brr[])
void sort3(unsigned long n, float ra[], float rb[], float rc[])
void spctrm(FILE *fp, float p[], int m, int k, int ovrlap)
void spear(float data1[], float data2[], unsigned long n, float *d, float *zd,
void sphbes(int n, float x, float *sj, float *sy, float *sjp, float *syp)
void splie2(float x1a[], float x2a[], float **ya, int m, int n, float **y2a)
void splin2(float x1a[], float x2a[], float **ya, float **y2a, int m, int n,
void spline(float x[], float y[], int n, float yp1, float ypn, float y2[])
void splint(float xa[], float ya[], float y2a[], int n, float x, float *y)
void spread(float y, float yy[], unsigned long n, float x, int m)
void sprsax(float sa[], unsigned long ija[], float x[], float b[],
void sprsin(float **a, int n, float thresh, unsigned long nmax, float sa[],
void sprspm(float sa[], unsigned long ija[], float sb[], unsigned long ijb[],
void sprstm(float sa[], unsigned long ija[], float sb[], unsigned long ijb[],
void sprstp(float sa[], unsigned long ija[], float sb[], unsigned long ijb[])
void sprstx(float sa[], unsigned long ija[], float x[], float b[],
void stifbs(float y[], float dydx[], int nv, float *xx, float htry, float eps,
void stiff(float y[], float dydx[], int n, float *x, float htry, float eps,
void stoerm(float y[], float d2y[], int nv, float xs, float htot, int nstep,
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,
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,
void tqli(float d[], float e[], int n, float **z)
void tred2(float **a, int n, float d[], float e[])
void tridag(float a[], float b[], float c[], float r[], float u[],
void trnspt(int iorder[], int ncity, int n[])
void ttest(float data1[], unsigned long n1, float data2[], unsigned long n2,
void tutest(float data1[], unsigned long n1, float data2[], unsigned long n2,
void twofft(float data1[], float data2[], float fft1[], float fft2[],
void usrfun(float *x,int n,float *fvec,float **fjac);
void vander(double x[], double w[], double q[], int n)
void vegas(float regn[], int ndim, float (*fxn)(float [], float), int init,
void voltra(int n, int m, float t0, float h, float *t, float **f,
void wt1(float a[], unsigned long n, int isign,
void wtn(float a[], unsigned long nn[], int ndim, int isign,
void wwghts(float wghts[], int n, float h,
void zbrak(float (*fx)(float), float x1, float x2, int n, float xb1[],
void zrhqr(float a[], int m, float rtr[], float rti[])
void zroots(fcomplex a[], int m, fcomplex roots[], int polish)
