#estimación de la probabilidad de aparición de dos caras 
#en tres lanzamientos de una moneda

#m<-10000
cat("Introduce el tamaño de la muestra")
m<-scan()
C=0
for (k in 1:m) {U=runif(3,0,1)
	X=sum(floor(2*U))#parte entera
	#X<-sum(round(2*U,0)) #redondea;si ponemos cero decimales=parte entera
	if(X==2) C=C+1
}
p<-C/m
cat("la probabilidad estimada es:",p)

#show(p) #muestra el resultado