[step:Compare the volumes of disjoint balls around the net points]
For each $u \in N$, define the open Euclidean ball
\begin{align*}
B\left(u,\frac{\varepsilon}{2}\right) := \left\{y \in \mathbb{R}^d : |y-u| < \frac{\varepsilon}{2}\right\}.
\end{align*}
If $u,v \in N$ are distinct, then
\begin{align*}
B\left(u,\frac{\varepsilon}{2}\right) \cap B\left(v,\frac{\varepsilon}{2}\right) = \varnothing.
\end{align*}
Indeed, if $y$ belonged to the intersection, the triangle inequality would give
\begin{align*}
|u-v| \leq |u-y| + |y-v| < \frac{\varepsilon}{2}+\frac{\varepsilon}{2} = \varepsilon,
\end{align*}
contradicting strict $\varepsilon$-separation.
Moreover, since $|u|=1$ for every $u \in N$, each point $y \in B(u,\varepsilon/2)$ satisfies
\begin{align*}
|y| \leq |y-u|+|u| < \frac{\varepsilon}{2}+1.
\end{align*}
Thus
\begin{align*}
\bigcup_{u \in N} B\left(u,\frac{\varepsilon}{2}\right)
\subset
B\left(0,1+\frac{\varepsilon}{2}\right).
\end{align*}
Let $\mathcal{L}^d$ denote the $d$-dimensional [Lebesgue measure](/page/Lebesgue%20Measure) on $\mathbb{R}^d$, and define
\begin{align*}
\alpha_d := \mathcal{L}^d(B(0,1)).
\end{align*}
Thus $\alpha_d$ is the $d$-dimensional Lebesgue measure of the open unit ball in $\mathbb{R}^d$. By translation and scaling invariance of $\mathcal{L}^d$,
\begin{align*}
\mathcal{L}^d\left(B\left(u,\frac{\varepsilon}{2}\right)\right)
= \alpha_d \left(\frac{\varepsilon}{2}\right)^d
\end{align*}
for every $u \in N$, and
\begin{align*}
\mathcal{L}^d\left(B\left(0,1+\frac{\varepsilon}{2}\right)\right)
= \alpha_d \left(1+\frac{\varepsilon}{2}\right)^d.
\end{align*}
Since the balls $B(u,\varepsilon/2)$ are pairwise disjoint and contained in $B(0,1+\varepsilon/2)$, finite additivity for pairwise disjoint measurable sets gives, for every finite subset $F \subset N$,
\begin{align*}
|F|\,\alpha_d\left(\frac{\varepsilon}{2}\right)^d
=
\mathcal{L}^d\left(\bigcup_{u \in F} B\left(u,\frac{\varepsilon}{2}\right)\right)
\leq
\alpha_d\left(1+\frac{\varepsilon}{2}\right)^d.
\end{align*}
Because $\alpha_d>0$, this implies
\begin{align*}
|F|
\leq
\left(\frac{1+\varepsilon/2}{\varepsilon/2}\right)^d
=
\left(1+\frac{2}{\varepsilon}\right)^d.
\end{align*}
The right-hand side is finite, so $N$ cannot contain arbitrarily large finite subsets. Therefore $N$ is finite, and taking $F=N$ yields
\begin{align*}
|N| \leq \left(1+\frac{2}{\varepsilon}\right)^d.
\end{align*}
[/step]