sample.n.fit function

Sample from a Gaussian process and fit an emulator to the points