|
%I #33 Sep 08 2022 08:45:11
%S 1,2,3,6,8,3,9,8,4,4,6,3,8,7,8,5,1,0,1,8,9,0,6,6,0,8,7,6,1,4,2,1,2,3,
%T 2,5,2,2,1,1,1,7,6,6,2,1,2,3,5,8,8,5,8,7,3,7,1,0,7,1,6,7,2,6,7,1,5,9,
%U 0,4,2,7,4,0,0,9,2,5,8,8,1,9,1,0,7,7,8,3,8,2,6,1,3,0,6,3,9,9,3,5,7,5,9,1
%N Decimal expansion of Feller's beta coin-tossing constant.
%H G. C. Greubel, <a href="/A086254/b086254.txt">Table of n, a(n) for n = 1..10000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Run.html">Run</a>
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Feller%27s_coin-tossing_constants">Feller's coin-tossing constants</a>
%F Equals (22 + (847 - 33*sqrt(33))^(1/3) + (11*(77 + 3*sqrt(33)))^(1/3))/33. - _Vaclav Kotesovec_, Oct 14 2018
%F Positive real root of 11*x^3 - 22*x^2 + 12*x - 2. - _Peter Luschny_, Oct 14 2018
%F Equals 2/(5 - A058265^2). - _Jon Maiga_, Nov 25 2018
%e 1.2368398446387851018906608761421232....
%p evalf[120](solve(11*x^3-22*x^2+12*x-2=0,x)[1]); # _Muniru A Asiru_, Nov 25 2018
%t alpha = Root[1-x+(x/2)^4, x, 1]; beta = (2-alpha)/(4-3*alpha); RealDigits[beta, 10, 102] // First (* _Jean-François Alcover_, Jun 03 2014 *)
%o (PARI) default(realprecision, 100); (22 + (847 - 33*sqrt(33))^(1/3) + (11*(77 + 3*sqrt(33)))^(1/3))/33 \\ _G. C. Greubel_, Nov 25 2018
%o (Magma) SetDefaultRealField(RealField(100)); (22 + (847 - 33*Sqrt(33))^(1/3) + (11*(77 + 3*Sqrt(33)))^(1/3))/33; // _G. C. Greubel_, Nov 25 2018
%o (Sage) numerical_approx((22 + (847 - 33*sqrt(33))^(1/3) + (11*(77 + 3*sqrt(33)))^(1/3))/33, digits=100) # _G. C. Greubel_, Nov 25 2018
%Y Cf. A086253, A058265.
%K nonn,cons
%O 1,2
%A _Eric W. Weisstein_, Jul 13 2003
|
