gibbs_sampling function

Markov chain Monte Carlo simulation for the probit model