Lemma数学ロードマップ

01 微分積分学 · 第 5 章

リーマン積分

目安 14〜19 時間定理など 16演習 10 問実験 1

この章の目標

  • ダルブーの上和・下和によってリーマン積分を定義し、可積分性の判定条件を使える
  • 連続関数・単調関数が可積分であることを証明できる
  • 微分積分学の基本定理を証明し、置換積分・部分積分・積分剰余のテイラーの定理を導ける
  • 広義積分の収束を判定でき、ガンマ関数・ベータ関数の基本的な関係式を導ける
  • ウォリスの公式とスターリングの公式を理解し、級数の積分判定法を使える

前提:第1章〜第4章

5.1 積分を定義し直す理由

高校では、ff の原始関数 FF(F′=fF' = f)を使って ∫abf(x) dx=F(b)−F(a)\int_a^b f(x)\ dx = F(b) - F(a) と定義し、それが「面積」を表すことを図で説明した。この定義にはいくつか問題がある。

  • 原始関数が式で書けない関数(例えば e−x2e^{-x^2})の積分は、そもそも存在するのか。
  • 不連続な関数(例えば階段関数)の積分はどう定義するのか。原始関数をもたない関数もある(第4章のダルブーの定理により、符号関数はどんな関数の導関数でもない)。
  • 「面積」とは何か。面積を定義しないまま積分を面積と呼ぶのは循環している。

そこで大学では、区間を細かく分けて長方形の面積の和で近似するという考え方で積分を直接定義し、「積分 = 原始関数の差」は定理(微分積分学の基本定理)として証明する。本章では、ダルブーによる上和・下和を使った定式化を採用する。

5.2 上和・下和とリーマン積分

以下、f ⁣:[a,b]→Rf\colon [a, b] \to \mathbb{R} は有界な関数とする。

定義 5.1 (1) 有限個の点 a=x0<x1<⋯<xn=ba = x_0 < x_1 < \cdots < x_n = b の集合 P={x0,…,xn}P = \lbrace x_0, \dots, x_n \rbrace を [a,b][a, b] の分割 (partition) という。Δxi=xi−xi−1\Delta x_i = x_i - x_{i-1} とし、∣P∣=max⁡iΔxi\lvert P \rvert = \max_i \Delta x_i を分割の幅 (mesh) という。分割 Q⊃PQ \supset P を PP の細分 (refinement) という。 (2) Mi=sup⁡[xi−1,xi]fM_i = \sup_{[x_{i-1}, x_i]} f、mi=inf⁡[xi−1,xi]fm_i = \inf_{[x_{i-1}, x_i]} f とおき、

U(f,P)=∑i=1nMiΔxi,L(f,P)=∑i=1nmiΔxiU(f, P) = \sum_{i=1}^{n} M_i \Delta x_i, \qquad L(f, P) = \sum_{i=1}^{n} m_i \Delta x_i

をそれぞれ上和 (upper sum)、下和 (lower sum) という。

上和はグラフを上から覆う長方形の面積の和、下和は下から埋める長方形の面積の和である。

補題 5.2 (1) QQ が PP の細分ならば L(f,P)≤L(f,Q)≤U(f,Q)≤U(f,P)L(f, P) \leq L(f, Q) \leq U(f, Q) \leq U(f, P)。 (2) 任意の分割 P,QP, Q について L(f,P)≤U(f,Q)L(f, P) \leq U(f, Q)。

証明. (1) PP に 1 点 y∈(xi−1,xi)y \in (x_{i-1}, x_i) を加えた場合を示せば、帰納法で一般の場合が従う。[xi−1,y][x_{i-1}, y] と [y,xi][y, x_i] 上の上限はどちらも MiM_i 以下なので、UU の第 ii 項 MiΔxiM_i \Delta x_i は sup⁡[xi−1,y]f⋅(y−xi−1)+sup⁡[y,xi]f⋅(xi−y)\sup_{[x_{i-1}, y]} f \cdot (y - x_{i-1}) + \sup_{[y, x_i]} f \cdot (x_i - y) 以下のものに置き換わる。よって UU は増えない。LL も同様。(2) P∪QP \cup Q は両方の細分なので L(f,P)≤L(f,P∪Q)≤U(f,P∪Q)≤U(f,Q)L(f, P) \leq L(f, P \cup Q) \leq U(f, P \cup Q) \leq U(f, Q)。□\square

補題 5.2 (2) より、上和全体は下に有界、下和全体は上に有界なので、次の上限・下限が存在する。

定義 5.3(リーマン積分, Riemann integral)

∫ab‾f=inf⁡PU(f,P),∫ab‾f=sup⁡PL(f,P)\overline{\int_a^b} f = \inf_P U(f, P), \qquad \underline{\int_a^b} f = \sup_P L(f, P)

をそれぞれ上積分・下積分という。補題 5.2 (2) より常に ∫ab‾f≤∫ab‾f\underline{\int_a^b} f \leq \overline{\int_a^b} f である。両者が等しいとき、ff は [a,b][a, b] でリーマン積分可能(または可積分, integrable)であるといい、その共通の値を ∫abf(x) dx\int_a^b f(x)\ dx と書く。

例 5.4 (1) 定数関数 f=cf = c は U(f,P)=L(f,P)=c(b−a)U(f, P) = L(f, P) = c(b - a) なので可積分で、∫abc dx=c(b−a)\int_a^b c\ dx = c(b-a)。 (2) ディリクレ関数 f=1Qf = \mathbf{1}_{\mathbb{Q}} は [0,1][0, 1] で可積分でない。どの小区間にも有理数と無理数がある(定理 1.23、系 1.24)ので、すべての分割で Mi=1M_i = 1、mi=0m_i = 0、したがって上積分は 11、下積分は 00 である。

実験リーマン和と上和・下和分割を細かくすると上和と下和が挟み撃ちで積分に近づく様子を観察します。

この実験は JavaScript を有効にすると動きます。

5.3 可積分性の条件

定理 5.5(リーマンの判定条件)有界関数 ff が [a,b][a, b] で可積分であるための必要十分条件は、任意の ε>0\varepsilon > 0 に対して U(f,P)−L(f,P)<εU(f, P) - L(f, P) < \varepsilon となる分割 PP が存在することである。

証明. (十分性) 任意の ε\varepsilon について 0≤∫‾f−∫‾f≤U(f,P)−L(f,P)<ε0 \leq \overline{\int} f - \underline{\int} f \leq U(f, P) - L(f, P) < \varepsilon なので、両者は等しい。(必要性) 共通の値を II とすると、下積分・上積分の定義より L(f,P1)>I−ε2L(f, P_1) > I - \frac{\varepsilon}{2}、U(f,P2)<I+ε2U(f, P_2) < I + \frac{\varepsilon}{2} となる分割 P1,P2P_1, P_2 がある。P=P1∪P2P = P_1 \cup P_2 とすれば補題 5.2 より U(f,P)−L(f,P)≤U(f,P2)−L(f,P1)<εU(f, P) - L(f, P) \leq U(f, P_2) - L(f, P_1) < \varepsilon。□\square

高校の区分求積法に対応するのが次の定理である。分割 PP と各小区間の点 ξi∈[xi−1,xi]\xi_i \in [x_{i-1}, x_i] に対して S(f,P,ξ)=∑i=1nf(ξi)ΔxiS(f, P, \xi) = \sum_{i=1}^{n} f(\xi_i)\Delta x_i をリーマン和 (Riemann sum) という。

定理 5.6 有界関数 ff について、次は同値である。

  1. ff は [a,b][a, b] で可積分で ∫abf=I\int_a^b f = I。
  2. 任意の ε>0\varepsilon > 0 に対して δ>0\delta > 0 が存在し、∣P∣<δ\lvert P \rvert < \delta を満たす任意の分割 PP と任意の点 ξi\xi_i の選び方について ∣S(f,P,ξ)−I∣<ε\lvert S(f, P, \xi) - I \rvert < \varepsilon。

証明. (1)⇒(2):∣f∣≤M\lvert f \rvert \leq M とする。ε>0\varepsilon > 0 に対し U(f,P0)−L(f,P0)<εU(f, P_0) - L(f, P_0) < \varepsilon となる分割 P0P_0 をとり、その内部の分点の個数を kk とする。∣P∣<δ\lvert P \rvert < \delta の分割 PP に P0P_0 の点を 1 個加えると、その点を含む PP の小区間は高々 1 個(分点と一致すれば 00 個)なので、UU は高々 2Mδ2M\delta しか減らない。よって U(f,P)≤U(f,P∪P0)+2Mkδ≤U(f,P0)+2MkδU(f, P) \leq U(f, P \cup P_0) + 2Mk\delta \leq U(f, P_0) + 2Mk\delta。同様に L(f,P)≥L(f,P0)−2MkδL(f, P) \geq L(f, P_0) - 2Mk\delta。したがって δ=ε4Mk+1\delta = \frac{\varepsilon}{4Mk + 1} とすれば U(f,P)−L(f,P)<2εU(f, P) - L(f, P) < 2\varepsilon。S(f,P,ξ)S(f, P, \xi) と II はともに [L(f,P),U(f,P)][L(f, P), U(f, P)] に属するので、∣S(f,P,ξ)−I∣<2ε\lvert S(f, P, \xi) - I \rvert < 2\varepsilon。

(2)⇒(1):ε>0\varepsilon > 0 に対する δ\delta をとり、∣P∣<δ\lvert P \rvert < \delta の分割 PP を一つ固定する。各 ii で f(ξi)>Mi−εb−af(\xi_i) > M_i - \frac{\varepsilon}{b - a} となる ξi\xi_i を選ぶと U(f,P)<S(f,P,ξ)+ε<I+2εU(f, P) < S(f, P, \xi) + \varepsilon < I + 2\varepsilon。同様に L(f,P)>I−2εL(f, P) > I - 2\varepsilon。よって U−L<4εU - L < 4\varepsilon であり、定理 5.5 より可積分で、積分値は II である。□\square

特に ff が可積分なら、等分割 xi=a+i(b−a)nx_i = a + \frac{i(b-a)}{n} に対して

lim⁡n→∞b−an∑i=1nf(a+i(b−a)n)=∫abf(x) dx\lim_{n \to \infty} \frac{b - a}{n}\sum_{i=1}^{n} f\left(a + \frac{i(b-a)}{n}\right) = \int_a^b f(x)\,dx

が成り立つ(区分求積法)。

5.4 可積分な関数

定理 5.7 [a,b][a, b] 上の連続関数は可積分である。

証明. ff はハイネ–カントールの定理(定理 3.25)により一様連続である。ε>0\varepsilon > 0 に対し、∣x−y∣<δ⇒∣f(x)−f(y)∣<εb−a\lvert x - y \rvert < \delta \Rightarrow \lvert f(x) - f(y) \rvert < \frac{\varepsilon}{b - a} となる δ\delta をとり、∣P∣<δ\lvert P \rvert < \delta の分割 PP をとる。各小区間で ff は最大値・最小値をとる(定理 3.21)ので Mi−mi=f(si)−f(ti)<εb−aM_i - m_i = f(s_i) - f(t_i) < \frac{\varepsilon}{b - a}(si,ti∈[xi−1,xi]s_i, t_i \in [x_{i-1}, x_i])。よって U(f,P)−L(f,P)<εb−a∑Δxi=εU(f, P) - L(f, P) < \frac{\varepsilon}{b-a}\sum \Delta x_i = \varepsilon。□\square

定理 5.8 [a,b][a, b] 上の単調関数は可積分である。

証明. ff を単調増加とし、nn 等分割 PnP_n をとると Mi=f(xi)M_i = f(x_i)、mi=f(xi−1)m_i = f(x_{i-1}) なので

U(f,Pn)−L(f,Pn)=∑i=1n(f(xi)−f(xi−1))b−an=(f(b)−f(a))(b−a)nU(f, P_n) - L(f, P_n) = \sum_{i=1}^{n} \bigl(f(x_i) - f(x_{i-1})\bigr)\frac{b - a}{n} = \frac{(f(b) - f(a))(b - a)}{n}

n→∞n \to \infty で 00 に収束するので、定理 5.5 より可積分。□\square

単調関数は不連続点を可算無限個もちうる(問題 3.10)が、それでも可積分である。可積分性は連続性よりずっと弱い条件である。

例 5.9 トマエ関数(例 3.17 (3))は [0,1][0, 1] で可積分で、積分は 00 である。すべての小区間に無理数があるので L(f,P)=0L(f, P) = 0。ε>0\varepsilon > 0 に対し、f(x)≥ε2f(x) \geq \frac{\varepsilon}{2} となる x∈[0,1]x \in [0,1](分母が 2ε\frac{2}{\varepsilon} 以下の有理数)は有限個で、その個数を NN とする。∣P∣<δ\lvert P \rvert < \delta の分割では、それらの点を含む小区間は高々 2N2N 個でその上で Mi≤1M_i \leq 1、残りの小区間では Mi<ε2M_i < \frac{\varepsilon}{2} なので、U(f,P)≤2Nδ+ε2U(f, P) \leq 2N\delta + \frac{\varepsilon}{2}。δ=ε4N\delta = \frac{\varepsilon}{4N} とすれば U(f,P)≤εU(f, P) \leq \varepsilon。この関数は有理数の点(稠密な集合)すべてで不連続である。

補足

有界関数がリーマン積分可能であるための必要十分条件は、その不連続点全体が「測度 00」(任意の ε>0\varepsilon > 0 に対して、長さの総和が ε\varepsilon 未満の可算個の開区間で覆える)であることが知られている(ルベーグの判定条件)。証明は 06-measure-integration 第3章 の内容と合わせて学ぶとよい。

5.5 積分の性質

定理 5.10 f,gf, g を [a,b][a, b] 上の可積分関数、c∈Rc \in \mathbb{R} とする。

  1. (線形性)f+gf + g、cfcf は可積分で、∫(f+g)=∫f+∫g\int (f + g) = \int f + \int g、∫cf=c∫f\int cf = c\int f。
  2. (単調性)f≤gf \leq g ならば ∫f≤∫g\int f \leq \int g。
  3. ∣f∣\lvert f \rvert は可積分で、∣∫abf∣≤∫ab∣f∣\left\lvert \int_a^b f \right\rvert \leq \int_a^b \lvert f \rvert。
  4. 積 fgfg は可積分である。
  5. (区間についての加法性)a<c<ba < c < b のとき、ff が [a,b][a, b] で可積分であることと [a,c][a, c]、[c,b][c, b] の両方で可積分であることは同値で、∫abf=∫acf+∫cbf\int_a^b f = \int_a^c f + \int_c^b f。

証明. (1) 各小区間で sup⁡(f+g)≤sup⁡f+sup⁡g\sup(f + g) \leq \sup f + \sup g、inf⁡(f+g)≥inf⁡f+inf⁡g\inf(f + g) \geq \inf f + \inf g なので L(f,P)+L(g,P)≤L(f+g,P)≤U(f+g,P)≤U(f,P)+U(g,P)L(f, P) + L(g, P) \leq L(f + g, P) \leq U(f + g, P) \leq U(f, P) + U(g, P)。U(f,P)−L(f,P)U(f, P) - L(f, P) と U(g,P)−L(g,P)U(g, P) - L(g, P) を同時に小さくする分割(共通細分)をとれば、定理 5.5 より f+gf + g は可積分で、積分値は ∫f+∫g\int f + \int g と ε\varepsilon 以内の差しかない。cfcf は c≥0c \geq 0 なら U(cf,P)=cU(f,P)U(cf, P) = cU(f, P) などから、c<0c < 0 なら U(cf,P)=cL(f,P)U(cf, P) = cL(f, P) から従う。 (2) g−f≥0g - f \geq 0 なので L(g−f,P)≥0L(g - f, P) \geq 0、よって ∫(g−f)≥0\int (g - f) \geq 0 で (1) を使う。 (3) 各小区間で sup⁡∣f∣−inf⁡∣f∣≤sup⁡f−inf⁡f\sup\lvert f \rvert - \inf \lvert f \rvert \leq \sup f - \inf f(∣∣f(s)∣−∣f(t)∣∣≤∣f(s)−f(t)∣\bigl\lvert \lvert f(s) \rvert - \lvert f(t) \rvert \bigr\rvert \leq \lvert f(s) - f(t) \rvert による)なので ∣f∣\lvert f \rvert も定理 5.5 の条件を満たす。不等式は −∣f∣≤f≤∣f∣-\lvert f \rvert \leq f \leq \lvert f \rvert と (2) から。 (4) ∣f∣≤M\lvert f \rvert \leq M なら各小区間で sup⁡f2−inf⁡f2≤2M(sup⁡∣f∣−inf⁡∣f∣)\sup f^2 - \inf f^2 \leq 2M(\sup\lvert f \rvert - \inf\lvert f \rvert) なので、(3) より f2f^2 は可積分。fg=14((f+g)2−(f−g)2)fg = \frac{1}{4}\bigl((f+g)^2 - (f-g)^2\bigr) と (1) による。 (5) [a,b][a, b] の分割に点 cc を加えても U−LU - L は増えないので、cc を含む分割だけを考えてよい。そのような分割は [a,c][a, c] と [c,b][c, b] の分割の組と対応し、上和・下和は和に分かれる。これと定理 5.5 から従う。□\square

a>ba > b のとき ∫abf=−∫baf\int_a^b f = -\int_b^a f、また ∫aaf=0\int_a^a f = 0 と約束すると、(5) の等式は a,b,ca, b, c の大小によらず成り立つ。

定理 5.11(積分の平均値定理)ff が [a,b][a, b] で連続、gg が可積分で g≥0g \geq 0 ならば、

∫abf(x)g(x) dx=f(c)∫abg(x) dx\int_a^b f(x)g(x)\,dx = f(c)\int_a^b g(x)\,dx

となる c∈[a,b]c \in [a, b] が存在する。特に g=1g = 1 として ∫abf=f(c)(b−a)\int_a^b f = f(c)(b - a)。

証明. m=min⁡fm = \min f、M=max⁡fM = \max f とすると mg≤fg≤Mgmg \leq fg \leq Mg なので m∫g≤∫fg≤M∫gm\int g \leq \int fg \leq M\int g。∫g=0\int g = 0 なら ∫fg=0\int fg = 0 で任意の cc でよい。∫g>0\int g > 0 なら ∫fg∫g∈[m,M]\frac{\int fg}{\int g} \in [m, M] であり、中間値の定理より f(c)f(c) がこの値に等しい cc がある。□\square

5.6 微分積分学の基本定理

定理 5.12(微分積分学の基本定理, fundamental theorem of calculus)ff を [a,b][a, b] で可積分とし、F(x)=∫axf(t) dtF(x) = \int_a^x f(t)\ dt とおく。

  1. FF は [a,b][a, b] でリプシッツ連続である。
  2. ff が x0∈[a,b]x_0 \in [a, b] で連続ならば、FF は x0x_0 で微分可能で F′(x0)=f(x0)F'(x_0) = f(x_0)。

証明. (1) ∣f∣≤M\lvert f \rvert \leq M なら ∣F(x)−F(y)∣=∣∫yxf∣≤M∣x−y∣\lvert F(x) - F(y) \rvert = \left\lvert \int_y^x f \right\rvert \leq M\lvert x - y \rvert。 (2) ε>0\varepsilon > 0 に対し、∣t−x0∣<δ⇒∣f(t)−f(x0)∣<ε\lvert t - x_0 \rvert < \delta \Rightarrow \lvert f(t) - f(x_0) \rvert < \varepsilon となる δ\delta をとる。0<∣h∣<δ0 < \lvert h \rvert < \delta(x0+h∈[a,b]x_0 + h \in [a, b])なら

∣F(x0+h)−F(x0)h−f(x0)∣=∣1h∫x0x0+h(f(t)−f(x0))dt∣≤ε\left\lvert \frac{F(x_0 + h) - F(x_0)}{h} - f(x_0) \right\rvert = \left\lvert \frac{1}{h}\int_{x_0}^{x_0 + h}\bigl(f(t) - f(x_0)\bigr)dt \right\rvert \leq \varepsilon

である。□\square

定理 5.13(微分積分学の基本定理・その 2)FF が [a,b][a, b] で微分可能(端点では片側微分)で、F′=fF' = f が [a,b][a, b] で可積分ならば

∫abf(x) dx=F(b)−F(a)\int_a^b f(x)\,dx = F(b) - F(a)

証明. 任意の分割 PP について、各小区間に平均値の定理を使うと

F(b)−F(a)=∑i=1n(F(xi)−F(xi−1))=∑i=1nf(ξi)Δxi(ξi∈(xi−1,xi))F(b) - F(a) = \sum_{i=1}^{n}\bigl(F(x_i) - F(x_{i-1})\bigr) = \sum_{i=1}^{n} f(\xi_i)\Delta x_i \quad (\xi_i \in (x_{i-1}, x_i))

右辺は L(f,P)L(f, P) と U(f,P)U(f, P) の間にある。PP は任意なので ∫‾f≤F(b)−F(a)≤∫‾f\underline{\int} f \leq F(b) - F(a) \leq \overline{\int} f であり、ff は可積分なので両端は ∫abf\int_a^b f に等しい。□\square

ff が連続なら、定理 5.12 より F(x)=∫axfF(x) = \int_a^x f は ff の原始関数であり、任意の原始関数 GG は FF と定数の差しかない(系 4.14)ので、∫abf=G(b)−G(a)\int_a^b f = G(b) - G(a)。こうして高校の「定義」は定理として回復される。以後、[G]ab=G(b)−G(a)[G]_a^b = G(b) - G(a) と書く。

注意

(1) 可積分関数は原始関数をもつとは限らない。[−1,1][-1, 1] 上の符号関数は(単調なので)可積分だが、ダルブーの定理により原始関数をもたない。 (2) FF が至るところ微分可能でも、F′F' が可積分とは限らない。F(x)=x2sin⁡1x2F(x) = x^2\sin\frac{1}{x^2}(x≠0x \neq 0)、F(0)=0F(0) = 0 は [0,1][0, 1] で微分可能だが、F′(x)=2xsin⁡1x2−2xcos⁡1x2F'(x) = 2x\sin\frac{1}{x^2} - \frac{2}{x}\cos\frac{1}{x^2} は有界でないのでリーマン積分可能でない。定理 5.13 で「F′F' が可積分」を仮定したのはこのためである。

5.7 置換積分と部分積分

定理 5.14 (1)(部分積分, integration by parts)f,gf, g が [a,b][a, b] で C1C^1 級ならば

∫abf(x)g′(x) dx=[f(x)g(x)]ab−∫abf′(x)g(x) dx\int_a^b f(x)g'(x)\,dx = \bigl[f(x)g(x)\bigr]_a^b - \int_a^b f'(x)g(x)\,dx

(2)(置換積分, substitution)φ ⁣:[α,β]→R\varphi\colon [\alpha, \beta] \to \mathbb{R} が C1C^1 級で、ff が φ([α,β])\varphi([\alpha, \beta]) を含む区間で連続ならば

∫φ(α)φ(β)f(x) dx=∫αβf(φ(t))φ′(t) dt\int_{\varphi(\alpha)}^{\varphi(\beta)} f(x)\,dx = \int_\alpha^\beta f(\varphi(t))\varphi'(t)\,dt

証明. (1) (fg)′=f′g+fg′(fg)' = f'g + fg' は連続なので、定理 5.13 より ∫ab(f′g+fg′)=[fg]ab\int_a^b (f'g + fg') = [fg]_a^b。(2) FF を ff の原始関数とすると、合成関数の微分より (F∘φ)′=f(φ)φ′(F \circ \varphi)' = f(\varphi)\varphi' は連続なので、右辺は F(φ(β))−F(φ(α))F(\varphi(\beta)) - F(\varphi(\alpha)) に等しく、これは左辺に等しい。□\square

置換積分で φ\varphi の単調性は不要である。多変数の変数変換(第9章)では単射性が必要になるので対照的である。

例 5.15 (1) ∫1elog⁡x dx=[xlog⁡x]1e−∫1ex⋅1x dx=e−(e−1)=1\int_1^e \log x\ dx = [x\log x]_1^e - \int_1^e x \cdot \frac{1}{x}\ dx = e - (e - 1) = 1。 (2) ∫01xex2dx\int_0^1 xe^{x^2}dx:φ(t)=t2\varphi(t) = t^2 として 12∫01eu du=e−12\frac{1}{2}\int_0^1 e^{u}\ du = \frac{e - 1}{2}。 (3) (三角関数を認めて)∫−111−x2 dx\int_{-1}^{1}\sqrt{1 - x^2}\ dx:x=sin⁡tx = \sin t(t∈[−π2,π2]t \in [-\frac{\pi}{2}, \frac{\pi}{2}])として ∫−π/2π/2cos⁡2t dt=π2\int_{-\pi/2}^{\pi/2}\cos^2 t\ dt = \frac{\pi}{2}。単位円の面積の半分である。

例 5.16(ウォリス積分)In=∫0π/2sin⁡nx dxI_n = \int_0^{\pi/2}\sin^n x\ dx とおく。I0=π2I_0 = \frac{\pi}{2}、I1=1I_1 = 1 であり、n≥2n \geq 2 で部分積分により

In=[−sin⁡n−1xcos⁡x]0π/2+(n−1)∫0π/2sin⁡n−2xcos⁡2x dx=(n−1)(In−2−In)I_n = \bigl[-\sin^{n-1}x\cos x\bigr]_0^{\pi/2} + (n - 1)\int_0^{\pi/2}\sin^{n-2}x\cos^2 x\,dx = (n-1)(I_{n-2} - I_n)

よって In=n−1nIn−2I_n = \frac{n-1}{n}I_{n-2}。これから

I2m=(2m−1)(2m−3)⋯12m(2m−2)⋯2⋅π2,I2m+1=2m(2m−2)⋯2(2m+1)(2m−1)⋯3=22m(m!)2(2m+1)!I_{2m} = \frac{(2m-1)(2m-3)\cdots 1}{2m(2m-2)\cdots 2}\cdot\frac{\pi}{2}, \qquad I_{2m+1} = \frac{2m(2m-2)\cdots 2}{(2m+1)(2m-1)\cdots 3} = \frac{2^{2m}(m!)^2}{(2m+1)!}

5.8 積分剰余のテイラーの定理

定理 5.17(テイラーの定理・積分剰余)ff が区間 II で Cn+1C^{n+1} 級で、a,x∈Ia, x \in I ならば

f(x)=∑k=0nf(k)(a)k!(x−a)k+1n!∫ax(x−t)nf(n+1)(t) dtf(x) = \sum_{k=0}^{n}\frac{f^{(k)}(a)}{k!}(x - a)^k + \frac{1}{n!}\int_a^x (x - t)^n f^{(n+1)}(t)\,dt

証明. nn についての帰納法。n=0n = 0 は f(x)=f(a)+∫axf′(t) dtf(x) = f(a) + \int_a^x f'(t)\ dt(定理 5.13)。n−1n - 1 で成り立つとし、その剰余項を u(t)=f(n)(t)u(t) = f^{(n)}(t)、v(t)=−(x−t)nn!v(t) = -\frac{(x - t)^n}{n!}(v′(t)=(x−t)n−1(n−1)!v'(t) = \frac{(x-t)^{n-1}}{(n-1)!})として部分積分すると

