import numpy.random as rd
def Y(n,p):
    X = rd.binomial(n,p)
    if X == 0 or X ==n :
        y = rd.randint(1,n)
    else : 
        y = X 
    return (y)

N = 10000
n = 9
p = 1/3

for i in range (0,N):
    val = Y(n,p)
    print(val)

Embed on website

To embed this project on your website, copy the following code and paste it into your website's HTML: