add gibbs_sampling.m