Feigenbaum constant calculation: Difference between revisions

Content deleted Content added

Inline

@@ Line 287: / Line 287: @@
 4.66920098
 4.66920537</pre>
+=={{header|Phix}}==
+{{trans|Ring}}
+<lang Phix>constant maxIt = 13,
+        maxItJ = 10
+atom a1 = 1.0,
+     a2 = 0.0,
+     d1 = 3.2
+puts(1," i d\n")
+for i=2 to maxIt do
+     atom a = a1 + (a1 - a2) / d1
+     for j=1 to maxItJ do
+          atom x = 0, y = 0
+          for k=1 to power(2,i) do
+               y = 1 - 2*y*x
+               x = a - x*x
+          end for
+          a = a - x/y
+     end for
+     atom d = (a1-a2)/(a-a1)
+     printf(1,"%2d %.8f\n",{i,d})
+     d1 = d
+     a2 = a1
+     a1 = a
+end for</lang>
+{{out}}
+<pre>
+ i d
+3.21851142
+4.38567760
+4.60094928
+4.65513050
+4.66611195
+4.66854858
+4.66906066
+4.66917155
+4.66919515
+4.66920026
+4.66920098
+4.66920537
+</pre>
 =={{header|Racket}}==