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 );
double betanc ( double x, double a, double b, double lambda, int *ifault );