bm_gaussianprocess function

Fit Gaussian Process models