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)
To embed this project on your website, copy the following code and paste it into your website's HTML: