STK: a Small (Matlab/Octave) Toolbox for Kriging
 STK_EXAMPLE_KB03  Ordinary kriging in 2D

 An anisotropic Matern covariance function is used for the Gaussian Process
 (GP) prior. The parameters of this covariance function (variance, regularity
 and ranges) are estimated using the Restricted Maximum Likelihood (ReML)
 method.

 The mean function of the GP prior is assumed to be constant and unknown. This
 default choice can be overridden by means of the model.lm property.

 The function is sampled on a space-filling Latin Hypercube design, and the
 data is assumed to be noiseless.