double alogam ( double x, int *ifault );
void beta_cdf_values ( int *n_data, double *a, double *b, double *x, 
  double *fx );
void gamma_log_values ( int *n_data, double *x, double *fx );
double mdbeta ( double x, double p, double q, int *ier );
void timestamp ( );