do_MCMC function

do_MCMC conducts an MCMC using Gibbs within Metropolis