How to calculate the Feigenbaum constant to high precision? I am looking for a reference on how to calculate the Feigenbaum constant $\delta$ to high precision. It seems that naive methods do not work well because they lead to solving very high degree polynomials.
I have searched Math Reviews, Googled around, and the best I found was Appendix A of Keith Brigg's thesis. Wikipedia gives a link which contains $\delta$ to high precision, but does not explain how it was computed.
 A: As far as I know O.E. Lanford III was the first to compute the Feigenbaum constant around 1980 using interval arithmetics. That is to say, he represented real numbers not by floating point arithmetics but rather using intervals containing the sought number. If I remember correctly he had at least fifty decimals at some point (he told me orally, so I might be completely wrong).
If you enter Author="Lanford" and Anywhere="Feigenbaum" in MathSciNet, you should find the relevant articles and some general remarks on using computers in rigorous proofs (I don't know how to post links here).
A: One idea for computing the Feigenbaum constant is as follows: $\delta$ is the unique expanding eigenvalue for period-doubling renormalization at its fixed point $F$. So take any real-analytic 1-parameter family of univalent maps that is transverse to the stable manifold of renormalization (for example, the real quadratic family). Then the iterated renormalizations of this 1-parameter family will converge to the unstable manifold of renormalization. The unstable manifold is of course mapped to itself by renormalization and the derivative at the fixed point will be $\delta$. 
It should be possible to implement this numerically, by keeping track of the power series for a real-analytic family of real-analytic univalent map, and applying the renormalization by replacing $f_\lambda$ with $f_\lambda \circ f_\lambda$ (and suitably rescaling). Because renormalization acts as a contraction on these 1-parameter families, this procedure should be computationally stable, and provide a number of digits proportional to the number of times one renormalizes. 
One classical reference for the stable/unstable manifold picture for period doubling renormalization is Iterated Maps on the Interval as Dynamical Systems by Collet and Eckmann. 
A: David Broadhurst computed some 15 years ago those constants to high precision (in particular $\delta$ to 1018 digits) using Keith Briggs' method
http://www.plouffe.fr/simon/constants/feigenbaum.txt
Note: the link provided by Andrej is not valid anymore, sorry if my answer is just a duplicate (today Wikipedia gives a link to OEIS.org for a value with about 103 digits).
