fix a bug in covDiff.m

Kss = sigma_f^2 * ones(n+nd, 1);