Add adaptive Gauss–Kronrod quadrature

[axkr]

The "Laguerre" and "Hermite" implementations of the GaussIntegratorFactory, doesn't seem to have "flexible changeable Infinity bounds", if used in a scripting language.

The below Java implementation is similar to what Julia supports in it's quadgk function:

• https://github.com/JuliaMath/QuadGK.jl
• https://juliamath.github.io/QuadGK.jl/stable/
• https://en.wikipedia.org/wiki/Gauss%E2%80%93Kronrod_quadrature_formula

I found this Java implementation (Apache license) on Github which accepts Infinity bounds.
Example:

• https://github.com/jonathanschilling/AdaptiveQuadrature/blob/b51098476a528a612abe34a358d0473757300874/src/test/java/de/labathome/tests/TestAdaptiveQuadrature.java#L60

It could also be implemented for "arbitrary-precision arithmetic" (i.e. FieldElement).