source("gibbs.R")

windows(width=10,height=7)
par(ask=TRUE)

gibbs.plot(.5,100)
gibbs.plot(.9,200)

par( mfrow=c(1,2),mar=c(3.6,3.6,0.1,0.1),lwd=2)
mhlocal(1000,20)
mhlocal(10000,5)