Update mcmc.m