adaptive 1d numerical Gauss–Kronrod integration in Julia - View it on GitHub