bivariate generating function