#====================
# Figure 7.1.R
#====================
outcomes=function(n,p) {
  u=runif(n)
  x=numeric(n)
  for (i in 1:n) {
    if (u[i]<=p) x[i]=1 else x[i]=0
  }
  return(x)
}

n=30
p=.26
num.sets=100
bat.ave=numeric(num.sets)
for (i in 1:num.sets) {
  bat.ave[i]=sum(outcomes(n,p))/n
}

hist(bat.ave)
stem(bat.ave)