#Examples for muCS2D
muCS2D(.5)
tseq<-seq(0.01,5,by=.05)
ltseq<-length(tseq)
mu<-vector()
for (i in 1:ltseq)
{
mu<-c(mu,muCS2D(tseq[i]))
}
plot(tseq, mu,type="l",xlab="t",ylab=expression(mu(t)),lty=1,xlim=range(tseq))
#Examples for asyvarCS2D
asyvarCS2D(.5)
tseq<-seq(0.01,10,by=.05)
ltseq<-length(tseq)
asyvar<-vector()
for (i in 1:ltseq)
{
asyvar<-c(asyvar,asyvarCS2D(tseq[i]))
}
oldpar <- par(no.readonly = TRUE)
par(mar=c(5,5,4,2))
plot(tseq, asyvar,type="l",xlab="t",ylab=expression(paste(sigma^2,"(t)")),lty=1,xlim=range(tseq))
par(oldpar)
Run the code above in your browser using DataLab