bernouli.py

def perform_bernoulli_trials(n, p):
    #Perform n Bernoulli trials with success probability p and return number of successes.
    # Initialize number of successes: n_success
    n_success = 0
# Perform trials
    for i in range(n):
        # Choose random number between zero and one: random_number
         random_number = np.random.random()
# If less than p, it's a success so add one to n_success
         if random_number<p:
            n_success+=1
return n_success