%I #16 May 21 2015 09:26:53
%S 0,1,1,48,1,9182736455463727,4,1,2,1,3,1,
%T 16413861141941053151166388889231063606316227031696978138434,9,9,2,1,
%U 3,8,1,1,19,1,1,2,3,1,7,1,1,4,1,1,1,3,3,1,1,2,2,2,1,5,2,1,1,1,1,5,1,78,1,21,1,1,5,3,2
%N Continued fraction expansion of constant A071873.
%C The next term is very large indeed (it has at least 400000 digits).
%C The constant A071873 is the sixth (of 10) decimal selvage numbers x that satisfy: x = Sum_{n>=0} ( floor(10*n*x) (mod 10) ) / 10^n.
%e x = 0.50505050505050505051616161616161616161627272727272...
%t Clear[a]; Do[a[n] = If[OddQ[n], 5, 0], {n, 1, n0=4}]; a[_] = 0; digits = 10^(n0-1); Do[a[n] = Mod[Floor[10*n*Sum[a[k]/10^k, {k, 1, n}]], 10], {n, n0+1, digits}]; x = FromDigits[{Table[a[n], {n, 1, digits}], 0}, 10] // N[#, digits+1]&; ContinuedFraction[x] (* _Jean-François Alcover_, May 12 2015 *)
%Y Cf. A071873.
%K nonn,cofr
%O 0,4
%A _Paul D. Hanna_, May 11 2015
|