Document the code and run 200 iterations of Gibbs (does not work yet)