double alngam ( double xvalue, int &ifault );
void gamma_log_values ( int &n_data, double &x, double &fx );
double lngamma ( double z, int &ier );