##### R-code for making Figure 4.4. It reads data from
##### "example41.dat" generated by the R-code "example41.r"

dat <- matrix(scan("example41.dat"),ncol=5,byrow=T)

xbar <- numeric()
nxbar = rep(0,20)

for(i in 1:20){
   xbar[i] <- mean(dat[i,])
   nxbar[i] <- xbar[i]*sqrt(5)
}

k=0.25

cn = rep(0,20)
cn[1] = max(0,nxbar[1]-k)
for(i in 2:20){
   cn[i] <- max(0,cn[i-1]+nxbar[i]-k)
}

ii <- seq(1,20)

postscript("fig44.ps",width=4.5,height=4.5,horizontal=F)

plot(ii,cn,type="o",lty=1,pch=16,xlab="n",
     ylab=expression(C[n]^{"+"}),mgp=c(2,1,0),xlim=c(0,21), 
     ylim=c(0,6.5),cex=0.8)
lines(ii,rep(5.597,20),lty=2,cex=0.8)

graphics.off()