float gammq(float a, float x);
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 gaussj(float **a, int n, float **b, int m);
void covsrt(float **covar, int ma, int ia[], int mfit);
void fpoly(float x, float p[], int np);
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 gser(float *gamser, float a, float x, float *gln);
void gcf(float *gammcf, float a, float x, float *gln);
float gammln(float xx);
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 spldrv(float xa[], float ya[], float y2a[], int n, float x, float *dy);
float deg_to_rad(float angle);