Solution
= Solution
The estimator is an affine transformation of the <multivariate normal distribution>. Its mean is $\beta$, and direct covariance calculation gives
$$
\boxed{
\widehat\beta\sim
N_p\left(\beta,
\sigma^2(X^T\Sigma_0^{-1}X)^{-1}\right)}.
$$