  Multivariate Poisson distributions have numerous applications. Fast
computation of these distributions, holding constant a fixed set of linear
combinations of these variables, has been explored by Sontag and Zeilberger.
This elaborates on their work.
