pred_gp function

Gaussian Process prediction