integrate_gk function

Adaptive Gauss-Kronrod Quadrature for multiple limits