Add Gauss Kronrod integrator

Very basic, stolen from my light bending code