def phi1(p):
a = 2
n = 0
diff = 2
while diff >=p:
a2 = a
a = 1 + 1.0/a
diff = abs(a2-a)
n = n+1
return [n,a]
print phi1(0.001)
print phi1(0.000001)
print phi1(0.000000001)
xxxxxxxxxx
from math import*
def phi2(p):
b = 2
n = 0
diff = 2
while diff >=p:
b2 = b
b = sqrt(1.0+b)
diff = abs(b2-b)
n = n+1
return [n,b]
print phi2(0.001)
print phi2(0.000001)
print phi2(0.000000001)
xxxxxxxxxx
from math import*
def phi3(p):
c = 2
n = 0
diff = 2
while diff >=p:
c2 = c
c = (c*c+1.0)/(2*c-1.0)
diff = abs(c2-c)
n = n+1
return [n,c]
print phi3(0.001)
print phi3(0.000001)
print phi3(0.000000001)
xxxxxxxxxx
def fibonacci(n):
v = 1
u = 1
for i in range(1,n):
w = v
v = u
u = u+w
return u
def ratio(n):
return (1.0*fibonacci(n+1))/fibonacci(n)
def convergence(p):
diff = 2
n = 0
while diff >= p:
diff = abs( ratio(n+1)-ratio(n) )
n = n+1
return [n,ratio(n)]
print convergence(0.001)
print convergence(0.000001)
print convergence(0.000000001)