Add oscillatory integral quadratures

[MikaelSlevinsky]

According to R. Piessens and M. Branders, On the computation of Fourier transforms of singular functions, J. Comp. Appl. Math., 43:159-169, 1992, the modified Chebyshev moments of the complex exponential (multiplied by algebraic endpoint singularities) satisfy a many-term recurrence relation that is neither stable in the forward nor backward directions. How to proceed?

Combination of moment-based quadrature with fast orthogonal polynomial transforms to change bases is now irresistible: for example, the modified Legendre moments of the complex exponential are proportional to spherical Bessel functions that satisfy a three-term recurrence relation with well-known stability properties. By changing bases, we would simplify the recurrence and secure the stability by design.

[dlfivefifty]

Should probably go here

https://github.com/JuliaApproximation/OscillatoryIntegrals.jl

[daanhb]

For reference, see also https://github.com/AndrewGibbs/NSDpackage (this is a Matlab package and not related to moment-based quadrature)