— Function: double gsl_sf_psi_n (int m, double x)
— Function: int gsl_sf_psi_n_e (int m, double x, gsl_sf_result * result)

These routines compute the polygamma function \psi^{(m)}(x) for m >= 0, x > 0.