まず、∣sint∣ の周期性が π であることから、積分区間を [0,π] に縮小し、積分を2倍します。 ∫02πe−∣sint∣dt=2∫0πe−∣sint∣dt 次に、区間 [0,π] で sint≥0 であるため、∣sint∣=sint となります。 したがって、
2∫0πe−∣sint∣dt=2∫0πe−sintdt 次に、t=π−u と変数変換すると、dt=−du となり、積分範囲は u:π→0 となります。 ∫0πe−sintdt=∫π0e−sin(π−u)(−du)=∫0πe−sin(π−u)du sin(π−u)=sinu であるため、 ∫0πe−sin(π−u)du=∫0πe−sinudu これは元の積分と同じです。したがって、この変数変換は直接的な計算には役立ちません。
ベッセル関数を使うことを考えます。
∫0πezcosθdθ=πI0(z) ここで、I0(z) は第一種0次の変形ベッセル関数です。 e−sint を展開することを考えます。 e−sint=∑n=0∞n!(−sint)n ここで、
∫0πsinntdt={02nπ(n/2n)n が奇数のときn が偶数のとき ∫0πe−sintdt=∫0π∑n=0∞n!(−sint)ndt=∑n=0∞n!(−1)n∫0π(sint)ndt ∫0πe−sintdt=∑k=0∞(2k)!(−1)2k∫0π(sint)2kdt=∑k=0∞(2k)!122kπ(k2k)=π∑k=0∞(2k)!122k(k!)2(2k)!=π∑k=0∞4k(k!)21=πI0(1) したがって、
∫02πe−∣sint∣dt=2∫0πe−sintdt=2πI0(1) I0(1)≈1.26606587752 2πI0(1)≈2π(1.26606587752)≈7.95774715459