|
%I #10 Sep 08 2022 08:45:42
%S 1,4,2,9,4,7,4,4,6,7,7,0,9,1,8,5,8,1,9,0,4,5,5,9,7,0,9,0,1,1,8,6,3,1,
%T 7,6,2,2,1,1,5,2,8,2,2,5,9,0,3,7,2,8,5,3,5,0,1,3,3,3,6,7,0,6,5,7,8,2,
%U 5,3,5,4,2,8,2,8,9,3,0,8,4,8,7,3,5,5,8,3,8,7,9,8,2,6,4,8,8,1,4,7,8,0,2,5,8
%N Decimal expansion of (33+8*sqrt(2))/31.
%C lim_{n -> infinity} b(n)/b(n-1) = (33+8*sqrt(2))/31 for n mod 3 = {1, 2}, b = A118674.
%C lim_{n -> infinity} b(n)/b(n-1) = (33+8*sqrt(2))/31 for n mod 3 = {0, 2}, b = A157646.
%H G. C. Greubel, <a href="/A157647/b157647.txt">Table of n, a(n) for n = 1..10000</a>
%F (33+8*sqrt(2))/31 = (8+sqrt(2))/(8-sqrt(2)).
%e (33+8*sqrt(2))/31 = 1.42947446770918581904...
%p with(MmaTranslator[Mma]): Digits:=100:
%p RealDigits(evalf((33+8*sqrt(2))/31))[1]; # _Muniru A Asiru_, Mar 31 2018
%t RealDigits[(33+8*Sqrt[2])/31,10,120][[1]] (* _Harvey P. Dale_, Mar 02 2014 *)
%o (PARI) (33+8*sqrt(2))/31 \\ _G. C. Greubel_, Mar 30 2018
%o (Magma) (33+8*Sqrt(2))/31; // _G. C. Greubel_, Mar 30 2018
%Y A118674, A157646, A002193 (decimal expansion of sqrt(2)), A157648 (decimal expansion of (1539+850*sqrt(2))/31^2).
%K cons,nonn
%O 1,2
%A _Klaus Brockhaus_, Mar 11 2009
|
