##### R-code for making Figure 8.11 about the weight function

pdf("weight.pdf",width=3.5,height=2.8)

par(mar=c(4,4,1,3))

lam=0.2
k=0.5

u=seq(0,1,0.01)

w=numeric()
for(i in 1:101){
   if (u[i] <= k) {w[i]=lam}
   else {w[i]=min(1,lam+(u[i]/k-1))}
}

plot(u,w,type="l",lty=1,xlab="u",ylab=expression(W(paste(u,";",lambda,",",kappa))),
     xlim=c(0,1), ylim=c(0,1), mgp=c(2,1,0), xaxt="n",yaxt="n",cex=1)
axis(1,at=c(0,0.5,0.9),labels=c(0,expression(kappa),expression(paste((2-lambda),kappa))),cex=1.0)
axis(2,at=c(0,0.2,1),labels=c(0,expression(lambda),1),cex=1.0)

graphics.off()