def mean(values):
    return sum(values) / float(len(values))

def variance(values, mean):
    return sum((x-mean)**2 for x in values)

def covariance(x, mean_x, y, mean_y):
    return sum((x[i]-mean_x)*(y[i]-mean_y) for i in range(len(x)))