import numpy as np

def Covariance(x, y):
    xbar, ybar = x.mean(), y.mean()
    return np.sum((x - xbar)*(y - ybar))/(len(x) - 1)

variance-covariance example from: https://stattrek.com/matrix-algebra/covariance-matrix.aspx