Glasserのマスター定理 〜積分計算への活用と定理の証明〜
はじめに
Glasserのマスター定理(Glasser's master theorem)と呼ばれる定理について解説している日本語の資料が私が探した限りネット上に見当たらなかったので、自力で頑張って証明しました。
Glasserのマスター定理 とは
$a, b_1, b_2, b_3, \ldots, b_{n-1}, b_n \in \mathbb{R}, \quad a_1, a_2, a_3, \ldots, a_{n-1}, a_n \in \mathbb{R^+}$
を満たす任意の$a, b_1, b_2, b_3, \ldots, b_{n-1}, b_n, a_1, a_2, a_3, \ldots, a_{n-1}, a_n$について、
$$
\int_{-\infty}^{\infty} f(x + a - \frac{a_1}{x + b_1} - \frac{a_2}{x + b_2} - \frac{a_3}{x + b_3} - \dotsb - \frac{a_{n-1}}{x + b_{n-1}} - \frac{a_n}{x + b_n}) \mathrm{d}x = \int_{-\infty}^{\infty} f(x) \mathrm{d}x
$$
が成り立つ。
今回はこれを証明していきます。
この式だけ見てもピンと来ないと思われるかもしれませんが、この定理を使ってみれば、この定理がいかに便利かおわかりいただけると思います。
まずは使ってみよう！
皆さんにGlasserのマスター定理に慣れ、この定理の有用性をご理解いただくために、まずはこの定理を用いて以下の積分を解いてみましょう。
$$ I = \int_{-\infty}^{\infty} \exp\left\{-\left(x + 3 - \frac{1}{x-4} - \frac{3}{x-1}\right)^2\right\} \mathrm{d}x $$
解法
Glasserのマスター定理より、この関数の$x + 3 - \frac{1}{x-4} - \frac{3}{x-1}$の部分は丸ごと$x$に置き換えられます。
よって、
$$
I = \int_{-\infty}^{\infty} \exp\left(-x^2\right) \mathrm{d}x
$$
これはガウス積分の形なので、
$$
I = \sqrtπ
$$
Glasserのマスター定理によって、あれほど複雑だった被積分関数の一部分を文字一つに置き換えて簡単にすることができました。
では、もう一つ例題を解いてみましょう。
$$
I = \int_{0}^{\infty} \frac{x^2}{x^4 - x^2 + 1} \mathrm{d}x
$$
まず、Glasserのマスター定理を使わずに解くとどのようになるのか見てみましょう。
解法1
\begin{align} I & = \int_{0}^{\infty} \frac{x^2}{x^4 - x^2 + 1} \mathrm{d}x \\ & = \int_{0}^{\infty} \frac{x^2}{(x^2 + 1)^2 - 3x^2} \mathrm{d}x \\ & = \int_{0}^{\infty} \frac{x^2}{(x^2 - \sqrt3 x + 1)(x^2 + \sqrt3 x + 1)} \mathrm{d}x \\ & = \frac{1}{2\sqrt3}\int_{0}^{\infty} \left(\frac{x}{x^2 - \sqrt3 x + 1} - \frac{x}{x^2 + \sqrt3 x + 1}\right) \mathrm{d}x \\ & = \frac{1}{4\sqrt3}\int_{0}^{\infty} \left(\frac{2x - \sqrt3}{x^2 - \sqrt3 x + 1} + \frac{\sqrt3}{x^2 - \sqrt3 x + 1} - \frac{2x + \sqrt3}{x^2 + \sqrt3 x + 1} + \frac{\sqrt3}{x^2 + \sqrt3 x + 1}\right) \mathrm{d}x \\ & = \frac{1}{4\sqrt3}\int_{0}^{\infty} \left\{\left(\ln \left| x^2 - \sqrt3 x + 1 \right|\right)' + \frac{\sqrt3}{\left(x - \frac{\sqrt3}{2}\right)^2 + \frac{1}{4}} - \left(\ln\left| x^2 + \sqrt3 x + 1 \right|\right)' + \frac{\sqrt3}{\left(x + \frac{\sqrt3}{2}\right)^2 + \frac{1}{4}}\right\} \mathrm{d}x \\ & = \frac{1}{4\sqrt3}\left[\ln \left| x^2 - \sqrt3 x + 1 \right| + 2\sqrt3\tan^{-1}\left\{2\left(x - \frac{\sqrt3}{2}\right)\right\} - \ln\left| x^2 + \sqrt3 x + 1 \right| + 2\sqrt3\tan^{-1}\left\{2\left(x + \frac{\sqrt3}{2}\right)\right\} \right]_0^\infty \\ & = \frac{1}{4\sqrt3}\left[\ln \left| \frac{x^2 - \sqrt3 x + 1}{x^2 + \sqrt3 x + 1} \right| + 2\sqrt3\tan^{-1}\left\{2\left(x - \frac{\sqrt3}{2}\right)\right\} + 2\sqrt3\tan^{-1}\left\{2\left(x + \frac{\sqrt3}{2}\right)\right\} \right]_0^\infty \\ & = \frac{1}{4\sqrt3}\left(2\sqrt3π - 0 \right) \\ & = \frac{π}{2} \end{align}
部分分数分解をしなければならなかったり、複雑な関数の極限が必要だったりと、ミスを誘う要素が多いですね。
では、Glasserのマスター定理を使用するとどうなるのかみてみましょう。
解法2
まず、Glasserのマスター定理を使用できる形まで、積分範囲や非積分関数を変形します。
\begin{align}
I & = \int_{0}^{\infty} \frac{x^2}{x^4 - x^2 + 1} \mathrm{d}x \\
& = \frac{1}{2}\int_{-\infty}^{\infty} \frac{x^2}{x^4 - x^2 + 1} \mathrm{d}x \\
& = \frac{1}{2}\int_{-\infty}^{\infty} \frac{\mathrm{d}x}{x^2 - 1 + \frac{1}{x^2}} \\
& = \frac{1}{2}\int_{-\infty}^{\infty} \frac{\mathrm{d}x}{\left(x - \frac{1}{x}\right)^2 + 1} \\
\end{align}
Glasserのマスター定理より、$x - \frac{1}{x}$の部分は丸ごと$x$に置き換えられるので、
\begin{align}
I & = \frac{1}{2}\int_{-\infty}^{\infty} \frac{\mathrm{d}x}{x^2 + 1} \\
& = \frac{1}{2}\left[\tan^{-1}x\right]_{-\infty}^\infty \\
& = \frac{π}{2}
\end{align}
このように、あんなにも複雑な操作を必要としていた定積分でも、Glasserのマスター定理を使用すると瞬殺できてしまうことがあるのです。
Glasserのマスター定理の主張や定理の有用性をご理解いただけましたでしょうか？
Glasserのマスター定理を使用する時は
$$a, b_1, b_2, b_3, \ldots, b_{n-1}, b_n \in \mathbb{R}, \quad a_1, a_2, a_3, \ldots, a_{n-1}, a_n \in \mathbb{R^+}$$
を満たしていないと使用できないことに留意しましょう。特に
$$
a_1, a_2, a_3, \ldots, a_{n-1}, a_n \in \mathbb{R^+}
$$
という条件は重要です。例えば、
$$
\int_{-\infty}^{\infty} \frac{\mathrm{d}x}{\left(x + \frac{1}{x}\right)^2 + 1}
$$
という積分の$x + \frac{1}{x}$の部分は、条件$a_1 \in \mathbb{R^+}$ を満たさないので$x$に置き換えることはできません。
Glasserのマスター定理の魅力をお伝えできたところで、次からはいよいよ証明に移っていきます。
証明のアイデア
いきなり定理の式を見せられて解答の方針が浮かぶ人はまずいないでしょう。
なので、まず$n = 1$のときについて考えて、その時の方針を$n$を一般化した時にも使えないか考えてみます。
$n = 1$のとき
すべての$a, v \in \mathbb{R}, u \in \mathbb{R^+}$について、
$$
\int_{-\infty}^{\infty} f(x + a - \frac{u}{x + v}) \mathrm{d}x = \int_{-\infty}^{\infty} f(x) \mathrm{d}x
$$
まず
$$
I = \int_{-\infty}^{\infty} f(x + a - \frac{u}{x + v}) \mathrm{d}x
$$
とおきます。
最終的な形が$ \int_{-\infty}^{\infty} f(x) \mathrm{d}x $なので、$f$の引数$x + a - \frac{u}{x + v}$を丸ごと一つの文字に帰着させられるような置換をすることを考えます。しかし、
$$
t = x + a - \frac{u}{x + v}
$$
と置換しようとすると、
\begin{align}
t(x + v) & = (x + a)(x + v) - u \\
& = x^2 + (a + v)x + av - u \\
\end{align}
$$
x^2 + (a + v - t)x + (a - t)v - u = 0
$$
\begin{align}
x & = \frac{-a - v + t \pm \sqrt{(a + v - t)^2 - 4\{(a - t)v - 4u\}}}{2} \\
& = \frac{-a - v + t \pm \sqrt{(a + v - t)^2 - 4(a - t)v + 16u}}{2}
\end{align}
となり、$x$が二つ出てきてしまいます。
(ちなみに、このとき$ u < 0 $ だと、ルートの中身が負になるような$t$が存在してしまうので、$u \in \mathbb{R^+}$でなければこの定理は成立しません。)
そこで、積分範囲を二つに分割することでこの重複を解消することを考えます。
ひとまず
$$
α(t) = \frac{-a - v + t - \sqrt{(a + v - t)^2 - 4(a - t)v + 16u}}{2}
$$
$$
β(t) = \frac{-a - v + t + \sqrt{(a + v - t)^2 - 4(a - t)v + 16u}}{2}
$$
と置いてみましょう。つまり、$x$について解いた時の式は、
$$
x = α(t), β(t)
$$
と書き換えられます。
ここで、$y = x + a - \frac{u}{x + v}$のグラフを見てみましょう。
![!FORMULA[43][-1057619717][0](!FORMULA[44][-521360394][0], !FORMULA[45][936347630][0], !FORMULA[46][1393275497][0]のとき)](https://firebasestorage.googleapis.com/v0/b/mathlog-361213.appspot.com/o/uploads%2Fmathdown%2Fl3kYdPNtrCo5cNuTP1Ze.png?alt=media)
$y = x + a - \frac{u}{x + v}$($a = -5$, $v = -2$, $u = \frac{5}{2}$のとき)
関数$y = x + a - \frac{u}{x + v}$が二つの曲線で構成されていることがわかります。
また、$y=α(t)$, $y = β(t)$のグラフを見てみましょう。
![!HTML[0][480702014][0], !HTML[1][-1912271112][0] (!FORMULA[52][-521360394][0], !FORMULA[53][936347630][0], !FORMULA[54][1393275497][0]のとき)](https://firebasestorage.googleapis.com/v0/b/mathlog-361213.appspot.com/o/uploads%2Fmathdown%2Fky4bnzchdGUlkpWlNu97.png?alt=media)
$y=α(t)$, $y = β(t)$ ($a = -5$, $v = -2$, $u = \frac{5}{2}$のとき)
これを見ると、図2のグラフと図1のグラフは直線$y=x$を軸に対称になっています。($α(t)$, $β(t)$は$t = x + a - \frac{u}{x + v}$をxについて解いた結果なので、当然と言えば当然ですね。)
そして、図1のグラフを構成していた曲線のうちの一つが$α$, もう一つが$β$になっているのもわかるでしょう。
図1のグラフにおいて二つの曲線の境界は$x = -v$なので、関数$α$,関数$β$の値域は、それぞれ$(-\infty, -v)$, $(-v, \infty)$です。そしてこれらの関数は連続で単調増加です。
ここまで来れば、積分区間を$(-\infty, -v)$, $(-v, \infty)$の二つに分割するという発想が最も自然であることはお分かりいただけるでしょう。
\begin{align}
I & = \int_{-\infty}^{\infty} f(x + a - \frac{u}{x + v}) \mathrm{d}x \\
& = \int_{-\infty}^{-v} f(x + a - \frac{u}{x + v}) \mathrm{d}x
+ \int_{-v}^{\infty} f(x + a - \frac{u}{x + v}) \mathrm{d}x \\
\end{align}
右辺の1項目の積分については$x = α(t)$、2項目の積分については$x = β(t)$と置換すると、
\begin{align}
I & = \int_{-\infty}^{\infty} f(t) \cdot α'(t) \mathrm{d}t
+ \int_{-\infty}^{\infty} f(t) \cdot β'(t) \mathrm{d}t \\
& = \int_{-\infty}^{\infty} f(t) \{α'(t) + β'(t)\} \mathrm{d}t
\end{align}
ここで、$α(t)$, $β(t)$は二次方程式$ x^2 + (a + v - t)x + (a - t)v - u = 0 $の解なので、解と係数の関係から、
\begin{align}
α(t) + β(t) & = -(a + v - t) \\
& = t - a - v
\end{align}
であり。ここで両辺を$t$について微分すると、
$$
α'(t) + β'(t) = 1
$$
がわかります。よって、
$$
I = \int_{-\infty}^{\infty} f(t) \mathrm{d}t
$$
となり、定理が導けました。
他のnのとき
$n = 1$のときと同じく、$f$の引数を丸ごと一つの文字に帰着させられるような置換をすることを考えます。
ただ、やはり単純に
$t = x + a - \frac{a_1}{x + b_1} - \frac{a_2}{x + b_2} - \frac{a_3}{x + b_3} - \dotsb - \frac{a_{n-1}}{x + b_{n-1}} - \frac{a_n}{x + b_n}$
と置換するのは厳しい。
$x$について整理すると$x$についての$n + 1$次方程式になり、それを解くと$x = (tの関数)$の形が$n + 1$通り出てきてしまいます。
そこで、$n = 1$の時と同じように、$x = (tの関数)$の形が1通りに定まるようにうまく積分範囲を分割することを考えます。
ただどのように積分範囲を分割すれば良いのでしょうか？
ここで、例として$ y = x + 10 - \frac{100}{x - 40} - \frac{100}{x + 30} - \frac{100}{x + 10} - \frac{100}{x - 20}$という関数のグラフを見てみましょう。
![!FORMULA[90][1895137698][0]](https://firebasestorage.googleapis.com/v0/b/mathlog-361213.appspot.com/o/uploads%2Fmathdown%2FfhO5X9mjPMHIlVDuA5RS.png?alt=media)
$ y = x + 10 - \frac{100}{x - 40} - \frac{100}{x + 30} - \frac{100}{x + 10} - \frac{100}{x - 20}$
この関数は、5つの曲線で構成されているのがお分かりいただけますでしょうか？
$n = 1$のときは、2つの曲線で構成されていました。なので積分範囲は二つに分割するだけで十分でした。しかし、5つの曲線で構成されているとなれば、5つに積分範囲を分割する必要があるでしょう。
また、$n = 1$のときは曲線と曲線の境界で分割していました。なので、この場合も曲線と曲線の境界で分割するのが自然そうです。
別の例も見てみましょう。以下は$y = x + 3 - \frac{80}{x - 50} - \frac{80}{x + 45} - \frac{50}{x - 11}$のグラフです。
![!FORMULA[94][1361036709][0]](https://firebasestorage.googleapis.com/v0/b/mathlog-361213.appspot.com/o/uploads%2Fmathdown%2FSBgGvKQTQDYSaxz4dCxE.png?alt=media)
$y = x + 3 - \frac{80}{x - 50} - \frac{80}{x + 45} - \frac{50}{x - 11}$
これもやはり4つの曲線で構成されていますので、4つに積分範囲を分割する必要がありそうです。
ここまでくれば、曲線と曲線の境目で積分区間を$n+1$区間に分割すれば、分割したそれぞれの積分において$x$を$t$を用いて一通りで表現できるのではないかという推測が立つでしょう。
本記事では、このような発想をもってGlasserのマスター定理を証明していきます。
証明
$b_1 > b_2 > b_3 > \dotsb > b_{n-1} > b_n$とする。
また、
$$g(x) = x + a - \frac{a_1}{x + b_1} - \frac{a_2}{x + b_2} - \frac{a_3}{x + b_3} - \dotsb - \frac{a_{n-1}}{x + b_{n-1}} - \frac{a_n}{x + b_n}$$
$$I = \int_{-\infty}^{\infty} f(x + a - \frac{a_1}{x + b_1} - \frac{a_2}{x + b_2} - \frac{a_3}{x + b_3} - \dotsb - \frac{a_{n-1}}{x + b_{n-1}} - \frac{a_n}{x + b_n}) \mathrm{d}x$$
とおく。すると、
$$
I = \int_{-\infty}^{\infty} f(g(x)) \mathrm{d}x
$$
となる。
任意の実数$t$で
\begin{align}
t & = g(α(t)) \tag{1}
\end{align}
となるような実関数$α$を考える。(以下、関数の引数を省略する場合がある。)
$$
t = α(t) + a - \frac{a_1}{α(t) + b_1} - \frac{a_2}{α(t) + b_2} - \frac{a_3}{α(t) + b_3} - \dotsb - \frac{a_{n-1}}{α(t) + b_{n-1}} - \frac{a_n}{α(t) + b_n}
$$
上の式を分母を払うと、
\begin{align}
t(α & + b_1)(α + b_2)(α + b_3) \dotsb (α + b_{n-1})(α + b_n) \\
= & (α + a)(α + b_1)(α + b_2)(α + b_3) \dotsb (α + b_{n-1})(α + b_n) \\
& - a_1(α + b_2)(α + b_3)(α + b_4) \dotsb (α + b_{n-1})(α + b_n) \\
& - a_2(α + b_1)(α + b_3)(α + b_4) \dotsb (α + b_{n-1})(α + b_n) \\
& - a_3(α + b_1)(α + b_2)(α + b_4) \dotsb (α + b_{n-1})(α + b_n) \\
& \quad \vdots \\
& - a_{n-1}(α + b_1)(α + b_2)(α + b_3) \dotsb (α + b_{n-2})(α + b_n) \\
& - a_n(α + b_1)(α + b_2)(α + b_3) \dotsb (α + b_{n-1})(α + b_n) \\
\end{align}
\begin{align}
\therefore \, & (α + a - t)(α + b_1)(α + b_2)(α + b_3) \cdots (α + b_{n-1})(α + b_n) \\
& - a_1(α + b_2)(α + b_3)(α + b_4) \dotsb (α + b_{n-1})(α + b_n) \\
& - a_2(α + b_1)(α + b_3)(α + b_4) \dotsb (α + b_{n-1})(α + b_n) \\
& - a_3(α + b_1)(α + b_2)(α + b_4) \dotsb (α + b_{n-1})(α + b_n) \\
& \quad \vdots \\
& - a_{n-1}(α + b_1)(α + b_2)(α + b_3) \dotsb (α + b_{n-2})(α + b_n) \\
& - a_n(α + b_1)(α + b_2)(α + b_3) \dotsb (α + b_{n-1})(α + b_n) \\
& = 0 \tag{2}
\end{align}
となり、$α(t)$に関する$n + 1$次方程式になる。よって、式$(1)$を満たす関数$α$は、高々$n + 1$個しか存在しない。
また、
\begin{align} \lim_{x \to -\infty} g(x) = -\infty \, & \lim_{x \to -b_1 - 0} g(x) = \infty, \\ \lim_{x \to -b_1 + 0} g(x) = -\infty \, & \lim_{x \to -b_2 - 0} g(x) = \infty, \\ \lim_{x \to -b_2 + 0} g(x) = -\infty \, & \lim_{x \to -b_3 - 0} g(x) = \infty, \\ \vdots \qquad \qquad \qquad \qquad \\ \lim_{x \to -b_{n-1} + 0} g(x) = -\infty \, & \lim_{x \to -b_n - 0} g(x) = \infty, \\ \lim_{x \to -b_n + 0} g(x) = -\infty \, & \lim_{x \to \infty} g(x) = \infty \qquad \qquad \qquad(★) \\ \end{align}
でかつ関数$ g $は$ (-\infty, -b_1), (-b_1, -b_2), (-b_2, -b_3), \ldots, (-b_{n-1}, -b_n), (-b_n, \infty) $の$n + 1$区間において単調増加で連続である。
よって、関数$g$の定義域を$ (-\infty, -b_1), (-b_1, -b_2), (-b_2, -b_3), \ldots, (-b_{n-1}, -b_n), (-b_n, \infty) $に限定した関数をそれぞれ$g_1, g_2, g_3, \ldots, g_n, g_{n+1}$とおくと、$1 ≤ k ≤ n + 1$を満たすすべての$k \in \mathbb{Z}$について、$g_k$は全単射で連続でかつ単調増加である。
ここで、$g_k$の逆関数を$g_k^{-1}$とおく。
式$(1)$において、$α(t) \in (-\infty, -b_1)$のとき、
$$
t = g_1(α(t))
$$
\begin{align}
g_1^{-1}(t) & = g_1^{-1}(g_1(α(t))) \\
& = α(t) \\
\end{align}
$$
\therefore α = g_1^{-1}
$$
同様に、$α(t) \in (-b_1, -b_2)$のとき、
$$
α = g_2^{-1}
$$
$α(t) \in (-b_2, -b_3)$のとき、
$$
α = g_3^{-1}
$$
というように、$ (-\infty, -b_1), (-b_1, -b_2), (-b_2, -b_3), \ldots, (-b_{n-1}, -b_n), (-b_n, \infty) $のすべての区間について同じことが言えるので、
$$α = g_1^{-1}, g_2^{-1}, g_3^{-1}, \ldots, g_n^{-1}, g_{n+1}^{-1}$$
は式$(1)$を満たす。
そして式$(1)$を満たす$α$は高々$n + 1$個しか存在しないので、これら以外に式$(1)$を満たす$α$は存在しない。そしてこれらは$(2)$の方程式の解のうちすべてである。
よって式$(2)$より、解と係数の関係から、
\begin{align}
g_1^{-1}(t) + g_2^{-1}(t) + g_3^{-1}(t) + \dotsb + g_n^{-1}(t) + g_{n + 1}^{-1}(t) & = -\left\{(a - t) + b_1 + b_2 + b_3 + \dotsb + b_{n-1} + b_n \right\} \\
& = t - a - b_1 - b_2 - b_3 - \dotsb - b_{n-1} - b_n \\
\end{align}
両辺を微分して、
$$
{g_1^{-1}}'(t) + {g_2^{-1}}'(t) + {g_3^{-1}}'(t) + \dotsb + {g_n^{-1}}'(t) + {g_{n + 1}^{-1}}'(t) = 1 \tag{3}
$$
ここで、$I$の積分範囲を下のように分割する。
\begin{align}
I = & \int_{-\infty}^{-b_1} f(g(x)) \mathrm{d}x \\
& + \int_{-b_1}^{-b_2} f(g(x)) \mathrm{d}x \\
& + \int_{-b_2}^{-b_3} f(g(x)) \mathrm{d}x \\
& \qquad \vdots \\
& + \int_{-b_{n-1}}^{-b_n} f(g(x)) \mathrm{d}x \\
& + \int_{-b_n}^{\infty} f(g(x)) \mathrm{d}x \\
\end{align}
そして、右辺の項のうち$k$項目を$I_k$とおく。つまり、
$$
I = I_1 + I_2 + I_3 + \dotsb + I_n + I_{n + 1}
$$
となる。
すべての$I_k$について、$x = g_k^{-1}(t)$と置換すると、$(★)$より、すべての$I_k$で積分範囲が$-\infty$から$\infty$になるので、
\begin{align}
I_k & = \int_{-\infty}^{\infty} f(g(g_k^{-1}(t))) \cdot {g_k^{-1}}'(t) \mathrm{d}t \\
& = \int_{-\infty}^{\infty} f(t) {g_k^{-1}}'(t) \mathrm{d}t
\end{align}
\begin{align}
\therefore I = & \int_{-\infty}^{\infty} f(t) {g_1^{-1}}'(t) \mathrm{d}t \\
& + \int_{-\infty}^{\infty} f(t) {g_2^{-1}}'(t) \mathrm{d}t \\
& + \int_{-\infty}^{\infty} f(t) {g_3^{-1}}'(t) \mathrm{d}t \\
& \qquad \vdots \\
& + \int_{-\infty}^{\infty} f(t) {g_n^{-1}}'(t) \mathrm{d}t \\
& + \int_{-\infty}^{\infty} f(t) {g_{n + 1}^{-1}}'(t) \mathrm{d}t \\
= & \int_{-\infty}^{\infty} f(t) \left\{{g_1^{-1}}'(t) + {g_2^{-1}}'(t) + {g_3^{-1}}'(t) + \dotsb + {g_n^{-1}}'(t) + {g_{n + 1}^{-1}}'(t)\right\} \mathrm{d}t \\
= & \int_{-\infty}^{\infty} f(t) \mathrm{d}t \qquad \left(\, \because (3) \, \right) \\
= & \int_{-\infty}^{\infty} f(x) \mathrm{d}x \\
\end{align}
よって定理は示された。
おわりに
個人的には日本語資料がないことが信じられないほど有用な定理だと思っています。なんで日本語資料ないんや...
この記事によってあなたの積分ライフがより豊かになることを願っています。
証明めっちゃ大変だった...
あと、$\tan$、$\cot$の部分分数展開と組み合わせると三角関数の積分にもこの定理が活用できることがあるらしいよ...(ﾎﾞｿｯ