Module QDAGS in IMSLSIntegrate a function (which may have endpoint singularities).
CALL QDAGS (F, A, B, ERRABS, ERRREL, RESULT, ERREST)| H2a1a1 | Automatic 1-D finite interval quadrature (user need only specify required accuracy), integrand available via user-defined procedure |