Problem: Let $X=[X_1,X_2]^T$ be a real-valued Gaussian random vector with mean and co-variance
$$\mu:=[1, 2]^T\,\,\,,\,\,\, \Lambda := \begin{pmatrix}2&-1\\-1&3\end{pmatrix}$$
Find the probability density $Y = AX + B$, where $A := [1, 1]$ and $B := 1$
What I did: I acquired the probability density of $X$, which is $$f_X(x)=\frac{1}{2\sqrt{5}\pi}\exp\left\{-\frac{3x_1^2+2x_2^2-10x_1-10x_2+2x_1x_2+15}{10}\right \}$$ (hopefully this form is correct). Since I was asked to find the probability density of $Y=X_1+X_2+1$, I tried to transform the $3x_1^2+2x_2^2-10x_1-10x_2+2x_1x_2+15$ part into a polynomial including multiple $(x_1+x_2+1)$ terms. Even though this sounds fair, it is not easy to find such a polynomial.