 Solution

For any $r_1,\dots ,r_n\in \mathbb{R}$, the linearity from part (a) givesalmost surely, and the right side has a normal distribution. This is the defining linear-combination criterion for a multivariate normal distribution; singular covariance matrices are allowed.

Passing to the limit in the inner products of the partial sums givesThe passage to the limit is justified by Cauchy-Schwarz inequality and $L^2$ convergence. Equivalently, the vector's characteristic function isThus both joint normality and the complete covariance matrix follow from the Hilbert-space inner product.