PG2PEtri(r=1.5)
PG2PEtri(r=1.4999999999)
PG2PEtri(r=1.5) / PG2PEtri(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,PG2PEtri(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,PG2PEtri(1.50))),pch=".",cex=3)
Run the code above in your browser using DataLab