I'm not quite sure how to calculate this integral:
$$I=\int_0^{2\pi}dx \int_0^{\pi} e^{\sin y(\cos x-\sin x)}\sin y\,dy$$
With Mathematica I see the result is $2\pi \sqrt 2 \sinh \sqrt 2$.
Is there any special technique for calculating this integral?
Any help will be appreciated.
Thank you for these brilliant solutions using Bessel function or complex analysis.
After looking into this integral for a long time I found another "primary" solution:
Rewrite it as
$$I=\int_0^{2\pi}d\theta \int_0^{\pi} e^{\sin \varphi(\cos \theta-\sin \theta)}\sin \varphi\,d\varphi$$
By polar coordinates it becomes
$$I=\int_{u^2+v^2+w^2=1} e^{u-v}\,dS$$
since $dS=\sin\varphi d\theta d\varphi$. Then choose an orthogonal transform such that $z=\frac{u-v}{\sqrt 2}$, we get
$$I=\int_{x^2+y^2+z^2=1} e^{\sqrt 2 z}\,dS$$
Going back to polar coordinates it becomes
$$I=\int_0^{2\pi}d\theta \int_0^{\pi} e^{\sqrt 2 \cos\varphi}\sin \varphi\,d\varphi$$
And finally we get
$$I=2\pi \int_{-1}^1 e^{\sqrt 2 t}\,dt$$
I see a solution using the modified Bessel function (not in the result, however) $$I_0(z)=\sum_{n=0}^\infty\frac{1}{n!^2}\left(\frac{z}{2}\right)^{2n}=\frac1\pi\int_0^\pi e^{z\cos t}\,dt.$$ Since $\cos x-\sin x=\sqrt2\cos(x+\pi/4)$ being integrated over its period, we may drop the $\pi/4$: $$\int_0^{2\pi}e^{r(\cos x-\sin x)}\,dx=\int_0^{2\pi}e^{r\sqrt2\cos x}\,dx=2\pi I_0(r\sqrt2),$$ and using the power series, we can compute \begin{align*} \int_0^\pi I_0(a\sin y)\sin y\,dy &=\sum_{n=0}^\infty\frac{a^{2n}}{2^{2n}n!^2}\int_0^\pi\sin^{2n+1}y\,dy \\&=\sum_{n=0}^\infty\frac{a^{2n}}{2^{2n}n!^2}\mathrm{B}\left(\frac12,n+1\right) \\&=\sum_{n=0}^\infty\frac{a^{2n}}{2^{2n} n!^2}\frac{\sqrt\pi\cdot n!}{2^{-n-1}(2n+1)!!\sqrt\pi} \\&=2\sum_{n=0}^\infty\frac{a^{2n}}{(2n+1)!}=\frac2a\,\sinh a. \end{align*} This gives the result you've got.