Feigenbaum constant calculation: Difference between revisions
Content added Content deleted
(Added Wren) |
|||
Line 1,441: | Line 1,441: | ||
12 4.66920098 |
12 4.66920098 |
||
13 4.66920537</pre> |
13 4.66920537</pre> |
||
=={{header|Wren}}== |
|||
{{trans|Ring}} |
|||
{{libheader|Wren-fmt}} |
|||
<lang ecmascript>import "/fmt" for Fmt |
|||
var feigenbaum = Fn.new { |
|||
var maxIt = 13 |
|||
var maxItJ = 10 |
|||
var a1 = 1 |
|||
var a2 = 0 |
|||
var d1 = 3.2 |
|||
System.print(" i d") |
|||
for (i in 2..maxIt) { |
|||
var a = a1 + (a1 - a2)/d1 |
|||
for (j in 1..maxItJ) { |
|||
var x = 0 |
|||
var y = 0 |
|||
for (k in 1..(1<<i)) { |
|||
y = 1 - 2*y*x |
|||
x = a - x*x |
|||
} |
|||
a = a - x/y |
|||
} |
|||
var d = (a1 - a2)/(a - a1) |
|||
System.print("%(Fmt.d(2, i)) %(Fmt.f(0, d, 8))") |
|||
d1 = d |
|||
a2 = a1 |
|||
a1 = a |
|||
} |
|||
} |
|||
feigenbaum.call()</lang> |
|||
{{out}} |
|||
<pre> |
|||
i d |
|||
2 3.21851142 |
|||
3 4.38567760 |
|||
4 4.60094928 |
|||
5 4.65513050 |
|||
6 4.66611195 |
|||
7 4.66854858 |
|||
8 4.66906066 |
|||
9 4.66917155 |
|||
10 4.66919515 |
|||
11 4.66920026 |
|||
12 4.66920098 |
|||
13 4.66920537 |
|||
</pre> |
|||
=={{header|zkl}}== |
=={{header|zkl}}== |