I did not write this - only the wrapping for use with fpm. Useful for integration of oscillatory functions. I was looking for a replacement for quadosc in mpmath, after migrating to fortran A test is included to calculate the normal stress in a layer under a circular load - a bessel function.
I’d like to add this method to the modernized quadpack i’m working on. The license is a little iffy though. Maybe we can email him and ask for a more standard BSD/MIT license.