covariance model. currently implemented:
matern: with covparms (var,range,smoothness)
esqe: exponential + squared exp with covparms (var1,range1,var2,range2)
If covmodel is a function it has to be able to take a distance matrix
and return a vector with distances which is of length k.