Lanzamiento de una moneda (distribución binomial)

Se lanza una moneda 500 veces. Obtener, exactamente y mediante simulación, la probabilidad de que el número de caras no difiera de 250 en más de 10.

#INICIO -------------------------
rm(list=ls(all=TRUE))

#Exactamente
pbinom(260,500,0.5)-pbinom(239,500,0.5)

#Mediante simulación
n<-1000000 #n es el número de pruebas
A<-replicate(n,abs(rbinom(1,500,0.5)-250))
f<-function(x) if(x<=10) 1 else 0
for (i in 1:n)
A[i]<-f(A[i])
sum(A)/n
#FIN -------------------------