def units_covariance(N_alpha, N_beta):
    """Compute the covariance of E_{lpha eta}(s) and E_{lpha eta}(t)"""
    if N_alpha == 0 or N_beta == 0:
        return float('nan')
    else:
        return (1 - 1 / N_alpha) * (1 - 1 / N_beta) - 1
