généralisation : (programme python)
def suites_moyennes( a, b ) :
""" a, b dans ]0 ; 1[
u1 = a u0 + (1-a) v0 ( ne peut pas sortir de [u0 ; v0] )
v1 = (1-b) u0 + b v0
d = u - v : geometrique de raison (a+b-1)
t = (1-b) u + (1-a) v : constante
L = [ (1-b) u0 + (1-a) v0 ] / (2-a-b)
"""
# controle
if a <= 0 or a >= 1 or b <= 0 or b >= 1 :
print "Erreur : a et b doivent dans ] 0 ; 1 ["
sys.exit(1)
u0 = 0
v0 = 1
d0 = u0 - v0
t0 = (1-b) * u0 + (1-a) * v0
n = 0
print "(u0, v0, d0, t0) =", (u0, v0, d0, t0)
print "L =", t0 / (2-a-b)
for i in range(30) :
u1 = a * u0 + (1-a) * v0
v1 = (1-b) * u0 + b * v0
d1 = u1 - v1
t1 = (1-b) * u1 + (1-a) * v1
n += 1
print "(u1, v1, d1, t1) =", (u1, v1, d1, t1)
u0 = u1
v0 = v1
return