gibbs_discrete function

Gibbs sampling of a bivariate discrete distribution