RandomFields package on CRAN from ?PrintModelList: * 'fractalB' (fractal Brownian motion) gamma(x) = x^a The parameter a is in (0,2]. (Implemented for up to three dimensions) cheers Ben