def covariance(x, y):
    n = len(x)
    return Ch4.dot(de_mean(x), de_mean(y))/(n-1)