ObtainMultBinaryDist function

Generating a multivariate Bernoulli joint-distribution