import numpy.random as rd

def simulation(n) :
    urne = [0,1]
    rouge = 1
    for k in range(n) :
        r = rd.randint(len(urne))
        if urne[r] == 0 :
            urne.append(0)
            rouge += 1
        else :
            urne.append(1)
    return rouge
    
def moyenne(n,N) :
    L = [0]*(n+3)
    for k in range(N) :
        L[simulation(n)] += 1
    return [i/N for i in L]