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

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