# \donttest{
#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))
# }
# \donttest{
#Examples for asy.varPE2D
asy.varPE2D(1.2)
rseq<-seq(1.01,5,by=.05)
lrseq<-length(rseq)
avar<-vector()
for (i in 1:lrseq)
{
avar<-c(avar,asy.varPE2D(rseq[i]))
}
oldpar <- 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