Feigenbaum constant calculation: Difference between revisions
Content deleted Content added
Amended template for converting this to a draft task. |
Added Kotlin |
||
| Line 1: | Line 1: | ||
{{draft task}} |
{{draft task}} |
||
Calculate the Feigenbaum constant. See details: [https://en.wikipedia.org/wiki/Feigenbaum_constants Feigenbaum constant] |
Calculate the Feigenbaum constant. See details: [https://en.wikipedia.org/wiki/Feigenbaum_constants Feigenbaum constant] |
||
=={{header|Kotlin}}== |
|||
{{trans|Ring}} |
|||
<lang scala>// Version 1.2.40 |
|||
fun feigenbaum() { |
|||
val maxIt = 13 |
|||
val maxItJ = 10 |
|||
var a1 = 1.0 |
|||
var a2 = 0.0 |
|||
var d1 = 3.2 |
|||
println(" i d") |
|||
for (i in 2..maxIt) { |
|||
var a = a1 + (a1 - a2) / d1 |
|||
for (j in 1..maxItJ) { |
|||
var x = 0.0 |
|||
var y = 0.0 |
|||
for (k in 1..(1 shl i)) { |
|||
y = 1.0 - 2.0 * y * x |
|||
x = a - x * x |
|||
} |
|||
a -= x / y |
|||
} |
|||
val d = (a1 - a2) / (a - a1) |
|||
println("%2d %.8f".format(i,d)) |
|||
d1 = d |
|||
a2 = a1 |
|||
a1 = a |
|||
} |
|||
} |
|||
fun main(args: Array<String>) { |
|||
feigenbaum() |
|||
}</lang> |
|||
{{output}} |
|||
<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|Ring}}== |
=={{header|Ring}}== |
||