The Museum of HP Calculators

HP Forum Archive 14

[ Return to Index | Top of Index ]

Feigenbaum-number
Message #1 Posted by Tizedes Csaba [Hungary] on 28 May 2004, 6:09 p.m.

Hello,

know anybody, what was the program, what used by Feigenbaum originally(!), to calculate this number on his HP65?

Thanks for answering!

Csaba

      
Re: Feigenbaum-number
Message #2 Posted by Julian Miranda (Spain) on 30 May 2004, 1:22 p.m.,
in response to message #1 by Tizedes Csaba [Hungary]

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)=rx[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)=rsin PIx(t). Where PI is 3.14159... as you may supposse.

I hope that can help you a little.

      
Re: Feigenbaum-number
Message #3 Posted by Valentin Albillo on 31 May 2004, 5:44 a.m.,
in response to message #1 by Tizedes Csaba [Hungary]

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:

Period Doubling Ad Infinitum

Best regards from V.

      
Re: Julian and Valentin - Bifurcation diagram on 28C
Message #4 Posted by Tizedes Csaba [Hungary] on 31 May 2004, 6:19 p.m.,
in response to message #1 by Tizedes Csaba [Hungary]

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


[ Return to Index | Top of Index ]

Go back to the main exhibit hall