Verify integration of a sine function

[yungman]

I use part of HS-Scientist derivation in another post thanks to his detail derivation. I want to solving $\int_0^{\pi} \sin(x sin(\theta)) d\theta$


\int_0^{\pi} \sin^m(\theta) d\theta=-\frac{1}{m}\left[sin^{m-1}(x)cos(\theta)\right]_0^{\pi}+\frac{m-1}{m} \int sin^{m-2}(\theta) d\theta= \frac{m-1}{m} \int sin^{m-2}(\theta) d\theta
\hbox{As }\;\left[sin^{m-1}(x)cos(\theta)\right]_0^{\pi}=0


Let $m=2n+1$
\Rightarrow\; \int_0^{\pi} \sin^m(\theta) d\theta=\int_0^{\pi} \sin^{2n+1}(\theta) d\theta=\frac{2n}{2n+1}\int_0^{\pi} \sin^{2n-1}\theta d\theta=\frac{2\cdot 4\cdot 6\cdot \cdot \cdot \cdot (2n)}{3\cdot 5\cdot 7\cdot \cdot \cdot (2n+1))}\int_0^{\pi}\sin\theta d\theta
\Rightarrow\; \int_0^{\pi} \sin^{2n+1}(\theta) d\theta=\frac{2\cdot 4\cdot 6\cdot \cdot \cdot \cdot (2n)}{3\cdot 5\cdot 7\cdot \cdot \cdot (2n+1))}\int_0^{\pi}\sin\theta d\theta =2\frac{2\cdot 4\cdot 6\cdot \cdot \cdot \cdot (2n)}{3\cdot 5\cdot 7\cdot \cdot \cdot (2n+1))}
\sin x =\sum_0^{\infty}\frac {(-1)^n x^{2n+1}}{(2n+1)!}


\Rightarrow\;\int_0^\pi \sin(x sin(\theta)) d\theta=\sum_0^{\infty} \frac {(-1)^n x^{2n+1}}{(2n+1)!} \int_0^{\pi}\sin^{2n+1}\theta d\theta\;=\;\left[\sum_0^{\infty} \frac {(-1)^n x^{2n+1}}{(2n+1)!}\right]\left[2\frac{2\cdot 4\cdot 6\cdot \cdot \cdot \cdot (2n)}{3\cdot 5\cdot 7\cdot \cdot \cdot (2n+1))}\right]

\Rightarrow\;\int_0^\pi \sin(x sin(\theta)) d\theta \;=\;\sum_0^{\infty} \frac {(-1)^n 2x^{2n+1}}{(2n+1)!} \frac{2\cdot 4\cdot 6\cdot \cdot \cdot \cdot (2n)}{3\cdot 5\cdot 7\cdot \cdot \cdot (2n+1))} \;=\;\;\sum_0^{\infty} \frac {(-1)^n 2x^{2n+1}}{3^2 \cdot 5^2\cdot 7^2\cdot \cdot \cdot \cdot (2n+1)^2}


I want to verify I am correct in my derivation.

Thanks

 

[Ray Vickson]

I use part of HS-Scientist derivation in another post thanks to his detail derivation. I want to solving $\int_0^{\pi} \sin(x sin(\theta)) d\theta$


\int_0^{\pi} \sin^m(\theta) d\theta=-\frac{1}{m}\left[sin^{m-1}(x)cos(\theta)\right]_0^{\pi}+\frac{m-1}{m} \int sin^{m-2}(\theta) d\theta= \frac{m-1}{m} \int sin^{m-2}(\theta) d\theta
\hbox{As }\;\left[sin^{m-1}(x)cos(\theta)\right]_0^{\pi}=0


