# load energies averaged over each run
e=read.csv("eavg.csv")
es<-e+1100
# compute average energy across samples
esm<-rowMeans(es)
esmat<-as.matrix(es)

postscript(file="eavg_acf.eps",horizontal=FALSE,width=5,height=5)

a<-acf(esmat[31:120,1],plot=FALSE)
plot(a,main="")