Compute a Distribution Function for Rationals in the Unit Interval 

Order the rational numbers between zero and one by their binary representation.

show complete Wolfram Language inputhide 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]=

The cumulative distribution functions of rationals in this ordering is the MinkowskiQuestionMark function.

In[3]:=

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

show complete Wolfram Language inputhide 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]=