xsinxの積分・xcosxの積分

xsinxの積分

$x\sin x$ の積分公式です。導出することは難しくありませんが、できれば覚えてしまって、少しでも時間を節約したいところです：
 $$\begin{array}{}\text{(1)}& \int x\sin xdx=\sin x-x\cos x+C\end{array}$$
【証明】$x\sin x$ の積分は部分積分によっても求められますが、$x\cos x$ を微分するほうがより簡単です：
 $$\begin{array}{}\text{(2)}& (x\cos x{)}^{\prime }=\cos x-x\sin x=(\sin x{)}^{\prime }–x\sin x\end{array}$$
右辺に $x\sin x$ が現れています。この式を整理すると
 $$\begin{array}{}\text{(3)}& x\sin x=(\sin x{)}^{\prime }–(x\cos x{)}^{\prime }\end{array}$$
両辺を積分すると、
 $$\begin{array}{}\text{(4)}& \int x\sin xdx=\sin x-x\cos x+C\end{array}$$
三角関数の変数が $mx$ になっているときは $t=mx$ と変換します：
 $$\begin{array}{}\text{(5)}& \begin{array}{rl}\int x\sin mxdx& =\frac{1}{m^2}\int t\sin tdt\\ & =\frac{1}{m^2}(\sin t-t\cos t)+C\\ & =\frac{1}{m^2}\sin mx-\frac{x}{m}\cos mx+C\end{array}\end{array}$$
余力があれば、この $m$ を含んだ公式を覚えてしまいましょう。

【Excel】xsinxのグラフ

Excel で描いた $f(x)=x\sin x$ のグラフを載せておきます。
 
(≫ グラフデータの入った Workbook をダウンロード)
 

 
奇関数×奇関数なので、$x\sin x$ は偶関数です。
$x\sin x$ は定積分に特徴があります。公式 (1) より、
 $$\begin{array}{}\text{(6)}& {\int }_0^{\pi }x\sin xdx={[\sin x-x\cos x]}_0^{\pi }=\pi \end{array}$$
を得ます。これは曲線 $y=x\sin x$ と線分 $[0,\pi ]$ で囲まれる領域の面積が円周率 $\pi$ (すなわち半径 $1$ の円の面積) に等しいことを意味します。

【VBA】xsinx の積分を計算するマクロ

$x\sin x$ を積分する VBA マクロを載せておきます。積分公式 (1) を使えば、以下のように簡単に実装できます。

'[VBA] xsinxの積分を返す関数
Function XSINX_INTEG(a As Double, b As Double)
    Dim val_a As Double, val_b As Double
    val_a = Sin(a) - a * Cos(a)
    val_b = Sin(b) - b * Cos(b)
    XSINX_INTEG = val_b - val_a
End Function

引数 a には積分の下端、引数 b には積分の上端を渡します。たとえば、ワークシートに「=XSINX_INTEG(0,PI())」と入力すると、円周率に近い値「3.141592654 …」が返ります。
 
シンプソンの公式を使って積分値を得るマクロは以下のようになります。

'xsinxの値を返すマクロ
Function XSINX(x)
    XSINX = x * Sin(x)
End Function

'シンプソンの公式でxsinxを数値積分するマクロ
Function XSINX_SIMPSON(a As Double, b As Double, n As Long)

    Dim h As Double, s As Double
    Dim k As Integer

    s = 0
    h = (b - a) / (2 * n)

    For k = 0 To n - 1
        s = s + XSINX(a + 2 * k * h) _
            + 4 * XSINX(a + (2 * k + 1) * h) _
            + XSINX(a + (2 * k + 2) * h)
    Next k

    XSINX_SIMPSON = s * h / 3

End Function

引数 a には積分の下端、引数 b には積分の上端を渡します。n は分割数で、大きな値を入れるほど精度が上がりますが処理時間も増えます。たとえば、ワークシートに「=XSINX_SIMPSON(0,PI(),100)」を入力すると、円周率に近い値「3.141592655 …」が返ります。

xcosxの積分

$x\cos x$ の不定積分公式です。
 $$\begin{array}{}\text{(7)}& \int x\cos xdx=\cos x+x\sin x+C\end{array}$$
【証明】 $x\sin x$ を微分すると、
 $$\begin{array}{}\text{(8)}& x\sin x^{\prime }=\sin x+x\cos x=(-\cos x{)}^{\prime }+x\cos x\end{array}$$
となるので、式を整理すると、
 $$\begin{array}{}\text{(9)}& x\cos x=(\cos x{)}^{\prime }+(x\sin x{)}^{\prime }\end{array}$$
両辺を積分すると、$x\cos x$ の積分を得ます：
 $$\begin{array}{}\text{(10)}& \int x\cos xdx=\cos x+x\sin x+C\end{array}$$
三角関数の変数が $mx$ になっているときは $t=mx$ と変換します：
 $$\begin{array}{}\text{(11)}& \begin{array}{rl}\int x\cos mxdx& =\frac{1}{m^2}\int t\cos tdt\\ & =\frac{1}{m^2}(\cos t+t\sin t)+C\\ & =\frac{1}{m^2}\cos mx+\frac{x}{m}\sin mx+C\end{array}\end{array}$$

【Excel】xcosx のグラフ

下の図は Excel で描いた $f(x)=x\cos x$ のグラフです：
 

 
奇関数×偶関数なので、$x\cos x$ は奇関数です。

【VBA】xcosxの積分を計算するマクロ

公式 (10) を使って $x\cos x$ の積分を計算するマクロを載せておきます。

