hw4_gibbs.m