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