Module DQC25S in CMLIBEstimates integral of function with algebraic-logarithmic singularities using 25 point Clenshaw-Curtis formula and gives error estimate.
CALL DQC25S(F, A, B, BL, BR, ALFA, BETA, RI, RJ, RG, RH, RESULT, ABSERR, RESASC) INTEGR, NEV)| 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 |