#1902: implemented neural network covariance function plus test case (comparison with GPML) for Gaussian processes