Compute uma função de distribuição para racionais no intervalo unitário 

Ordene os números racionais entre zero e um por sua representação binária.

mostre o input completo de Wolfram Languageoculte o input

In[1]:=

X RankedRational[n_Integer?Positive] := Block[{dig = IntegerDigits[n, 2], avec, cf}, avec = Reverse[Pick[Range[Length[dig] - 1, 0, -1], dig, 1]]; cf = ConstantArray[0, Length[avec] + 1]; cf[[3 ;; -1]] = Differences[avec]; cf[[-1]] += 1; cf[[2]] = avec[[1]] + 1; cf // FromContinuedFraction ]

In[2]:=

X Table[RankedRational[n], {n, 1, 16}]

Out[2]=

A função de distribuição cumulativa de racionais nesta ordenação é a função MinkowskiQuestionMark.

In[3]:=

X cdf[x_] = CDF[EmpiricalDistribution[Table[RankedRational[n], {n, 1024}]], x];

mostre o input completo de Wolfram Languageoculte o input

In[4]:=

X Plot[{MinkowskiQuestionMark[x], cdf[x] }, {x, 0, 1}, PlotStyle -> {Red, Directive[Thick, Black, Dashing[.04]]}, PlotTheme -> "Detailed", AspectRatio -> 1, ImageSize -> Medium]

Out[4]=