Added adaptive step sizing to both Simpson and Boole integration. We can no longer support integrals over infinite regions; we need to modify refinedPoints to keep increasing the region.
Added support for infinite domains, and replaced the various named rules with a single "Integrator" function (currently using Boole's rule).
Fixed an incorrect test case.