%r


repeticionesExperimento=10000 #sería el tamaño de muestreo para cada v.a.
cantidadPersonas=10000 #sería la cantidad de v.a.

resultado <- 0

for (i in 1:cantidadPersonas) {
	resultadoPersona <-  rbinom(repeticionesExperimento, 2, 0.5) #muestreo una binomial muchas veces
	resultado <- (resultado + resultadoPersona) #voy sumando los resultados: X_1 + X_2 + X_3 + ...
}
resultadoPromedio <- resultado/cantidadPersonas #calculo el promedio de la suma: sería el X_raya
hist(resultadoPromedio, breaks=50, probability=TRUE)
curve( dnorm(x, mean=mean(resultadoPromedio), sd=sd(resultadoPromedio)), add=TRUE)