Podemos escribir su integral como
$$I=\int_0^\infty \int_{\alpha x}^\infty \int_{\beta x}^\infty \frac{\cos u \cos v}{uv } du\: dv\: dx,$$
o, lo que es lo mismo,
$$\int_0^\infty \int_{x}^\infty \int_{x}^\infty \frac{\cos \beta u \cos \alpha v}{uv } du\: dv\: dx.$$
Ahora la región de integración es $$\{(x,u,v): x<u, x<v\}$$
que, para un conjunto de medida cero, podemos escribir como distinto de la unión de las regiones $$\{(x,u,v): x<u<v\}$$ and $$\{(x,u,v): x<v<u\}.$$
Ahora, por ejemplo, para la segunda región que hemos
$$\int_0^\infty \int_{x}^\infty \int_{v}^\infty \square\: du\: dv\: dx = \int_0^\infty \int_{0}^u \int_{0}^v \square\: dx\: dv\: du $$
y para nosotros, esto da
$$ \int_0^\infty \int_{0}^u \int_{0}^v \frac{\cos \beta u \cos \alpha v}{uv } dx\: dv\: du = \int_0^\infty \int_{0}^u \frac{\cos \beta u \cos \alpha v}{u } dv\: du = \int_0^\infty \frac{\cos \beta u \sin \alpha u}{\alpha u } du$$
Cambiar los roles de $u$, $v$, y los roles de las $\alpha$$\beta$, y la adición de la resultante de dos integrales, obtenemos que
$$I(\alpha, \beta)=\int_0^\infty \frac{\beta \cos \beta t \sin \alpha t + \alpha \sin \beta t \cos \alpha t}{\alpha \beta t } dt.$$
Alguien puede ser capaz de tomar de allí. El parecido con la integral del seno me sugiere que la adaptación de uno de los métodos utilizados para evaluar $\int_0^\infty \frac{\sin t}{t} dt$ puede trabajar.
Edit: gracias a Didier Piau, aquí es el paso final de la solución (cita directa de su post)
A partir de la penúltima expresión en Bruno la solución, es decir,
$$I(\alpha, \beta)=\int_0^\infty \frac{\beta \cos \beta t \sin \alpha t + \alpha \sin \beta t \cos \alpha t}{\alpha \beta t } \mathrm dt,$$
vamos a utilizar las relaciones trigonométricas
$$
2\cos \beta t \sin \alpha t =\sin(\alpha+\beta)t+\sin(\alpha\beta)t,\quad
2\sin \beta t \cos \alpha t =\sin(\alpha+\beta)t+\sin(\beta\alpha)t,
$$
y el hecho de que para cada $\gamma\ne0$, el cambio de las variables de $t\to\gamma t$ rendimientos
$$
\int_0^\infty \frac{\sin \gamma t}{t} \mathrm dt=\text{sgn}(\gamma)\int_0^\infty \frac{\sen t}{t} \mathrm dt=\text{sgn}(\gamma)\frac{\pi}2,
$$
donde $\text{sgn}(\gamma)$ $+1$ si $\gamma\gt0$, $-1$ si $\gamma\lt0$, e $0$ si $\gamma=0$. Este rendimientos
$$
I(\alpha,\beta)=\frac\pi{4\alpha}(1+\text{sgn}(\alpha-\beta))+\frac\pi{4\beta}(1+\text{sgn}(\beta-\alpha)).
$$
Esta expresión es simétrica en $(\alpha,\beta)$, como debe ser. Si $\alpha\gt\beta$, el segundo término es cero y la primera es $\pi/(2\alpha)=\pi/(2\max(\alpha,\beta))$. Por último, si $\alpha=\beta$, ambos términos son $\pi/(4\alpha)=\pi/(4\beta)$ por lo tanto la suma es $\pi/(2\alpha)=\pi/(2\beta)$. Esto demuestra la fórmula deseada.