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