# \donttest{
#Examples for muPEend.int
muPEend.int(1.2)
rseq<-seq(1.01,5,by=.1)
lrseq<-length(rseq)
mu.end<-vector()
for (i in 1:lrseq)
{
mu.end<-c(mu.end,muPEend.int(rseq[i]))
}
plot(rseq, mu.end,type="l",
ylab=expression(paste(mu,"(r)")),xlab="r",lty=1,xlim=range(rseq),ylim=c(0,1))
# }
# \donttest{
#Examples for asy.varPEend.int
asy.varPEend.int(1.2)
rseq<-seq(1.01,5,by=.1)
lrseq<-length(rseq)
var.end<-vector()
for (i in 1:lrseq)
{
var.end<-c(var.end,asy.varPEend.int(rseq[i]))
}
oldpar <- par(mar=c(5,5,4,2))
plot(rseq, var.end,type="l",
xlab="r",ylab=expression(paste(sigma^2,"(r)")),lty=1,xlim=range(rseq))
par(oldpar)
# }
Run the code above in your browser using DataLab