Let $m=2n+1$
\Rightarrow\; \int_0^{\pi} \sin^m(\theta) d\theta=\int_0^{\pi} \sin^{2n+1}(\theta) d\theta=\frac{2n}{2n+1}\int_0^{\pi} \sin^{2n-1}\theta d\theta=\frac{2\cdot 4\cdot 6\cdot \cdot \cdot \cdot (2n)}{3\cdot 5\cdot 7\cdot \cdot \cdot (2n+1))}\int_0^{\pi}\sin\theta d\theta
\Rightarrow\; \int_0^{\pi} \sin^{2n+1}(\theta) d\theta=\frac{2\cdot 4\cdot 6\cdot \cdot \cdot \cdot (2n)}{3\cdot 5\cdot 7\cdot \cdot \cdot (2n+1))}\int_0^{\pi}\sin\theta d\theta =\frac{2\cdot 4\cdot 6\cdot \cdot \cdot \cdot (2n)}{3\cdot 5\cdot 7\cdot \cdot \cdot (2n+1))}\pi
\sin x =\sum_0^{\infty}\frac {(-1)^n x^{2n+1}}{(2n+1)!}


\Rightarrow\;\int_0^\pi \sin(x sin(\theta)) d\theta=\sum_0^{\infty} \frac {(-1)^n x^{2n+1}}{(2n+1)!} \int_0^{\pi}\sin^{2n+1}\theta d\theta\;=\;\left[\sum_0^{\infty} \frac {(-1)^n x^{2n+1}}{(2n+1)!}\right]\left[\frac{2\cdot 4\cdot 6\cdot \cdot \cdot \cdot (2n)}{3\cdot 5\cdot 7\cdot \cdot \cdot (2n+1))}\pi\right]

\Rightarrow\;\int_0^\pi \sin(x sin(\theta)) d\theta \;=\;\sum_0^{\infty} \frac {(-1)^n x^{2n+1}}{(2n+1)!} \frac{2\cdot 4\cdot 6\cdot \cdot \cdot \cdot (2n)}{3\cdot 5\cdot 7\cdot \cdot \cdot (2n+1))}\pi \;=\;\;\sum_0^{\infty} \frac {(-1)^n x^{2n+1}}{2^2 \cdot 4^2\cdot 6^2\cdot \cdot \cdot \cdot (2n)^2}\pi


I want to verify I am correct in my derivation.

Thanks

Here is what I get using Maple 11:
f:=sin(x*sin(t));
f := sin(x sin(t)) <--- echo of input
J:=int(f,t =0..Pi);
J := Pi StruveH(0, x) <--- the answer

The function 'StruveH' is non-elementary and is related to solutions of non-homogeneous Bessel differential equations.

We can look at the series expansion:

series(J,x=0,10); <--- get the first 10 terms of the series expn about x = 0:

J = 2*x-(2/9)*x^3+(2/225)*x^5-(2/11025)*x^7+(2/893025)*x^9+O(x^10)

Note that the terms do NOT contain π, so this does not match your solution.

 

[yungman]

Thanks for the reply. I double checked, I made two mistakes.

1) my in cut and paste in the very last step was wrong. I since corrected the original post already.

2)My other mistake was
\int_0^{\pi}\sin\theta d\theta=-\cos\theta|_0^{\pi} =-[\cos(\pi)-cos(0)] =-[-1-1] =2\;\hbox {not }\;\piI did the substitution. I got

n=0\;\Rightarrow \;\frac {2 x}{1^2}=2 x
n=1\;\Rightarrow \;-\frac {2 x^3}{3^2}=-\frac {2 x^3}{9}
n=2\;\Rightarrow \;\frac {2 x^5}{3^2\cdot 5^2}=\frac {2 x^5}{225}
n=3\;\Rightarrow \;-\frac {2 x^7}{3^2\cdot 5^2\cdot 7^2}=-\frac {2 x^7}{11025}
n=4\;\Rightarrow \;\frac {2 x^9}{3^2\cdot 5^2\cdot 7^2\cdot 9^2}=\frac {2 x^9} {893925}

Now it fits your result.

Thanks