IF = (N,p) ->
    [bI,bS] = [0,N]
    somme = 0
    until somme >= 0.025
        somme = somme + binomiale(N,p,bI)
        bI = bI + 1
    bI = bI - 1
    somme = 0
    until somme >= 0.025
        somme = somme + binomiale(N,p,bS)
        bS = bS - 1
    bS = bS + 1
    [bI,bS]