Let $Z\stackrel{d}{=}\mu + \xi A U$ where $\xi$ is a non specified (one dimensional) random variable with $\mathbb{E}(\xi^2)=n$ independent of $U$, $U$ is uniformly distributed on the $n$-dimensional sphere $\mathbb{S^{n-1}}$ and A is a $n\times n$ matrix. $Z$ and $\mu$ are vectors of length $n$.
Denote with $A_i$ the i-th row of $A$ and with $Z_i, \mu_i, U_i$ the i-th entry.
Also I know that $\mathbb{E}(U_i)=0$ for all $i$ and the covariance matrix of U equals $\frac{1}{n}\mathbb{I}$, where $\mathbb{I}$ denotes the identity matrix.
Now I want to prove that the convariance of $Z$ equals $A A^T$.
It is obvious that $\mathbb{E}(Z_i)=\mu_i$. Now I calculate $\mathbb{E}(Z_iZ_j)$ and $\mathbb{E}(Z_i^2)$. For this reason we have
$$\mathbb{E}(Z_iZ_j)=\mathbb{E}((\mu_i+\xi A_i U)(\mu_j+\xi A_j U))$$ or
$$\mathbb{E}(Z_i^2)=\mathbb{E}((\mu_i+\xi A_i U)(\mu_i+\xi A_i U))$$
but now I'm stuck what to do next.
Of course I can solve the two brackets but then I get something like
$$\mathbb{E}(\xi A_i U\xi A_i U)$$ I can't just switch the vectors or can I??