from random import * def urnePolya(r,b,etapes): assert(r+b>=1) urne = ['r']*r+['b']*b for i in range(etapes): urne.append(choice(urne)) print urne.count('r')/len(urne)