1(n−1)!∫ax(x−t)n−1f(n)(t) dt=[−(x−t)nn!f(n)(t)]ax+1n!∫ax(x−t)nf(n+1)(t) dt\frac{1}{(n-1)!}\int_a^x (x - t)^{n-1}f^{(n)}(t)\,dt = \left[-\frac{(x-t)^n}{n!}f^{(n)}(t)\right]_a^x + \frac{1}{n!}\int_a^x (x-t)^n f^{(n+1)}(t)\,dt

第 1 項は f(n)(a)n!(x−a)n\frac{f^{(n)}(a)}{n!}(x - a)^n なので、nn の場合が得られる。□\square

x>ax > a のとき (x−t)n≥0(x - t)^n \geq 0 なので、積分の平均値定理(定理 5.11)より剰余項は f(n+1)(c)n!∫ax(x−t)ndt=f(n+1)(c)(n+1)!(x−a)n+1\frac{f^{(n+1)}(c)}{n!}\int_a^x (x-t)^n dt = \frac{f^{(n+1)}(c)}{(n+1)!}(x - a)^{n+1} となり、ラグランジュ剰余が再び得られる(x<ax < a も同様)。積分剰余はラグランジュ剰余より精密な評価ができることがある。

例 5.18 −1<x≤1-1 < x \leq 1 で log⁡(1+x)=∑k=1∞(−1)k−1kxk\log(1 + x) = \sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}x^k。0≤x≤10 \leq x \leq 1 は例 4.22 (2) で示した。−1<x<0-1 < x < 0 のとき、f(n+1)(t)=(−1)nn!(1+t)n+1f^{(n+1)}(t) = \frac{(-1)^n n!}{(1+t)^{n+1}} より剰余項の絶対値は ∫x0(t−x)n(1+t)n+1dt\int_x^0 \frac{(t - x)^n}{(1 + t)^{n+1}}dt。x≤t≤0x \leq t \leq 0 では t(1+x)≤0t(1 + x) \leq 0 より t−x≤∣x∣(1+t)t - x \leq \lvert x \rvert(1 + t) なので、これは

