Programs for numerical integration of piece-wise smooth functions

A series of fortran programs for the numerical integration of piece-wise smooth functions in the one- and manydimensional spaces is described. The programs use Newton-Cotes quadrature formulas, make accuracy checks at each step, find efficiently the places of degraded smoothness (narrow peaks, splittings, breakpoints), optimize the step size and the degree of quadrature formula, assure the closeness of the actual accuracy of the result the asked one, and have improved reliability