  A procedure to include the uncertainty on the background estimate for upper
limit calculations using Poissonian sampling is presented for the case where a
Gaussian assumption on the uncertainty can be made. Under that hypothesis an
analytic expression of the likelihood is derived which can be written in terms
of polynomials defined by recursion. This expression may lead to a significant
speed up of computing applications that extract the upper limits using Toy
Monte Carlo.
