# NOT RUN {
n=10 #number of observations
Ndisc=1000; t=1/Ndisc #functions represented at grid points 1/Ndisc, 2/Ndisc,...,1
X=matrix(rep(0,Ndisc*n),nrow=n)
for(i in 1:n){
x=rnorm(Ndisc,mean=0,sd=1)
X[i,]=cumsum(x*sqrt(t))
}
InfKern(X,2,2)
# }
Run the code above in your browser using DataLab