##### R-code for making Figure 4.11

x=c(8,  17, 10,  8, 11, 14,  8, 16, 8, 13,   
    21, 15, 20, 18, 17)

k=14.524
h=4.0

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

ii <- seq(1,length(x))

postscript("fig411.ps",width=4.5,height=4.5,horizontal=FALSE)

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

graphics.off()