Evaluating $\int_{-\infty}^\infty\frac{\cos(2x)}{x^2+4}\:\mathrm{d}x$ As stated in the title, I want to evaluate the integral
$$I=\int_{-\infty}^\infty\frac{\cos(2x)}{x^2+4}\:\mathrm{d}x$$
I'm pretty sure it evaluates to
$$\frac{\pi}{2e^4}$$
But I'm not sure how to evaluate it.
I have read an Instagram post where 3 different methods are provided for proving that
\begin{equation}
I(t)=\int_{-\infty}^\infty\frac{\cos(tx)}{x^2+1}\:\mathrm{d}x=\frac{\pi}{e^t}
\end{equation}
and I think similar logic can be applied here, but I am not sure how yet.
One of the methods mentioned in the post uses laplace transform to prove it but it's a little bit long. I'm wondering if there's any elegant method for evaluating $I$
I encountered this integral when I tried to solve this integral from one of the members of the Instagram math community.
$$\omega=\int_0^{\infty}\frac{x^2-4}{x^2+4}\:\frac{\sin 2x}{x}\mathrm{d}x$$
I first split the integral, used a property of laplace transform and some properties of the sine integral then used integration by parts and got to here
$$\omega=\frac{\pi}{2}-\left(2\int_{-\infty}^\infty\frac{\cos(2s)}{s^2+4}\:\mathrm{d}s+\pi\right)$$
Thank you so much for your help and attention! (BTW I'm not so proficient in complex analysis so I would prefer a solution without one :P)

Too long for a comment, just reducing this case, to the case you know from instagram and some notes.
For $t \in \mathbb R,a > 0$ let:$$I(t,a) = \int_{-\infty}^\infty \frac{\cos(tx)}{x^2+a^2}dx $$
Note that it converges for every $t \in \mathbb R,a> 0$. Taking substitution $x=ay$, $dx=ady$ we get:
$$ I(t,a) = \int_{-\infty}^\infty \frac{\cos(aty)}{a^2y^2+a^2}(ady) = \frac{1}{a} \int_{-\infty}^\infty \frac{\cos(aty)}{y^2+1}dy = \frac{1}{a} \cdot I(ta,1) $$
So it boils down to evaluate $$I(s) := I(s,1) = \int_{-\infty}^\infty \frac{\cos(sx)}{x^2+1}dx$$
There are many ways of calculating it. Probably the easiest one might be either complex analysis or noticing it is almost the characteristic function of Cauchy distribution (it does not require $"$complex analysis$"$ (even though there are complex numbers under integral sign) to calculate, however it would be a long road for you if you're not familiar with notion of characteristic function and inverse formula for them. It can be calculated via taking derivative and manipulations, however one need to be a bit careful with showing we can go with derivative under integral sign, since $\frac{d}{ds}(\frac{\cos(sx)}{x^2+1}) = -\frac{x\sin(xs)}{x^2+1}$ and integral of the latter does not converge when treated as improper lebesgue integral on the whole line (so dominated convergence theorem cannot be applied straightforward). However, it does converge when treated as improper Riemann integral or limit of proper Lebesgue integrals, so it actually makes sense. We'll proceed by a bit different way. It is amazing what a substitution can do:
Let $s>0$ and take substitution: $y=sx, dy = sdx$, then:
$$ I(s) = \int_{-\infty}^\infty \frac{s\cos(y)}{s^2+y^2}dy $$
Derivative of function under integral (with respect to $s$) yields $\frac{\cos(y)(s^2+y^2) - 2s^2\cos(y)}{(s^2+y^2)^2} = \frac{cos(y)(y^2-s^2)}{(s^2+y^2)^2}$, which is integrable on the whole line, treated as Lebesgue improper integral, so dominated convergence theorem allows us to go with derivative under integral. Taking integral one more time (justification is the same) we get:
$$ \frac{d^2}{ds^2} I(s) = \int_{-\infty}^\infty \cos(y) \cdot (\frac{d^2}{ds^2} \frac{s}{y^2+s^2}) dy = - \int_{-\infty}^\infty \cos(y) \cdot (\frac{d^2}{dy^2} \frac{s}{y^2+s^2})dy$$
Integrating by parts gives us:
$$ \frac{d^2}{ds^2}I(s) = -\frac{d}{dy}(\frac{s}{y^2+s^2})\cos(y)|_{-\infty}^\infty -\int_{-\infty}^\infty \sin(y) \frac{d}{dy}(\frac{s}{y^2+s^2})dy $$
I'll leave calculations so that boundary terms goes to zero. One time again:
$$ \frac{d^2}{ds^2}I(s) = - \sin(y)\frac{s}{y^2+s^2}|_{-\infty}^\infty + \int_{-\infty}^\infty \frac{s\cos(y)}{y^2+s^2}dy = I(s)$$
Hence the general solution is $I(s) = Ae^s + Be^{-s}$ for some constants $A,B$. We can find them by letting $s \to \infty$ and $s \to 0^+$.
Indeed, back to first form of integral, by dominated convergence $\lim_{s \to 0^+} I(s) =  \lim_{s \to 0^+} \frac{\cos(sx)}{x^2+1}dx = \pi$, so $A + B = \pi$.
To justify limit as $s \to \infty$ we use integration by parts with $cos(sx)$ and $\frac{1}{x^2+1}$, getting:
$$ \lim_{s \to \infty}I(s) = \lim_{s \to \infty} \int_{-\infty}^\infty \frac{\cos(sx)}{x^2+1}dx = \lim_{s \to \infty} \int_{-\infty}^\infty \frac{2x\sin(sx)}{s(x^2+1)^2}dx $$ which tends to zero, since we can bound $|\sin(sx)| \le 1$, and we have something that tends to zero left.
But if $Ae^s + Be^{-s} \to 0$ as $s \to \infty$, then $A=0$. So $B=\pi$.
And we get $I(s) = \pi e^{-s}$ for $s > 0$ and by symetry and easy calculation for $s=0$, we get for any $s \in \mathbb R$:
$$I(s) = \pi e^{-|s|}$$
This means that $$ I(t,a) = \frac{1}{a} I(ta,1) = \frac{\pi}{a}e^{-|ta|}$$
So your integral is equal indeed, to $\frac{\pi}{2e^4}$
