# 留数定理と定積分への応用：孤立特異点の分類から実積分の計算まで

> 孤立特異点を除去可能・極・真性の三つに分類し、留数を定義して留数定理を証明する。さらに有理関数の広義積分、三角関数の周期積分、ジョルダンの補題を使うフーリエ型積分を、輪郭の選び方から計算例まで具体的に扱う。
> https://rikai.mugen-giken.com/mathematics/complex-analysis/residue-theorem

## 0. この記事の要点

- 正則関数を孤立特異点 $a$ のまわりでローラン展開したとき、$(z-a)^{-1}$ の係数 $c_{-1}$ だけが $a$ を囲む閉曲線に沿う積分に生き残ります。この $c_{-1}$ を留数と呼びます。
- 孤立特異点は主要部の形で三つに分かれます。除去可能特異点は「近くで有界」、$m$ 位の極は「$(z-a)^m f(z)$ が $0$ でない有限の極限を持つ」、真性特異点は「どんな小さい近傍の像も複素平面で稠密」と、いずれも極限の言葉に翻訳できます。
- 留数定理は「単純閉曲線に沿う積分 $= 2\pi i \times$（内部の留数の和）」です。証明は、多重連結領域に対するコーシーの積分定理で積分路を各特異点まわりの小円に分解し、小円上でローラン級数を項別積分するだけです。
- 実積分への応用は三段構えです。(1) 被積分関数を複素平面に延長して積分路を閉じる、(2) 付け足した部分の寄与が $0$ に消えることを不等式で示す、(3) 内部の留数を足す。
- 型ごとに使う輪郭が決まっています。有理関数の $(-\infty,\infty)$ 積分は上半平面の半円、周期 $2\pi$ の三角関数の積分は単位円、$e^{i\lambda x}$ を含む積分はジョルダンの補題つきの半円、実軸上に単純極があるときは極を小円でくぼませて主値を取ります。

## 1. 動機：実積分を複素平面に持ち上げる

微分積分学で $\int_{-\infty}^{\infty}\frac{dx}{x^2+1}$ を計算するのは簡単です。原始関数 $\arctan x$ が分かっているので、$\lim_{R\to\infty}[\arctan x]_{-R}^{R} = \pi$ で終わります。ところが分母を $x^4+1$ に変えただけで景色が一変します。実数の範囲では $x^4+1 = (x^2+\sqrt{2}x+1)(x^2-\sqrt{2}x+1)$ と因数分解し、部分分数に分け、四つの項をそれぞれ $\arctan$ と $\log$ で積分し、最後に極限を取る、という長い計算になります。$\int_0^{\infty}\frac{\sin x}{x}\,dx$ に至っては、初等関数で書ける原始関数がそもそも存在しないので、パラメータ微分などの技巧なしには実数の世界で手が出ません。

この章で扱う留数定理は、この種の積分を次の手続きに置き換えます。被積分関数を複素平面に延長し、積分路を閉曲線になるように閉じ、閉曲線の内部にある特異点で「留数」という一つの数を計算して足す。それだけです。$\int_{-\infty}^{\infty}\frac{dx}{x^4+1}$ は、後で見るように二つの留数の足し算だけで $\pi/\sqrt{2}$ と分かります。

なぜそんなことができるのでしょうか。出発点は [コーシーの積分定理](/mathematics/complex-analysis/cauchy-integral-theorem)（<Ref to="mathematics/complex-analysis/cauchy-integral-theorem#thm-cauchy-green" />）です。単連結領域で正則な関数は閉曲線に沿う積分が $0$ になります。つまり、閉曲線に沿う積分が $0$ でないとしたら、その値は**内部にある特異点だけが原因**です。では特異点はどれだけ寄与するのでしょうか。中心 $a$ の小さい円周 $|z-a| = r$ 上で $(z-a)^n$ を積分してみます。$n \ne -1$ のとき $(z-a)^n$ は $\mathbb{C}\setminus\{a\}$ 全体で一価な原始関数 $\frac{(z-a)^{n+1}}{n+1}$ を持つので、閉曲線に沿う積分は $0$ です。ところが $n = -1$ のときだけ、原始関数の候補は $\log(z-a)$ という多価関数しかなく、円周を一周すると偏角が $2\pi$ 増えて、積分は $2\pi i$ を返します。

したがって、ローラン展開 $f(z) = \sum_{n=-\infty}^{\infty}c_n(z-a)^n$ を円周上で項別に積分すると、無限個の項のうち $c_{-1}(z-a)^{-1}$ ただ一つだけが生き残ります。この「取り残される数」がフランス語で résidu（残りもの）、日本語で留数です。コーシーは 1826 年の『Exercices de mathématiques』でこの語を導入し、実積分の計算に系統的に用いました。この章は、その計算術を最初から最後まで具体的に組み立てます。

<div data-gated data-pagefind-ignore>

## 2. 準備：記号とすでに証明した道具

**記号**：領域とは連結開集合のことです。$D(a,r) = \{z \in \mathbb{C} : |z-a| < r\}$、$D'(a,r) = \{z : 0 < |z-a| < r\}$（中心を抜いた円板、除外近傍）と書きます。曲線はすべて区分的に $C^1$ 級とし、$\int_\gamma f(z)\,dz = \int_\alpha^\beta f(\gamma(t))\gamma'(t)\,dt$ で定義します。円周 $|z-a| = r$ には断らない限り正の向き（反時計回り）を与え、閉曲線に沿う積分は $\oint$ と書きます。

**長さと最大値による評価（ML 不等式）**：曲線 $\gamma$ の長さを $L(\gamma)$ とすると

$$
\left|\int_\gamma f(z)\,dz\right| \le \int_\gamma |f(z)|\,|dz| \le L(\gamma)\cdot\max_{z \in \gamma}|f(z)|
$$

が成り立ちます（<Ref to="mathematics/complex-analysis/cauchy-integral-theorem#lem-ml" />）。この章の評価はほとんどすべてこの不等式に帰着します。

**ローラン展開**：[正則関数の強力な性質](/mathematics/complex-analysis/properties-of-holomorphic-functions) の <Ref to="mathematics/complex-analysis/properties-of-holomorphic-functions#thm-laurent" /> で示したとおり、$f$ が円環 $r < |z-a| < R$ で正則ならば、そこで

$$
f(z) = \sum_{n=-\infty}^{\infty}c_n (z-a)^n,
\qquad
c_n = \frac{1}{2\pi i}\oint_{|\zeta-a|=\rho}\frac{f(\zeta)}{(\zeta-a)^{n+1}}\,d\zeta
\quad (r < \rho < R)
$$

と一意的に展開でき、この級数は円環の任意のコンパクト部分集合上で一様収束します。係数 $c_n$ は $\rho$ の取り方によりません。

**多重連結領域に対するコーシーの積分定理**：$C$ を正の向きの単純閉曲線、$\gamma_1,\dots,\gamma_n$ を $C$ の内部にあり互いに他の外部にある正の向きの単純閉曲線とします。$f$ が「$C$ の内側かつすべての $\gamma_j$ の外側」という閉領域を含むある開集合で正則ならば

$$
\oint_C f(z)\,dz = \sum_{j=1}^{n}\oint_{\gamma_j} f(z)\,dz
$$

が成り立ちます。証明は前章の <Ref to="mathematics/complex-analysis/cauchy-integral-theorem#thm-deformation" /> と同じ要領で、$C$ と各 $\gamma_j$ を細い「橋」（切断線）でつないで単連結な部分に分け、各橋を往復する積分が符号違いで相殺することを使います。留数定理はこの定理の言い換えにすぎません。

## 3. 孤立特異点の分類

<Definition id="def-isolated-singularity" title="孤立特異点">
$a \in \mathbb{C}$ とする。ある $R > 0$ が存在して $f$ が除外近傍 $D'(a,R) = \{z : 0 < |z-a| < R\}$ で正則であるとき、$a$ を $f$ の**孤立特異点**と呼ぶ。$f$ は $a$ 自身では定義されていなくてもよいし、定義されていても正則でなくてよい。
</Definition>

孤立特異点のまわりでは $f$ は円環 $0 < |z-a| < R$ で正則ですから、§2 のローラン展開が使えます。展開の負べきの部分が特異点の「悪さ」を測る量です。

<Definition id="def-classification" title="主要部と特異点の三分類">
$a$ を $f$ の孤立特異点とし、$D'(a,R)$ でのローラン展開を $f(z) = \sum_{n=-\infty}^{\infty}c_n(z-a)^n$ とする。負べきの部分

$$
P_a(z) = \sum_{n \le -1} c_n (z-a)^n
$$

を $a$ における $f$ の**主要部**と呼ぶ。主要部の形により、$a$ を次のように分類する。

- $c_n = 0$（すべての $n \le -1$）、すなわち主要部が $0$ のとき、$a$ を**除去可能特異点**という。
- ある正の整数 $m$ について $c_{-m} \ne 0$ かつ $c_n = 0$（すべての $n < -m$）のとき、$a$ を **$m$ 位の極**という。$m = 1$ のときは**単純極**ともいう。
- $c_n \ne 0$ となる $n \le -1$ が無限個あるとき、$a$ を**真性特異点**という。

ローラン展開は一意なので、この分類は展開の取り方によらず定まる。
</Definition>

<Figure caption="孤立特異点の三分類。主要部の項数を数えるだけで型が決まり、それぞれに極限による判定条件が対応します。">
<Mermaid code={`flowchart TD
  A["孤立特異点 a のまわりでローラン展開する"] --> B["主要部（負べきの項の全体）を見る"]
  B --> C["負べきの項が無い"]
  B --> D["負べきが有限個で、最低次が -m"]
  B --> E["負べきが無限個ある"]
  C --> C2["除去可能特異点：a まで正則に延長できる。近くでの有界性で判定できる"]
  D --> D2["m 位の極：z が a に近づくと絶対値が無限大に発散する"]
  E --> E2["真性特異点：どんな近傍の像も複素平面で稠密になる"]`} />