∫x0∣x∣n1+tdt=∣x∣n(−log⁡(1+x))→0(n→∞)\int_x^0 \frac{\lvert x \rvert^n}{1 + t}dt = \lvert x \rvert^n\bigl(-\log(1 + x)\bigr) \to 0 \quad (n \to \infty)

以下である(ラグランジュ剰余では xx が −1-1 に近いとき評価がうまくいかない)。

5.9 広義積分

区間が無限であったり、関数が非有界であったりすると、リーマン積分は定義できない。そのような場合は、有界閉区間上の積分の極限として積分を定める。

定義 5.19(広義積分, improper integral)−∞<a<b≤∞-\infty < a < b \leq \infty とし、ff は [a,b)[a, b) で定義され、任意の c∈(a,b)c \in (a, b) について [a,c][a, c] で可積分とする。極限

∫abf(x) dx:=lim⁡c→b−0∫acf(x) dx\int_a^b f(x)\,dx := \lim_{c \to b - 0}\int_a^c f(x)\,dx

(b=∞b = \infty なら c→∞c \to \infty)が存在するとき、広義積分は収束するという。左端が問題になる場合も同様に定め、両端が問題になる場合は内部の点で分けて、両方が収束するとき収束するという。∫ab∣f∣\int_a^b \lvert f \rvert が収束するとき、絶対収束するという。

