import numpy as np
import numpy.random as rd

def Y(p,n) :
    Y = [0]*(n+1)
    for _ in range(int(1e3)) :
        z = 0
        for _ in range(p) :
            k = rd.randint(1,n)
            if k > z :
                z += 1
            else :
                z -= 1
        Y[z] += 1
    return [y*1e-3 for y in Y]

n = 10
for p in [100,101,120,151] :
    print('p = ',p,'Y = ',Y(p,n))