I think that there are two ways to tackle this problem by exploiting the analyticity of the integrand ecosθcos(sinθ). We can calculate it in the real domain by Taylor series or in the complex domain by Cauchy's integral formula.
Solution by Taylor Series
First, let us introduce the function
f(r)=∫02πercosθcos(rsinθ)dθ, so we have the equation
I=f(1) . By Taylor expansion, we have
f(r)=n=0∑∞n!f(n)(0)rn. It is easy to show that all the coefficients of this Taylor expansion vanish except for the first one. Namely, f(0)=∫02πe0cosθcos(0sinθ)dθ=2π. The n-th derivatives of f at 0 must be in the form of
f(n)(0)=∫02πk even ∧0≤k≤n∑(−1)2k(kn)cosn−kθsinkθdθ=∫02πcos(nθ)dθ=0. Hence, we have f(r)=2π=I.
Solution by Cauchy's Integral Formula
First, let us introduce the same function
f(r)=∫02πercosθcos(rsinθ)dθ. The integrand can be converted to a complex exponential,
ercosθcos(rsinθ)=ℜ{ercosθcos(rsinθ)+isin(rsinθ)}=ℜ{ercosθ+irsinθ}=ℜ{ereiθ}. Cauchy's Integral formula at 0 is
g(0)=2πi1∮γzg(z)dz. If we take g(z)=ez, and γ:reiθ, we have
g(0)=2πi1∫02πreiθg(reiθ)ireiθdθ=2π1∫02πg(reiθ)dθ. Note that we have
ℜ{g(reiθ)}=ℜ{ereiθ}=ercosθcos(rsinθ). Hence,
g(0)g(0)f(r)I=e0=1=2π1f(r)=2π=f(1)=2π.