x = 0.1;
xp = 0.2;
y=covr(x,xp);

fprintf(1,'y = %.10f\n',y)
