float_inc.jl

function not_inc(x)
    Float64(x) == Float64(x+1)
end

function perc_inc(e, n)
    s = 0
    for i in 0:n
        s += Int(not_inc(2^e + i))
    end
    1-s/n
end

# julia> perc_inc(53, 2^20)
# 0.4999990463256836

# julia> perc_inc(54, 2^20)
# 0.2499990463256836

# julia> perc_inc(55, 2^20)
# 0.1249990463256836