Después de un cambio de variables xn=y (suponemos n>0) la generalización de la integral lee
F(n)=1n∫∞01ysin(y)cos(y)+cosh(y)=F(1)n
por lo que es suficiente para evaluar la integral para n=1.
Esto puede hacerse mediante el teorema de los residuos. El uso de la paridad escribimos
F(1)=12ℑ[P∫Reiyy(cos(y)+cosh(y))]
el integrando muy bien converge en la mitad superior del plano complejo y tiene polos en z0=0 zk,±=(±1+i)(π2+2πk) en este dominio.
Hoy día no es demasiado difícil mostrar que Res(zk+)+Res(zk−)=0 cualquier k∈N.
De las consideraciones anteriores (el factor de πi =2πi/2 es debido al hecho de que la singularidad es en el contorno de integración), llegamos a la
F(1)=12ℑ[πiRes(z0)]=π4
o
F(n)=F(1)n=π4n
poner a n=π4 obtenemos I=1 como se espera de OP consideraciones