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