def Bienaymé(n,p):
    bI = int(p*n-racine(n))+1
    bS = int(p*n+racine(n))
    return 1-sum( loi_binomiale(n,p)[bI:bS+1] )