Binomial-distributed Random Generator

É a distribuição de probabilidade discreta de um número específico de sucessos em uma sequência de n experimentos independentes, cada um fazendo uma pergunta sim-não e cada um com seu próprio resultado de valor booleano: sucesso com probabilidade p ou fracasso com probabilidade [math]q = 1 - p[/math].

Exemplo em Python

Exemplo com 1 milhão de sorteios

def pseudo_binomial(n=100, p=0.5, size=1):
    """
    Binomial generator from uniform generator
    """
    binom = []
    for _ in range(size):
    
        # Sets seed based on the decimal portion of the current system clock
        t = time.perf_counter()
        seed = int(10**9*float(str(t-int(t))[0:]))
        
        U = pseudo_uniform(seed=seed, size=n)
        Y = (U <= p).astype(int)
        binom.append(np.sum(Y))
    
    return binom