DE quadrature module

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 :slight_smile: A test is included to calculate the normal stress in a layer under a circular load - a bessel function.


See also: David H. Bailey’s routines · Issue #13 · jacobwilliams/quadpack (

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.