##### R-code for making Figure 5.1

n=seq(1,30)

sigma005 = rep(0,length(n))
sigma01  = rep(0,length(n))
sigma02  = rep(0,length(n))

for(i in 1:length(n)){
   sigma005[i] = (0.05/(2-0.05))*(1-(1-0.05)^(2*i))
   sigma01[i]  = (0.1/(2-0.1))*(1-(1-0.1)^(2*i))
   sigma02[i]  = (0.2/(2-0.2))*(1-(1-0.2)^(2*i))
}

postscript("fig51.ps",width=5,height=5,horizontal=F)

plot(n,sigma02,type="l",lty=1,xlab="n",yaxt="n",
     ylab=expression(sigma[E[n]]^2),mgp=c(2,1,0),xlim=c(0,31), 
     ylim=c(0,0.12),cex=0.8)
lines(n,sigma01,lty=2,cex=0.8)
lines(n,sigma005,lty=3,cex=0.8)
axis(2,at=c(0,0.04,0.08,0.12),labels=c("0.00","0.04","0.08","0.12"),cex=0.8)
legend(20,0.1,legend=c("lambda=0.2","lambda=0.1","lambda=0.05"),
       lty=1:3,cex=0.7)

graphics.off()