1 /* Compare the integral of f(x) with the integral of |f(x)|
2 to determine if f(x) covers both positive and negative values */
5 test_positivity (double result, double resabs);
8 test_positivity (double result, double resabs)
10 int status = (fabs (result) >= (1 - 50 * GSL_DBL_EPSILON) * resabs);