</Figure>

分類の定義はローラン係数を全部調べよと言っているように見えますが、実際にはもっと安上がりな判定法があります。まず除去可能特異点です。

<Theorem id="thm-riemann-removable" title="リーマンの除去可能特異点定理">
$a$ を $f$ の孤立特異点とし、$f$ は $D'(a,R)$ で正則とする。次の 4 条件は同値である。

- **(i)** $a$ は $f$ の除去可能特異点である。
- **(ii)** $f(a)$ の値を適当に定めると、$f$ は $D(a,R)$ 全体で正則になる。
- **(iii)** ある $0 < r \le R$ と $M > 0$ があって、$0 < |z-a| < r$ で $|f(z)| \le M$ が成り立つ。
- **(iv)** $\lim_{z \to a}(z-a)f(z) = 0$。
</Theorem>

<Proof of="thm-riemann-removable">
**(i) $\Rightarrow$ (ii)**：主要部が $0$ なので $D'(a,R)$ で $f(z) = \sum_{n \ge 0}c_n(z-a)^n$ です。右辺のべき級数は $D'(a,R)$ で収束するので収束半径は $R$ 以上であり、$D(a,R)$ 全体で正則関数を定めます。$f(a) := c_0$ と定めれば、$f$ はこのべき級数と $D(a,R)$ 全体で一致し、正則になります。

**(ii) $\Rightarrow$ (iii)**：正則関数は連続なので、$r = R/2$ とすればコンパクト集合 $\overline{D(a,r)}$ 上で $|f|$ は最大値 $M$ を取ります。

**(iii) $\Rightarrow$ (iv)**：$|(z-a)f(z)| \le M|z-a| \to 0$（$z \to a$）です。

**(iv) $\Rightarrow$ (i)**：$0 < \rho < R$ に対し $\eta(\rho) := \max_{|z-a| = \rho}|(z-a)f(z)|$ とおくと、仮定より $\eta(\rho) \to 0$（$\rho \to 0$）です。§2 のローラン係数の式と ML 不等式から、$n \le -1$ に対して

$$
|c_n| \le \frac{1}{2\pi}\cdot 2\pi\rho \cdot \frac{\max_{|z-a|=\rho}|f(z)|}{\rho^{\,n+1}}
= \rho^{-n}\max_{|z-a|=\rho}|f(z)|
= \rho^{-n-1}\,\eta(\rho)
$$

が成り立ちます。ここで $n \le -1$ より $-n-1 \ge 0$ なので、$\rho < 1$ のとき $\rho^{-n-1} \le 1$ です。よって $|c_n| \le \eta(\rho) \to 0$ となり、$c_n$ は $\rho$ によらない定数ですから $c_n = 0$ です。これがすべての $n \le -1$ で成り立つので主要部は $0$ です。
</Proof>

条件 (iii) が重要です。「近くで有界」という、ローラン係数を一つも計算しなくてよい条件だけで、特異点が消せると分かります。次は極です。

<Proposition id="prop-pole-characterization" title="極の特徴づけ">
$a$ を $f$ の孤立特異点、$m$ を正の整数とする。

**(1)** 次の 3 条件は同値である。

- **(i)** $a$ は $f$ の $m$ 位の極である。
- **(ii)** 極限 $\lim_{z \to a}(z-a)^m f(z)$ が存在し、その値は $0$ でない。
- **(iii)** ある $r > 0$ と $D(a,r)$ で正則な関数 $g$ で $g(a) \ne 0$ を満たすものがあり、$0 < |z-a| < r$ で $f(z) = \dfrac{g(z)}{(z-a)^m}$ と書ける。

**(2)** $a$ が（何らかの位数の）極であるための必要十分条件は $\lim_{z \to a}|f(z)| = \infty$ である。
</Proposition>

<Proof of="prop-pole-characterization">
**(1) (i) $\Rightarrow$ (iii)**：$m$ 位の極なら $D'(a,R)$ で $f(z) = \sum_{n \ge -m}c_n(z-a)^n$ かつ $c_{-m} \ne 0$ です。両辺に $(z-a)^m$ を掛けると

$$
(z-a)^m f(z) = \sum_{n \ge -m} c_n (z-a)^{n+m} = \sum_{k \ge 0} c_{k-m}(z-a)^k
$$

となり、右辺は $D'(a,R)$ で収束するべき級数ですから収束半径は $R$ 以上で、$D(a,R)$ で正則関数 $g$ を定めます。$g(a) = c_{-m} \ne 0$ であり、$0 < |z-a| < R$ で $f(z) = g(z)/(z-a)^m$ です。

**(iii) $\Rightarrow$ (ii)**：$(z-a)^m f(z) = g(z) \to g(a) \ne 0$（$z \to a$、$g$ の連続性）です。

**(ii) $\Rightarrow$ (i)**：$h(z) := (z-a)^m f(z)$ とおくと、$h$ は $D'(a,R)$ で正則で $z \to a$ のとき有限の極限 $L \ne 0$ を持つので、$a$ の近くで有界です。<Ref to="thm-riemann-removable" /> の (iii) $\Rightarrow$ (i),(ii) より $a$ は $h$ の除去可能特異点で、$h(a) = L$ と定めれば $h$ は $D(a,R)$ で正則になります。$h$ をテイラー展開して $h(z) = \sum_{k \ge 0}b_k(z-a)^k$（$b_0 = L \ne 0$）とすると、$0 < |z-a| < R$ で

$$
f(z) = \frac{h(z)}{(z-a)^m} = \sum_{k \ge 0} b_k (z-a)^{k-m}
$$

となります。これは $f$ のローラン展開そのもの（一意性より）で、最低次の項は $b_0(z-a)^{-m}$、$b_0 \ne 0$ ですから $a$ は $m$ 位の極です。

**(2) 十分性**：$a$ が $m$ 位の極なら (iii) より $f = g/(z-a)^m$、$g(a) \ne 0$ なので、$|f(z)| = |g(z)|\,|z-a|^{-m} \to |g(a)| \cdot \infty = \infty$ です。

**必要性**：$|f(z)| \to \infty$ とします。ある $r > 0$ があって $0 < |z-a| < r$ で $|f(z)| \ge 1$、特に $f(z) \ne 0$ です。そこで $u := 1/f$ とおくと $u$ は $D'(a,r)$ で正則、$|u| \le 1$ と有界なので、<Ref to="thm-riemann-removable" /> より $a$ は $u$ の除去可能特異点で、$u$ は $D(a,r)$ で正則に延長されます。$|f| \to \infty$ より $u(a) = 0$ です。$u$ は $D'(a,r)$ で $0$ にならないので $u \not\equiv 0$ であり、<Ref to="mathematics/complex-analysis/properties-of-holomorphic-functions#lem-zero-isolated" /> で見たとおり、恒等的に $0$ でない正則関数の零点は有限の位数を持ちます。その位数を $m$ とすれば $u(z) = (z-a)^m v(z)$、$v$ は $D(a,r')$ で正則、$v(a) \ne 0$ と書けます。よって $0 < |z-a| < r'$ で $f(z) = \dfrac{1/v(z)}{(z-a)^m}$ となり、$1/v$ は $a$ の近くで正則（$v(a)\ne 0$ と連続性より近傍で $v \ne 0$）かつ $a$ で $0$ でないので、(1)(iii) から $a$ は $m$ 位の極です。
</Proof>

残るのが真性特異点です。<Ref to="thm-riemann-removable" /> と <Ref to="prop-pole-characterization" /> により、$a$ が真性特異点であることは「$z \to a$ で $f(z)$ が有限の極限も持たず、$|f(z)| \to \infty$ でもない」ことと同値です。極限が存在しないだけでなく、値の散らばり方は極端です。

<Theorem id="thm-casorati-weierstrass" title="カゾラーティ・ワイエルシュトラスの定理">
$a$ を $f$ の真性特異点とし、$f$ は $D'(a,R)$ で正則とする。このとき任意の $0 < r \le R$ に対して像 $f(D'(a,r))$ は $\mathbb{C}$ で稠密である。すなわち、任意の $w \in \mathbb{C}$ と任意の $\varepsilon > 0$ に対し、$0 < |z-a| < r$ かつ $|f(z) - w| < \varepsilon$ を満たす $z$ が存在する。
</Theorem>

<Proof of="thm-casorati-weierstrass">
稠密でないと仮定します。すると、ある $w \in \mathbb{C}$ と $\delta > 0$ があって、$0 < |z-a| < r$ のすべての $z$ で $|f(z) - w| \ge \delta$ となります。そこで

$$
g(z) := \frac{1}{f(z) - w} \qquad (0 < |z-a| < r)
$$

とおくと、分母は $0$ にならないので $g$ は $D'(a,r)$ で正則で、$|g| \le 1/\delta$ と有界です。<Ref to="thm-riemann-removable" /> より $a$ は $g$ の除去可能特異点で、$g$ は $D(a,r)$ で正則に延長されます。

$g(a) \ne 0$ の場合：$f(z) = w + 1/g(z)$ は $a$ の近くで有界なので、再び <Ref to="thm-riemann-removable" /> により $a$ は $f$ の除去可能特異点です。

$g(a) = 0$ の場合：$g$ は $D'(a,r)$ で $0$ にならないので $g \not\equiv 0$ であり、$a$ は $g$ の有限位数 $m$ の零点です。<Ref to="prop-pole-characterization" /> の (1)(iii) の形に書けば $f(z) - w$ は $a$ で $m$ 位の極を持ち、定数 $w$ を足しても主要部は変わらないので $f$ も $a$ で $m$ 位の極を持ちます。

どちらの場合も $a$ が真性特異点であることに反します。
</Proof>

<Example id="ex-classify" title="特異点の分類の実例">
**(a)** $f(z) = \dfrac{\sin z}{z}$ の $z = 0$。$\sin z = z - \frac{z^3}{3!} + \frac{z^5}{5!} - \cdots$ を $z$ で割ると $f(z) = 1 - \frac{z^2}{6} + \frac{z^4}{120} - \cdots$ となり、負べきの項がありません。よって $0$ は除去可能特異点で、$f(0) := 1$ と定めれば整関数になります。