例 5.20 (1) ∫01x−p dx\int_0^1 x^{-p}\ dx は p<1p < 1 で収束して値は 11−p\frac{1}{1-p}、p≥1p \geq 1 で発散する。p≠1p \neq 1 なら ∫c1x−pdx=1−c1−p1−p\int_c^1 x^{-p}dx = \frac{1 - c^{1-p}}{1 - p} で、c→+0c \to +0 のとき c1−pc^{1-p} は p<1p < 1 なら 00 に、p>1p > 1 なら ∞\infty に近づく。p=1p = 1 なら −log⁡c→∞-\log c \to \infty。 (2) ∫1∞x−p dx\int_1^\infty x^{-p}\ dx は p>1p > 1 で収束して値は 1p−1\frac{1}{p - 1}、p≤1p \leq 1 で発散する。 (3) ∫0∞e−xdx=lim⁡c→∞(1−e−c)=1\int_0^\infty e^{-x}dx = \lim_{c \to \infty}(1 - e^{-c}) = 1。

定理 5.21 (1)(コーシーの判定法)∫abf\int_a^b f(bb で広義)が収束するための必要十分条件は、任意の ε>0\varepsilon > 0 に対して c0<bc_0 < b が存在し、c0<c<c′<bc_0 < c < c' < b ならば ∣∫cc′f∣<ε\left\lvert \int_c^{c'} f \right\rvert < \varepsilon となることである。 (2) 絶対収束する広義積分は収束する。 (3)(比較判定法)[a,b)[a, b) で ∣f∣≤g\lvert f \rvert \leq g であり ∫abg\int_a^b g が収束すれば、∫abf\int_a^b f は絶対収束する。

