Let me remind first the construction of Brownian motion.
Fix a vector $x \in \mathbb{R}^n$ and define $p(t,x,y) := (2\pi t )^{-\frac{n}{2}} \cdot \exp{\left( - \frac{|x-y|^2}{2t} \right)},$ for $y \in \mathbb{R}^n$, $t >0$.
Then for $0 \leq t_1 \leq t_2 \leq \ldots \leq t_k$ we define a measure $\nu_{t_1, \ldots, t_k}$ on $\mathbb{R}^{nk}$ by $\nu_{t_1, \ldots, t_k}(F_1 \times \ldots \times F_k)= \int \limits_{F_1 \times \ldots \times F_k}^{} p(t_1, x, x_1)\prod_{j=1}^{k-1}p(t_{j+1}-t_j,x_{j}, x_{j+1}) dx_1\ldots dx_k,$ where the following conventions are used $dy=dy_1\ldots dy_k$ for Lebesgue measure and $p(0, x, y)dy=\delta_x(y)$ (the unit point mass at $x$).
By Kolmogorov's extension theorem applied to the probability measures $\nu_{t_1, \ldots, t_k}$ (which easily satisfies all the assumtions of that theorem) there exists a probability space $(\Omega, \mathcal{F}, P^x)$ and a stochastic process $\{B_t\}_{t \geq 0}$ on $\Omega$ such that the finite distributions of $B_t$ are given by $ (\star) \ \ P^x(B_{t_1} \in F_1, \ldots, B_{t_k} \in F_k)= \int \limits_{F_1 \times \ldots \times F_k}^{} p(t_1, x, x_1)\prod_{j=1}^{k-1}p(t_{j+1}-t_j,x_{j}, x_{j+1}) dx_1\ldots dx_k.$
Problem
I want to show that for the random variable $Z = (B_{t_1}, \ldots, B_{t_k}) \in \mathbb{R}^{nk}$ there exist a vector $M \in \mathbb{R}^{nk}$ and a non-negative definite matrix $C \in \text{M}_{nk \times nk}(\mathbb{R})$ such that $E^x\left[\exp\left(i\left \right)\right]= \exp\left( -\frac{1}{2} \left
I was trying to calculate it explicitly by writing the left hand side by its definition and applying the ($\star$) formula, but the integral which I've got is not so nice and I cannot conclude what I want to.
This exercise is nothing else, but showing that $B_t$ is a Gaussian process (so it is not hard to guess what should be $M$ and $C$ in this case).
I hope you know some tricks how to compute such an integral in an easy way. Thanks in advance for any help.