We explore the combination of deterministic and Monte Carlo methods to
facilitate efficient automatic numerical computation of multidimensional
integrals with singular integrands. Two adaptive algorithms are presented that
employ recursion and are runtime and memory optimised, respectively. SINGINT, a
C implementation of the algorithms, is introduced and its utilisation in the
calculation of particle scattering amplitudes is exemplified.Comment: 15 pages, 1 figur