Jump to content

Feigenbaum constant calculation: Difference between revisions

(adding lambdatalk)
Line 952:
ReadChar
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}}==
Anonymous user
Cookies help us deliver our services. By using our services, you agree to our use of cookies.