I have the following expression that is part of a more complicated integral \begin{equation} \cos(\sqrt{x^2+2axy+y^2}) \end{equation}
For $a=\pm 1$ the expression reduces to sum and difference angle formulas.
I would like to express it, for any $-1\le a \le 1$, in terms of $\cos x$, $\sin x$, $\cos y$, $\sin y$, or more generally the sum of separable trigonometric products (e.g. trig($x$)trig($y$) ).
Expanding the cosine into a Taylor series, the square root disappears. One can then apply the binomial / multinomial theorem, and manipulate factorials, but it is unclear to me how to reduce the resulting series to sums and products of simple trig functions.
My final goal is to evaluate the average \begin{equation} \lim_{T\rightarrow\infty}\frac{1}{(2T)^2}\int^T_{-T}\int^T_{-T}\sin(x)\sin(y)\cos(\sqrt{x^2+2axy+y^2})\ dx\ dy \end{equation}
as a function of $a$.