Bernoulli

set.seed(8)
y=NULL
s=1
p=0.5
n=20
for(i in 1:500)
{
 x=rbinom(n,s,p) #size:numero de intentos
 x
 y[i]=sum(x)
}
y

hist(y, freq=FALSE)  #Para que no salga la frecuencia
E_y=n*p
V_y=n*p*(1-p)
val=(seq(min(y),max(y),length=1000))
lines(val, dnorm(val,E_y,sqrt(V_y)))

Z=(y-E_y)/sqrt(V_y)
hist(Z, freq=FALSE)
val2=(seq(-3,3,length=1000))
lines(val, dnorm(val2,0,1),lty=2)
legend(-3,0.6,legend="hist.normal",lty=2)