/*****************************************************************************/
/*                             FILE DCDFLIB.H                                */
/*****************************************************************************/

double algdiv (double *a,double *b);
double alnrel (double *a);
double apser  (double *a,double *b,double *x,double *eps);
double basym  (double *a,double *b,double *lambda,double *eps);
double bcorr  (double *a0,double *b0);
double betaln (double *a0,double *b0);
double bfrac  (double *a,double *b,double *x,double *y,double *lambda, double *eps);
void   bgrat  (double *a,double *b,double *x,double *y,double *w,double *eps,int *ierr);
double bpser  (double *a,double *b,double *x,double *eps);
double brcmp1 (int *mu,double *a,double *b,double *x,double *y);
double brcomp (double *a,double *b,double *x,double *y);
double bup    (double *a,double *b,double *x,double *y,int *n,double *eps);
void   cdff   (int *which,double *p,double *q,double *f,double *dfn,double *dfd,int *status,double *bound);
void   cumf   (double *f,double *dfn,double *dfd,double *cum,double *ccum);
void   dinvr  (int *status,double *x,double *fx,unsigned long *qleft,unsigned long *qhi);
void   dstinv (double *zsmall,double *zbig,double *zabsst,double *zrelst,double *zstpmu,double *zabsto,double *zrelto);
void   dzror  (int *status,double *x,double *fx,double *xlo,double *xhi,unsigned long *qleft,unsigned long *qhi);
void   dstzr  (double *zxlo,double *zxhi,double *zabstl,double *zreltl);
double erf1   (double *x);
double erfc1  (int *ind,double *x);
double esum   (int *mu,double *x);
double exparg (int *l);
double fpser  (double *a,double *b,double *x,double *eps);
double gam1   (double *a);
double gamln  (double *a);
double gamln1 (double *a);
void   grat1  (double *a,double *x,double *r,double *p,double *q,double *eps);
double gsumln (double *a,double *b);
double psi    (double *xx);
double rexp   (double *x);
double rlog1  (double *x);
double spmpar (int *i);

/************************************************************************
FIFDMAX1:
returns the maximum of two numbers a and b
************************************************************************/
double fifdmax1 (double a,double b);
/************************************************************************
FIFDMIN1:
returns the minimum of two numbers a and b
************************************************************************/
double fifdmin1 (double a,double b);
/************************************************************************
FIFDSIGN:
transfers the sign of the variable "sign" to the variable "mag"
************************************************************************/
double fifdsign (double mag,double sign);
/************************************************************************
FIFIDINT:
Truncates a double precision number to a long integer
************************************************************************/
long   fifidint (double a);
/************************************************************************
FTNSTOP:
Prints msg to standard error and then exits
************************************************************************/
void   ftnstop  (char* msg);