'xcosxの定積分の値を返すユーザー定義関数
Function XCOSX_INTEG(a As Double, b As Double)
    Dim val_a As Double, val_b As Double
    val_a = Sin(a) - a * Cos(a)
    val_b = Sin(b) - b * Cos(b)
    XCOSX_INTEG = val_b - val_a
End Function

たとえば、ワークシートの適当なセルに「=XCOSX_INTEG(0,2*PI())」を入力すると「-6.283185307」が返ります。

x^2 sinxの積分・x^2cosx の積分

公式を使って、$x^2\sin x$ の積分を求めることができます。$x^2\cos x$ を微分すると
 $$\begin{array}{}\text{(12)}& (x^2\cos x{)}^{\prime }=2x\cos x-x^2\sin x\end{array}$$
となって、右辺に $x^2\sin x$ の項が現れます。移項して整理すると
 $$\begin{array}{}\text{(13)}& x^2\sin x=2\cos xdx-(x^2\cos x{)}^{\prime }\end{array}$$
両辺を積分すると
 $$\begin{array}{}\text{(14)}& \int x^2\sin xdx=\int 2\cos xdx-x^2\cos x\end{array}$$
ここで公式 (10) を使って右辺第 1 項の積分を実行すると、
 $$\begin{array}{}\text{(15)}& \int x^2\sin xdx=2(x\sin x+\cos x)-x^2\cos x+C\end{array}$$
を得ます。整理すると
 $$\begin{array}{}\text{(16)}& \int x^2\sin xdx=(2-x^2)\cos x+2x\sin x+C\end{array}$$
となります。同様に $x^2\sin x$ を微分して公式 (1) を使うと、
 $$\begin{array}{}\text{(17)}& \int x^2\cos xdx=(x^2-2)\sin x+2x\cos x+C\end{array}$$
が得られます。

f(x)sinxの積分・f(x)cosxの積分

次の公式は無理に覚える必要はありませんが、どこかに控えておくと、証明に役立つことがあります。$f^n(x)$ は $n$ 階導関数を表す記号です。
 $$\begin{array}{}\text{(18)}& \int f(x)\sin xdx=\sum _{k=0}^{\mathrm{\infty }}(-1{)}^k[f^{(2k+1)}(x)\sin x-f^{(2k)}(x)\cos x]\end{array}$$
積分定数 $C$ は省略してあります。

【証明】単純に部分積分を実行します。
 $$\begin{array}{}\text{(19)}& \begin{array}{rl}I(x)=\int f(x)\sin xdx& =-f(x)\cos x+\int f^{(1)}(x)\cos xdx\\ & =-f(x)\cos x+f^{(1)}(x)\sin x-\int f^{(2)}(x)\sin xdx\end{array}\end{array}$$
最後の積分は $I(x)$ の $f(x)$ のところを $f”(x)$ に置き換えただけなので、あとは繰り返しになります。
 $$\begin{array}{}\text{(20)}& \begin{array}{rl}I(x)=& -f(x)\cos x+f^{(1)}(x)sinx\\ & -\{-f^{(2)}(x)\cos x+f^{(3)}(x)\sin x-\int f^{(4)}(x)\sin xdx\}\\ & =[f^{(1)}(x)-f^{(3)}(x)+f^{(5)}(x)+\cdots ]\sin x\\ & +[-f^{(0)}(x)+f^{(2)}(x)-f^{(4)}(x)+\cdots ]\cos x\\ & =\sum _{k=0}^{\mathrm{\infty }}[(-1{)}^k{f}^{(2k+1)}(x)\sin x+(-1{)}^{k+1}{f}^{(2k)}(x)\cos x]\\ & =\sum _{k=0}^{\mathrm{\infty }}(-1{)}^k[f^{(2k+1)}(x)\sin x-f^{(2k)}(x)\cos x]\end{array}\end{array}$$
和は $n$ 階導関数 $f^n(x)$ が $0$ になるまでとればよいので、$f(x)$ が三角関数のように微分に対して循環しない限りは有限項の級数となります。ですから実用的に考えると、
 $$\begin{array}{}\text{(21)}& \int f(x)\sin xdx=[f^{(1)}(x)-f^{(3)}(x)]\sin x+[f^{(2)}(x)-f^{(0)}(x)]\cos x\end{array}$$
ぐらいまでを覚えておけば充分です。$f(x)=x,f(x)=x^2$ について直ちに次の公式が得られます。
 $$\begin{array}{}\text{(22)}& \begin{array}{rl}& \int x\sin xdx=-x\cos x-\sin x+C\\ & \int x^2\sin xdx=(2-x^2)\cos x+2x\sin x+C\end{array}\end{array}$$　
同じようにして、
 $$\begin{array}{}\text{(23)}& \int f(x)\cos xdx=\sum _{k=0}^{\mathrm{\infty }}(-1{)}^k[f^{(2k)}(x)\sin x+f^{(2k+1)}(x)\cos x]\end{array}$$
という公式も得られます。$f^{(3)}(x)$ までの導関数を使って簡略的に書くと
 $$\begin{array}{}\text{(24)}& \int f(x)\cos xdx=[f^{(0)}(x)-f^{(2)}(x)]\sin x+[f^{(1)}(x)-f^{(3)}(x)]\cos x\end{array}$$
となります。たとえば、$f(x)=x,f(x)=x^2$ とすると
 $$\begin{array}{}\text{(25)}& \begin{array}{rl}& \int x\cos xdx=x\sin x+\cos x+C\\ & \int x^2\cos xdx=(x^2-2)\sin x+2x\cos x+C\end{array}\end{array}$$
という公式を得ます。