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