# Multivariate Normal Difference Distribution

Since the distribution of a difference of two normally distributed variates X and Y with means and variances $(\mu_x,\sigma_x^2)$ and $(\mu_y,\sigma_y^2)$ respectively is given by another normal distribution with mean $\mu_x-\mu_y$ and variance $\sigma_x^2+\sigma_y^2$ http://mathworld.wolfram.com/NormalDifferenceDistribution.html.

Assuming that the distributions are independent, can we find a distribution of a difference of two multivariate normally distributed variates?

If two multivariate normal random variables $X\sim\mathcal N(\mu_X,\Sigma_X)$ and $Y\sim\mathcal N(\mu_Y,\Sigma_Y)$ are independent (and of the same dimension), then their sum is still normal, and you can sum the mean and variance directly: $$X+Y\sim\mathcal N(\mu_X+\mu_Y,\Sigma_X+\Sigma_Y).$$
The easiest way to see this is with the expression of the characteristic function, as can be found on the Wikipedia article for example: $$\phi_X(u)=\exp(iu'\mu_X-u'\Sigma_Xu).$$
Now back to your original question, it suffices to see that $-Y$ has the same variance as $Y$, but with opposite mean to conclude: $$-Y\sim \mathcal N(-\mu_Y,\Sigma_Y),$$ so indeed $$X-Y\sim\mathcal N(\mu_X-\mu_Y,\Sigma_X+\Sigma_Y).$$