Datos de ejercicio:

P=0.1

Z=1.96

sd=0.05

N=1000

fP=function(P, Z, sd, N){
n0 =(P*(1-P)/(sd/Z)^2)
n =(n0/(1+(n0/N)))
return (n)
}

fP(0.1, 1.96, 0.05, 1000)
## [1] 121.4951