Complete Gibbs sampling part 2
Finished Gibb sampling code
test