[step:Approximate arbitrary measurable sets by elements of the cylinder algebra]Let $E, F \in \mathcal{B}$ and $\varepsilon > 0$. By Step 1, $\mathcal{A}$ is an algebra with $\sigma(\mathcal{A}) = \mathcal{B}$, and $\mu$ is a probability measure on $(\Omega, \mathcal{B})$. By the [Caratheodory Approximation Theorem](/theorems/???) (for every $G \in \sigma(\mathcal{A})$ and every $\delta > 0$ there exists $G' \in \mathcal{A}$ with $\mu(G \triangle G') < \delta$), there exist $E', F' \in \mathcal{A}$ with
\begin{align*}
\mu(E \triangle E') < \varepsilon, \qquad \mu(F \triangle F') < \varepsilon,
\end{align*}
where $E \triangle E' := (E \setminus E') \cup (E' \setminus E)$ denotes the [symmetric difference](/page/Symmetric%20Difference).
For every $N \in \mathbb{N}$ define the mixing error
\begin{align*}
D_N(E, F) := \left|\mu(\sigma^{-N}E \cap F) - \mu(E)\mu(F)\right|.
\end{align*}
*Intersection perturbation.* Because preimages commute with set operations and $\sigma$ preserves $\mu$,
\begin{align*}
\mu\bigl((\sigma^{-N}E) \triangle (\sigma^{-N}E')\bigr) = \mu(\sigma^{-N}(E \triangle E')) = \mu(E \triangle E') < \varepsilon.
\end{align*}
For any measurable $G, G', H, H'$ the inclusion
\begin{align*}
(G \cap H) \triangle (G' \cap H') \subseteq (G \triangle G') \cup (H \triangle H')
\end{align*}
holds (if $\omega \in G \cap H$ but $\omega \notin G' \cap H'$, then $\omega \notin G'$ or $\omega \notin H'$; the symmetric case is identical), so subadditivity of $\mu$ gives
\begin{align*}
\left|\mu(G \cap H) - \mu(G' \cap H')\right| \leq \mu(G \triangle G') + \mu(H \triangle H').
\end{align*}
Applying this with $G = \sigma^{-N}E$, $G' = \sigma^{-N}E'$, $H = F$, $H' = F'$,
\begin{align*}
\left|\mu(\sigma^{-N}E \cap F) - \mu(\sigma^{-N}E' \cap F')\right| < 2\varepsilon.
\end{align*}
*Product perturbation.* The bound $|\mu(G) - \mu(G')| \leq \mu(G \triangle G')$ follows from $G \setminus G', G' \setminus G \subseteq G \triangle G'$. Using $\mu(F), \mu(E') \leq 1$,
\begin{align*}
|\mu(E)\mu(F) - \mu(E')\mu(F')|
&\leq |\mu(E) - \mu(E')|\mu(F) + \mu(E')|\mu(F) - \mu(F')| \\
&\leq \mu(E \triangle E') + \mu(F \triangle F') < 2\varepsilon.
\end{align*}
Combining the two perturbation bounds by the triangle inequality,
\begin{align*}
D_N(E, F) \leq D_N(E', F') + 4\varepsilon.
\end{align*}
By Step 2 applied to $E', F' \in \mathcal{A}$, there exists $N_0 \in \mathbb{N}$ such that $D_N(E', F') = 0$ for every $N \geq N_0$. Therefore
\begin{align*}
\limsup_{N \to \infty} D_N(E, F) \leq 4\varepsilon.
\end{align*}
Since $\varepsilon > 0$ was arbitrary and $D_N(E, F) \geq 0$,
\begin{align*}
\lim_{N \to \infty} \left|\mu(\sigma^{-N}E \cap F) - \mu(E)\mu(F)\right| = 0.
\end{align*}[/step]