function result = get_wl_ub(N, M, eps)
    k1 = sqrt((M - 1) / N);
    k2 = log(M * N);
    k3 = sqrt(k2 - log(eps));
    k4 = sqrt(log(2 * M) - log(eps));
    
    delta_1 = 24 * k1 * k2 * (2 * k3 + 1);
    delta_2 = 1.5 * k1 * (2 * k4 + 1);

    numerator = N * (1/8 - delta_1 - delta_2) ^ 2;
    denominator = 81 * M * k2 * (1 + 2 * delta_2 / 3);
    result = numerator / denominator;
end