Create gibbs.m