##### R-code for making Figure 3.2

postscript("fig32.ps",width=4,height=4,horizontal=F)

k=seq(0,3,0.1)

m=5
beta5=pnorm(-k*sqrt(m)+3)-pnorm(-k*sqrt(m)-3)

m=10
beta10=pnorm(-k*sqrt(m)+3)-pnorm(-k*sqrt(m)-3)

m=20
beta20=pnorm(-k*sqrt(m)+3)-pnorm(-k*sqrt(m)-3)

plot(k,beta5,type="l",lty=1,xlab="k",ylab=expression(beta),
     mgp=c(2,1,0),xlim=c(0,3), ylim=c(0,1), cex=0.8)
lines(k,beta10,lty=2,cex=0.8)
lines(k,beta20,lty=3,cex=0.8)
legend(2,1,legend=c("m=5","m=10","m=20"),lty=1:3,cex=0.7)

graphics.off()