import random def urnesPolya(r,b,etapes): assert(r+b>=1) L=[r/(r+b)] for i in range(etapes): n = random.randint(1,r+b) if n<=r: r,b = r+1,b else: r,b = r,b+1 L.append(r/(r+b))