generate random number with GStools
This is not possible with GSTools, since your covariance model is unbound, which refers to an infinite variance of the random field. GSTools only supports random fields with finite mean and variance.
How can I generate random field with GStools for my correlation function: c(x)= (x)^(a-2)/(x^2 + 1)^(a/2) for differences values x and a=0.5?