Feigenbaum-number « Next Oldest | Next Newest »

 ▼ Tizedes Csaba [Hungary] Member Posts: 97 Threads: 25 Joined: Jan 1970 05-28-2004, 06:09 PM Hello, know anybody, what was the program, what used by Feigenbaum originally(!), to calculate this number on his HP65? Thanks for answering! Csaba ▼ Julián Miranda (Spain) Member Posts: 57 Threads: 3 Joined: Feb 2007 05-30-2004, 01:22 PM I don't have the program. From the book "Chaos. Making a New Science" by James Gleick, I can tell you the equations he was working with. First he worked with x(t+1)=r·x·[1-x(t)]. A quadratic equation. y=rx-rx^2. You start with a value for "r" and a value for "x" and the result is the next value for "x" if you keep iterating finally you get always the same value. With some pairs of values you end in an alternation between two values. He searched for the values of "r" that produced this duplication. Samples of values r=2.7 x=0.02 that ends in 0.6296. r=3.5 x1=0.850 x2=0.3828 x3=0.8269 x4=0.5009, in this case there is a double duplication so we end with four x values. Then he repeated the process with the following equation x(t+1)=r·sin PI·x(t). Where PI is 3.14159... as you may supposse. I hope that can help you a little. Valentin Albillo Posting Freak Posts: 1,755 Threads: 112 Joined: Jan 2005 05-31-2004, 05:44 AM I don't think the exact HP-65 program is available, but it can probably be reconstructed from the equations he worked with (see J. Miranda's post) and the very detailed (and to the point) info on this URL: Best regards from V. Tizedes Csaba [Hungary] Member Posts: 97 Threads: 25 Joined: Jan 1970 05-31-2004, 06:19 PM Thanks for your answers! I read that book some years ago, and I wrote a little program in PASCAL. The next prg is the HP28C version: ```FEIG: << (2.63,0) PMIN (4,1) PMAX CLLCD 2.63 4 FOR R .3 'X' STO 1 20 START R NEWX NEXT 1 12 START R NEWX R X R->C PIXEL NEXT .01 STEP >> NEWX: << -> R << R X * 1 X - * 'X' STO >> >> ``` Press 'FEIG' to start. Running time is about 12 minutes. (Slow, but nice... ;) ) The greatest problem, how to isolate the bifurcation points on the diagram... Csaba

 Possibly Related Threads... Thread Author Replies Views Last Post Re: Feigenbaum number Valentin Albillo 5 333 06-07-2004, 04:23 PM Last Post: Tizedes Csaba [Hungary] HP-65 helped discovering the Feigenbaum number Juergen (CH) 32 1,630 06-30-2003, 11:22 AM Last Post: Paul Brogger

Forum Jump: