I feel kind of silly to ask this question, but how am I supposed to pass the integrand to the QAWF routine, whose signature is integration_qawf(a::Real, epsabs::Real, limit::Integer)?
Other integration functions have the integrand as their first argument:
integration_cquad(f::Ref{gsl_function}, ...)
integration_qng(f::Ref{gsl_function},...)
michele