OFFSET
1,1
COMMENTS
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 = 4.
LINKS
Stanislav Sykora, Table of n, a(n) for n = 1..2000
SpectralCalc, Calculation of Blackbody Radiance, Appendix C.
Wikipedia, Planck's law
EXAMPLE
3.9206903948728863435608913526135362205256273712079845304011750...
MATHEMATICA
RealDigits[x/.FindRoot[x==4(1-Exp[-x]), {x, 3}, WorkingPrecision->120]] [[1]] (* Harvey P. Dale, May 08 2017 *)
PROG
(PARI) a4=solve(x=0.1, 10, x-4*(1-exp(-x))) \\ Use real precision in excess
CROSSREFS
KEYWORD
AUTHOR
Stanislav Sykora, Apr 01 2015
STATUS
approved