#========================
# Event simulator function
#========================
success=function(n,p) {
  u=runif(n)
  x=sum((u<=p)*1)
  return(x)
}