— Function: double gsl_sf_atanint (double x)
— Function: int gsl_sf_atanint_e (double x, gsl_sf_result * result)

These routines compute the Arctangent integral, which is defined as AtanInt(x) = \int_0^x dt \arctan(t)/t.