**(b)** $f(z) = \dfrac{1}{z^2(z-1)}$。$z = 0$ では $z^2 f(z) = \frac{1}{z-1} \to -1 \ne 0$ なので <Ref to="prop-pole-characterization" /> より $2$ 位の極、$z = 1$ では $(z-1)f(z) = \frac{1}{z^2} \to 1 \ne 0$ なので単純極です。

**(c)** $f(z) = \dfrac{1}{\sin z}$ の $z = k\pi$（$k \in \mathbb{Z}$）。$\sin(k\pi) = 0$ かつ $\frac{d}{dz}\sin z\big|_{z=k\pi} = \cos(k\pi) = (-1)^k \ne 0$ なので $k\pi$ は $\sin z$ の $1$ 位の零点、したがって

$$
(z-k\pi)f(z) = \frac{z-k\pi}{\sin z - \sin(k\pi)} \longrightarrow \frac{1}{\cos(k\pi)} = (-1)^k \ne 0
$$

（微分係数の定義そのもの）となり、$k\pi$ は単純極です。

**(d)** $f(z) = e^{1/z}$ の $z = 0$。$e^{w} = \sum_{n\ge0}w^n/n!$ に $w = 1/z$ を代入して

$$
e^{1/z} = \sum_{n \ge 0}\frac{1}{n!}z^{-n} = 1 + \frac{1}{z} + \frac{1}{2!z^2} + \cdots
\qquad (z \ne 0)
$$

が $0 < |z|$ で収束します。負べきの項が無限個あるので $0$ は真性特異点です。実際、値の散らばり方を直接確かめられます。$w \ne 0$ を任意に取り $w = |w|e^{i\theta}$ と書くと、$z_k = \dfrac{1}{\log|w| + i(\theta + 2\pi k)}$（$k \in \mathbb{N}$）は $k \to \infty$ で $0$ に収束し、$e^{1/z_k} = w$ です。つまり $0$ のどんな近傍でも、$f$ は $0$ 以外のすべての複素数を無限回取ります。<Ref to="thm-casorati-weierstrass" /> の「稠密」よりずっと強い状況です。
</Example>

<Remark id="rem-picard">
<Ref to="ex-classify" /> (d) の現象は一般に成り立ちます。ピカールの大定理によれば、真性特異点のどんな除外近傍でも、$f$ は高々 1 個の例外値を除くすべての複素数値を無限回取ります。証明はこの章の範囲を超えるので、Ahlfors『Complex Analysis』の第 8 章（Picard の定理を扱う章）を参照してください。
</Remark>

<Remark id="rem-non-isolated">
すべての特異点が孤立しているわけではありません。$f(z) = \dfrac{1}{\sin(1/z)}$ は $z = 1/(k\pi)$（$k$ は $0$ でない整数）で極を持ちますが、これらは $k \to \infty$ で $0$ に集積します。したがって $0$ のどんな除外近傍にも $f$ が正則でない点が含まれ、$0$ は <Ref to="def-isolated-singularity" /> の意味の孤立特異点ではありません。ローラン展開も留数も定義できないので、この章の道具は $0$ には使えません。
</Remark>

## 4. 留数とその計算法

<Definition id="def-residue" title="留数">
$a$ を $f$ の孤立特異点とし、$D'(a,R)$ でのローラン展開を $f(z) = \sum_{n=-\infty}^{\infty}c_n(z-a)^n$ とする。係数 $c_{-1}$ を $f$ の $a$ における**留数**と呼び、

$$
\operatorname{Res}(f, a) := c_{-1}
$$

と書く。ローラン係数の積分表示（§2）で $n = -1$ とすれば、$0 < \rho < R$ に対して

$$
\operatorname{Res}(f,a) = \frac{1}{2\pi i}\oint_{|z-a| = \rho} f(z)\,dz
$$

でもある。
</Definition>

定義そのものは「ローラン展開せよ」と言っていますが、極の場合には展開せずに留数を取り出す公式があります。実用上はこの公式だけで九割方の計算が片付きます。

<Proposition id="prop-residue-formulas" title="極における留数の計算法">
$a$ を $f$ の孤立特異点とする。

**(1)** $a$ が $f$ の高々 $m$ 位の極である（すなわちローラン展開で $n < -m$ の係数がすべて $0$ である）とき、$\varphi(z) := (z-a)^m f(z)$ は $a$ で正則に延長され、

$$
\operatorname{Res}(f,a) = \frac{1}{(m-1)!}\lim_{z \to a}\frac{d^{\,m-1}}{dz^{\,m-1}}\Bigl[(z-a)^m f(z)\Bigr].
$$

**(2)** 特に $a$ が単純極（$m=1$）のとき $\operatorname{Res}(f,a) = \lim_{z \to a}(z-a)f(z)$。

**(3)** $f = g/h$ とし、$g,h$ は $a$ の近傍で正則、$g(a) \ne 0$、$h(a) = 0$、$h'(a) \ne 0$ とする。このとき $a$ は $f$ の単純極で

