Try This
from math import sqrt from random import random def ecartfreq(n=100,p=1/2,k=2): sigma = sqrt(p*(1-p)) binf = p-k*sigma/sqrt(n) bsup = p+k*sigma/sqrt(n) N = 1000 interieur = 0 for i in range(N): S = 0 for j in range(n): if (random() < p): S = S+1 f = S/n if (binf < f) and (f < bsup): interieur = interieur + 1 F = interieur/N return F print(" ",ecartfreq()," ")
Run