def bernoulli(p):
    if random()<p:
        return(1)
    else:
        return(0)