if (FALSE) {
#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))
}
if (FALSE) {
#Examples for asyvarPEend.int
asyvarPEend.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,asyvarPEend.int(rseq[i]))
}
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))
}
Run the code above in your browser using DataLab