# \donttest{
Pdom.num2PEtri(r=1.5)
Pdom.num2PEtri(r=1.4999999999)
Pdom.num2PEtri(r=1.5) / Pdom.num2PEtri(r=1.4999999999)
rseq<-seq(1.01,1.49999999999,l=20) #try also l=100
lrseq<-length(rseq)
pg2<-vector()
for (i in 1:lrseq)
{
pg2<-c(pg2,Pdom.num2PEtri(rseq[i]))
}
plot(rseq, pg2,type="l",xlab="r",
ylab=expression(paste("P(", gamma, "=2)")),
lty=1,xlim=range(rseq)+c(0,.01),ylim=c(0,1))
points(rbind(c(1.50,Pdom.num2PEtri(1.50))),pch=".",cex=3)
# }
Run the code above in your browser using DataLab