Module DQC25F in CMLIBClenshaw-Curtis integration rule for function with cos or sin factor, also uses Gauss Kronrod formula.
CALL DQC25F(F, A, B, OMEGA, INTEGR, NRMOM, MAXP1, KSAVE, RESULT, ABSERR, NEVAL, RESABS, RESASC, MOMCOM, CHEBMO)| H2a2a2 | Nonautomatic 1-D finite interval quadrature (special integrand including weight functions, oscillating and singular integrands, principal value integrals, splines, etc.), integrand available via user-defined procedure |