def histo(N,n,p): L =[0]*(n+1) for k in range(N): ns = ech(n,p).count(1) L[ns] = L[ns] + 1 for k in range(n+1): fill_rect(2*k,200-L[k],2,L[k],bleu)