Estimating the Feigenbaum Constant from a One-Parameter Scaling Law
Requires a Wolfram Notebook System
Interact on desktop, mobile and cloud with the free Wolfram CDF Player or other Wolfram Language products.
Requires a Wolfram Notebook System
Edit on desktop, mobile and cloud with any Wolfram Language product.
Mitchell J. Feigenbaum's one-parameter scaling law for one-dimensional iterative maps with -unimodality is given by[more]
for any ,
where 1. is the order of the period-doubling pitchfork bifurcation; 2. is the control parameter of the iterative maps; 3. is the super-stable parameter value  for each bifurcation order (e.g. for period 1, for period 2, for period 4, for period 8, and so on); 4. is a constant; 5. is the Feigenbaum constant as a function of . For , .
On the left is the plot of versus . For any , the filled-in blue circles within the fitting interval fall nicely on a straight line with the slope , where is the uncertainty, indicating the above scaling law. The uncertainty on can be found from a standard linear regression analysis.
On the right is the plot of versus . The accuracy of the estimated value of can be increased by decreasing the length of the fitting interval, which is displayed in the dropdown menu button.[less]
1. The test map used in this Demonstration, [8, 9], generalizes the well-known logistic map [3–7].
2. The program for this Demonstration presents a large collection of super-stable parameter values () for period-doubling periodic attractors within and , which were obtained using a high-precision Newton algorithm  with fixed precision 50 . These values are accurate up to 45 decimal places, so you can use them for your own research or study.
3. Keith Briggs calculated precise values of for in his 1997 thesis . For even , the estimated values of in this Demonstration are consistent with his calculation. But, for odd , the estimated values are slightly different from his calculation. It would be nice if someone could explain this. Refer to the following table for comparison:
Here D and K refers to for Delbourgo and Kenny's calculation ; Ch, Cv, and R for Christiansen, Cvitanović, and Rugh's calculation ; B for Briggs' calculation; and M for the calculation in this Demonstration for the fitting interval .
4. Of particular interest is that the following limits exist:
 In mathematics, specifically bifurcation theory, the Feigenbaum constants are two mathematical constants that both express ratios in a bifurcation diagram for a nonlinear map. They are named after the mathematician/physicist Mitchell J. Feigenbaum. The first Feigenbaum constant is the limiting ratio of each bifurcation interval to the next between every period doubling of iterative maps with degree-2 unimodality. The second Feigenbaum constant is the ratio between the width of a tine and the width of one of its two subtines (except the tine closest to the fold). Both numbers are believed to be transcendental, although they have not been proven to be so. See the related Wikipedia article, "Feigenbaum Constants." See also [3-9].
 M. J. Feigenbaum, "Quantitative Universality for a Class of Non-Linear Transformations," Journal of Statistical Physics, 19, 1978 pp. 25–52.
 M. J. Feigenbaum, "The Universal Metric Properties of Nonlinear Transformations," Journal of Statistical Physics, 21, 1979 pp. 669–706.
 K. T. Alligood, T. D. Sauer, and J. A. Yorke, Chaos: An Introduction to Dynamical Systems, New York: Springer, 1996.
 S. H. Strogatz, Nonlinear Dynamics and Chaos, Studies in Nonlinearity, New York: Perseus Books Publishing, 1994.
 S. Wolfram, A New Kind of Science, Champaign, IL: Wolfram Media, 2002.
 K.-J. Moon, "Reducible Expansions and Related Sharp Crossovers in Feigenbaum's Renormalization Field," Chaos: An Interdisciplinary Journal of Nonlinear Sciences, 18, 2008 pp. 023104.
 K.-J. Moon, "Erratum: Reducible Expansions and Related Sharp Crossovers in Feigenbaum's Renormalization Field," Chaos: An Interdisciplinary Journal of Nonlinear Sciences, 20, 2010 pp. 049902.
 In numerical analysis, the Newton method (also known as the Newton–Raphson method), named after Issac Newton and Joseph Raphson, is a method for finding successively better approximations to the roots of a real-valued function. See the related Wikipedia article, "Newton's Method."
 In computing, a fixed-point number representation is a real data type for a number that has a fixed number of digits below the decimal point. Fixed-point number representation can be compared to the more complicated (and more computationally demanding) floating-point number representation. See the related Wikipedia article, "Fixed-Point Arithmetic."
 K. Briggs, "Feigenbaum Scaling in Discrete Dynamical Systems," Ph.D. thesis, Department of Mathematics, Melbourne University, Australia, 1997. keithbriggs.info/documents/Keith_Briggs_PhD.pdf.
 R. Delbourgo, B. G. Kenny, "Universality Relations," Physical Review A, 33, 1986 pp. 3293–3302.
 F. Christiansen, P. Cvitanović, and H. H. Rugh, "The Spectrum of the Period-Doubling Operator in Terms of Cycles," Journal of Physics A, 23, 1990 pp. 713–717.
 P. Collet, J.-P. Eckmann, and O. E. Lanford, "Universal Properties of Maps on an Interval," Communications in Mathematical Physics, 76, 1980 pp. 211–254.
 J.-P. Eckmann and P. Witter, Computer Methods and Borel Summability Applied to Feigenbaum's Equations, New York: Springer–Verlag, 1985.
 J. P. van der Weele, H. W. Capel, and R. Kluiving, "On the Scaling Factors and ," Physics Letters A, 119, 1986 pp. 15–20.
 J. P. van der Weele, H. W. Capel, and R. Kluiving, "Period-Doubling in Maps with a Maximum of order ," Physica A, 145, 1987 pp. 425–460.