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 xinbta ( double p, double q, double beta, double alpha, int *ifault );