$$
\operatorname{Res}(f,a) = \frac{g(a)}{h'(a)}.
$$
</Proposition>

<Proof of="prop-residue-formulas">
**(1)**：仮定より $D'(a,R)$ で $f(z) = \sum_{n \ge -m}c_n(z-a)^n$ なので、

$$
\varphi(z) = (z-a)^m f(z) = \sum_{k \ge 0}c_{k-m}(z-a)^k
$$

であり、右辺は $D'(a,R)$ で収束するべき級数、したがって収束半径は $R$ 以上で $D(a,R)$ で正則関数を定めます（<Ref to="prop-pole-characterization" /> の証明と同じ議論です）。この級数の $(z-a)^{m-1}$ の係数は $k = m-1$ の項、すなわち $c_{-1}$ です。一方、正則関数のテイラー係数の公式（<Ref to="mathematics/complex-analysis/properties-of-holomorphic-functions#thm-taylor" />）から $(z-a)^{m-1}$ の係数は $\varphi^{(m-1)}(a)/(m-1)!$ です。両者を等しいとおけば結論を得ます。$\varphi$ は $a$ で正則なので $\varphi^{(m-1)}$ は $a$ で連続で、$\varphi^{(m-1)}(a) = \lim_{z\to a}\varphi^{(m-1)}(z)$ と極限で書けます。

**(2)**：(1) で $m = 1$ とすると $0! = 1$、$0$ 階微分は恒等写像なので、$\operatorname{Res}(f,a) = \lim_{z\to a}(z-a)f(z)$ です。

**(3)**：$h(a) = 0$ かつ $h'(a) \ne 0$ より $a$ は $h$ の $1$ 位の零点で、$a$ の近傍で $h(z) = (z-a)h_1(z)$、$h_1$ は正則で $h_1(a) = h'(a) \ne 0$ と書けます。よって $f(z) = \dfrac{g(z)/h_1(z)}{z-a}$ であり、$g/h_1$ は $a$ で正則、$a$ での値は $g(a)/h'(a) \ne 0$ ですから、<Ref to="prop-pole-characterization" /> (1)(iii) より $a$ は単純極です。(2) を使って

$$
\operatorname{Res}(f,a) = \lim_{z\to a}(z-a)f(z) = \lim_{z\to a}\frac{g(z)}{h_1(z)} = \frac{g(a)}{h'(a)}
$$

を得ます。
</Proof>

<Aside type="tip">
(1) は「**高々** $m$ 位」で使えるのが便利な点です。実際の位数が $m$ より小さくても、$n < -m$ の係数が $0$ である以上、証明はそのまま通ります。位数が確定できないときは大きめの $m$ を取っておけば安全です（そのぶん微分の回数は増えます）。
</Aside>

<Example id="ex-residue-computations" title="留数の計算例">
**(a)** $f(z) = \dfrac{z}{(z-1)(z+2)^2}$。$z=1$ は単純極なので <Ref to="prop-residue-formulas" /> (2) より

$$
\operatorname{Res}(f,1) = \lim_{z\to1}\frac{z}{(z+2)^2} = \frac{1}{9}.
$$

$z=-2$ は $2$ 位の極なので (1) で $m=2$ とし、

$$
\operatorname{Res}(f,-2) = \lim_{z\to-2}\frac{d}{dz}\left[\frac{z}{z-1}\right]
= \lim_{z\to-2}\frac{(z-1) - z}{(z-1)^2} = \frac{-1}{9}.
$$

二つの留数の和が $0$ になっているのは偶然ではありません。分母の次数が分子より $2$ 以上大きい有理関数では、後の <Ref to="lem-arc-vanishing" /> (i) と <Ref to="thm-residue" /> を大きな円周に適用すると留数の総和が $0$ と分かります。検算に使えます。

**(b)** $f(z) = \dfrac{e^z}{z^3}$。$m=3$ として

$$
\operatorname{Res}(f,0) = \frac{1}{2!}\lim_{z\to0}\frac{d^2}{dz^2}\bigl[e^z\bigr] = \frac{1}{2}.
$$

ローラン展開 $e^z/z^3 = z^{-3} + z^{-2} + \frac{1}{2}z^{-1} + \frac{1}{6} + \cdots$ の $z^{-1}$ の係数と一致します。

**(c)** $f(z) = \cot z = \dfrac{\cos z}{\sin z}$ の $z = k\pi$。$g = \cos$、$h = \sin$ とすると $g(k\pi) = (-1)^k \ne 0$、$h(k\pi) = 0$、$h'(k\pi) = \cos(k\pi) = (-1)^k \ne 0$ なので <Ref to="prop-residue-formulas" /> (3) より

$$
\operatorname{Res}(\cot, k\pi) = \frac{\cos(k\pi)}{\cos(k\pi)} = 1.
$$

すべての極で留数が $1$ です。

**(d)** $f(z) = z^2 e^{1/z}$ の $z = 0$。真性特異点なので (1)〜(3) は使えません。定義に戻ってローラン展開します。

$$
z^2 e^{1/z} = z^2\sum_{n\ge0}\frac{1}{n!}z^{-n} = \sum_{n \ge 0}\frac{1}{n!}z^{2-n}
$$

で、$z^{-1}$ が現れるのは $2-n = -1$、すなわち $n = 3$ のときですから $\operatorname{Res}(f,0) = \frac{1}{3!} = \frac{1}{6}$ です。留数は極の専売特許ではありません。
</Example>

## 5. 留数定理

<Theorem id="thm-residue" title="留数定理">
$D \subset \mathbb{C}$ を単連結領域とする。$C$ を $D$ に含まれる区分的に $C^1$ 級の単純閉曲線とし、正の向き（内部を左に見る向き）を与える。$C$ の内部（ジョルダンの曲線定理でいう有界成分）を $\operatorname{int}C$ と書く。相異なる有限個の点 $a_1,\dots,a_n \in \operatorname{int}C$ を取り、$f$ は $D \setminus \{a_1,\dots,a_n\}$ で正則とする。このとき

$$
\oint_C f(z)\,dz = 2\pi i \sum_{j=1}^{n}\operatorname{Res}(f,a_j).
$$
</Theorem>

<Proof of="thm-residue">
$D$ は単連結なので $\operatorname{int}C \subset D$ です（単連結領域内の閉曲線の内部は領域に含まれます）。したがって各 $a_j$ は $f$ の孤立特異点で、留数が定義されています。

**小円を取る**。$\varepsilon > 0$ を

$$
\varepsilon < \tfrac{1}{2}\min_{j \ne k}|a_j - a_k|,
\qquad
\varepsilon < \min_{j}\operatorname{dist}(a_j, \partial(\operatorname{int}C))
$$

を満たすように取ります（$a_j$ は有限個で、$\operatorname{int}C$ は開集合なので、このような $\varepsilon$ が存在します）。$\gamma_j$ を正の向きの円周 $|z - a_j| = \varepsilon$ とすると、閉円板 $\overline{D(a_j,\varepsilon)}$ は互いに交わらず、すべて $\operatorname{int}C$ に含まれます。

**積分路を分解する**。$f$ は「$C$ の内側かつすべての $\gamma_j$ の外側」の閉領域を含む開集合 $D \setminus \{a_1,\dots,a_n\}$ で正則です。§2 の多重連結領域に対するコーシーの積分定理より

$$
\oint_C f(z)\,dz = \sum_{j=1}^{n}\oint_{\gamma_j}f(z)\,dz .
$$

**小円上で項別積分する**。$j$ を固定します。$f$ は $D'(a_j, R_j)$ でローラン展開 $f(z) = \sum_{k=-\infty}^{\infty}c^{(j)}_k(z-a_j)^k$ を持ちます。必要なら $\varepsilon$ をさらに小さくして $\varepsilon < R_j$ としておきます。円周 $\gamma_j$ は円環 $D'(a_j,R_j)$ のコンパクト部分集合なので、この級数は $\gamma_j$ 上で一様収束します。一様収束する連続関数列は長さ有限の曲線上で項別積分できる（<Ref to="mathematics/complex-analysis/properties-of-holomorphic-functions#lem-termwise-integration" />）ので

$$
\oint_{\gamma_j} f(z)\,dz = \sum_{k=-\infty}^{\infty}c^{(j)}_k \oint_{\gamma_j}(z-a_j)^k\,dz .
$$

右辺の積分を計算します。$z = a_j + \varepsilon e^{i\theta}$（$0 \le \theta \le 2\pi$）と置くと $dz = i\varepsilon e^{i\theta}d\theta$ なので

$$
\oint_{\gamma_j}(z-a_j)^k dz
= \int_0^{2\pi}\varepsilon^k e^{ik\theta}\cdot i\varepsilon e^{i\theta}\,d\theta
= i\,\varepsilon^{\,k+1}\int_0^{2\pi}e^{i(k+1)\theta}\,d\theta .
$$

$k \ne -1$ のとき $\int_0^{2\pi}e^{i(k+1)\theta}d\theta = \left[\dfrac{e^{i(k+1)\theta}}{i(k+1)}\right]_0^{2\pi} = 0$（$e^{2\pi i(k+1)} = 1$）です。$k = -1$ のときは被積分関数が $1$ なので $2\pi$ となり、積分は $i \cdot \varepsilon^0 \cdot 2\pi = 2\pi i$ です。よって

$$
\oint_{\gamma_j}f(z)\,dz = 2\pi i\, c^{(j)}_{-1} = 2\pi i \operatorname{Res}(f,a_j)
$$

となり（最後の等号は <Ref to="def-residue" />）、これを $j$ について足せば結論を得ます。
</Proof>

<Remark id="rem-orientation">
向きは結果の符号を決めるので常に確認してください。$C$ を逆向きに回れば右辺は $-2\pi i \sum \operatorname{Res}$ になります。また、特異点が $C$ の**上**にある場合はこの定理は適用できません。実軸上に極がある実積分を扱う §8 では、輪郭を変形して極を避ける必要が出てきます。
</Remark>

ここから先は応用です。実積分を計算するときの手順はいつも同じで、次の三つに分かれます。

| 実積分の型 | 使う輪郭 | 余分な部分を消す道具 |
|---|---|---|
| $\int_{-\infty}^{\infty}P(x)/Q(x)\,dx$（$\deg Q \ge \deg P + 2$） | 上半平面の半円 | <Ref to="lem-arc-vanishing" /> (i) |
| $\int_0^{2\pi}R(\cos\theta,\sin\theta)\,d\theta$ | 単位円 $\lvert z\rvert = 1$ | 不要（輪郭が最初から閉じている） |
| $\int_{-\infty}^{\infty}f(x)e^{i\lambda x}dx$（$\lambda > 0$） | 上半平面の半円 | <Ref to="lem-arc-vanishing" /> (ii)（ジョルダンの補題） |
| 実軸上に単純極があるとき | 極を小円でくぼませた輪郭 | <Ref to="lem-indentation" /> と主値 |
| $\int_0^{\infty}x^{a-1}g(x)\,dx$（$a$ は非整数） | 鍵穴（Appendix） | 分枝の指定と大小円の評価 |

## 6. 有理関数の広義積分

<Figure caption="有理関数の広義積分に使う輪郭。実軸上の線分と上半平面の半円弧をつなぎ、上半平面にある極（×）だけを内部に取り込みます。半径 R を大きくすると弧の寄与が消えます。">
<svg viewBox="0 0 440 235" width="100%" role="img" aria-label="上半平面の半円輪郭と、その内部にある二つの極">
  <line x1="20" y1="185" x2="425" y2="185" stroke="currentColor" stroke-width="1" opacity="0.45" />
  <line x1="220" y1="222" x2="220" y2="18" stroke="currentColor" stroke-width="1" opacity="0.45" />
  <path d="M 60 185 A 160 160 0 0 1 380 185" fill="none" stroke="var(--sl-color-accent)" stroke-width="2.5" />
  <line x1="60" y1="185" x2="380" y2="185" stroke="var(--sl-color-accent)" stroke-width="2.5" />
  <polygon points="152,185 138,180 138,190" fill="var(--sl-color-accent)" />
  <polygon points="212,25 226,20 226,30" fill="var(--sl-color-accent)" />
  <line x1="165" y1="110" x2="175" y2="120" stroke="currentColor" stroke-width="1.8" />
  <line x1="165" y1="120" x2="175" y2="110" stroke="currentColor" stroke-width="1.8" />
  <line x1="267" y1="80" x2="277" y2="90" stroke="currentColor" stroke-width="1.8" />
  <line x1="267" y1="90" x2="277" y2="80" stroke="currentColor" stroke-width="1.8" />
  <text x="160" y="107" font-size="13" fill="currentColor" text-anchor="end">a₁</text>
  <text x="262" y="77" font-size="13" fill="currentColor" text-anchor="end">a₂</text>
  <text x="60" y="203" font-size="13" fill="currentColor" text-anchor="middle">−R</text>
  <text x="380" y="203" font-size="13" fill="currentColor" text-anchor="middle">R</text>
  <text x="230" y="201" font-size="13" fill="currentColor">O</text>
  <text x="345" y="58" font-size="14" fill="var(--sl-color-accent)">Γ<tspan font-size="10" dy="4">R</tspan></text>
  <text x="120" y="176" font-size="13" fill="var(--sl-color-accent)">実軸上の線分</text>
</svg>
</Figure>

輪郭を閉じるために付け足した半円弧の寄与は、半径を大きくすると消えてほしいわけです。それを保証するのが次の補題です。(i) と (ii) の仮定の違いに注意してください。(ii) のほうが $f$ に対する要求がゆるく、そのぶん $e^{i\lambda z}$ という因子が必要です。

<Lemma id="lem-arc-vanishing" title="大円弧の寄与が消える条件">
$R_0 > 0$ とし、$\Gamma_R$ を上半平面の半円弧 $z = Re^{i\theta}$（$0 \le \theta \le \pi$）とする。$f$ は $\{z : |z| \ge R_0,\ \operatorname{Im}z \ge 0\}$ で連続とし、$M(R) := \max_{0 \le \theta \le \pi}|f(Re^{i\theta})|$ とおく。

**(i)** $\lim_{R\to\infty} R\,M(R) = 0$ ならば $\displaystyle\lim_{R\to\infty}\int_{\Gamma_R}f(z)\,dz = 0$。

**(ii)（ジョルダンの補題）** $\lambda > 0$ とする。$R > R_0$ のとき

$$
\left|\int_{\Gamma_R}f(z)e^{i\lambda z}\,dz\right| \le \frac{\pi}{\lambda}M(R)
$$

が成り立つ。特に $\lim_{R\to\infty}M(R) = 0$ ならば $\displaystyle\lim_{R\to\infty}\int_{\Gamma_R}f(z)e^{i\lambda z}dz = 0$。
</Lemma>

<Proof of="lem-arc-vanishing">
**(i)**：$\Gamma_R$ の長さは $\pi R$ なので、§2 の ML 不等式より

$$
\left|\int_{\Gamma_R}f(z)\,dz\right| \le \pi R \cdot M(R) \longrightarrow 0 .
$$

**(ii)**：$z = Re^{i\theta} = R\cos\theta + iR\sin\theta$ に対し

$$
|e^{i\lambda z}| = |e^{i\lambda R\cos\theta}|\cdot|e^{-\lambda R\sin\theta}| = e^{-\lambda R\sin\theta}
$$

です（第一の因子は絶対値 $1$）。$|dz| = R\,d\theta$ に注意して

$$
\left|\int_{\Gamma_R}f(z)e^{i\lambda z}dz\right|
\le \int_0^{\pi}M(R)\,e^{-\lambda R \sin\theta}\,R\,d\theta
= 2RM(R)\int_0^{\pi/2}e^{-\lambda R\sin\theta}\,d\theta
$$

となります（$\sin(\pi-\theta) = \sin\theta$ による対称性で $[0,\pi/2]$ の積分の 2 倍にしました）。

ここでジョルダンの不等式 $\sin\theta \ge \dfrac{2\theta}{\pi}$（$0 \le \theta \le \pi/2$）を使います。証明しておきます。$\psi(\theta) = \sin\theta/\theta$（$\theta > 0$）とおくと $\psi'(\theta) = \dfrac{\theta\cos\theta - \sin\theta}{\theta^2}$ です。分子を $u(\theta) = \theta\cos\theta - \sin\theta$ とすると $u(0) = 0$、$u'(\theta) = \cos\theta - \theta\sin\theta - \cos\theta = -\theta\sin\theta < 0$（$0 < \theta \le \pi/2$）なので $u < 0$、したがって $\psi$ は $(0,\pi/2]$ で狭義単調減少です。よって $\psi(\theta) \ge \psi(\pi/2) = \dfrac{1}{\pi/2} = \dfrac{2}{\pi}$、すなわち $\sin\theta \ge 2\theta/\pi$ です。

これを使うと $e^{-\lambda R\sin\theta} \le e^{-2\lambda R\theta/\pi}$ なので

$$
\int_0^{\pi/2}e^{-\lambda R\sin\theta}d\theta
\le \int_0^{\pi/2}e^{-2\lambda R\theta/\pi}d\theta
= \frac{\pi}{2\lambda R}\left(1 - e^{-\lambda R}\right)
\le \frac{\pi}{2\lambda R}
$$

となり、これを代入して $2RM(R)\cdot\dfrac{\pi}{2\lambda R} = \dfrac{\pi}{\lambda}M(R)$ を得ます。
</Proof>

<Proposition id="prop-rational-integral" title="有理関数の広義積分">
$P, Q$ を複素係数の多項式とし、次を仮定する。

- $\deg Q \ge \deg P + 2$。
- $Q$ は実軸上に零点を持たない。

このとき $f = P/Q$ に対し $\int_{-\infty}^{\infty}f(x)\,dx$ は絶対収束し、

$$
\int_{-\infty}^{\infty}\frac{P(x)}{Q(x)}\,dx = 2\pi i \sum_{\operatorname{Im}a > 0}\operatorname{Res}(f, a)
$$

が成り立つ。和は上半平面にある $Q$ のすべての零点 $a$ にわたる。
</Proposition>

<Proof of="prop-rational-integral">
$p = \deg P$、$q = \deg Q$、最高次係数を $\alpha, \beta$ とします。$|z|\to\infty$ のとき $z^{\,q-p}P(z)/Q(z) \to \alpha/\beta$ なので、ある $R_1 \ge 1$ と $K > 0$ があって

$$
|z| \ge R_1 \implies |f(z)| \le K|z|^{\,p-q} \le K|z|^{-2}
$$

です（最後の不等号に $q - p \ge 2$ と $|z| \ge R_1 \ge 1$ を使いました）。したがって $\int_{|x| \ge R_1}|f(x)|dx \le 2\int_{R_1}^{\infty}Kx^{-2}dx < \infty$ であり、$[-R_1,R_1]$ 上では $f$ が連続（$Q$ は実軸上で $0$ にならない）なので、広義積分は絶対収束します。また $R \ge R_1$ に対し $R\,M(R) \le K/R \to 0$ なので <Ref to="lem-arc-vanishing" /> (i) の仮定が満たされます。

$Q$ の零点は有限個です。$R > \max(R_1, \max_a|a|)$ を取り、$C_R$ を線分 $[-R,R]$ と半円弧 $\Gamma_R$ をつないだ正の向きの単純閉曲線とします。$C_R$ とその内部は単連結領域 $\mathbb{C}$ に含まれ、$f$ は $Q$ の零点を除いて正則で、$C_R$ の内部にある零点はちょうど上半平面の零点全部です（実軸上には零点がなく、$R$ の取り方から下半平面の零点は内部に入りません）。<Ref to="thm-residue" /> より

$$
\int_{-R}^{R}f(x)\,dx + \int_{\Gamma_R}f(z)\,dz = 2\pi i \sum_{\operatorname{Im}a > 0}\operatorname{Res}(f,a)
$$

です。右辺は $R$ によりません。$R \to \infty$ とすると、左辺第 1 項は絶対収束性より $\int_{-\infty}^{\infty}f(x)dx$ に収束し、第 2 項は <Ref to="lem-arc-vanishing" /> (i) より $0$ に収束します。
</Proof>

<Example id="ex-x4-plus-1" title="分母が 4 次の場合">
$\displaystyle\int_{-\infty}^{\infty}\frac{dx}{x^4+1}$ を計算します。$P = 1$、$Q = z^4+1$ で $\deg Q - \deg P = 4 \ge 2$、$Q(x) = x^4 + 1 \ge 1 > 0$ より実軸上に零点はないので <Ref to="prop-rational-integral" /> が使えます。

$z^4 = -1 = e^{i\pi}$ の解は $z_k = e^{i\pi(2k+1)/4}$（$k = 0,1,2,3$）、すなわち $e^{i\pi/4},\ e^{3i\pi/4},\ e^{5i\pi/4},\ e^{7i\pi/4}$ です。上半平面にあるのは偏角が $\pi/4$ と $3\pi/4$ の二つです。いずれも $Q$ の単純零点（$Q' = 4z^3$ は $0$ にならない）なので、<Ref to="prop-residue-formulas" /> (3) より、$a$ をそのどちらかとして

$$
\operatorname{Res}(f,a) = \frac{1}{4a^3} = \frac{a}{4a^4} = \frac{a}{4\cdot(-1)} = -\frac{a}{4}
$$

です（$a^4 = -1$ を使って $a^3$ を $a$ に直しました。この書き換えで計算がぐっと楽になります）。したがって

$$
\sum \operatorname{Res} = -\frac{1}{4}\left(e^{i\pi/4} + e^{3i\pi/4}\right)
= -\frac{1}{4}\left(\frac{\sqrt2}{2}(1+i) + \frac{\sqrt2}{2}(-1+i)\right)
= -\frac{1}{4}\cdot\sqrt{2}\,i
$$

となり、

$$
\int_{-\infty}^{\infty}\frac{dx}{x^4+1} = 2\pi i \cdot\left(-\frac{\sqrt2 i}{4}\right) = \frac{2\pi\sqrt2}{4} = \frac{\pi}{\sqrt2} \approx 2.2214 .
$$

被積分関数は正なので答えが正の実数になったのは整合的です。また $|x| \ge 1$ では $\frac{1}{x^4+1} \le \frac{1}{x^2+1}$、$|x| \le 1$ では $\frac{1}{x^4+1} \le 1$ なので、積分値は $2 + \int_{|x|\ge1}\frac{dx}{x^2+1} = 2 + \frac{\pi}{2} \approx 3.57$ 以下と分かり、$2.22$ という値と矛盾しません。
</Example>

## 7. 三角関数の周期積分：単位円を回る

$\theta$ について周期 $2\pi$ の三角有理式の積分 $\int_0^{2\pi}R(\cos\theta,\sin\theta)\,d\theta$ は、実積分が最初から閉じた区間を一周しているので、輪郭を付け足す必要すらありません。$z = e^{i\theta}$ と置けば、これはそのまま単位円周に沿う複素積分になります。オイラーの公式から

$$
\cos\theta = \frac{z + z^{-1}}{2},\qquad
\sin\theta = \frac{z - z^{-1}}{2i},\qquad
d\theta = \frac{dz}{iz}
$$

です。最後の関係は、単位円周を $\gamma(\theta) = e^{i\theta}$（$0 \le \theta \le 2\pi$）で径数付けると $\gamma'(\theta) = ie^{i\theta} = iz$ なので、任意の連続関数 $F$ に対し

$$
\oint_{|z|=1}\frac{F(z)}{iz}\,dz = \int_0^{2\pi}\frac{F(e^{i\theta})}{ie^{i\theta}}\,ie^{i\theta}\,d\theta = \int_0^{2\pi}F(e^{i\theta})\,d\theta
$$

が成り立つ、という意味です。

<Proposition id="prop-trig-integral" title="a + b cosθ 型の積分">
$a, b$ を実数とし $a > |b| \ge 0$ とする。このとき $a + b\cos\theta \ge a - |b| > 0$ で被積分関数は連続であり、

$$
\int_0^{2\pi}\frac{d\theta}{a + b\cos\theta} = \frac{2\pi}{\sqrt{a^2-b^2}} .
$$
</Proposition>

<Proof of="prop-trig-integral">
$b = 0$ のときは左辺 $= 2\pi/a$、右辺 $= 2\pi/\sqrt{a^2} = 2\pi/a$（$a > 0$）で成立します。$b < 0$ のときは $\theta \mapsto \theta + \pi$ と置換すると $\cos(\theta+\pi) = -\cos\theta$ で、被積分関数は周期 $2\pi$ なので積分区間は $[0,2\pi]$ のままでよく、$b$ を $-b$ に取り替えた式に帰着します。右辺は $b$ の符号によらないので、以下 $b > 0$ とします。

$z = e^{i\theta}$ と置く上の対応により

$$
I = \oint_{|z|=1}\frac{1}{a + \frac{b}{2}\left(z + \frac1z\right)}\cdot\frac{dz}{iz}
= \oint_{|z|=1}\frac{2\,dz}{i\left(bz^2 + 2az + b\right)} .
$$

分母の 2 次式 $bz^2+2az+b$ の零点は

$$
z_{\pm} = \frac{-a \pm \sqrt{a^2-b^2}}{b}
$$

で、$a > b > 0$ より $a^2 - b^2 > 0$ なのでどちらも実数、かつ相異なります。解と係数の関係から $z_+z_- = b/b = 1$ です。

$|z_+| < 1$ を確かめます。$z_+ = \dfrac{-a+\sqrt{a^2-b^2}}{b}$ は $\sqrt{a^2-b^2} < a$ より負で、$|z_+| = \dfrac{a - \sqrt{a^2-b^2}}{b}$ です。$|z_+| < 1 \iff a - b < \sqrt{a^2-b^2}$ であり、両辺とも正（$a>b$）なので 2 乗して $\ (a-b)^2 < (a-b)(a+b) \iff a-b < a+b \iff b > 0$。これは仮定です。よって $|z_+|<1$、そして $|z_-| = 1/|z_+| > 1$ です。特に単位円周上に零点はありません。

したがって単位円の内部にある特異点は単純極 $z_+$ だけです。$g(z) = \dfrac{2}{i(bz^2+2az+b)}$ とすると、<Ref to="prop-residue-formulas" /> (2) より

$$
\operatorname{Res}(g, z_+) = \lim_{z\to z_+}\frac{2(z-z_+)}{i\,b(z-z_+)(z-z_-)}
= \frac{2}{i\,b(z_+-z_-)}
= \frac{2}{i\,b\cdot\frac{2\sqrt{a^2-b^2}}{b}}
= \frac{1}{i\sqrt{a^2-b^2}} .
$$

<Ref to="thm-residue" /> より $I = 2\pi i \cdot \dfrac{1}{i\sqrt{a^2-b^2}} = \dfrac{2\pi}{\sqrt{a^2-b^2}}$ です。
</Proof>

たとえば $\displaystyle\int_0^{2\pi}\frac{d\theta}{5+4\cos\theta} = \frac{2\pi}{\sqrt{25-16}} = \frac{2\pi}{3}$ です。次は分子に $\cos 2\theta$ が入る場合で、$0$ に多重極が現れます。

<Example id="ex-cos2theta" title="分子に cos 2θ がある積分">
$\displaystyle I = \int_0^{2\pi}\frac{\cos2\theta}{5-4\cos\theta}\,d\theta$ を計算します。$5 - 4\cos\theta \ge 1 > 0$ なので被積分関数は連続です。$z = e^{i\theta}$ とすると $\cos2\theta = \frac{z^2+z^{-2}}{2}$、$5-4\cos\theta = 5 - 2(z+z^{-1}) = \dfrac{-2z^2+5z-2}{z}$ なので

$$
I = \oint_{|z|=1}\frac{(z^2+z^{-2})/2}{(-2z^2+5z-2)/z}\cdot\frac{dz}{iz}
= \oint_{|z|=1}\frac{z^4+1}{2z^2\,i\,(-2z^2+5z-2)}\,dz
= \frac{i}{2}\oint_{|z|=1}\frac{z^4+1}{z^2(2z-1)(z-2)}\,dz
$$

です（$-2z^2+5z-2 = -(2z-1)(z-2)$ と $\frac{1}{-2i} = \frac{i}{2}$ を使いました）。被積分関数を $g$ とすると、単位円の内部の特異点は $z=0$（$2$ 位の極）と $z=1/2$（単純極）で、$z=2$ は外部です。

$z = 1/2$：<Ref to="prop-residue-formulas" /> (2) より

$$
\operatorname{Res}(g,\tfrac12) = \frac{z^4+1}{z^2\cdot 2(z-2)}\bigg|_{z=1/2}
= \frac{\frac{1}{16}+1}{\frac14 \cdot 2\left(\frac12-2\right)}
= \frac{17/16}{-3/4} = -\frac{17}{12}.
$$

$z = 0$：$m=2$ として $\varphi(z) = z^2 g(z) = \dfrac{z^4+1}{2z^2-5z+2}$ とおくと

$$
\varphi'(z) = \frac{4z^3(2z^2-5z+2) - (z^4+1)(4z-5)}{(2z^2-5z+2)^2},
\qquad
\varphi'(0) = \frac{0 - 1\cdot(-5)}{4} = \frac{5}{4},
$$

なので $\operatorname{Res}(g,0) = \varphi'(0)/1! = 5/4$ です。念のためローラン展開でも確かめます。$\dfrac{1}{2-5z+2z^2} = \dfrac{1}{2}\cdot\dfrac{1}{1 - \frac{5z-2z^2}{2}} = \dfrac12\left(1 + \frac{5z}{2} + O(z^2)\right) = \frac12 + \frac{5}{4}z + O(z^2)$ なので、$(z^4+1)$ を掛けて $z^2$ で割ると $g(z) = \dfrac{1}{2z^2} + \dfrac{5}{4z} + \cdots$ となり、$z^{-1}$ の係数は $5/4$ で一致します。

したがって留数の和は $-\frac{17}{12} + \frac{5}{4} = -\frac{17}{12}+\frac{15}{12} = -\frac16$ で、<Ref to="thm-residue" /> より

$$
I = \frac{i}{2}\cdot 2\pi i \cdot\left(-\frac16\right) = -\pi\cdot\left(-\frac16\right) = \frac{\pi}{6}.
$$
</Example>

## 8. 指数因子と実軸上の極

### 8.1. フーリエ型の積分

$\int_{-\infty}^{\infty}f(x)\cos\lambda x\,dx$ のような積分では、$\cos\lambda z$ をそのまま複素平面に延長してはいけません。上半平面では $\cos(iy) = \cosh y$ が $y \to \infty$ で指数的に発散するからです。代わりに $e^{i\lambda z}$ を使います。これは上半平面で $|e^{i\lambda z}| = e^{-\lambda\operatorname{Im}z} \le 1$（$\lambda>0$）と抑えられ、しかも実軸上では $\cos\lambda x + i\sin\lambda x$ なので、最後に実部と虚部を取れば目的の積分が得られます。

<Aside type="caution">
上半平面で有界なのは $e^{i\lambda z}$（$\lambda > 0$）です。$\lambda < 0$ のときは下半平面（時計回りの輪郭、留数の和に $-2\pi i$ が掛かる）を使ってください。$\cos\lambda z$ や $\sin\lambda z$ を輪郭に載せると、円弧上の評価が破綻します。
</Aside>

<Example id="ex-fourier-cos" title="ローレンツ型のフーリエ変換">
$\displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{x^2+1}\,dx$ を求めます。$f(z) = \dfrac{1}{z^2+1}$、$\lambda = 1$ として $F(z) = f(z)e^{iz}$ を <Ref to="ex-x4-plus-1" /> と同じ半円輪郭 $C_R$（$R>1$）で積分します。$F$ の特異点は $z = \pm i$ の単純極で、上半平面にあるのは $z = i$ だけです。<Ref to="prop-residue-formulas" /> (3) を $g(z) = e^{iz}$、$h(z)=z^2+1$、$h'(z) = 2z$ に適用して

$$
\operatorname{Res}(F, i) = \frac{e^{i\cdot i}}{2i} = \frac{e^{-1}}{2i}.
$$

<Ref to="thm-residue" /> より $\displaystyle\int_{-R}^{R}\frac{e^{ix}}{x^2+1}dx + \int_{\Gamma_R}F(z)\,dz = 2\pi i\cdot\frac{e^{-1}}{2i} = \frac{\pi}{e}$ です。

円弧の評価には <Ref to="lem-arc-vanishing" /> (ii) を使います。$|z| = R > 1$ のとき $|f(z)| = \dfrac{1}{|z^2+1|} \le \dfrac{1}{R^2-1}$（三角不等式 $|z^2+1| \ge |z|^2 - 1$）なので $M(R) \le \dfrac{1}{R^2-1}\to0$ であり、$\left|\int_{\Gamma_R}F\right| \le \dfrac{\pi}{1}\cdot\dfrac{1}{R^2-1}\to0$ です。よって

$$
\int_{-\infty}^{\infty}\frac{e^{ix}}{x^2+1}\,dx = \frac{\pi}{e} .
$$

実部を取ると $\displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{x^2+1}dx = \frac{\pi}{e}$、虚部を取ると $\displaystyle\int_{-\infty}^{\infty}\frac{\sin x}{x^2+1}dx = 0$ です。後者は被積分関数が奇関数であることからも分かり、検算になります。なお $\dfrac{\cos x}{x^2+1}$ の絶対値は $\dfrac{1}{x^2+1}$ 以下なので、この広義積分は絶対収束しています。
</Example>

### 8.2. 実軸上の単純極と主値

被積分関数が実軸上に極を持つと、広義積分は普通の意味では発散します。それでも「対称に切り取れば」極限が存在することがあります。

<Definition id="def-principal-value" title="コーシーの主値">
$c \in \mathbb{R}$ とし、$f$ は $\mathbb{R}\setminus\{c\}$ で連続とする。極限

$$
\mathrm{p.v.}\int_{-\infty}^{\infty}f(x)\,dx
:= \lim_{\substack{R \to \infty \\ \varepsilon \to +0}}
\left(\int_{-R}^{c-\varepsilon}f(x)\,dx + \int_{c+\varepsilon}^{R}f(x)\,dx\right)
$$

が存在するとき、その値を $f$ の**コーシーの主値**と呼ぶ。左右を同じ $\varepsilon$ で切り落とすのが要点である。通常の広義積分が存在するときは、主値も存在して同じ値になる。
</Definition>

輪郭のほうは、実軸上の極を半径 $\varepsilon$ の小円弧で迂回して避けます。その小円弧の寄与を測るのが次の補題です。「一周すれば $2\pi i \times$ 留数」だったものが、「角度 $\beta-\alpha$ だけ回れば $i(\beta-\alpha)\times$ 留数」になる、と憶えると自然です。

<Lemma id="lem-indentation" title="単純極まわりの小円弧（くぼみ補題）">
$c$ を $f$ の単純極とし、$\rho = \operatorname{Res}(f,c)$ とおく。実数 $\alpha, \beta$ に対し、$\gamma_\varepsilon$ を $z = c + \varepsilon e^{i\theta}$（$\theta$ が $\alpha$ から $\beta$ まで動く）で定まる円弧とする。このとき

$$
\lim_{\varepsilon \to +0}\int_{\gamma_\varepsilon}f(z)\,dz = i(\beta - \alpha)\rho .
$$

（$\beta > \alpha$ なら反時計回り、$\beta < \alpha$ なら時計回りである。）
</Lemma>

<Proof of="lem-indentation">
$c$ は単純極なので、ローラン展開は $f(z) = \dfrac{\rho}{z-c} + \sum_{n \ge 0}c_n(z-c)^n$ の形です。第 2 項を $g(z)$ とおくと、$g$ は $D(c,R)$ で正則なべき級数の和ですから、ある $r>0$ と $M>0$ に対し $|z-c| \le r$ で $|g(z)| \le M$ です（連続関数のコンパクト集合上の有界性）。

主要部の寄与は $\varepsilon$ によらず計算できます。$z = c+\varepsilon e^{i\theta}$、$dz = i\varepsilon e^{i\theta}d\theta$ より

$$
\int_{\gamma_\varepsilon}\frac{\rho}{z-c}\,dz
= \rho\int_\alpha^\beta \frac{i\varepsilon e^{i\theta}}{\varepsilon e^{i\theta}}\,d\theta
= i\rho(\beta-\alpha).
$$

残りは ML 不等式で消えます。$\gamma_\varepsilon$ の長さは $|\beta-\alpha|\varepsilon$ なので、$\varepsilon \le r$ のとき

$$
\left|\int_{\gamma_\varepsilon}g(z)\,dz\right| \le |\beta-\alpha|\,\varepsilon\, M \longrightarrow 0
\qquad (\varepsilon \to +0).
$$

両者を足して結論を得ます。
</Proof>

<Example id="ex-dirichlet" title="ディリクレ積分">
$\displaystyle\int_0^{\infty}\frac{\sin x}{x}\,dx = \frac{\pi}{2}$ を示します。$f(z) = \dfrac{e^{iz}}{z}$ とし、$0 < \varepsilon < R$ に対して次の四つをつないだ輪郭 $C$ を考えます。線分 $[\varepsilon, R]$、大円弧 $\Gamma_R$（反時計回りに $0$ から $\pi$ へ）、線分 $[-R,-\varepsilon]$、小円弧 $\gamma_\varepsilon$（$\theta$ が $\pi$ から $0$ へ、すなわち時計回り）。

$C$ とその内部は、複素平面から半直線 $\{iy : y \le 0\}$ を除いた単連結領域 $D$ に含まれ、$f$ は $D$ で正則です（$0$ は $D$ に属しません）。よって <Ref to="thm-residue" /> を特異点 $0$ 個の場合に使って（これはコーシーの積分定理そのものです）

$$
\int_\varepsilon^R \frac{e^{ix}}{x}dx + \int_{\Gamma_R}f\,dz + \int_{-R}^{-\varepsilon}\frac{e^{ix}}{x}dx + \int_{\gamma_\varepsilon}f\,dz = 0 .
$$

第 3 項は $x = -t$ と置換すると $\displaystyle\int_{-R}^{-\varepsilon}\frac{e^{ix}}{x}dx = \int_{R}^{\varepsilon}\frac{e^{-it}}{-t}(-dt) = -\int_\varepsilon^R\frac{e^{-it}}{t}dt$ なので、第 1 項と合わせて

$$
\int_\varepsilon^R\frac{e^{ix}-e^{-ix}}{x}\,dx = 2i\int_\varepsilon^R\frac{\sin x}{x}\,dx
$$

になります。第 2 項は <Ref to="lem-arc-vanishing" /> (ii) で $f(z) = 1/z$、$\lambda = 1$ とすると $M(R) = 1/R \to 0$ なので $0$ に収束します。第 4 項は、$0$ が $f$ の単純極で $\operatorname{Res}(f,0) = \lim_{z\to0}z\cdot\frac{e^{iz}}{z} = 1$ ですから、<Ref to="lem-indentation" /> を $\alpha = \pi$、$\beta = 0$ として適用し $i(0-\pi)\cdot1 = -i\pi$ に収束します。

以上を合わせて $\varepsilon \to +0$、$R \to \infty$ とすると

$$
2i\int_0^{\infty}\frac{\sin x}{x}\,dx - i\pi = 0,
\qquad\text{すなわち}\qquad
\int_0^{\infty}\frac{\sin x}{x}\,dx = \frac{\pi}{2}.
$$

同じ計算を実部・虚部に分けて読むと $\mathrm{p.v.}\displaystyle\int_{-\infty}^{\infty}\frac{e^{ix}}{x}dx = i\pi$ であり、実部から $\mathrm{p.v.}\displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{x}dx = 0$（奇関数なので当然です）、虚部から上の結果が出ます。

注意すべき点が二つあります。第一に、$\dfrac{\sin x}{x}$ は $x=0$ で除去可能特異点を持つ（<Ref to="ex-classify" /> (a)）ので $\varepsilon \to 0$ の極限は問題なく、$\varepsilon$ を入れたのは $e^{iz}/z$ を扱うための便宜です。第二に、この広義積分は絶対収束しません（<Ref to="mathematics/calculus/borwein-integrals#prop-conditional" />）。実際 $\int_0^{\infty}\left|\frac{\sin x}{x}\right|dx = \infty$ であり、値は $R \to \infty$ の極限としてのみ意味を持ちます。実解析の側からの別証明は <Ref to="mathematics/calculus/borwein-integrals#thm-dirichlet" /> にあり、この繊細さが表に出る有名な例として [ボールウェイン積分](/mathematics/calculus/borwein-integrals) があります。
</Example>

## 9. 演習

<Exercise id="exr-residues" difficulty="易">
$f(z) = \dfrac{e^z}{z^2(z-1)}$ のすべての孤立特異点を求め、それぞれの位数と留数を求めてください。

<Solution>
特異点は分母の零点 $z=0$（$2$ 位）と $z=1$（$1$ 位）です。$e^z$ はこれらの点で $0$ にならないので、位数は分母のとおりです。実際 $z^2f(z) = \frac{e^z}{z-1} \to -1 \ne 0$、$(z-1)f(z) = \frac{e^z}{z^2}\to e \ne 0$ なので <Ref to="prop-pole-characterization" /> より $0$ は $2$ 位の極、$1$ は単純極です。

$z=1$：<Ref to="prop-residue-formulas" /> (2) より $\operatorname{Res}(f,1) = \lim_{z\to1}\frac{e^z}{z^2} = e$。

$z=0$：同 (1) で $m=2$ とし、$\varphi(z) = \dfrac{e^z}{z-1}$ に対し

$$
\varphi'(z) = \frac{e^z(z-1) - e^z}{(z-1)^2} = \frac{e^z(z-2)}{(z-1)^2},
\qquad
\operatorname{Res}(f,0) = \varphi'(0) = \frac{1\cdot(-2)}{1} = -2 .
$$

したがって留数の和は $e - 2$ です。
</Solution>
</Exercise>

<Exercise id="exr-even-rational" difficulty="標準">
$\displaystyle\int_0^{\infty}\frac{x^2}{(x^2+1)(x^2+4)}\,dx$ を留数定理で計算し、部分分数分解による初等的な計算と一致することを確かめてください。

<Solution>
被積分関数は偶関数なので $\displaystyle\int_0^{\infty} = \frac12\int_{-\infty}^{\infty}$ です。$f(z) = \dfrac{z^2}{(z^2+1)(z^2+4)}$ は $\deg$ の差が $2$、分母は実軸上で $0$ にならないので <Ref to="prop-rational-integral" /> が使えます。上半平面の極は単純極 $z=i$ と $z=2i$ です。

$$
\operatorname{Res}(f,i) = \frac{z^2}{(z+i)(z^2+4)}\bigg|_{z=i} = \frac{-1}{2i\cdot3} = \frac{-1}{6i} = \frac{i}{6},
$$

$$
\operatorname{Res}(f,2i) = \frac{z^2}{(z^2+1)(z+2i)}\bigg|_{z=2i} = \frac{-4}{(-3)\cdot4i} = \frac{1}{3i} = -\frac{i}{3}.
$$

（$\frac1i = -i$ を使いました。）和は $\frac{i}{6}-\frac{i}{3} = -\frac{i}{6}$ なので

$$
\int_{-\infty}^{\infty}f(x)\,dx = 2\pi i\left(-\frac{i}{6}\right) = \frac{\pi}{3},
\qquad
\int_0^{\infty}f(x)\,dx = \frac{\pi}{6}.
$$

初等的な確認：$\dfrac{x^2}{(x^2+1)(x^2+4)} = \dfrac{A}{x^2+1}+\dfrac{B}{x^2+4}$ とおくと $x^2 = A(x^2+4)+B(x^2+1)$ より $A+B=1$、$4A+B=0$、すなわち $A = -\frac13$、$B=\frac43$ です。$\int_0^{\infty}\frac{dx}{x^2+a^2} = \frac{\pi}{2a}$（$a>0$）を使うと

$$
-\frac13\cdot\frac{\pi}{2} + \frac43\cdot\frac{\pi}{4} = -\frac{\pi}{6}+\frac{\pi}{3} = \frac{\pi}{6}
$$

で一致します。
</Solution>
</Exercise>

<Exercise id="exr-trig-square" difficulty="標準">
$\displaystyle\int_0^{2\pi}\frac{d\theta}{(2+\cos\theta)^2}$ を計算してください。

<Solution>
$z = e^{i\theta}$ と置くと $2+\cos\theta = 2 + \dfrac{z+z^{-1}}{2} = \dfrac{z^2+4z+1}{2z}$ なので

$$
I = \oint_{|z|=1}\frac{4z^2}{(z^2+4z+1)^2}\cdot\frac{dz}{iz}
= \frac{4}{i}\oint_{|z|=1}\frac{z\,dz}{(z^2+4z+1)^2}.
$$

$z^2+4z+1 = 0$ の解は $z_0 = -2+\sqrt3 \approx -0.268$ と $z_1 = -2-\sqrt3 \approx -3.73$ で、単位円の内部にあるのは $z_0$ だけです（$z_0z_1 = 1$ なので一方が内側、他方が外側）。分母は $(z-z_0)^2(z-z_1)^2$ なので $z_0$ は $2$ 位の極です。$m=2$ として <Ref to="prop-residue-formulas" /> (1) を使います。$\varphi(z) = \dfrac{z}{(z-z_1)^2}$ とすると

$$
\varphi'(z) = \frac{(z-z_1)^2 - z\cdot2(z-z_1)}{(z-z_1)^4} = \frac{(z-z_1)-2z}{(z-z_1)^3} = \frac{-z-z_1}{(z-z_1)^3}.
$$

$z_0 - z_1 = 2\sqrt3$、$-z_0-z_1 = -(-2+\sqrt3)-(-2-\sqrt3) = 4$ なので

$$
\operatorname{Res} = \varphi'(z_0) = \frac{4}{(2\sqrt3)^3} = \frac{4}{24\sqrt3} = \frac{1}{6\sqrt3}.
$$

<Ref to="thm-residue" /> より

$$
I = \frac{4}{i}\cdot2\pi i\cdot\frac{1}{6\sqrt3} = \frac{8\pi}{6\sqrt3} = \frac{4\pi}{3\sqrt3} = \frac{4\sqrt3\,\pi}{9} \approx 2.418 .
$$

検算：<Ref to="prop-trig-integral" /> の式 $\int_0^{2\pi}\frac{d\theta}{a+b\cos\theta} = 2\pi(a^2-b^2)^{-1/2}$ を $a$ で微分すると、左辺は $-\int_0^{2\pi}\frac{d\theta}{(a+b\cos\theta)^2}$、右辺は $-2\pi a(a^2-b^2)^{-3/2}$ です。$a=2$、$b=1$ とすると $\frac{4\pi}{3^{3/2}} = \frac{4\pi}{3\sqrt3}$ で一致します。
</Solution>
</Exercise>

<Exercise id="exr-jordan" difficulty="難">
$a > 0$ とするとき $\displaystyle\int_0^{\infty}\frac{x\sin x}{x^2+a^2}\,dx = \frac{\pi}{2}e^{-a}$ を示してください。$\dfrac{z}{z^2+a^2}$ には <Ref to="lem-arc-vanishing" /> (i) が使えないことに注意してください。

<Solution>
$f(z) = \dfrac{z}{z^2+a^2}$、$F(z) = f(z)e^{iz}$ とし、半径 $R > a$ の半円輪郭 $C_R$ を取ります。$F$ の上半平面の特異点は単純極 $z = ia$ だけで、<Ref to="prop-residue-formulas" /> (3) を $g(z) = ze^{iz}$、$h(z) = z^2+a^2$、$h'(z)=2z$ に適用して

$$
\operatorname{Res}(F, ia) = \frac{ia\,e^{i\cdot ia}}{2ia} = \frac{e^{-a}}{2}.
$$

<Ref to="thm-residue" /> より $\displaystyle\int_{-R}^{R}\frac{xe^{ix}}{x^2+a^2}dx + \int_{\Gamma_R}F\,dz = 2\pi i\cdot\frac{e^{-a}}{2} = i\pi e^{-a}$ です。

円弧の評価：$|z|=R>a$ のとき $|f(z)| \le \dfrac{R}{R^2-a^2}$ なので $M(R) \to 0$ ですが、$R\,M(R) \to 1 \ne 0$ なので <Ref to="lem-arc-vanishing" /> (i) は使えません。指数因子を持つ (ii)（ジョルダンの補題）なら $\lambda = 1$ として

$$
\left|\int_{\Gamma_R}F(z)\,dz\right| \le \pi M(R) \le \frac{\pi R}{R^2-a^2}\longrightarrow 0
$$

と結論できます。よって $\displaystyle\lim_{R\to\infty}\int_{-R}^{R}\frac{xe^{ix}}{x^2+a^2}dx = i\pi e^{-a}$。虚部を取ると

$$
\int_{-\infty}^{\infty}\frac{x\sin x}{x^2+a^2}\,dx = \pi e^{-a}
$$

であり、被積分関数は偶関数（奇関数どうしの積）なので $\displaystyle\int_0^{\infty} = \frac{\pi}{2}e^{-a}$ です。なおこの積分は絶対収束せず、$R \to \infty$ の対称極限として存在します。
</Solution>
</Exercise>

## 参考文献

- L. V. Ahlfors, *Complex Analysis*, 3rd ed., McGraw-Hill, 1979 — 第 4 章に留数計算と定積分への応用がまとまっています。真性特異点の値の分布（ピカールの定理）は第 8 章です。
- E. M. Stein and R. Shakarchi, *Complex Analysis* (Princeton Lectures in Analysis II), Princeton University Press, 2003 — 第 3 章「Meromorphic Functions and the Logarithm」が留数公式と本記事で扱った型の実積分を扱っています。
- W. Rudin, *Real and Complex Analysis*, 3rd ed., McGraw-Hill, 1987 — 第 10 章に孤立特異点の分類と、巻き数を使った一般形の留数定理があります。
- 神保道夫『複素関数入門』岩波書店、2003 — 日本語で読める入門書。留数と実積分への応用を丁寧に扱っています。
- A.-L. Cauchy, *Exercices de mathématiques*, 1826 — 「résidu（留数）」という語と計算法が導入された原典です。

## Appendix: 分数べきを含む積分と鍵穴輪郭

**扱う積分**。$0 < a < 1$ に対して

$$
\int_0^{\infty}\frac{x^{a-1}}{1+x}\,dx = \frac{\pi}{\sin \pi a}
$$

を示します。$x^{a-1}$ は $a$ が整数でないと一価な複素関数に延長できないので、これまでの輪郭は使えません。切断（分枝の切れ目）を正の実軸に置き、その両岸を往復する「鍵穴輪郭」を使います。

**分枝を決める**。$\operatorname{Log}z = \log|z| + i\arg z$、$\arg z \in (0, 2\pi)$ と定め、$z^{a-1} := e^{(a-1)\operatorname{Log}z}$ とします。これは $\mathbb{C}$ から非負実軸を除いた領域で正則です。$f(z) = \dfrac{z^{a-1}}{1+z}$ の特異点は $z = -1 = e^{i\pi}$ の単純極だけです。

**輪郭**。$0 < \varepsilon < 1 < R$ とし、切断の上岸を $\varepsilon$ から $R$ まで進み、大円 $|z|=R$ を反時計回りに一周し、下岸を $R$ から $\varepsilon$ まで戻り、小円 $|z|=\varepsilon$ を時計回りに一周して閉じます。厳密には上岸・下岸を実軸から $\delta$ だけずらした線分に取り、最後に $\delta \to 0$ とします。この閉曲線の内部にある特異点は $z=-1$ だけです。

**円周上の評価**。大円上では $|z^{a-1}| = R^{a-1}$、$|1+z| \ge R-1$ なので、ML 不等式より寄与は $2\pi R\cdot\dfrac{R^{a-1}}{R-1} \sim 2\pi R^{a-1} \to 0$（$a<1$）です。小円上では $2\pi\varepsilon\cdot\dfrac{\varepsilon^{a-1}}{1-\varepsilon} \sim 2\pi\varepsilon^{a}\to 0$（$a>0$）です。$0 < a < 1$ という仮定は、この二つの評価を同時に成立させるためにあります。

**両岸の差**。上岸では $\arg z \to 0$ なので $z^{a-1} \to x^{a-1}$、下岸では $\arg z \to 2\pi$ なので $z^{a-1}\to x^{a-1}e^{2\pi i(a-1)} = x^{a-1}e^{2\pi ia}$ です。下岸は逆向きに進むので、二つの線分の寄与の和は $I := \int_0^{\infty}\frac{x^{a-1}}{1+x}dx$ を使って $\left(1 - e^{2\pi i a}\right)I$ と書けます。切断があるおかげで両岸が相殺せず、求めたい積分が残るのがこの輪郭の要点です。

**留数と結論**。$\operatorname{Res}(f,-1) = (-1)^{a-1} = e^{(a-1)i\pi} = -e^{i\pi a}$ なので、<Ref to="thm-residue" /> より $\left(1-e^{2\pi ia}\right)I = 2\pi i\left(-e^{i\pi a}\right)$、すなわち

$$
I = \frac{-2\pi i\,e^{i\pi a}}{1 - e^{2\pi i a}}
= \frac{-2\pi i}{e^{-i\pi a} - e^{i\pi a}}
= \frac{-2\pi i}{-2i\sin\pi a}
= \frac{\pi}{\sin\pi a}.
$$

**検算**。$a = 1/2$ とすると $I = \pi/\sin(\pi/2) = \pi$ です。一方 $x = t^2$ と置換すれば

$$
\int_0^{\infty}\frac{dx}{\sqrt{x}(1+x)} = \int_0^{\infty}\frac{2t\,dt}{t(1+t^2)} = 2\bigl[\arctan t\bigr]_0^{\infty} = \pi
$$

で一致します。この積分はベータ関数・ガンマ関数の相反公式 $\Gamma(a)\Gamma(1-a) = \pi/\sin\pi a$ の複素解析による証明にもなっています。


</div>
