[step:Verify joint Gaussianity and the covariance formula $\mathbb{E}[B_s B_t] = s \wedge t$]
**Joint Gaussianity.** For any $0 \leq t_1 \leq \cdots \leq t_n$ and $c_1, \ldots, c_n \in \mathbb{R}$, the linear combination $\sum_{k=1}^n c_k B_{t_k} = \sum_{k=1}^n c_k I(\mathbb{1}_{[0,t_k]}) = I\!\left(\sum_{k=1}^n c_k \mathbb{1}_{[0,t_k]}\right)$, using linearity of $I$. Since $\sum_{k=1}^n c_k \mathbb{1}_{[0,t_k]} \in L^2(\mathbb{R}_+)$, [Theorem 2067](/theorems/2067) gives that $I\!\left(\sum c_k \mathbb{1}_{[0,t_k]}\right) \sim N(0, \|\sum c_k \mathbb{1}_{[0,t_k]}\|_{L^2}^2)$. Hence every linear combination of $(B_{t_1}, \ldots, B_{t_n})$ is a centered Gaussian, so the vector is jointly Gaussian.
**Covariance.** For $s, t \geq 0$, the isometry property of [Theorem 2067](/theorems/2067) gives:
\begin{align*}
\mathbb{E}[B_s B_t] = \mathbb{E}[I(\mathbb{1}_{[0,s]}) \cdot I(\mathbb{1}_{[0,t]})] = (\mathbb{1}_{[0,s]}, \mathbb{1}_{[0,t]})_{L^2(\mathbb{R}_+)}.
\end{align*}
Computing the inner product:
\begin{align*}
(\mathbb{1}_{[0,s]}, \mathbb{1}_{[0,t]})_{L^2(\mathbb{R}_+)} = \int_{\mathbb{R}_+} \mathbb{1}_{[0,s]}(u) \, \mathbb{1}_{[0,t]}(u) \, d\mathcal{L}^1(u) = \int_{\mathbb{R}_+} \mathbb{1}_{[0,s] \cap [0,t]}(u) \, d\mathcal{L}^1(u) = \mathcal{L}^1([0, s \wedge t]) = s \wedge t.
\end{align*}
Therefore $\mathbb{E}[B_s B_t] = s \wedge t$.
This completes the verification that $(B_t)_{t \geq 0}$ satisfies properties (1)-(4) of Brownian motion, establishing that it is a Brownian motion in all respects except possibly the continuity of sample paths.
[/step]