We present an algorithm which allows a fast numerical computation of Feldma
n-Cousins confidence intervals for Poisson processes, even when the number
of background events is relatively large. This algorithm incorporates an ap
propriate treatment of the singularities that arise as a consequence of the
discreteness of the variable. (C) 2000 Elsevier Science B.V. All rights re
served.