23 lines
439 B
Python
23 lines
439 B
Python
max_it = 13
|
|
max_it_j = 10
|
|
a1 = 1.0
|
|
a2 = 0.0
|
|
d1 = 3.2
|
|
a = 0.0
|
|
|
|
print " i d"
|
|
for i in range(2, max_it + 1):
|
|
a = a1 + (a1 - a2) / d1
|
|
for j in range(1, max_it_j + 1):
|
|
x = 0.0
|
|
y = 0.0
|
|
for k in range(1, (1 << i) + 1):
|
|
y = 1.0 - 2.0 * y * x
|
|
x = a - x * x
|
|
a = a - x / y
|
|
d = (a1 - a2) / (a - a1)
|
|
print("{0:2d} {1:.8f}".format(i, d))
|
|
d1 = d
|
|
a2 = a1
|
|
a1 = a
|