I'm reading through Andrew Ng Lecture Notes for CS229 and he makes the statement that, for a random variable $Z \in \mathbb{R}^{n}$,
\begin{align} Cov(Z) &= E[(Z - E[Z])(Z - E[Z])^{T}]\\ &= E[ZZ^{T}] - E[Z](E[Z])^{T} \end{align}
I've started with: \begin{align} Cov(Z) &= E[(Z - E[Z])(Z - E[Z])^{T}]\\ &= E[(ZZ^{T}) - 2(E[Z]Z^T) + E[Z](E[Z])^{T}] \end{align}
But I'm not sure how to continue.
Just derive as an extension of the single variable case. Derive Variance(X), where X is a random variable, as a function of E(X), you get E($X^2$) - $(E(X))^2$. Just extend to the multivariate case, while remembering some of the other hints listed in above comments.