from random import random import matplotlib.pyplot as plt def echantillon(p, n): L = [] for i in range(n): if random()<p: L.append(1) else: L.append(0) return L def fluct_echant(nb_echant, p, n): L_freq = [] for i in range(nb_echant): eff = 0 for i in range(n): if random()<p: eff += 1 freq = eff/n L_freq.append(freq) L_num = list(range(1,nb_echant+1)) plt.plot([1,nb_echant],[p,p], color='red') plt.scatter(L_num, L_freq, color='blue') plt.show()