##### R-code for making Figure 4.6.

delta=seq(0,2,0.1)
n=length(delta)

ARL025 = rep(0,n)
ARL05 = rep(0,n)
ARL075 = rep(0,n)

for(i in 1:n){
   ARL025[i]=xcusum.arl(k=0.25,h=5.597,mu=delta[i])
   ARL05[i] =xcusum.arl(k=0.5,h=3.502,mu=delta[i])
   ARL075[i]=xcusum.arl(k=0.75,h=2.481,mu=delta[i])
}

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

plot(delta,log(ARL025),type="l",lty=1,xlab=expression(delta),
     yaxt="n",
     ylab=expression(ARL[1]),mgp=c(2,1,0),xlim=c(0,2), 
     ylim=c(0,5.5),cex=0.8)
lines(delta,log(ARL05),lty=2,cex=0.8)
lines(delta,log(ARL075),lty=3,cex=0.8)
axis(2,at=c(log(1),log(5),log(15),log(50),log(200)),
     labels=c("0","5","15","50","200"),cex=0.8)
legend(1.4,5.5,legend=c("k=0.25","k=0.5","k=0.75"),lty=1:3,cex=0.7)

graphics.off()