blob: 143e6c3eb6c6c44d90829e6df546181065cb7ae0 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
#!/usr/bin/bc
scale = 20
phi = (1 + sqrt(5)) / 2
loglimit = 1000 - 1 + l(sqrt(5))/l(10)
scale = 0
n = 1
min = -1
max = -1
while (1) {
print n, " (", min, ", ", max, ")\n"
scale = 20
if (l(phi^n - (1 - phi)^n)/l(10) > loglimit) {
max = n
} else {
min = n
}
scale = 0
if (max - min == 1) {
max
halt
}
if (max == -1) {
n *= 2
} else {
n = (max + min) / 2
}
}
|