%I #19 Aug 22 2018 03:08:39
%S 1,5,9,3,6,2,4,2,6,0,0,4,0,0,4,0,0,9,2,3,2,3,0,4,1,8,7,5,8,7,5,1,6,0,
%T 2,4,1,7,8,9,0,0,2,4,2,4,8,1,8,8,5,9,3,6,4,9,9,9,5,0,4,5,1,1,6,9,6,0,
%U 8,4,9,8,4,8,1,6,1,8,7,9,5,0,2,3,2,7,4,9,9,2,7,6,6,1,8,4,4,0,7,1,4,1,7,0,6
%N Decimal expansion of the positive solution to x = 2*(1-exp(-x)).
%C Each of the positive solutions to x = q*(1-exp(-x)) obtained for q = 2, 3, 4, and 5, appears in several formulas pertinent to Planck's black-body radiation law. For a given q, the solution can be also written as q+W(-q/exp(q)), where W is the Lambert function. Here q = 2.
%C The constant appears in asymptotic formula for A007820. - _Vladimir Reshetnikov_, Oct 10 2016
%H Stanislav Sykora, <a href="/A256500/b256500.txt">Table of n, a(n) for n = 1..2000</a>
%H V. Kotesovec, <a href="https://oeis.org/wiki/User:Vaclav_Kotesovec">Non-attacking chess pieces</a>, 6ed, 2013, p. 249.
%H SpectralCalc, <a href="http://www.spectralcalc.com/blackbody/blackbody.html">Calculation of Blackbody Radiance</a>, Appendix C.
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Planck%27s_law">Planck's law</a>
%e 1.5936242600400400923230418758751602417890024248188593649995...
%t RealDigits[2 + LambertW[-2 Exp[-2]], 10, 100][[1]] (* _Vladimir Reshetnikov_, Oct 10 2016 *)
%o (PARI) a2=solve(x=0.1,10,x-2*(1-exp(-x))) \\ Use real precision in excess
%Y Cf. A194567 (q=3), A256501 (q=4), A256502 (q=5).
%Y Cf. A191236, A217905.
%K nonn,cons
%O 1,2
%A _Stanislav Sykora_, Mar 31 2015
|