Mathematica Stack Exchange is a question and answer site for users of Mathematica. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I am working on a problem to find the Feigenbaum constant. In this problem I have to get $f^{2^n}(1/2)$ with $f(x) = kx(1-x)$. I use Nest to get $f^{2^n}(1/2)$ but it takes Mathematica forever to run with $n>5$ but my goal is to get to $n=14$.

Does any body know any function or trick to get $f^{2^n}(1/2)$ with large value of $n$?

Thank you very much.

share|improve this question
Ha! another stupid decision from SO 10K users – Dr. belisarius Apr 29 '13 at 14:38
You do know that iterating the logistic equation eventually results in a periodic sequence, no? Find the period, and then you can extrapolate on what the $2^n$-th iterate ought to be. – J. M. Apr 29 '13 at 14:43
@belisarius Try k=2.9 or 3.1 or 3.4. You happen to be choosing a value of k right at a bifurcation, where the convergence to the fixed point is absurdly slow. But, of course, the iteration need not always converge to a periodic orbit; that's the point behind chaos. – Mark McClure Apr 29 '13 at 15:01
@belisarius No, I don't think you can analytically find the attractive orbit in terms of k in a simple way. What I think J.M. is suggesting is extrapolate from when the bifurcations occur. In fact, the existence of the Feigenbaum constant is really an assertion that such an extrapolation works. – Mark McClure Apr 29 '13 at 15:10
@thewanderer I assume that you mean when n>5, but note that for n=5, your polynomial has degree 2^2^5. and it's simply not feasible to solve this using simple techniques. You should look into some of the papers of Keith Briggs to learn his extrapolation techniques: – Mark McClure Apr 29 '13 at 21:26

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.