double alngam ( double xvalue, int *ifault );
void beta_inc_values ( int *n_data, double *a, double *b, double *x, 
  double *fx );
double betain ( double x, double p, double q, double beta, int *ifault );
double r8_abs ( double x );
double r8_max ( double x, double y );
double xinbta ( double p, double q, double beta, double alpha, int *ifault );
void timestamp ( void );