dfn = max(
        np.round(2 * dfd * dfd * (dfd - 2) /
                 (v * (dfd - 4) * (dfd - 2) ** 2 - 2 * dfd * dfd)), 1)