File: feigenbaum.en

package info (click to toggle)
euler 1.61.0-12
  • links: PTS, VCS
  • area: main
  • in suites: bookworm, bullseye, forky, sid, trixie
  • size: 5,164 kB
  • sloc: ansic: 24,761; sh: 8,314; makefile: 141; cpp: 47; php: 1
file content (15 lines) | stat: -rw-r--r-- 432 bytes parent folder | download | duplicates (8) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
% Try the following notebook, which produces the famous
% chaotic behaviour, which was investigated by Feigenbaum.
>function f(x,l)
$return l*x*(1-x);
$endfunction
% This function is iterated, starting from 0.5.
>l=3:0.005:4;
>x=0.5*ones(size(l));
% The niterate function is given an additional argument
% l.
>h=niterate("f",x,200;l); h=h[100:200];
>style("m.");
% xmark is a bit slow for that many points.
>xmark(l,h); wait(20);
>