#**********************Plot Pre Processing***************************************
m<-250		#number of points to plot
step<-round(nsim/m)
tplot<-array(0,dim=c(m,1))
xplot<-array(1,dim=c(m,1))
temp<-cumsum(that)/c(1:nsim)
tplot[1]<-temp[1];xplot<-1;
for(i in 1:m) {tplot[i]<-temp[i*step];xplot[i]<-i*step}
par(mfrow=c(1,2))
plot(xplot,tplot,type="l",xlim=c(0,nsim),ylim=c(3.6,3.8),lwd=2,xlab="",ylab="")
hist(that,col="grey",breaks=25,main="")