def covariance(points, m):
    """Returns the covariance of a list of complex points"""
    covar = 0.0
    for point in points:
        covar += (point.real - m.real) * (point.imag - m.imag)
    return covar