証明. (1) Φ(c)=∫acf\Phi(c) = \int_a^c f に命題 3.12 を適用すればよい(Φ(c′)−Φ(c)=∫cc′f\Phi(c') - \Phi(c) = \int_c^{c'} f)。(2) ∣∫cc′f∣≤∫cc′∣f∣\left\lvert \int_c^{c'} f \right\rvert \leq \int_c^{c'}\lvert f \rvert と (1)。(3) ∫cc′∣f∣≤∫cc′g\int_c^{c'}\lvert f \rvert \leq \int_c^{c'} g と (1)。□\square

実用上は、例 5.20 と比較して次のように判定することが多い:x→∞x \to \infty で f(x)=O(x−p)f(x) = O(x^{-p})(p>1p > 1)なら ∫∞f\int^\infty f は絶対収束し、x→+0x \to +0 で f(x)=O(x−p)f(x) = O(x^{-p})(p<1p < 1)なら ∫0f\int_0 f は絶対収束する。

例 5.22(条件収束する広義積分)∫0∞sin⁡xx dx\int_0^\infty \frac{\sin x}{x}\ dx は収束するが、絶対収束しない。x→0x \to 0 で sin⁡xx→1\frac{\sin x}{x} \to 1 なので [0,1][0, 1] では問題ない(00 での値を 11 とすれば連続)。c>1c > 1 で部分積分すると

∫1csin⁡xx dx=[−cos⁡xx]1c−∫1ccos⁡xx2 dx\int_1^c \frac{\sin x}{x}\,dx = \left[-\frac{\cos x}{x}\right]_1^c - \int_1^c \frac{\cos x}{x^2}\,dx

で、cos⁡cc→0\frac{\cos c}{c} \to 0、∣cos⁡xx2∣≤1x2\left\lvert \frac{\cos x}{x^2} \right\rvert \leq \frac{1}{x^2} より第 2 項は絶対収束するので、c→∞c \to \infty で収束する。一方 ∫kπ(k+1)π∣sin⁡x∣xdx≥1(k+1)π∫kπ(k+1)π∣sin⁡x∣dx=2(k+1)π\int_{k\pi}^{(k+1)\pi}\frac{\lvert \sin x \rvert}{x}dx \geq \frac{1}{(k+1)\pi}\int_{k\pi}^{(k+1)\pi}\lvert \sin x \rvert dx = \frac{2}{(k+1)\pi} で、∑1k+1\sum \frac{1}{k+1} は発散するので絶対収束しない。値が π2\frac{\pi}{2} であることは 05-complex-analysis 第5章 の留数定理などで示される。

5.10 ガンマ関数とベータ関数

定義 5.23 s>0s > 0、p>0p > 0、q>0q > 0 に対して

Γ(s)=∫0∞xs−1e−x dx,B(p,q)=∫01xp−1(1−x)q−1 dx\Gamma(s) = \int_0^\infty x^{s-1}e^{-x}\,dx, \qquad B(p, q) = \int_0^1 x^{p-1}(1 - x)^{q-1}\,dx

と定め、Γ\Gamma をガンマ関数 (gamma function)、BB をベータ関数 (beta function) という。

収束の確認:Γ(s)\Gamma(s) は x→+0x \to +0 で xs−1e−x≤xs−1x^{s-1}e^{-x} \leq x^{s-1}(s−1>−1s - 1 > -1)、x→∞x \to \infty で xs+1e−x→0x^{s+1}e^{-x} \to 0(問題 3.8)より xs−1e−x=O(x−2)x^{s-1}e^{-x} = O(x^{-2}) なので、両端で絶対収束する。B(p,q)B(p,q) は x→+0x \to +0 で O(xp−1)O(x^{p-1})、x→1−0x \to 1 - 0 で O((1−x)q−1)O((1-x)^{q-1}) なので収束する。

定理 5.24 (1) Γ(s+1)=sΓ(s)\Gamma(s + 1) = s\Gamma(s)(s>0s > 0)、Γ(1)=1\Gamma(1) = 1。特に Γ(n+1)=n!\Gamma(n + 1) = n!(n∈Z≥0n \in \mathbb{Z}_{\geq 0})。 (2) B(p,q)=B(q,p)B(p, q) = B(q, p)、B(p,q)=2∫0π/2sin⁡2p−1θcos⁡2q−1θ dθB(p, q) = 2\int_0^{\pi/2}\sin^{2p-1}\theta\cos^{2q-1}\theta\ d\theta。 (3) B(p+1,q)=pp+qB(p,q)B(p + 1, q) = \frac{p}{p + q}B(p, q)。 (4) 正の整数 m,nm, n について B(m,n)=(m−1)!(n−1)!(m+n−1)!=Γ(m)Γ(n)Γ(m+n)B(m, n) = \frac{(m-1)!(n-1)!}{(m+n-1)!} = \frac{\Gamma(m)\Gamma(n)}{\Gamma(m+n)}。

証明. 広義積分での部分積分・置換積分は、有界閉区間 [ε,R][\varepsilon, R] や [ε,1−ε][\varepsilon, 1 - \varepsilon] で行ってから極限をとればよい。

(1) ∫εRxse−xdx=[−xse−x]εR+s∫εRxs−1e−xdx\int_\varepsilon^R x^s e^{-x}dx = \bigl[-x^s e^{-x}\bigr]_\varepsilon^R + s\int_\varepsilon^R x^{s-1}e^{-x}dx で、ε→0\varepsilon \to 0、R→∞R \to \infty とすると境界項は 00 に収束する。Γ(1)=∫0∞e−xdx=1\Gamma(1) = \int_0^\infty e^{-x}dx = 1。

(2) 前半は x=1−yx = 1 - y、後半は x=sin⁡2θx = \sin^2\theta(dx=2sin⁡θcos⁡θ dθdx = 2\sin\theta\cos\theta\ d\theta)と置換する。

(3) 部分積分で B(p,q+1)=[xpp(1−x)q]01+qp∫01xp(1−x)q−1dx=qpB(p+1,q)B(p, q + 1) = \left[\frac{x^p}{p}(1-x)^q\right]_0^1 + \frac{q}{p}\int_0^1 x^p(1-x)^{q-1}dx = \frac{q}{p}B(p + 1, q)。また xp−1(1−x)q−1=xp−1(1−x)q−1(x+(1−x))x^{p-1}(1-x)^{q-1} = x^{p-1}(1-x)^{q-1}\bigl(x + (1 - x)\bigr) より B(p,q)=B(p+1,q)+B(p,q+1)=(1+qp)B(p+1,q)B(p, q) = B(p + 1, q) + B(p, q + 1) = \left(1 + \frac{q}{p}\right)B(p+1, q)。

(4) B(1,n)=∫01(1−x)n−1dx=1nB(1, n) = \int_0^1(1-x)^{n-1}dx = \frac{1}{n} と (3) から、mm についての帰納法で B(m,n)=m−1m+n−1B(m−1,n)=(m−1)!(n−1)!(m+n−1)!B(m, n) = \frac{m-1}{m+n-1}B(m-1, n) = \frac{(m-1)!(n-1)!}{(m+n-1)!}。□\square

(4) の関係式は、実はすべての p,q>0p, q > 0 で成り立つ:

B(p,q)=Γ(p)Γ(q)Γ(p+q)B(p, q) = \frac{\Gamma(p)\Gamma(q)}{\Gamma(p + q)}

この証明には重積分の変数変換を使うので、第9章で行う。そこでは Γ(12)=π\Gamma\left(\frac{1}{2}\right) = \sqrt{\pi} も示す(本章の問題 5.10 でも別の方法で示す)。

