# NOT RUN {
# Load FracKrigingR library
library(FracKrigingR)
#Fractional Brownian vector field
K = 10; H = 0.5; p = 2
#Generate coordinates
X<-matrix(0,ncol=p, nrow=K)
for(j in 1:p){
for(i in 1:K){
X[i,j] = rnorm(1, 0, 1)
}
}
FracMatrix(H, K, X)
# }
Run the code above in your browser using DataLab