# NOT RUN {
#Examples for muPEend1D
muPEend1D(1.2)
rseq<-seq(1.01,5,by=.05)
lrseq<-length(rseq)
mu.end<-vector()
for (i in 1:lrseq)
{
mu.end<-c(mu.end,muPEend1D(rseq[i]))
}
plot(rseq, mu.end,type="l",
ylab=expression(paste(mu,"(r)")),xlab="r",lty=1,xlim=range(rseq),ylim=c(0,1))
#Examples for asyvarPEend1D
asyvarPEend1D(1.2)
rseq<-seq(1.01,5,by=.05)
lrseq<-length(rseq)
var.end<-vector()
for (i in 1:lrseq)
{
var.end<-c(var.end,asyvarPEend1D(rseq[i]))
}
oldpar <- par(no.readonly = TRUE)
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