Feigenbaum constant calculation: Difference between revisions
Content added Content deleted
(adding lambdatalk) |
|||
| Line 952: | Line 952: | ||
ReadChar |
ReadChar |
||
END Feigenbaum.</lang> |
END Feigenbaum.</lang> |
||
=={{header|Nim}}== |
|||
{{trans|Kotlin}} |
|||
<lang Nim>import strformat |
|||
iterator feigenbaum(): tuple[n: int; δ: float] = |
|||
## Yield |
|||
const |
|||
MaxI = 13 |
|||
MaxJ = 10 |
|||
var |
|||
a1 = 1.0f |
|||
a2 = 0.0 |
|||
δ = 3.2 |
|||
for i in 2..MaxI: |
|||
var a = a1 + (a1 - a2) / δ |
|||
for j in 1..MaxJ: |
|||
var x, y = 0.0 |
|||
for _ in 1..(1 shl i): |
|||
y = 1 - 2 * y * x |
|||
x = a - x * x |
|||
a -= x / y |
|||
δ = (a1 - a2) / (a - a1) |
|||
a2 = a1 |
|||
a1 = a |
|||
yield (i, δ) |
|||
echo " i δ" |
|||
for n, δ in feigenbaum(): |
|||
echo fmt"{n:2d} {δ:.8f}"</lang> |
|||
{{out}} |
|||
<pre> i δ |
|||
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|Perl}}== |
=={{header|Perl}}== |
||