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