cov_holder function

Calculate the covariance for two time points with sparse sampling