updated Gibbs_sampling.m