5.11 ウォリスの公式とスターリングの公式

n!n! は nn とともに急激に増大する。その増え方を正確にとらえるのがスターリングの公式であり、その定数 2π\sqrt{2\pi} を決めるのにウォリス積分を使う。

命題 5.25(ウォリスの公式)例 5.16 の InI_n について nInIn−1=π2nI_nI_{n-1} = \frac{\pi}{2} であり、n In→π2\sqrt{n}\ I_n \to \sqrt{\frac{\pi}{2}}。また

lim⁡n→∞22n(n!)2(2n)!n=π\lim_{n \to \infty}\frac{2^{2n}(n!)^2}{(2n)!\sqrt{n}} = \sqrt{\pi}

証明. In=n−1nIn−2I_n = \frac{n-1}{n}I_{n-2} より nInIn−1=(n−1)In−1In−2=⋯=1⋅I1I0=π2nI_nI_{n-1} = (n-1)I_{n-1}I_{n-2} = \cdots = 1 \cdot I_1 I_0 = \frac{\pi}{2}。0≤sin⁡x≤10 \leq \sin x \leq 1 より (In)(I_n) は正の単調減少列なので n−1n=InIn−2≤InIn−1≤1\frac{n-1}{n} = \frac{I_n}{I_{n-2}} \leq \frac{I_n}{I_{n-1}} \leq 1、よって InIn−1→1\frac{I_n}{I_{n-1}} \to 1。nIn2=π2⋅InIn−1→π2nI_n^2 = \frac{\pi}{2}\cdot\frac{I_n}{I_{n-1}} \to \frac{\pi}{2}。最後の式は、例 5.16 の I2n+1=22n(n!)2(2n+1)!I_{2n+1} = \frac{2^{2n}(n!)^2}{(2n+1)!} より

22n(n!)2(2n)!n=(2n+1)I2n+1n=2n+1n2n+1⋅2n+1 I2n+1→2⋅π2=π\frac{2^{2n}(n!)^2}{(2n)!\sqrt{n}} = \frac{(2n+1)I_{2n+1}}{\sqrt{n}} = \frac{2n+1}{\sqrt{n}\sqrt{2n+1}}\cdot\sqrt{2n+1}\,I_{2n+1} \to \sqrt{2}\cdot\sqrt{\frac{\pi}{2}} = \sqrt{\pi}

である。□\square

