N = 1000
s2 = 45 
d = 500
alpha = 0.05

z = qnorm(1-alpha/2)
z

n1 = (d^2/(z^2*s2*N^2)+1/N)^{-1}
n = ceiling(n1)
n