double gamain ( double x, double p, int *ifault );
void gamma_inc_p_values ( int &n_data, double &a, double &x, double &fx );