#####   R-code for making Figure 5.7

n = seq(1,100)

w005 = rep(0,100)
w01  = rep(0,100)
w02  = rep(0,100)
w05  = rep(0,100)

for(i in 1:100){
   w005[i] = 0.05*(1-0.05)^(100-i)
   w01[i]  = 0.1*(1-0.1)^(100-i)
   w02[i]  = 0.2*(1-0.2)^(100-i)
   w05[i]  = 0.5*(1-0.5)^(100-i)
}

postscript("fig57.ps",width=6.5,height=6.5,horizontal=F)

par(mfrow=c(2,2), mar=c(4,4,2,2))

plot(n,w005,type="l",lty=1,pch=16,xlab="i",
     ylab=expression(w[i]),mgp=c(2,1,0),xlim=c(0,100), 
     ylim=c(0,0.05), cex=0.8)
title(xlab="(a)",cex=0.9)

plot(n,w01,type="l",lty=1,pch=16,xlab="i",
     ylab=expression(w[i]),mgp=c(2,1,0),xlim=c(0,100), 
     ylim=c(0,0.1), cex=0.8)
title(xlab="(b)",cex=0.9)

plot(n,w02,type="l",lty=1,pch=16,xlab="i",
     ylab=expression(w[i]),mgp=c(2,1,0),xlim=c(0,100), 
     ylim=c(0,0.2), cex=0.8)
title(xlab="(c)",cex=0.9)

plot(n,w05,type="l",lty=1,pch=16,xlab="i",
     ylab=expression(w[i]),mgp=c(2,1,0),xlim=c(0,100), 
     ylim=c(0,0.5), cex=0.8)
title(xlab="(d)",cex=0.9)

graphics.off()