29 lines
766 B
Groovy
29 lines
766 B
Groovy
class Feigenbaum {
|
|
static void main(String[] args) {
|
|
int max_it = 13
|
|
int max_it_j = 10
|
|
double a1 = 1.0
|
|
double a2 = 0.0
|
|
double d1 = 3.2
|
|
double a
|
|
|
|
println(" i d")
|
|
for (int i = 2; i <= max_it; i++) {
|
|
a = a1 + (a1 - a2) / d1
|
|
for (int j = 0; j < max_it_j; j++) {
|
|
double x = 0.0
|
|
double y = 0.0
|
|
for (int k = 0; k < 1 << i; k++) {
|
|
y = 1.0 - 2.0 * y * x
|
|
x = a - x * x
|
|
}
|
|
a -= x / y
|
|
}
|
|
double d = (a1 - a2) / (a - a1)
|
|
printf("%2d %.8f\n", i, d)
|
|
d1 = d
|
|
a2 = a1
|
|
a1 = a
|
|
}
|
|
}
|
|
}
|