The emulation properties of the GP are driven by the choice of the mean and covariance functions. In this work, we standarize the observations \(\mathbf {y}_n\)  such that they have a mean of zero; that is we set \(\mu (\theta ) = 0\)  in the prior assumption. On the other hand, we use a covariance function from the Matern class [1]}, [2]}, [3]}, given by
\(k(\theta , \theta ^{\prime }) = {\Gamma (\nu )}\left({\ell } \right)^\nu K_{\nu } \left({\ell } \right),\) 
