BivPoiss function

Bivariate Poisson distribution