train_gp function

Learning hyper-parameters of a Gaussian Process