[guided]Assume that $f$ is uniformly continuous. By definition, this means that for every tolerance $\varepsilon>0$ there is at least one positive radius $\eta>0$ such that the implication
\begin{align*}d_X(x,x')<\eta \implies d_Y(f(x),f(x'))<\varepsilon\end{align*}
holds for all $x,x'\in X$.
The point of the theorem is to turn these separate existence statements into one function of $\varepsilon$. To avoid making an arbitrary choice of one radius for every $\varepsilon$, we collect all radii that work. For each $\varepsilon>0$, define
\begin{align*}A_\varepsilon:=\{r\in(0,1]:\text{ for all }x,x'\in X,\ d_X(x,x')<r\implies d_Y(f(x),f(x'))<\varepsilon\}.\end{align*}
The restriction $r\le 1$ is harmless; it only makes the set bounded above. We verify that $A_\varepsilon$ is nonempty. Uniform continuity gives an admissible $\eta>0$ for this $\varepsilon$. Then
\begin{align*}\frac{\min\{\eta,1\}}{2} \in (0,1]\end{align*}
and whenever
\begin{align*}d_X(x,x')<\frac{\min\{\eta,1\}}{2}\end{align*}
we also have $d_X(x,x')<\eta$, so the uniform continuity implication gives $d_Y(f(x),f(x'))<\varepsilon$. Hence
\begin{align*}\frac{\min\{\eta,1\}}{2}\in A_\varepsilon\end{align*}.
Because $A_\varepsilon$ is a nonempty subset of $(0,1]$, its supremum exists as a real number. Define
\begin{align*}
s(\varepsilon):=\sup A_\varepsilon.
\end{align*}
The previous paragraph shows $s(\varepsilon)>0$, and the inclusion $A_\varepsilon\subset(0,1]$ gives $s(\varepsilon)\le 1$. Now define the candidate modulus $\delta:(0,\infty)\to(0,\infty)$ by
\begin{align*}
\delta(\varepsilon):=\frac{s(\varepsilon)}{2}.
\end{align*}
This is positive because $s(\varepsilon)>0$.
It remains to check that this particular value actually works. Fix $\varepsilon>0$. Since $\delta(\varepsilon)<s(\varepsilon)=\sup A_\varepsilon$, the definition of supremum gives some $r_\varepsilon\in A_\varepsilon$ with
\begin{align*}
\delta(\varepsilon)<r_\varepsilon.
\end{align*}
Now take arbitrary points $x,x'\in X$ satisfying $d_X(x,x')<\delta(\varepsilon)$. Then
\begin{align*}
d_X(x,x')<r_\varepsilon.
\end{align*}
Since $r_\varepsilon\in A_\varepsilon$, its defining property applies and yields
\begin{align*}
d_Y(f(x),f(x'))<\varepsilon.
\end{align*}
Thus the function $\delta$ satisfies the required epsilon-delta modulus condition.[/guided]