— Function: double gsl_sf_erfc (double x)
— Function: int gsl_sf_erfc_e (double x, gsl_sf_result * result)

These routines compute the complementary error function erfc(x) = 1 - erf(x) = (2/\sqrt(\pi)) \int_x^\infty \exp(-t^2).