double trigamma ( double x, int *ifault );
void trigamma_values ( int *n_data, double *x, double *fx );