定理 5.26(スターリングの公式, Stirling's formula)

lim⁡n→∞n!2πn (n/e)n=1\lim_{n \to \infty}\frac{n!}{\sqrt{2\pi n}\,(n/e)^n} = 1

さらに 1<n!2πn(n/e)n<e1/(12n)1 < \frac{n!}{\sqrt{2\pi n}(n/e)^n} < e^{1/(12n)} が成り立つ。

証明. an=n! ennn+1/2a_n = \frac{n!\ e^n}{n^{n + 1/2}} とおく。直接計算すると

log⁡anan+1=(n+12)log⁡(1+1n)−1\log\frac{a_n}{a_{n+1}} = \left(n + \frac{1}{2}\right)\log\left(1 + \frac{1}{n}\right) - 1

t=12n+1t = \frac{1}{2n + 1} とおくと 1+1n=1+t1−t1 + \frac{1}{n} = \frac{1 + t}{1 - t}、n+12=12tn + \frac{1}{2} = \frac{1}{2t} であり、例 5.18 の展開 log⁡1+t1−t=2(t+t33+t55+⋯ )\log\frac{1+t}{1-t} = 2\left(t + \frac{t^3}{3} + \frac{t^5}{5} + \cdots\right)(0<t<10 < t < 1)を使うと

log⁡anan+1=t23+t45+t67+⋯\log\frac{a_n}{a_{n+1}} = \frac{t^2}{3} + \frac{t^4}{5} + \frac{t^6}{7} + \cdots

これは正であり、t23(1+t2+t4+⋯ )=t23(1−t2)=112n(n+1)=112n−112(n+1)\frac{t^2}{3}(1 + t^2 + t^4 + \cdots) = \frac{t^2}{3(1 - t^2)} = \frac{1}{12n(n+1)} = \frac{1}{12n} - \frac{1}{12(n+1)} 以下である。よって (an)(a_n) は単調減少、(ane−1/(12n))\bigl(a_n e^{-1/(12n)}\bigr) は単調増加であり、an≥a1e−1/12>0a_n \geq a_1 e^{-1/12} > 0。単調収束定理より an→C>0a_n \to C > 0 で、C<an<Ce1/(12n)C < a_n < Ce^{1/(12n)}。

CC を決める。n!=annn+1/2e−nn! = a_n n^{n+1/2}e^{-n} をウォリスの公式に代入すると

22n(n!)2(2n)!n=22nan2n2n+1e−2na2n(2n)2n+1/2e−2nn=an22 a2n→C2\frac{2^{2n}(n!)^2}{(2n)!\sqrt{n}} = \frac{2^{2n}a_n^2 n^{2n+1}e^{-2n}}{a_{2n}(2n)^{2n+1/2}e^{-2n}\sqrt{n}} = \frac{a_n^2}{\sqrt{2}\,a_{2n}} \to \frac{C}{\sqrt{2}}

命題 5.25 よりこれは π\sqrt{\pi} なので C=2πC = \sqrt{2\pi}。□\square

数値例:10!=362880010! = 3628800 に対して 20π(10/e)10=3598695.6…\sqrt{20\pi}(10/e)^{10} = 3598695.6\ldots で、比は 1.00837…1.00837\ldots(e1/120=1.00837…e^{1/120} = 1.00837\ldots 未満)である。

5.12 級数の積分判定法

定理 5.27(積分判定法, integral test)f ⁣:[1,∞)→[0,∞)f\colon [1, \infty) \to [0, \infty) が単調減少ならば、∑n=1∞f(n)\sum_{n=1}^{\infty}f(n) と ∫1∞f(x) dx\int_1^\infty f(x)\ dx は収束・発散をともにする。さらに dn=∑k=1nf(k)−∫1nf(x) dxd_n = \sum_{k=1}^{n}f(k) - \int_1^n f(x)\ dx は単調減少で 0≤dn≤f(1)0 \leq d_n \leq f(1) を満たし、収束する。

証明. ff は単調なので各有界閉区間で可積分(定理 5.8)。k≤x≤k+1k \leq x \leq k + 1 で f(k+1)≤f(x)≤f(k)f(k+1) \leq f(x) \leq f(k) なので

f(k+1)≤∫kk+1f(x) dx≤f(k)f(k + 1) \leq \int_k^{k+1}f(x)\,dx \leq f(k)

k=1,…,n−1k = 1, \dots, n - 1 で足すと ∑k=2nf(k)≤∫1nf≤∑k=1n−1f(k)\sum_{k=2}^{n}f(k) \leq \int_1^n f \leq \sum_{k=1}^{n-1}f(k)。f≥0f \geq 0 なので部分和も ∫1cf\int_1^c f も cc について単調増加であり、一方が有界なら他方も有界である。また dn−dn+1=∫nn+1f−f(n+1)≥0d_n - d_{n+1} = \int_n^{n+1}f - f(n+1) \geq 0、dn=f(n)+∑k=1n−1(f(k)−∫kk+1f)≥0d_n = f(n) + \sum_{k=1}^{n-1}\left(f(k) - \int_k^{k+1}f\right) \geq 0、dn≤d1=f(1)d_n \leq d_1 = f(1)。□\square

例 5.28 (1) 実数 pp について ∑1np\sum \frac{1}{n^p} は p>1p > 1 で収束、p≤1p \leq 1 で発散する(例 5.20 (2)。例 2.36 の拡張)。 (2) ∑n≥21nlog⁡n\sum_{n \geq 2}\frac{1}{n\log n} は発散する:∫2cdxxlog⁡x=log⁡log⁡c−log⁡log⁡2→∞\int_2^c \frac{dx}{x\log x} = \log\log c - \log\log 2 \to \infty。 (3)(オイラーの定数)f(x)=1xf(x) = \frac{1}{x} とすると、dn=1+12+⋯+1n−log⁡nd_n = 1 + \frac{1}{2} + \cdots + \frac{1}{n} - \log n は収束する。その極限 γ=0.5772…\gamma = 0.5772\ldots をオイラーの定数 (Euler's constant) という。したがって Hn=log⁡n+γ+o(1)H_n = \log n + \gamma + o(1) であり、調和級数の部分和は対数の速さで発散する。

まとめ

  • リーマン積分は上和の下限(上積分)と下和の上限(下積分)が一致することで定義される。可積分性は「U(f,P)−L(f,P)U(f,P) - L(f,P) をいくらでも小さくできる」ことと同値であり、リーマン和の極限としても特徴づけられる。
  • 連続関数(ハイネ–カントールの定理による)と単調関数は可積分である。可積分性は連続性よりずっと弱い(トマエ関数)。
  • 微分積分学の基本定理:ff が連続なら ddx∫axf=f(x)\frac{d}{dx}\int_a^x f = f(x)。F′=fF' = f が可積分なら ∫abf=F(b)−F(a)\int_a^b f = F(b) - F(a)。高校の積分の定義はこの定理として回復される。
  • 置換積分・部分積分・積分剰余のテイラーの定理は基本定理から導かれる。
  • 広義積分は有界閉区間上の積分の極限であり、比較判定法と x−px^{-p} との比較で収束を判定する。絶対収束しないが収束する広義積分もある。
  • Γ(s+1)=sΓ(s)\Gamma(s+1) = s\Gamma(s)、Γ(n+1)=n!\Gamma(n+1) = n!、B(p+1,q)=pp+qB(p,q)B(p+1,q) = \frac{p}{p+q}B(p,q)。一般の B=ΓΓ/ΓB = \Gamma\Gamma/\Gamma は第9章で示す。
  • ウォリス積分から 22n(n!)2(2n)!n→π\frac{2^{2n}(n!)^2}{(2n)!\sqrt{n}} \to \sqrt{\pi} が得られ、スターリングの公式 n!∼2πn(n/e)nn! \sim \sqrt{2\pi n}(n/e)^n の定数が決まる。
  • 単調減少な正値関数では、級数と広義積分の収束・発散は一致する(積分判定法)。

演習問題

問題 5.1 ★ ∫01x2 dx=13\int_0^1 x^2\ dx = \frac{1}{3} を、定義(上和・下和)から直接示せ。

解答

nn 等分割 PnP_n では、x2x^2 が [0,1][0,1] で単調増加なので

U(f,Pn)=∑i=1ni2n2⋅1n=n(n+1)(2n+1)6n3,L(f,Pn)=∑i=1n(i−1)2n3=(n−1)n(2n−1)6n3U(f, P_n) = \sum_{i=1}^{n}\frac{i^2}{n^2}\cdot\frac{1}{n} = \frac{n(n+1)(2n+1)}{6n^3}, \qquad L(f, P_n) = \sum_{i=1}^{n}\frac{(i-1)^2}{n^3} = \frac{(n-1)n(2n-1)}{6n^3}

U−L=1n→0U - L = \frac{1}{n} \to 0 なので可積分で、∫‾≥L(f,Pn)→13\underline{\int} \geq L(f,P_n) \to \frac{1}{3}、∫‾≤U(f,Pn)→13\overline{\int} \leq U(f,P_n) \to \frac{1}{3} より積分は 13\frac{1}{3}。

問題 5.2 ★ [a,b][a, b] 上の有界関数で、不連続点が有限個しかないものは可積分であることを示せ。

解答

∣f∣≤M\lvert f \rvert \leq M とし、不連続点を c1<⋯<ckc_1 < \cdots < c_k とする。ε>0\varepsilon > 0 に対し、各 cjc_j を中心とする長さ η\eta の小区間 JjJ_j(互いに交わらないように η\eta を小さくとる)を考える。[a,b][a,b] から JjJ_j の内部を除いた部分は有限個の閉区間の和で、その上で ff は連続なので可積分であり、それぞれ U−L<ε2(k+1)U - L < \frac{\varepsilon}{2(k+1)} となる分割がとれる。これらと JjJ_j の端点を合わせた分割 PP では、JjJ_j 上の寄与は 2Mη2M\eta 以下なので U(f,P)−L(f,P)<ε2+2MkηU(f, P) - L(f, P) < \frac{\varepsilon}{2} + 2Mk\eta。η<ε4Mk\eta < \frac{\varepsilon}{4Mk} とすれば ε\varepsilon 未満。

問題 5.3 ★ 次の極限を求めよ。 (1) lim⁡n→∞∑k=1nnn2+k2\lim_{n \to \infty}\sum_{k=1}^{n}\frac{n}{n^2 + k^2} (2) lim⁡n→∞1n((n+1)(n+2)⋯(2n))1/n\lim_{n \to \infty}\frac{1}{n}\bigl((n+1)(n+2)\cdots(2n)\bigr)^{1/n}

解答

(1) 1n∑k=1n11+(k/n)2→∫01dx1+x2=arctan⁡1=π4\frac{1}{n}\sum_{k=1}^{n}\frac{1}{1 + (k/n)^2} \to \int_0^1\frac{dx}{1 + x^2} = \arctan 1 = \frac{\pi}{4}(定理 5.6)。

(2) 対数をとると 1n∑k=1nlog⁡(1+kn)→∫01log⁡(1+x) dx=[(1+x)log⁡(1+x)−x]01=2log⁡2−1\frac{1}{n}\sum_{k=1}^{n}\log\left(1 + \frac{k}{n}\right) \to \int_0^1\log(1+x)\ dx = [(1+x)\log(1+x) - x]_0^1 = 2\log 2 - 1。exp⁡\exp の連続性より極限は e2log⁡2−1=4ee^{2\log 2 - 1} = \frac{4}{e}。

問題 5.4 ★ (1) ddx∫xx2et2dt\frac{d}{dx}\int_x^{x^2}e^{t^2}dt を求めよ。(2) lim⁡x→01x3∫0xlog⁡(1+t2) dt\lim_{x \to 0}\frac{1}{x^3}\int_0^x\log(1 + t^2)\ dt を求めよ。

解答

(1) G(u)=∫0uet2dtG(u) = \int_0^u e^{t^2}dt とすると G′(u)=eu2G'(u) = e^{u^2}(定理 5.12)で、積分は G(x2)−G(x)G(x^2) - G(x)。合成関数の微分より 2xex4−ex22xe^{x^4} - e^{x^2}。

(2) ロピタルの定理 (i) と定理 5.12 より lim⁡log⁡(1+x2)3x2=13\lim\frac{\log(1+x^2)}{3x^2} = \frac{1}{3}(定理 3.37 (2))。

問題 5.5 ★★ ff が [a,b][a, b] で可積分、f≥0f \geq 0、∫abf=0\int_a^b f = 0 ならば、ff の連続点 x0x_0 では f(x0)=0f(x_0) = 0 であることを示せ。また、ff が 00 でない点をもちうることを例で示せ。

解答

連続点 x0x_0 で f(x0)>0f(x_0) > 0 なら、∣x−x0∣<δ\lvert x - x_0 \rvert < \delta で f(x)>f(x0)2f(x) > \frac{f(x_0)}{2} となる δ\delta があり、[a,b][a, b] と (x0−δ,x0+δ)(x_0 - \delta, x_0 + \delta) の共通部分を含む長さ ℓ>0\ell > 0 の閉区間 JJ がとれる。f≥0f \geq 0 と定理 5.10 より ∫abf≥∫Jf≥f(x0)2ℓ>0\int_a^b f \geq \int_J f \geq \frac{f(x_0)}{2}\ell > 0 となり矛盾。例:トマエ関数(例 5.9)は f≥0f \geq 0、∫f=0\int f = 0 だが有理数で正の値をとる。

問題 5.6 ★★ 次の広義積分の収束・発散を判定し、収束するものは値を求めよ。 (1) ∫01log⁡x dx\int_0^1 \log x\ dx (2) ∫0∞x2e−xdx\int_0^\infty x^2e^{-x}dx (3) ∫2∞dxx(log⁡x)p\int_2^\infty\frac{dx}{x(\log x)^p}(p∈Rp \in \mathbb{R}) (4) ∫01dxx(1−x)\int_0^1\frac{dx}{\sqrt{x(1-x)}}

解答

(1) ∫ε1log⁡x dx=[xlog⁡x−x]ε1=−1−εlog⁡ε+ε→−1\int_\varepsilon^1\log x\ dx = [x\log x - x]_\varepsilon^1 = -1 - \varepsilon\log\varepsilon + \varepsilon \to -1(例 4.17 (2))。値は −1-1。

(2) Γ(3)=2!=2\Gamma(3) = 2! = 2。

(3) u=log⁡xu = \log x と置換すると ∫log⁡2log⁡cduup\int_{\log 2}^{\log c}\frac{du}{u^p} なので、例 5.20 (2) より p>1p > 1 で収束(値 (log⁡2)1−pp−1\frac{(\log 2)^{1-p}}{p-1})、p≤1p \leq 1 で発散。

(4) B(12,12)=2∫0π/2dθ=πB\left(\frac{1}{2}, \frac{1}{2}\right) = 2\int_0^{\pi/2}d\theta = \pi(定理 5.24 (2))。

問題 5.7 ★ ベータ関数を用いて ∫01x3(1−x)4dx\int_0^1 x^3(1-x)^4dx と ∫0π/2sin⁡5θcos⁡3θ dθ\int_0^{\pi/2}\sin^5\theta\cos^3\theta\ d\theta を求めよ。

解答

∫01x3(1−x)4dx=B(4,5)=3! 4!8!=14440320=1280\int_0^1 x^3(1-x)^4dx = B(4, 5) = \frac{3!\ 4!}{8!} = \frac{144}{40320} = \frac{1}{280}。∫0π/2sin⁡5θcos⁡3θ dθ=12B(3,2)=12⋅2! 1!4!=124\int_0^{\pi/2}\sin^5\theta\cos^3\theta\ d\theta = \frac{1}{2}B(3, 2) = \frac{1}{2}\cdot\frac{2!\ 1!}{4!} = \frac{1}{24}。(検算:u=sin⁡θu = \sin\theta とすると ∫01u5(1−u2)du=16−18=124\int_0^1 u^5(1 - u^2)du = \frac{1}{6} - \frac{1}{8} = \frac{1}{24}。)

問題 5.8 ★★ ∑n≥21n(log⁡n)p\sum_{n \geq 2}\frac{1}{n(\log n)^p} が収束する pp の範囲を求めよ。また ∑n≥31nlog⁡nlog⁡log⁡n\sum_{n \geq 3}\frac{1}{n\log n\log\log n} の収束・発散を判定せよ。

解答

f(x)=1x(log⁡x)pf(x) = \frac{1}{x(\log x)^p} は p≥0p \geq 0 なら x≥2x \geq 2 で単調減少。積分判定法と問題 5.6 (3) より、p>1p > 1 で収束、0≤p≤10 \leq p \leq 1 で発散。p<0p < 0 なら項は 1n\frac{1}{n} 以上(n≥3n \geq 3)なので発散。よって p>1p > 1 のときに限り収束。第 2 の級数は ∫3cdxxlog⁡xlog⁡log⁡x=log⁡log⁡log⁡c−log⁡log⁡log⁡3→∞\int_3^c\frac{dx}{x\log x\log\log x} = \log\log\log c - \log\log\log 3 \to \infty なので発散する。

問題 5.9 ★★ ff が [0,1][0, 1] で連続ならば、lim⁡n→∞n∫01xnf(x) dx=f(1)\lim_{n \to \infty}n\int_0^1 x^nf(x)\ dx = f(1) であることを示せ。

解答

n∫01xndx=nn+1→1n\int_0^1 x^n dx = \frac{n}{n+1} \to 1 なので、n∫01xn(f(x)−f(1))dx→0n\int_0^1 x^n\bigl(f(x) - f(1)\bigr)dx \to 0 を示せばよい。∣f∣≤M\lvert f \rvert \leq M とし、ε>0\varepsilon > 0 に対して 1−δ≤x≤11 - \delta \leq x \leq 1 で ∣f(x)−f(1)∣<ε\lvert f(x) - f(1) \rvert < \varepsilon となる δ∈(0,1)\delta \in (0,1) をとる。

∣n∫01xn(f(x)−f(1))dx∣≤2Mn∫01−δxndx+εn∫1−δ1xndx≤2M(1−δ)n+1+ε\left\lvert n\int_0^1 x^n(f(x) - f(1))dx \right\rvert \leq 2Mn\int_0^{1-\delta}x^ndx + \varepsilon n\int_{1-\delta}^1 x^ndx \leq 2M(1-\delta)^{n+1} + \varepsilon

第 1 項は n→∞n \to \infty で 00 に収束するので、十分大きな nn で全体は 2ε2\varepsilon 未満。

問題 5.10 ★★★(ガウス積分)G=∫0∞e−x2dxG = \int_0^\infty e^{-x^2}dx が収束することを示し、G=π2G = \frac{\sqrt{\pi}}{2} を示せ。(ヒント:(1−x2)n≤e−nx2≤1(1+x2)n(1 - x^2)^n \leq e^{-nx^2} \leq \frac{1}{(1 + x^2)^n} とウォリス積分。)これから Γ(12)=π\Gamma\left(\frac{1}{2}\right) = \sqrt{\pi} を導け。

解答

収束:x≥1x \geq 1 で e−x2≤e−xe^{-x^2} \leq e^{-x}。不等式:1−u≤e−u1 - u \leq e^{-u} と eu≥1+ue^u \geq 1 + u(例 4.15)に u=x2u = x^2 を代入し nn 乗する(左は 0≤x≤10 \leq x \leq 1 で使う)。∫0∞e−nx2dx=Gn\int_0^\infty e^{-nx^2}dx = \frac{G}{\sqrt{n}}(x=y/nx = y/\sqrt{n})なので

∫01(1−x2)ndx≤∫01e−nx2dx≤Gn≤∫0∞dx(1+x2)n\int_0^1(1 - x^2)^ndx \leq \int_0^1 e^{-nx^2}dx \leq \frac{G}{\sqrt{n}} \leq \int_0^\infty\frac{dx}{(1 + x^2)^n}

左端は x=cos⁡θx = \cos\theta で I2n+1I_{2n+1}、右端は x=tan⁡θx = \tan\theta で ∫0π/2cos⁡2n−2θ dθ=I2n−2\int_0^{\pi/2}\cos^{2n-2}\theta\ d\theta = I_{2n-2}(cos⁡\cos の積分は θ↦π2−θ\theta \mapsto \frac{\pi}{2} - \theta で sin⁡\sin の積分に等しい)。よって nI2n+1≤G≤nI2n−2\sqrt{n}I_{2n+1} \leq G \leq \sqrt{n}I_{2n-2}。命題 5.25 の mIm→π/2\sqrt{m}I_m \to \sqrt{\pi/2} より

nI2n+1=n2n+12n+1 I2n+1→12π2=π2\sqrt{n}I_{2n+1} = \sqrt{\frac{n}{2n+1}}\sqrt{2n+1}\,I_{2n+1} \to \frac{1}{\sqrt{2}}\sqrt{\frac{\pi}{2}} = \frac{\sqrt{\pi}}{2}

同様に nI2n−2→π2\sqrt{n}I_{2n-2} \to \frac{\sqrt{\pi}}{2}。はさみうちにより G=π2G = \frac{\sqrt{\pi}}{2}。最後に x=y2x = y^2 と置換すると Γ(12)=∫0∞x−1/2e−xdx=2∫0∞e−y2dy=π\Gamma\left(\frac{1}{2}\right) = \int_0^\infty x^{-1/2}e^{-x}dx = 2\int_0^\infty e^{-y^2}dy = \sqrt{\pi}。

この章を読み終えたら

「読了」にすると学習記録とロードマップに反映されます。演習の自己採点もお忘れなく。

この章の誤りを報告GitHub で見る