adding in new integrator for knm using higher order newton-cotes methods and analytic solution for square aperture