double alogam ( double x, int *ifault );
double beta_noncentral_cdf ( double a, double b, double lambda, double x, 
  double error_max );
void beta_noncentral_cdf_values ( int *n_data, double *a, double *b, 
  double *lambda, double *x, double *fx );
double betain ( double x, double p, double q, double beta, int *ifault );
void timestamp ( );