from math import sqrt
from random import random

def epreuve(p):
    if random() < p:
        return 1
    else:
        return 0

def var_binom(n,p):
    somme = 0
    for k in range(n):
        somme = somme + epreuve(p)
    return somme

def simul10000(n,p):
    dehors = 0
    for k in range(10000):
        if abs(var_binom(n,p)-n*p)>=...:
            dehors = dehors + 1
    return ...