Mathematics · Book 5 · Bachelor Year 3

University Mathematics — Year 3

University Mathematics — Year 3 · Bachelor Year 3

10The Lebesgue Integral

Riemann’s integral slices the domain into small intervals; Lebesgue’s slices the range: to integrate ff, measure the sets {f>t}\{f > t\}. The change looks innocent and is revolutionary. Limits and integrals, forever quarreling in the Riemann theory (uniform convergence required!), are reconciled by three convergence theorems — monotone convergence, Fatou, dominated convergence — whose hypotheses are almost embarrassingly weak. This chapter constructs the integral over an arbitrary measure space (X,A,μ)(X, \mathcal A, \mu), proves the three theorems, settles the exact relationship with Riemann’s integral (a bounded function is Riemann-integrable iff it is continuous almost everywhere), and industrializes the differentiation of parameter-dependent integrals — the technique that the weekend problem uses to compute ∫0∞sin⁡xx  ⁣dx\int_0^\infty\frac{\sin x}x\,\dd x and ∫Re−x2 ⁣dx\int_\R \eu^{-x^2}\dd x.

10.1 Measurable functions

Definition 10.1

Let (X,A)(X, \mathcal A), (Y,B)(Y, \mathcal B) be measurable spaces. f ⁣:X→Yf \colon X \to Y is measurable if f−1(B)∈Af^{-1}(B) \in \mathcal A for every B∈BB \in \mathcal B. For real (or [−∞,+∞][-\infty,+\infty]-valued) functions, Y=RY = \R carries its Borel σ\sigma-algebra, and it suffices to check f−1((t, +∞))={f>t}∈Af^{-1}(\intoo t{+\infty}) = \{f > t\} \in \mathcal A for all t∈Rt \in \R: the good sets {B:f−1(B)∈A}\{B : f^{-1}(B) \in \mathcal A\} form a σ\sigma-algebra (preimages commute with set operations) containing the generating rays (Definition 9.2, Method 9.17).

Proposition 10.2

(a) Compositions of measurable maps are measurable; continuous maps are Borel-measurable. (b) If f,g ⁣:X→Rf, g \colon X \to \R are measurable, so are f+gf + g, fgfg, max⁡(f,g)\max(f,g), ∣f∣\abs f, λf\lambda f. (c) If (fn)(f_n) are measurable with values in [−∞,+∞][-\infty, +\infty], then sup⁡nfn\sup_nf_n, inf⁡nfn\inf_nf_n, lim sup⁡fn\limsup f_n, lim inf⁡fn\liminf f_n are measurable; if fn→ff_n \to f pointwise, ff is measurable.

Proof. (a) (g∘f)−1(B)=f−1(g−1(B))(g\circ f)^{-1}(B) = f^{-1}(g^{-1}(B)); continuity gives measurability via the generating open sets (Problem 9.1, question 10, in general form). (b) (f,g) ⁣:X→R2(f, g) \colon X \to \R^2 is measurable for the Borel σ\sigma-algebra of R2\R^2 — check on open boxes, which generate (R2\R^2’s opens are countable unions of rational boxes): (f,g)−1(U×V)=f−1(U)∩g−1(V)(f,g)^{-1}(U\times V) = f^{-1}(U)\cap g^{-1}(V) — and +,×,max⁡+, \times, \max are continuous R2→R\R^2 \to \R: compose. (c) {sup⁡fn>t}=⋃n{fn>t}\{\sup f_n > t\} = \bigcup_n\{f_n > t\}; inf⁡=−sup⁡(−)\inf = -\sup(-); lim sup⁡=inf⁡Nsup⁡n≥N\limsup = \inf_N\sup_{n \geq N}; a pointwise limit is its own lim sup⁡\limsup. ∎

Definition 10.3

A simple function is a measurable function with finitely many values: s=∑i=1nci 1Ais = \sum_{i=1}^n c_i\,\mathbf 1_{A_i}, Ai∈AA_i \in \mathcal A disjoint, ci≥0c_i \geq 0 (for the nonnegative theory). Its integral is

∫s  ⁣dμ=∑ici μ(Ai)∈[0,+∞]\int s\,\dd\mu = \sum_i c_i\,\mu(A_i) \in [0, +\infty]

(convention 0⋅∞=00\cdot\infty = 0); the value does not depend on the representation (refine two partitions).

Theorem 10.4 (Approximation by simple functions)

Every measurable f ⁣:X→[0,+∞]f \colon X \to [0, +\infty] is the pointwise limit of an increasing sequence of simple functions:

sn=∑k=1n2nk−12n 1{k−12n≤f<k2n}+n 1{f≥n}↗f.s_n = \sum_{k=1}^{n2^n} \frac{k-1}{2^n}\, \mathbf 1_{\{\frac{k-1}{2^n} \leq f < \frac k{2^n}\}} + n\,\mathbf 1_{\{f \geq n\}} \nearrow f .

Proof. Each sns_n is simple (the sets are preimages of Borel sets). Monotonicity: passing from nn to n+1n+1 splits each dyadic level in two and never decreases the assigned value (a point with k−12n≤f(x)<k2n\frac{k-1}{2^n} \leq f(x) < \frac k{2^n} gets either 2k−22n+1\frac{2k-2}{2^{n+1}} or 2k−12n+1\frac{2k-1}{2^{n+1}}, both ≥k−12n\geq \frac{k-1}{2^n}; the cap nn rises too). Convergence: if f(x)<∞f(x) < \infty, for n>f(x)n > f(x) we have f(x)−sn(x)≤2−nf(x) - s_n(x) \leq 2^{-n}; if f(x)=∞f(x) = \infty, sn(x)=n→∞s_n(x) = n \to \infty. ∎

10.2 The integral and the convergence theorems

Definition 10.5

For measurable f≥0f \geq 0:

∫f  ⁣dμ=sup⁡{∫s  ⁣dμ:s simple, 0≤s≤f}∈[0,+∞].\int f \,\dd\mu = \sup\Bigl\{\int s\,\dd\mu : s \text{ simple}, \ 0 \leq s \leq f\Bigr\} \in [0, +\infty].

It is monotone in ff by construction, and extends the simple case (for simple ff, the sup is attained at ff: comparison of simple integrals via common refinements).

Theorem 10.6 (Monotone convergence, Beppo Levi)

If 0≤fn↗f0 \leq f_n \nearrow f pointwise (measurable), then

∫fn  ⁣dμ↗∫f  ⁣dμ.\int f_n\,\dd\mu \nearrow \int f\,\dd\mu .

Proof. ff is measurable (Proposition 10.2(c)) and ∫fn\int f_n increases to some L≤∫fL \leq \int f (monotonicity). Conversely, fix a simple s=∑ci1Ai≤fs = \sum c_i\mathbf 1_{A_i} \leq f and θ∈(0,1)\theta \in (0,1); the sets En={fn≥θs}E_n = \{f_n \geq \theta s\} are measurable and increase to XX (where s(x)>0s(x) > 0: f(x)≥s(x)>θs(x)f(x) \geq s(x) > \theta s(x), so eventually fn(x)≥θs(x)f_n(x) \geq \theta s(x); where s(x)=0s(x) = 0: trivially). Then

∫fn≥∫Enθs  ⁣dμ=θ∑ici μ(Ai∩En)→n→∞θ∑ici μ(Ai)=θ∫s\int f_n \geq \int_{E_n}\theta s\,\dd\mu = \theta\sum_i c_i\,\mu(A_i \cap E_n) \xrightarrow[n\to\infty]{} \theta\sum_ic_i\,\mu(A_i) = \theta\int s

by continuity from below (Proposition 9.6(c)). So L≥θ∫sL \geq \theta\int s for all θ<1\theta < 1 and all simple s≤fs \leq f: L≥∫fL \geq \int f. ∎

Corollary 10.7

For measurable f,g≥0f, g \geq 0 and c≥0c \geq 0: ∫(f+g)=∫f+∫g\int(f + g) = \int f + \int g and ∫cf=c∫f\int cf = c\int f; for a series of nonnegative measurable functions, ∫∑nfn=∑n∫fn\int\sum_nf_n = \sum_n\int f_n.

Proof. For simple functions, additivity is a computation on a common refinement. In general take sn↗fs_n \nearrow f, tn↗gt_n \nearrow g (Theorem 10.4): sn+tn↗f+gs_n + t_n \nearrow f + g, and MCT passes additivity to the limit. The series statement is MCT applied to the partial sums. ∎

Theorem 10.8 (Fatou’s lemma)

For measurable fn≥0f_n \geq 0:

∫lim inf⁡nfn  ⁣dμ  ≤  lim inf⁡n∫fn  ⁣dμ.\int \liminf_n f_n \,\dd\mu \;\leq\; \liminf_n \int f_n\,\dd\mu .

Proof. Let gN=inf⁡n≥Nfng_N = \inf_{n\geq N}f_n: measurable, 0≤gN↗lim inf⁡fn0 \leq g_N \nearrow \liminf f_n, and gN≤fng_N \leq f_n for every n≥Nn \geq N, so ∫gN≤inf⁡n≥N∫fn\int g_N \leq \inf_{n \geq N}\int f_n. Apply MCT to the left side: ∫lim inf⁡fn=lim⁡N∫gN≤lim⁡Ninf⁡n≥N∫fn=lim inf⁡∫fn\int\liminf f_n = \lim_N\int g_N \leq \lim_N\inf_{n\geq N}\int f_n = \liminf\int f_n. ∎

Definition 10.9

A measurable f ⁣:X→Rf \colon X \to \R (or C\C) is integrable if ∫∣f∣  ⁣dμ<∞\int\abs f\,\dd\mu < \infty; then ∫f=∫f+−∫f−\int f = \int f^+ - \int f^- (positive and negative parts; real and imaginary parts in the complex case). The integral is linear on integrable functions (decompose and recombine positive parts; the complex case reduces to the real one) and satisfies ∣∫f∣≤∫∣f∣\abs{\int f} \leq \int\abs f (real case: ±∫f=∫(±f)≤∫∣f∣\pm\int f = \int(\pm f) \leq \int\abs f; complex case: multiply by a unimodular constant to make the integral real). A property holds almost everywhere (a.e.) if it fails only on a μ\mu-null set; modifying ff on a null set changes no integral (the difference is dominated by ∞⋅1N\infty\cdot\mathbf 1_N, of integral 00).

Theorem 10.10 (Dominated convergence)

Let fn→ff_n \to f a.e., with ∣fn∣≤g\abs{f_n} \leq g a.e. for a fixed integrable gg. Then ff is integrable and

∫fn  ⁣dμ⟶∫f  ⁣dμ,indeed∫∣fn−f∣  ⁣dμ→0.\int f_n\,\dd\mu \longrightarrow \int f\,\dd\mu, \qquad\text{indeed}\quad \int\abs{f_n - f}\,\dd\mu \to 0 .

Proof. Discard a null set to make the hypotheses pointwise. ∣f∣≤g\abs f \leq g: ff is integrable. The functions hn=2g−∣fn−f∣≥0h_n = 2g - \abs{f_n - f} \geq 0 satisfy lim inf⁡hn=2g\liminf h_n = 2g; Fatou gives

∫2g≤lim inf⁡∫(2g−∣fn−f∣)=∫2g−lim sup⁡∫∣fn−f∣,\int 2g \leq \liminf\int\bigl(2g - \abs{f_n - f}\bigr) = \int 2g - \limsup\int\abs{f_n - f},

so lim sup⁡∫∣fn−f∣≤0\limsup\int\abs{f_n - f} \leq 0 (the subtraction is legal: ∫2g<∞\int 2g < \infty). Finally ∣∫fn−∫f∣≤∫∣fn−f∣→0\abs{\int f_n - \int f} \leq \int\abs{f_n - f} \to 0. ∎

Method 10.11

Faced with lim⁡n∫fn\lim_n\int f_n: try, in order — (1) is the sequence monotone (or a series of nonnegative terms)? MCT, no integrability needed. (2) Is there a single integrable dominator g≥∣fn∣g \geq \abs{f_n}, found by crude bounds (“sup⁡n\sup_n” the estimates)? DCT. (3) No domination, no monotonicity? Fatou still bounds one side, and equality may genuinely fail: the escaping bump fn=n1(0, 1/n)f_n = n\mathbf 1_{\intoo0{1/n}} has ∫fn=1\int f_n = 1 but fn→0f_n \to 0 a.e. Domination is exactly what forbids mass from escaping to infinity, vertically or horizontally.

10.3 Riemann versus Lebesgue

Theorem 10.12 (Lebesgue’s criterion)

Let f ⁣:[a, b]→Rf \colon \intcc ab \to \R be bounded. Then ff is Riemann-integrable iff ff is continuous λ\lambda-almost everywhere; in that case ff is Lebesgue-integrable and the two integrals coincide.

Proof. For a subdivision σ=(a=x0<⋯<xN=b)\sigma = (a = x_0 < \dots < x_N = b), let LσL_\sigma and UσU_\sigma be the step functions equal, on each (xi−1, xi)\intoo{x_{i-1}}{x_i}, to mi=inf⁡[xi−1,xi]fm_i = \inf_{[x_{i-1}, x_i]}f and Mi=sup⁡M_i = \sup; the Darboux sums are their integrals (Riemann and Lebesgue agree on step functions, both giving ∑miΔxi\sum m_i\Delta x_i). Take a sequence of subdivisions σn\sigma_n, each refining the last, of mesh →0\to 0, with Darboux sums converging to the lower and upper Darboux integrals of ff. The refinements make LσnL_{\sigma_n} nondecreasing and UσnU_{\sigma_n} nonincreasing pointwise off the countable set DD of all division points; call the limits ℓ\ell and uu (measurable, Proposition 10.2). For x∉Dx \notin D, writing In(x)I_n(x) for the open σn\sigma_n-interval containing xx: ℓ(x)=sup⁡ninf⁡In(x)f\ell(x) = \sup_n\inf_{I_n(x)}f and u(x)=inf⁡nsup⁡In(x)fu(x) = \inf_n\sup_{I_n(x)}f; since the meshes shrink to 00, these are the lower and upper envelopes of ff at xx — u(x)−ℓ(x)u(x) - \ell(x) is the oscillation of ff at xx — so that ℓ(x)=u(x)\ell(x) = u(x) iff ff is continuous at xx. By MCT/DCT (bounded, finite interval):

∫[a, b]ℓ  ⁣dλ=lim⁡n∫Lσn=∫‾f,∫[a, b]u  ⁣dλ=∫‾f.\int_{\intcc ab}\ell\,\dd\lambda = \lim_n\int L_{\sigma_n} = \underline{\int}f, \qquad \int_{\intcc ab}u\,\dd\lambda = \overline{\int}f .

ff Riemann-integrable   ⟺  \iff ∫‾f=∫‾f\underline\int f = \overline\int f   ⟺  \iff ∫(u−ℓ)=0\int(u - \ell) = 0   ⟺  \iff u=ℓu = \ell a.e. (u−ℓ≥0u - \ell \geq 0; Exercise 10.5)   ⟺  \iff ff continuous a.e. In that case ℓ≤f≤u\ell \leq f \leq u with ℓ=u\ell = u a.e.: ff equals the measurable ℓ\ell a.e., hence is Lebesgue-measurable (completeness of λ\lambda) with ∫f  ⁣dλ=∫ℓ  ⁣dλ=∫‾f=∫abf\int f\,\dd\lambda = \int\ell\,\dd\lambda = \underline\int f = \int_a^bf. ∎

Example 10.13

1Q\mathbf 1_\Q is nowhere continuous: not Riemann-integrable — but Lebesgue-trivial: ∫1Q  ⁣dλ=λ(Q)=0\int\mathbf 1_\Q\,\dd\lambda = \lambda(\Q) = 0. Thomae’s function ( 1q\frac1q at rationals pq\frac pq, 00 elsewhere) is continuous exactly at the irrationals: Riemann-integrable with integral 00. And improper Riemann integrals are a different notion: ∫0∞sin⁡xx  ⁣dx\int_0^\infty\frac{\sin x}x\,\dd x converges as a limit of ∫0A\int_0^A (the weekend problem computes it =π2= \frac\pi2), but sin⁡xx∉L1((0, +∞))\frac{\sin x}x \notin L^1(\intoo0{+\infty}): the absolute integral diverges like the harmonic series (Exercise 10.6). Lebesgue’s theory trades conditional convergence for robust limit theorems.

10.4 Integrals with parameters

Throughout, (X,A,μ)(X, \mathcal A, \mu) is a measure space, TT a metric space (the parameter), and f ⁣:T×X→Cf \colon T \times X \to \C with f(t,⋅)f(t, \cdot) integrable for each tt; set F(t)=∫Xf(t,x)  ⁣dμ(x)F(t) = \int_X f(t, x)\,\dd\mu(x).

Theorem 10.14 (Continuity)

Suppose: t↦f(t,x)t \mapsto f(t,x) is continuous at t0t_0 for a.e. xx, and there is an integrable gg with ∣f(t,x)∣≤g(x)\abs{f(t,x)} \leq g(x) for all tt in a neighborhood of t0t_0 and a.e. xx. Then FF is continuous at t0t_0.

Proof. For any sequence tn→t0t_n \to t_0: f(tn,⋅)→f(t0,⋅)f(t_n, \cdot) \to f(t_0, \cdot) a.e., dominated by gg: DCT gives F(tn)→F(t0)F(t_n) \to F(t_0); sequential continuity suffices in metric spaces (Remark 6.8). ∎

Theorem 10.15 (Differentiation under the integral)

Let TT be an open interval of R\R. Suppose: for a.e. xx, t↦f(t,x)t \mapsto f(t,x) is differentiable on TT, with

∣∂f∂t(t,x)∣≤g(x)for all t∈T, a.e. x,\Bigl|\frac{\partial f}{\partial t}(t, x)\Bigr| \leq g(x) \quad \text{for all } t \in T,\ \text{a.e. } x,

gg integrable. Then FF is differentiable on TT with F′(t)=∫X∂f∂t(t,x)  ⁣dμ(x)F'(t) = \int_X \frac{\partial f}{\partial t}(t, x)\,\dd\mu(x).

Proof. Fix tt and hn→0h_n \to 0: the difference quotients

φn(x)=f(t+hn,x)−f(t,x)hn⟶∂f∂t(t,x)a.e.,\varphi_n(x) = \frac{f(t + h_n, x) - f(t, x)}{h_n} \longrightarrow \frac{\partial f}{\partial t}(t,x) \quad\text{a.e.},

and the mean value inequality bounds ∣φn(x)∣≤sup⁡s∣∂tf(s,x)∣≤g(x)\abs{\varphi_n(x)} \leq \sup_{s}\abs{\partial_tf(s,x)} \leq g(x): DCT applies, and F(t+hn)−F(t)hn=∫φn→∫∂tf(t,⋅)\frac{F(t + h_n) - F(t)}{h_n} = \int\varphi_n \to \int\partial_t f(t, \cdot). ∎

Example 10.16 (The Gamma function)

For t>0t > 0 let

Γ(t)=∫0+∞xt−1e−x  ⁣dx.\Gamma(t) = \int_0^{+\infty} x^{t-1}\eu^{-x}\,\dd x .

The integral converges: near 00, xt−1x^{t-1} is integrable (t>0t > 0); at infinity, xt−1e−x≤Ce−x/2x^{t-1}\eu^{-x} \leq C\eu^{-x/2}. Integration by parts (on [ε,A][\varepsilon, A], then limits via MCT) gives the functional equation Γ(t+1)=t Γ(t)\Gamma(t + 1) = t\,\Gamma(t), whence Γ(n+1)=n!\Gamma(n+1) = n!: the factorial interpolated. On every [a, b]⊆(0, +∞)\intcc ab \subseteq \intoo0{+\infty}, ∂t(xt−1e−x)=ln⁡x⋅xt−1e−x\partial_t\bigl(x^{t-1}\eu^{-x}\bigr) = \ln x\cdot x^{t-1}\eu^{-x} is dominated by ∣ln⁡x∣(xa−1+xb−1)e−x\abs{\ln x}(x^{a-1} + x^{b-1})\eu^{-x}, integrable: Γ\Gamma is C1\mathcal C^1, and by induction C∞\mathcal C^\infty, with Γ(k)(t)=∫0∞(ln⁡x)kxt−1e−x ⁣dx\Gamma^{(k)}(t) = \int_0^\infty(\ln x)^kx^{t-1}\eu^{-x}\dd x. The value Γ(12)=π\Gamma(\frac12) = \sqrt\pi is the Gaussian integral in disguise (Problem 10.1).

The integrand xx: the improper integral ∈t_0∈fty converges by alternating cancellation between the arches, but the areas | | of the arches behave like 2π k — a harmonic series: xx ∉ L1. Lebesgue integrability is absolute integrability.
The integrand sin⁡xx\frac{\sin x}x: the improper integral ∫0∞\int_0^\infty converges by alternating cancellation between the arches, but the areas ∣⋅∣\abs{\cdot} of the arches behave like 2πk\frac2{\pi k} — a harmonic series: sin⁡xx∉L1\frac{\sin x}x \notin L^1. Lebesgue integrability is absolute integrability.

10.5 Exercises

Exercise 10.1 ★

(a) Show that a monotone function R→R\R \to \R is Borel-measurable, and that a derivative (of an everywhere differentiable function) is Borel-measurable. (b) Show that f ⁣:X→Rf \colon X \to \R is measurable iff {f>q}∈A\{f > q\} \in \mathcal A for every rational qq.

Solution

Solution of Exercise 10.1.

(a) If ff is nondecreasing, {f>t}\{f > t\} is ∅\varnothing, R\R, or a ray (a, +∞)\intoo a{+\infty} / [a, +∞)\intco a{+\infty}: Borel in every case; nonincreasing likewise. A derivative: f′(x)=lim⁡nn(f(x+1n)−f(x))f'(x) = \lim_n n\bigl(f(x + \frac1n) - f(x)\bigr) is a pointwise limit of continuous (hence measurable) functions: Proposition 10.2(c).

(b) {f>t}=⋃q∈Q, q>t{f>q}\{f > t\} = \bigcup_{q \in \Q,\, q > t}\{f > q\}: if the rational levels are measurable, all levels are, and the rays generate B(R)\mathcal B(\R).

Exercise 10.2 ★

Compute, with full justification:

lim⁡n→∞∫0+∞cos⁡x(1+x/n)n  ⁣dx,lim⁡n→∞∫01n xn−11+x  ⁣dx.\lim_{n\to\infty}\int_0^{+\infty} \frac{\cos x}{(1 + x/n)^{n}}\,\dd x, \qquad \lim_{n\to\infty}\int_0^1 \frac{n\,x^{n-1}}{1 + x}\,\dd x .

(For the second: substitute u=xnu = x^n before dominating.)

Solution

Solution of Exercise 10.2.

First: (1+x/n)n↗ex(1 + x/n)^n \nearrow \eu^x for x≥0x \geq 0, so the integrand tends pointwise to e−xcos⁡x\eu^{-x}\cos x; for n≥2n \geq 2, (1+x/n)n≥(1+x/2)2(1 + x/n)^n \geq (1 + x/2)^2, giving the integrable dominator (1+x/2)−2(1 + x/2)^{-2}. DCT:

lim⁡n∫0∞cos⁡x(1+x/n)n ⁣dx=∫0∞e−xcos⁡x  ⁣dx=Re⁡∫0∞e−(1−i)x ⁣dx=Re⁡11−i=12.\lim_n\int_0^\infty\frac{\cos x}{(1 + x/n)^n}\dd x = \int_0^\infty \eu^{-x}\cos x\,\dd x = \operatorname{Re}\int_0^\infty\eu^{-(1 - \iu)x}\dd x = \operatorname{Re}\frac{1}{1 - \iu} = \frac12 .

Second: substitute u=xnu = x^n (a C1\mathcal C^1 bijection of (0, 1)\intoo01):

∫01nxn−11+x ⁣dx=∫01 ⁣du1+u1/n⟶∫01 ⁣du2=12,\int_0^1\frac{nx^{n-1}}{1 + x}\dd x = \int_0^1\frac{\dd u}{1 + u^{1/n}} \longrightarrow \int_0^1\frac{\dd u}{2} = \frac12,

by DCT: for u∈(0, 1)u \in \intoo01, u1/n→1u^{1/n} \to 1, and the integrand is bounded by 11 on a finite measure space.

Exercise 10.3 ★★

(a) Exhibit strict inequality in Fatou’s lemma. (b) Exhibit fn→0f_n \to 0 pointwise with ∫fn=1\int f_n = 1 in three ways: escape in height, in width, to infinity. Which single hypothesis of DCT does each violate? (c) Show that in Fatou’s lemma one cannot replace lim inf⁡\liminf by lim sup⁡\limsup on either side.

Solution

Solution of Exercise 10.3.

(a) fn=n 1(0, 1/n)f_n = n\,\mathbf 1_{\intoo0{1/n}}: lim inf⁡fn=0\liminf f_n = 0 pointwise, ∫fn=1\int f_n = 1: 0<10 < 1.

(b) Height: n1(0, 1/n)n\mathbf 1_{\intoo0{1/n}}; width: 1n1(0, n)\frac1n\mathbf 1_{\intoo0n}; translation: 1(n, n+1)\mathbf 1_{\intoo n{n+1}}. All tend to 00 pointwise with ∫=1\int = 1. In each case the domination hypothesis fails: sup⁡nfn\sup_nf_n is ≈1/x\approx 1/x near 00, ≈\approx a nonintegrable constant profile, 1(1, ∞)\mathbf 1_{\intoo1\infty}-like — never integrable.

(c) “∫lim sup⁡fn≥lim sup⁡∫fn\int\limsup f_n \geq \limsup\int f_n” fails for the translating bump: left side 00, right side 11. “lim sup⁡∫≤∫lim sup⁡\limsup\int \leq \int\limsup” is the same statement. And Fatou for lim sup⁡\limsup with ≤\leq reversed (“reverse Fatou”) requires a dominator — the same bump is the counterexample.

Exercise 10.4 ★★

(a) Show ∫0+∞xex−1  ⁣dx=∑n≥11n2=π26\displaystyle\int_0^{+\infty}\frac{x}{\eu^x - 1}\,\dd x = \sum_{n\geq1}\frac1{n^2} = \frac{\pi^2}6 (expand 1ex−1\frac1{\eu^x - 1} in a geometric series and integrate term by term — which theorem permits it?). (b) (Sophomore’s dream) Show ∫01x−x  ⁣dx=∑n≥1n−n\displaystyle\int_0^1 x^{-x}\,\dd x = \sum_{n\geq1}n^{-n}. (Write x−x=e−xln⁡x=∑k(−xln⁡x)kk!x^{-x} = \eu^{-x\ln x} = \sum_k\frac{(-x\ln x)^k}{k!} and compute ∫01(−xln⁡x)k ⁣dx\int_0^1(-x\ln x)^k\dd x by substituting x=e−u/(k+1)x = \eu^{-u/(k+1)}, recognizing Γ\Gamma.)

Solution

Solution of Exercise 10.4.

(a) For x>0x > 0: 1ex−1=e−x1−e−x=∑n≥1e−nx\frac1{\eu^x - 1} = \frac{\eu^{-x}}{1 - \eu^{-x}} = \sum_{n\geq1}\eu^{-nx}, so xex−1=∑n≥1xe−nx\frac{x}{\eu^x - 1} = \sum_{n\geq1}x\eu^{-nx}, a series of nonnegative measurable functions: Corollary 10.7 allows term-by-term integration:

∫0∞x  ⁣dxex−1=∑n≥1∫0∞xe−nx ⁣dx=∑n≥11n2=π26\int_0^\infty\frac{x\,\dd x}{\eu^x - 1} = \sum_{n\geq1}\int_0^\infty x\eu^{-nx}\dd x = \sum_{n\geq1}\frac1{n^2} = \frac{\pi^2}6

(∫0∞xe−nx ⁣dx=n−2\int_0^\infty x\eu^{-nx}\dd x = n^{-2} by parts; Basel from the Year 2 volume, or Exercise 13.5 to come).

(b) On (0, 1)\intoo01, −xln⁡x≥0-x\ln x \geq 0, so x−x=e−xln⁡x=∑k(−xln⁡x)kk!x^{-x} = \eu^{-x\ln x} = \sum_k\frac{(-x\ln x)^k}{k!} is a series of nonnegative terms: interchange again. Substituting x=e−u/(k+1)x = \eu^{-u/(k+1)}:

∫01(−xln⁡x)k ⁣dx=∫0∞(uk+1)ke−kuk+1  e−uk+1k+1  ⁣du=1(k+1)k+1∫0∞uke−u ⁣du=k!(k+1)k+1.\int_0^1(-x\ln x)^k\dd x = \int_0^\infty\Bigl(\frac{u}{k+1}\Bigr)^{k} \eu^{-\frac{ku}{k+1}}\;\frac{\eu^{-\frac u{k+1}}}{k+1}\,\dd u = \frac{1}{(k+1)^{k+1}}\int_0^\infty u^k\eu^{-u}\dd u = \frac{k!}{(k+1)^{k+1}} .

Hence ∫01x−x ⁣dx=∑k≥01(k+1)k+1=∑n≥1n−n\int_0^1x^{-x}\dd x = \sum_{k\geq0}\frac{1}{(k+1)^{k+1}} = \sum_{n\geq1}n^{-n}: the sophomore’s dream, rigorously.

Exercise 10.5 ★★

(a) Show that f≥0f \geq 0 measurable with ∫f  ⁣dμ=0\int f\,\dd\mu = 0 satisfies f=0f = 0 a.e. (Consider {f≥1/n}\{f \geq 1/n\} and Markov’s inequality: μ({f≥a})≤1a∫f\mu(\{f \geq a\}) \leq \frac1a\int f — prove it.) (b) Show that an integrable ff is finite a.e. (c) Show that if ∫Af  ⁣dμ=0\int_A f\,\dd\mu = 0 for every measurable AA, then f=0f = 0 a.e.

Solution

Solution of Exercise 10.5.

(a) Markov: a 1{f≥a}≤fa\,\mathbf 1_{\{f \geq a\}} \leq f, integrate: μ({f≥a})≤1a∫f\mu(\{f \geq a\}) \leq \frac1a\int f. If ∫f=0\int f = 0: μ({f≥1n})=0\mu(\{f \geq \frac1n\}) = 0 for every nn, and {f>0}=⋃n{f≥1n}\{f > 0\} = \bigcup_n\{f \geq \frac1n\} is null.

(b) μ({∣f∣=∞})≤μ({∣f∣≥n})≤1n∫∣f∣→0\mu(\{\abs f = \infty\}) \leq \mu(\{\abs f \geq n\}) \leq \frac1n\int\abs f \to 0.

(c) Take A={f>0}A = \{f > 0\}: ∫f+ ⁣dμ=∫Af  ⁣dμ=0\int f^+\dd\mu = \int_Af\,\dd\mu = 0, so f+=0f^+ = 0 a.e. by (a); likewise f−=0f^- = 0 a.e.

Exercise 10.6 ★★

(a) Apply Theorem 10.12 to decide Riemann integrability of: 1Q\mathbf 1_\Q; Thomae’s function; 1K\mathbf 1_K for KK a fat Cantor set (Exercise 9.5). (b) Show that ∫1+∞∣sin⁡xx∣  ⁣dx=+∞\int_1^{+\infty}\abs{\frac{\sin x}x}\,\dd x = +\infty, while lim⁡A→∞∫1Asin⁡xx  ⁣dx\lim_{A\to\infty}\int_1^A\frac{\sin x}x\,\dd x exists (integrate by parts): improper convergence without integrability.

Solution

Solution of Exercise 10.6.

(a) 1Q\mathbf 1_\Q: discontinuous everywhere, not Riemann-integrable (Theorem 10.12); its Lebesgue integral is λ(Q)=0\lambda(\Q) = 0. Thomae: continuous at every irrational (given ε\varepsilon, only finitely many rationals in [0, 1]\intcc01 have denominator ≤1/ε\leq 1/\varepsilon; avoid them by a small neighborhood), discontinuous at rationals (density of irrationals): continuous a.e., Riemann-integrable, integral 00 (it vanishes a.e.). 1K\mathbf 1_K, KK a fat Cantor set: the discontinuity set is ∂K=K\partial K = K (closed with empty interior), of measure 12>0\frac12 > 0: not Riemann-integrable — yet Lebesgue-integrable with integral λ(K)=12\lambda(K) = \frac12.

(b) ∫kπ(k+1)π∣sin⁡x∣x ⁣dx≥1(k+1)π∫kπ(k+1)π∣sin⁡x∣ ⁣dx=2(k+1)π\int_{k\pi}^{(k+1)\pi}\frac{\abs{\sin x}}x\dd x \geq \frac1{(k+1)\pi}\int_{k\pi}^{(k+1)\pi}\abs{\sin x}\dd x = \frac{2}{(k+1)\pi}: the series diverges. Convergence of the improper integral: for A>πA > \pi,

∫πAsin⁡xx ⁣dx=[−cos⁡xx]πA−∫πAcos⁡xx2 ⁣dx,\int_\pi^A\frac{\sin x}x\dd x = \Bigl[\frac{-\cos x}x\Bigr]_\pi^A - \int_\pi^A\frac{\cos x}{x^2}\dd x,

and both terms converge as A→∞A \to \infty (1x2\frac1{x^2} is integrable): conditional convergence without absolute integrability.

Exercise 10.7 ★★

Justify that F(t)=∫0+∞e−x2cos⁡(tx)  ⁣dxF(t) = \int_0^{+\infty}\eu^{-x^2}\cos(tx)\,\dd x is C1\mathcal C^1 on R\R and satisfies F′(t)=−t2F(t)F'(t) = -\frac t2F(t) (integrate by parts); deduce F(t)=F(0) e−t2/4F(t) = F(0)\,\eu^{-t^2/4}. (With F(0)=π2F(0) = \frac{\sqrt\pi}2 from the weekend problem: the Gaussian is essentially its own Fourier transform — Chapter 14 will systematize this.)

Solution

Solution of Exercise 10.7.

Domination: ∣∂t(e−x2cos⁡(tx))∣=∣xe−x2sin⁡(tx)∣≤xe−x2\abs{\partial_t(\eu^{-x^2}\cos(tx))} = \abs{x\eu^{-x^2}\sin(tx)} \leq x\eu^{-x^2}, integrable and independent of tt: Theorem 10.15 applies globally,

F′(t)=−∫0∞xe−x2sin⁡(tx)  ⁣dx=[12e−x2sin⁡(tx)]0∞−t2∫0∞e−x2cos⁡(tx) ⁣dx=−t2F(t)F'(t) = -\int_0^\infty x\eu^{-x^2}\sin(tx)\,\dd x = \Bigl[\tfrac12\eu^{-x^2}\sin(tx)\Bigr]_0^\infty - \frac t2\int_0^\infty\eu^{-x^2}\cos(tx)\dd x = -\frac t2F(t)

(integration by parts with  ⁣dv=xe−x2 ⁣dx\dd v = x\eu^{-x^2}\dd x). The linear ODE gives F(t)=F(0)e−t2/4F(t) = F(0)\eu^{-t^2/4}; with F(0)=π2F(0) = \frac{\sqrt\pi}2 (Problem 10.1), the Gaussian reproduces itself under this cosine transform.

Exercise 10.8 ★★★

(Frullani) Let 0<a<b0 < a < b. Show

∫0+∞e−ax−e−bxx  ⁣dx=ln⁡ba,\int_0^{+\infty}\frac{\eu^{-ax} - \eu^{-bx}}{x}\,\dd x = \ln\frac ba,

by writing the integrand as ∫abe−xt  ⁣dt\int_a^b \eu^{-xt}\,\dd t and justifying the interchange via the nonnegative theory (Corollary 10.7 in continuous form — anticipate Tonelli, or slice [a,b][a,b] into nn equal parts and pass to the limit).

Solution

Solution of Exercise 10.8.

The integral converges: near 00 the integrand tends to b−ab - a (bounded), and it decays like e−ax\eu^{-ax} at infinity. Fix aa and view I(b)=∫0∞e−ax−e−bxx ⁣dxI(b) = \int_0^\infty\frac{\eu^{-ax} - \eu^{-bx}}x\dd x as a function of b∈[a, +∞)b \in \intco a{+\infty}. For all b≥ab \geq a: ∣∂b(integrand)∣=e−bx≤e−ax\abs{\partial_b(\text{integrand})} = \eu^{-bx} \leq \eu^{-ax}, and ∫0∞e−ax ⁣dx=1a<∞\int_0^\infty\eu^{-ax}\dd x = \frac1a < \infty: an integrable dominator. So Theorem 10.15 gives I′(b)=∫0∞e−bx ⁣dx=1bI'(b) = \int_0^\infty\eu^{-bx}\dd x = \frac1b, and I(a)=0I(a) = 0:

I(b)=∫ab ⁣dtt=ln⁡ba.I(b) = \int_a^b\frac{\dd t}t = \ln\frac ba .

(Equivalently, the hint’s route: the integrand is ∫abe−xt ⁣dt≥0\int_a^b\eu^{-xt}\dd t \geq 0 and the interchange is the continuous analogue of Corollary 10.7, i.e. Tonelli — proved in Chapter 11; the parameter route stays within this chapter.)

Exercise 10.9 ★★

Let f≥0f \geq 0 be measurable on (X,A,μ)(X, \mathcal A, \mu). Show that ν(A)=∫Af  ⁣dμ\nu(A) = \int_A f\,\dd\mu defines a measure (density ff with respect to μ\mu), and that ∫g  ⁣dν=∫gf  ⁣dμ\int g\,\dd\nu = \int gf\,\dd\mu for all measurable g≥0g \geq 0 (prove it for indicators, then simple functions, then MCT — the standard machine).

Solution

Solution of Exercise 10.9.

ν(∅)=0\nu(\varnothing) = 0; for disjoint (An)(A_n), f1⨆An=∑nf1Anf\mathbf 1_{\bigsqcup A_n} = \sum_nf\mathbf 1_{A_n} (pointwise, all terms ≥0\geq 0), and Corollary 10.7 gives σ\sigma-additivity. The formula ∫g  ⁣dν=∫gf  ⁣dμ\int g\,\dd\nu = \int gf\,\dd\mu: for g=1Ag = \mathbf 1_A it is the definition of ν\nu; for simple gg, linearity; for g≥0g \geq 0 measurable, take simple sn↗gs_n \nearrow g (Theorem 10.4): snf↗gfs_nf \nearrow gf, and MCT on both sides passes to the limit. (This “indicator →\to simple →\to MCT” escalator is the standard machine of the theory.)

Exercise 10.10 ★★★

(A Weierstrass-style failure) Define f(t)=∫0+∞sin⁡(tx)x(1+x2)  ⁣dxf(t) = \int_0^{+\infty}\frac{\sin(tx)}{x(1 + x^2)}\,\dd x. (a) Show that ff is well defined and continuous on R\R, and C1\mathcal C^1 with f′(t)=∫0∞cos⁡(tx)1+x2 ⁣dxf'(t) = \int_0^\infty\frac{\cos(tx)}{1 + x^2}\dd x for every tt — but that differentiating again under the integral is illegitimate. (b) Admitting f′(t)=π2e−tf'(t) = \frac\pi2\eu^{-t} for t>0t > 0 (proved in Chapter 17), what is ∫0∞xsin⁡(tx)1+x2 ⁣dx\int_0^\infty\frac{x\sin(tx)}{1+x^2}\dd x for t>0t > 0, and why does its formula confirm the failure in (a)?

Solution

Solution of Exercise 10.10.

(a) ∣sin⁡(tx)∣≤∣t∣x\abs{\sin(tx)} \leq \abs tx gives ∣sin⁡(tx)x(1+x2)∣≤∣t∣1+x2\abs{\frac{\sin(tx)}{x(1+x^2)}} \leq \frac{\abs t}{1+x^2}: the integral converges, and on ∣t∣≤T\abs t \leq T the dominator T1+x2\frac{T}{1+x^2} yields continuity (Theorem 10.14). Differentiation: ∣∂t∣=∣cos⁡(tx)1+x2∣≤11+x2\abs{\partial_t} = \abs{\frac{\cos(tx)}{1+x^2}} \leq \frac1{1+x^2}, integrable: f′(t)=∫0∞cos⁡(tx)1+x2 ⁣dxf'(t) = \int_0^\infty\frac{\cos(tx)}{1+x^2}\dd x for all tt. A second differentiation would require integrating xsin⁡(tx)1+x2\frac{x\sin(tx)}{1 + x^2}, whose absolute value behaves like ∣sin⁡(tx)∣x\frac{\abs{\sin(tx)}} x at infinity: not integrable — no dominator exists and Theorem 10.15 cannot be applied again.

(b) Admitting f′(t)=π2e−tf'(t) = \frac\pi2\eu^{-t} for t>0t > 0: by oddness of ff, f′f' is even, so f′(t)=π2e−∣t∣f'(t) = \frac\pi2\eu^{-\abs t} — which is not differentiable at 00: ff is C1\mathcal C^1 but not C2\mathcal C^2, confirming that the blocked second differentiation was not a technical accident. For t>0t > 0 the improper integral ∫0∞xsin⁡(tx)1+x2 ⁣dx\int_0^\infty\frac{x\sin(tx)}{1+x^2}\dd x equals −f′′(t)=π2e−t-f''(t) = \frac\pi2\eu^{-t} (differentiating the admitted formula where it is legitimate, i.e. on (0, ∞)\intoo0\infty) — an improper, non-Lebesgue value.

Exercise 10.11 ★★

(Scheffé’s lemma) Let fn,f≥0f_n, f \geq 0 be integrable with fn→ff_n \to f a.e. and ∫fn→∫f\int f_n \to \int f. (a) Show that ∫∣fn−f∣→0\int\abs{f_n - f} \to 0. (Apply dominated convergence to gn=(f−fn)+≤fg_n = (f - f_n)^+ \leq f, and write ∫∣fn−f∣=2∫gn−∫(f−fn)\int\abs{f_n - f} = 2\int g_n - \int(f - f_n).) (b) Show by example that the hypothesis ∫fn→∫f\int f_n \to \int f cannot be dropped (a sliding or concentrating bump), and that the conclusion fails for signed fnf_n without absolute-value control: fn=n1(0, 1/n]−n1(−1/n, 0]f_n = n\mathbf 1_{\intoc0{1/n}} - n\mathbf 1_{\intoc{-1/n}0} has fn→0f_n \to 0 a.e., ∫fn=0→0\int f_n = 0 \to 0, yet ∫∣fn∣=2\int\abs{f_n} = 2. (c) Application (densities): if probability densities pn→pp_n \to p a.e., then automatically ∫∣pn−p∣→0\int\abs{p_n - p} \to 0: pointwise convergence of densities is L1L^1 convergence — a convergence upgrade for free.

Solution

Solution of Exercise 10.11.

(a) Let gn=(f−fn)+g_n = (f - f_n)^+: then 0≤gn≤f0 \leq g_n \leq f (positivity of fnf_n), gn→0g_n \to 0 a.e., and ff is an integrable dominator: ∫gn→0\int g_n \to 0 (DCT). Since ∣fn−f∣=2(f−fn)+−(f−fn)\abs{f_n - f} = 2(f - f_n)^+ - (f - f_n),

∫∣fn−f∣=2∫gn−(∫f−∫fn)⟶0+0.\int\abs{f_n - f} = 2\int g_n - \Bigl(\int f - \int f_n\Bigr) \longrightarrow 0 + 0 .

(b) The sliding bump fn=1[n, n+1]f_n = \mathbf 1_{\intcc n{n+1}} has fn→0f_n \to 0 a.e. and ∫fn=1↛0\int f_n = 1 \not\to 0: without the convergence of integrals, L1L^1 convergence fails (and so does the hypothesis). The signed example: fn→0f_n \to 0 at every x≠0x \neq 0, ∫fn=0\int f_n = 0, but ∫∣fn∣=2\int\abs{f_n} = 2: for signed sequences the theorem is genuinely about ∣fn∣\abs{f_n}-type control, and positivity was used exactly in gn≤fg_n \leq f.

(c) Densities satisfy ∫pn=1=∫p\int p_n = 1 = \int p: the hypothesis of (a) is automatic, so pn→pp_n \to p a.e. forces ∥pn−p∥L1→0\norm{p_n - p}_{L^1} \to 0 — and hence convergence of the probabilities ∫Apn→∫Ap\int_Ap_n \to \int_Ap uniformly over all measurable AA (∣∫A(pn−p)∣≤∥pn−p∥1\abs{\int_A(p_n - p)} \leq \norm{p_n - p}_1): Scheffé turns pointwise convergence of densities into total-variation convergence of laws.

Exercise 10.12 ★★

Classical limits, with full justification via MCT/DCT:

(a) lim⁡n→∞∫0n(1−xn)nex/2  ⁣dx,(b) lim⁡n→∞∫01n xn−11+x  ⁣dx,\text{(a)}\ \lim_{n\to\infty}\int_0^n\Bigl(1 - \frac xn\Bigr)^n\eu^{x/2}\,\dd x, \qquad \text{(b)}\ \lim_{n\to\infty}\int_0^1\frac{n\,x^{n-1}}{1 + x}\,\dd x,
(c) lim⁡n→∞∫0∞ ⁣dx(1+x/n)n x1/n.\text{(c)}\ \lim_{n\to\infty}\int_0^\infty \frac{\dd x}{(1 + x/n)^n\,x^{1/n}} .

(For (a): (1−x/n)n↗e−x(1 - x/n)^n \nearrow \eu^{-x} for fixed xx — prove the monotonicity via log⁡\log; for (b), integrate by parts or substitute x=u1/nx = u^{1/n} and identify a boundary concentration; for (c), find an integrable dominator valid for all n≥2n \geq 2 by splitting at x=1x = 1.)

Solution

Solution of Exercise 10.12.

(a) On (0, n)\intoo0n, φn(x)=nlog⁡(1−xn)\varphi_n(x) = n\log(1 - \frac xn) increases in nn to −x-x (the map t↦log⁡(1−xt)tt \mapsto \frac{\log(1 - xt)}{t} decreases as t=1n↓0t = \frac1n \downarrow 0; or expand: φn′≥0\varphi_{n}' \geq 0 in nn via log⁡(1−u)+u1−u≥0\log(1-u) + \frac{u}{1-u} \geq 0). So (1−xn)nex/21x<n↗e−x/2(1 - \frac xn)^n\eu^{x/2}\mathbf 1_{x<n} \nearrow \eu^{-x/2}, and MCT gives

lim⁡n∫0n(1−xn)nex/2 ⁣dx=∫0∞e−x/2 ⁣dx=2.\lim_n\int_0^n\Bigl(1 - \frac xn\Bigr)^n\eu^{x/2}\dd x = \int_0^\infty\eu^{-x/2}\dd x = 2 .

(b) Substitute u=xnu = x^n (x=u1/nx = u^{1/n}, nxn−1 ⁣dx= ⁣dun x^{n-1}\dd x = \dd u):

∫01nxn−11+x ⁣dx=∫01 ⁣du1+u1/n.\int_0^1\frac{nx^{n-1}}{1 + x}\dd x = \int_0^1\frac{\dd u}{1 + u^{1/n}} .

For u∈(0, 1)u \in \intoo01: u1/n→1u^{1/n} \to 1, so the integrand tends to 12\frac12, dominated by 11: the limit is 12\frac12 (DCT). (The mass of nxn−1nx^{n-1} concentrates at x=1x = 1, where 11+x=12\frac1{1+x} = \frac12: the substitution makes the concentration visible.)

(c) Pointwise, (1+x/n)n↗ex(1 + x/n)^n \nearrow \eu^x and x1/n→1x^{1/n} \to 1 (x>0x > 0): the integrand tends to e−x\eu^{-x}. Dominator for n≥2n \geq 2: on (0, 1]\intoc01, x−1/n≤x−1/2x^{-1/n} \leq x^{-1/2} and (1+x/n)−n≤1(1 + x/n)^{-n} \leq 1: bound x−1/2x^{-1/2}, integrable; on (1, ∞)\intoo1\infty, x−1/n≤1x^{-1/n} \leq 1 and (1+x/n)n≥1+(n2)x2n2≥1+x24(1 + x/n)^n \geq 1 + \binom n2\frac{x^2}{n^2} \geq 1 + \frac{x^2}4: bound 44+x2\frac{4}{4 + x^2}, integrable. DCT:

lim⁡n∫0∞ ⁣dx(1+x/n)nx1/n=∫0∞e−x ⁣dx=1.\lim_n\int_0^\infty\frac{\dd x}{(1 + x/n)^nx^{1/n}} = \int_0^\infty\eu^{-x}\dd x = 1 .

10.6 Problem: two celebrated integrals

Problem 10.1

Weekend problem — the Gaussian integral and Dirichlet’s integral, by parameters alone

Two integrals rule applied analysis:

G=∫−∞+∞e−x2 ⁣dx=π,D=∫0+∞sin⁡xx  ⁣dx=π2G = \int_{-\infty}^{+\infty}\eu^{-x^2}\dd x = \sqrt\pi, \qquad D = \int_0^{+\infty}\frac{\sin x}{x}\,\dd x = \frac\pi2

(the second as an improper integral, Example 10.13). We prove both using only this chapter’s tools.

Part I — The Gaussian. For t≥0t \geq 0 set

A(t)=(∫0te−x2 ⁣dx)2,B(t)=∫01e−t2(1+x2)1+x2  ⁣dx.A(t) = \Bigl(\int_0^t\eu^{-x^2}\dd x\Bigr)^{2}, \qquad B(t) = \int_0^1\frac{\eu^{-t^2(1 + x^2)}}{1 + x^2}\,\dd x .
  1. Justify that AA and BB are C1\mathcal C^1 on (0, +∞)\intoo0{+\infty} and compute A′A' and B′B'; show A′(t)+B′(t)=0A'(t) + B'(t) = 0. (In B′B', substitute u=txu = tx.)
  2. Compute A(0)+B(0)A(0) + B(0) and lim⁡t→+∞(A+B)(t)\lim_{t\to+\infty}(A + B)(t) — justify the limit under the integral in BB.
  3. Conclude ∫0∞e−x2 ⁣dx=π2\int_0^\infty \eu^{-x^2}\dd x = \frac{\sqrt\pi}2, hence G=πG = \sqrt\pi, and deduce Γ(12)=π\Gamma(\tfrac12) = \sqrt\pi (substitute x=u2x = u^2 in Γ(12)\Gamma(\frac12)).

Part II — Dirichlet’s integral. For t≥0t \geq 0 set

F(t)=∫0+∞e−tx sin⁡xx  ⁣dx.F(t) = \int_0^{+\infty}\eu^{-tx}\,\frac{\sin x}{x}\,\dd x .
  1. Show that the integral defining F(t)F(t) converges for every t>0t > 0 as a Lebesgue integral, and for t=0t = 0 as an improper integral; show that D=lim⁡A→∞∫0Asin⁡xx ⁣dxD = \lim_{A\to\infty}\int_0^A\frac{\sin x}x\dd x exists (integrate by parts on [π,A][\pi, A]).
  2. Show that FF is C1\mathcal C^1 on (0, +∞)\intoo0{+\infty} with

    F′(t)=−∫0+∞e−txsin⁡x  ⁣dx=−11+t2F'(t) = -\int_0^{+\infty}\eu^{-tx}\sin x\,\dd x = -\frac{1}{1 + t^2}

    (domination on [t0,∞)[t_0, \infty) for each t0>0t_0 > 0; the last integral by two integrations by parts or complex exponentials).

  3. Show F(t)→0F(t) \to 0 as t→+∞t \to +\infty, and deduce F(t)=π2−arctan⁡tF(t) = \frac\pi2 - \arctan t on (0, +∞)\intoo0{+\infty}.
  4. The delicate point: D=lim⁡t→0+F(t)D = \lim_{t\to0^+}F(t). Prove it by uniform control of the tail: for 0≤t≤10 \leq t \leq 1 and A≥πA \geq \pi, integrate by parts to show

    ∣∫A+∞e−txsin⁡xx  ⁣dx∣≤CA\Bigl|\int_A^{+\infty}\eu^{-tx}\frac{\sin x}x\,\dd x\Bigr| \leq \frac{C}{A}

    with CC independent of tt (differentiate e−txx\frac{\eu^{-tx}}x and bound ∣cos⁡∣\abs{\cos} by 11; note te−tx≤1/x⋅(txe−tx)t\eu^{-tx} \leq 1/x\cdot(tx\eu^{-tx}) with sup⁡u≥0ue−u<1\sup_{u\geq0}u\eu^{-u} < 1); then split F(t)−DF(t) - D into [0,A][0, A] (where DCT applies as t→0t \to 0) and [A,∞)[A, \infty).

  5. Conclude: D=π2D = \frac\pi2.

Part III — Dividends.

  1. Compute ∫0+∞sin⁡2xx2  ⁣dx\int_0^{+\infty}\frac{\sin^2x}{x^2}\,\dd x (integrate by parts and reduce to DD via sin⁡2x=2sin⁡xcos⁡x\sin 2x = 2\sin x\cos x).
  2. Compute ∫0+∞1−cos⁡xx2  ⁣dx\int_0^{+\infty}\frac{1 - \cos x}{x^2}\,\dd x, and check the consistency of the two results.
  3. For a>0a > 0, compute ∫0+∞sin⁡(ax)x ⁣dx\int_0^{+\infty}\frac{\sin(ax)}x\dd x and ∫−∞+∞e−ax2 ⁣dx\int_{-\infty}^{+\infty}\eu^{-ax^2}\dd x, and record the scaling rules (they will be the workhorses of Chapter 14).
  4. Explain precisely why DD could not have been treated by DCT directly at t=0t = 0 (no integrable dominator on [0,1]×[0,∞)[0,1]\times[0,\infty)), and why the tail-splitting of question 7 is the honest substitute — this pattern (“uniform integrability of tails”) recurs throughout analysis.

Part IV — The Gamma function according to Bohr and Mollerup. The function Γ\Gamma (Example 10.16) satisfies Γ(1)=1\Gamma(1) = 1 and Γ(x+1)=xΓ(x)\Gamma(x+1) = x\Gamma(x) — but so do infinitely many other functions (multiply by any 11-periodic wobble). One convexity condition pins Γ\Gamma down uniquely, and its deeper identities then fall out mechanically. A positive function ff on an interval is log-convex if log⁡f\log f is convex.

  1. Show that log-convex implies convex, that products of log-convex functions and their compositions with affine maps are log-convex, and — via the two-function Hölder inequality ∫∣uv∣≤(∫∣u∣p)1/p(∫∣v∣q)1/q\int\abs{uv} \leq \bigl(\int\abs u^p\bigr)^{1/p}\bigl(\int\abs v^q\bigr)^{1/q}, proved directly from Young’s inequality — that Γ\Gamma is log-convex on (0, ∞)\intoo0\infty.
  2. (Slope lemma) Let gg be convex on (0, ∞)\intoo0\infty with g(n+1)−g(n)=log⁡ng(n+1) - g(n) = \log n for every integer n≥1n \geq 1. For x∈(0, 1]x \in \intoc01 and n≥2n \geq 2, compare the slopes of gg over [n−1,n][n-1, n], [n,n+x][n, n+x] and [n,n+1][n, n+1], and deduce

    xlog⁡(n−1)  ≤  g(n+x)−g(n)  ≤  xlog⁡n.x\log(n-1) \;\leq\; g(n + x) - g(n) \;\leq\; x\log n .
  3. (Bohr–Mollerup) Let f>0f > 0 satisfy f(1)=1f(1) = 1, f(x+1)=xf(x)f(x+1) = xf(x), and log⁡f\log f convex. Unwinding the recursion into f(n+x)=x(x+1)⋯(x+n−1) f(x)f(n + x) = x(x+1)\cdots(x + n - 1)\,f(x) and f(n)=(n−1)!f(n) = (n-1)!, deduce from question 14 that for x∈(0, 1]x \in \intoc01

    f(x)=lim⁡n→∞n! nxx(x+1)⋯(x+n):f(x) = \lim_{n\to\infty} \frac{n!\,n^x}{x(x+1)\cdots(x+n)} :

    ff is unique, hence f=Γf = \Gamma, and Gauss’s limit formula holds (extend to all x>0x > 0 by the recursion).

  4. Define the Beta function B(x,y)=∫01tx−1(1−t)y−1  ⁣dtB(x, y) = \int_0^1t^{x-1}(1-t)^{y-1}\,\dd t (x,y>0x, y > 0). Prove convergence, the recursion B(x+1,y)=xx+y B(x,y)B(x+1, y) = \frac{x}{x+y}\,B(x, y) (integrate by parts), and B(1,y)=1yB(1, y) = \frac1y.
  5. Show that x↦B(x,y)x \mapsto B(x, y) is log-convex (Hölder again), and apply Bohr–Mollerup to

    f(x)=B(x,y) Γ(x+y)Γ(y)f(x) = \frac{B(x, y)\,\Gamma(x + y)}{\Gamma(y)}

    to conclude Euler’s formula: B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x, y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} — no double integrals anywhere.

  6. Compute B(12,12)B(\frac12, \frac12) directly (substitute t=sin⁡2θt = \sin^2\theta) and deduce Γ(12)=π\Gamma(\frac12) = \sqrt\pi: the Gaussian integral of Part I, recovered by pure convexity. Compare the two proofs in one sentence each.
  7. (Legendre duplication) Show that

    g(x)=2x−1π Γ(x2)Γ(x+12)g(x) = \frac{2^{x-1}}{\sqrt\pi}\, \Gamma\Bigl(\frac x2\Bigr) \Gamma\Bigl(\frac{x+1}2\Bigr)

    satisfies the three Bohr–Mollerup hypotheses, and conclude g=Γg = \Gamma, i.e. Γ(2z)=22z−1π Γ(z) Γ(z+12)\Gamma(2z) = \frac{2^{2z-1}}{\sqrt\pi}\,\Gamma(z)\,\Gamma(z + \tfrac12) for all z>0z > 0.

  8. Deduce the closed form Γ(n+12)=(2n)!4n n!π\Gamma\bigl(n + \tfrac12\bigr) = \dfrac{(2n)!}{4^n\,n!}\sqrt\pi, and prove, by the slope lemma applied to log⁡Γ\log\Gamma around large integers, the asymptotics

    Γ(n+12)Γ(n) n⟶1.\frac{\Gamma(n + \frac12)}{\Gamma(n)\,\sqrt n} \longrightarrow 1 .
  9. Combine the last two questions into the central binomial asymptotics

    (2nn)∼4nπn,\binom{2n}{n} \sim \frac{4^n}{\sqrt{\pi n}},

    and verify numerically for n=10n = 10 ((2010)=184756\binom{20}{10} = 184756, against 410/10π≈1870794^{10}/\sqrt{10\pi} \approx 187079: ratio ≈0.988\approx 0.988).

  10. (Synthesis) The constant π\sqrt\pi has now appeared as the Gaussian integral (Part I), as B(12,12)B(\frac12, \frac12) (question 18), and inside duplication (question 19); the central binomial estimate anticipates both Stirling (Chapter 11’s weekend problem) and de Moivre–Laplace. Map the connections: which statements are equivalent to which, and what does each technique — differentiation under the integral versus convexity — contribute that the other cannot?

Part V — Three more dividends.

  1. (Wallis, by Beta) For p>−1p > -1, substitute t=sin⁡2θt = \sin^2\theta to show

    Wp=∫0π/2sin⁡pθ  ⁣dθ=12 B(p+12,12),W_p = \int_0^{\pi/2}\sin^p\theta\,\dd\theta = \frac12\,B\Bigl(\frac{p+1}2, \frac12\Bigr),

    and deduce from the Beta recursion (question 16) that Wn+2=n+1n+2 WnW_{n+2} = \frac{n+1}{n+2}\,W_n for integers n≥0n \geq 0. Compute W2nW_{2n} and W2n+1W_{2n+1} in closed form, show W2n+1/W2n→1W_{2n+1}/W_{2n} \to 1 by squeezing, and conclude with Wallis’s product

    π2=lim⁡n→∞∏k=1n4k24k2−1.\frac\pi2 = \lim_{n\to\infty} \prod_{k=1}^{n}\frac{4k^2}{4k^2 - 1} .
  2. (The Gaussian meets a frequency) For b∈Rb \in \R set

    Φ(b)=∫−∞+∞e−x2cos⁡(2bx)  ⁣dx.\Phi(b) = \int_{-\infty}^{+\infty} \eu^{-x^2}\cos(2bx)\,\dd x .

    Show that Φ\Phi is C1\mathcal C^1 on R\R, that an integration by parts yields the differential equation Φ′(b)=−2b Φ(b)\Phi'(b) = -2b\,\Phi(b), and conclude

    Φ(b)=π e−b2:\Phi(b) = \sqrt\pi\,\eu^{-b^2} :

    the Gaussian reproduces itself under this transform — the single identity on which Chapter 14 will run.

  3. (Frullani’s integral) For 0<a<b0 < a < b, show that

    ∫0+∞e−ax−e−bxx  ⁣dx=log⁡ba,\int_0^{+\infty} \frac{\eu^{-ax} - \eu^{-bx}}{x}\,\dd x = \log\frac ba ,

    by differentiating in the parameter aa (justify the domination on every [a0, +∞)\intco{a_0}{+\infty}, a0>0a_0 > 0, and identify the constant by letting a→ba \to b). Where exactly does the integrand need its removable singularity at x=0x = 0?

Solution

Solution of Problem 10.1.

1. AA is C1\mathcal C^1 by the fundamental theorem of calculus and the chain rule: A′(t)=2e−t2∫0te−x2 ⁣dxA'(t) = 2\eu^{-t^2}\int_0^t\eu^{-x^2}\dd x. For BB: ∂t(e−t2(1+x2)1+x2)=−2t e−t2(1+x2)\partial_t\bigl(\frac{\eu^{-t^2(1+x^2)}}{1+x^2}\bigr) = -2t\,\eu^{-t^2(1+x^2)}, continuous and bounded on [t0,T]×[0, 1][t_0, T] \times \intcc01 for any 0<t0<T0 < t_0 < T (bounded domination on a finite measure space suffices): BB is C1\mathcal C^1 on (0, +∞)\intoo0{+\infty} with

B′(t)=−2t∫01e−t2(1+x2) ⁣dx=−2e−t2∫01t e−(tx)2 ⁣dx=−2e−t2∫0te−u2 ⁣du=−A′(t)B'(t) = -2t\int_0^1 \eu^{-t^2(1+x^2)}\dd x = -2\eu^{-t^2}\int_0^1 t\,\eu^{-(tx)^2}\dd x = -2\eu^{-t^2}\int_0^t\eu^{-u^2}\dd u = -A'(t)

(substitution u=txu = tx).

2. A(0)=0A(0) = 0 and B(0)=∫01 ⁣dx1+x2=π4B(0) = \int_0^1\frac{\dd x}{1+x^2} = \frac\pi4. A+BA + B has zero derivative on (0, +∞)\intoo0{+\infty} and is continuous at 00 (BB by domination 11+x2\frac1{1+x^2} and Theorem 10.14): A+B≡π4A + B \equiv \frac\pi4. As t→+∞t\to+\infty: 0≤B(t)≤e−t2→00 \leq B(t) \leq \eu^{-t^2} \to 0, and A(t)→(∫0∞e−x2 ⁣dx)2A(t) \to \bigl(\int_0^\infty\eu^{-x^2}\dd x\bigr)^2 (MCT or just monotone convergence of the inner integral).

3. Hence (∫0∞e−x2 ⁣dx)2=π4\bigl(\int_0^\infty\eu^{-x^2}\dd x\bigr)^2 = \frac\pi4: ∫0∞e−x2 ⁣dx=π2\int_0^\infty\eu^{-x^2}\dd x = \frac{\sqrt\pi}2, and by evenness G=πG = \sqrt\pi. Also Γ(12)=∫0∞x−1/2e−x ⁣dx=x=u22∫0∞e−u2 ⁣du=π\Gamma(\frac12) = \int_0^\infty x^{-1/2}\eu^{-x}\dd x \overset{x = u^2}{=} 2\int_0^\infty\eu^{-u^2}\dd u = \sqrt\pi.

4. For t>0t > 0: ∣e−txsin⁡xx∣≤e−tx\abs{\eu^{-tx}\frac{\sin x}x} \leq \eu^{-tx}, integrable. For t=0t = 0, improper convergence: on [π,A][\pi, A],

∫πAsin⁡xx ⁣dx=[−cos⁡xx]πA−∫πAcos⁡xx2 ⁣dx,\int_\pi^A\frac{\sin x}x\dd x = \Bigl[-\frac{\cos x}x\Bigr]_\pi^A - \int_\pi^A\frac{\cos x}{x^2}\dd x,

both terms convergent as A→∞A \to \infty; near 00 the integrand extends continuously by 11.

5. On [t0,+∞)[t_0, +\infty) (t0>0t_0 > 0): ∣∂t(e−txsin⁡xx)∣=e−tx∣sin⁡x∣≤e−t0x\abs{\partial_t\bigl(\eu^{-tx}\tfrac{\sin x}x\bigr)} = \eu^{-tx}\abs{\sin x} \leq \eu^{-t_0x}, integrable: Theorem 10.15 applies on every such interval, so on all of (0, +∞)\intoo0{+\infty}:

F′(t)=−∫0∞e−txsin⁡x  ⁣dx=−Im⁡∫0∞e−(t−i)x ⁣dx=−Im⁡1t−i=−11+t2.F'(t) = -\int_0^\infty\eu^{-tx}\sin x\,\dd x = -\operatorname{Im}\int_0^\infty\eu^{-(t - \iu)x}\dd x = -\operatorname{Im}\frac{1}{t - \iu} = -\frac{1}{1 + t^2}.

6. ∣F(t)∣≤∫0∞e−tx ⁣dx=1t→0\abs{F(t)} \leq \int_0^\infty\eu^{-tx}\dd x = \frac1t \to 0. Integrating F′=−11+t2F' = -\frac1{1+t^2}: F(t)=C−arctan⁡tF(t) = C - \arctan t, and t→∞t \to \infty forces C=π2C = \frac\pi2: F(t)=π2−arctan⁡tF(t) = \frac\pi2 - \arctan t on (0, +∞)\intoo0{+\infty}.

7. Integrate by parts on [A,R][A, R] with sin⁡x=(−cos⁡x)′\sin x = (-\cos x)' and let R→∞R \to \infty:

∫A∞e−txsin⁡xx ⁣dx=cos⁡A  e−tAA−∫A∞cos⁡x (tx+1x2)e−tx ⁣dx.\int_A^{\infty}\eu^{-tx}\frac{\sin x}x\dd x = \frac{\cos A\;\eu^{-tA}}{A} - \int_A^\infty \cos x\,\Bigl(\frac tx + \frac1{x^2}\Bigr)\eu^{-tx}\dd x .

Bounding ∣cos⁡∣≤1\abs{\cos} \leq 1: the first term is ≤1A\leq \frac1A; the integral is at most ∫A∞te−tx ⁣dxx+∫A∞ ⁣dxx2≤1A∫0∞te−tx ⁣dx+1A=2A\int_A^\infty t\eu^{-tx}\frac{\dd x}x + \int_A^\infty\frac{\dd x}{x^2} \leq \frac1A\int_0^\infty t\eu^{-tx}\dd x + \frac1A = \frac2A. Total: ≤3A\leq \frac 3A, uniformly for t∈[0,1]t \in [0, 1] (the case t=0t = 0 included). Now

∣F(t)−D∣≤∣∫0A(e−tx−1) sin⁡xx  ⁣dx∣+6A.\abs{F(t) - D} \leq \Bigl|\int_0^A(\eu^{-tx} - 1)\,\frac{\sin x}x\,\dd x\Bigr| + \frac6A .

On [0,A][0, A]: ∣e−tx−1∣≤tx≤tA\abs{\eu^{-tx} - 1} \leq tx \leq tA and ∣sin⁡xx∣≤1\abs{\frac{\sin x}x} \leq 1, so the first term is at most tA2tA^2. Choose AA with 6A<ε\frac6A < \varepsilon, then t<ε/A2t < \varepsilon/A^2: ∣F(t)−D∣<2ε\abs{F(t) - D} < 2\varepsilon.

8. Therefore D=lim⁡t→0+F(t)=lim⁡t→0+(π2−arctan⁡t)=π2D = \lim_{t\to0^+}F(t) = \lim_{t\to0^+}\bigl(\frac\pi2 - \arctan t\bigr) = \frac\pi2.

9. By parts (u=sin⁡2xu = \sin^2x, v′=x−2v' = x^{-2}):

∫0∞sin⁡2xx2 ⁣dx=[−sin⁡2xx]0∞+∫0∞2sin⁡xcos⁡xx ⁣dx=∫0∞sin⁡2xx ⁣dx=D=π2\int_0^\infty\frac{\sin^2x}{x^2}\dd x = \Bigl[-\frac{\sin^2x}{x}\Bigr]_0^\infty + \int_0^\infty\frac{2\sin x\cos x}{x}\dd x = \int_0^\infty\frac{\sin 2x}{x}\dd x = D = \frac\pi2

(substitute u=2xu = 2x in the last step; boundary terms vanish: sin⁡2x/x→0\sin^2 x/x \to 0 at both ends).

10. By parts (u=1−cos⁡xu = 1 - \cos x, v′=x−2v' = x^{-2}): ∫0∞1−cos⁡xx2 ⁣dx=∫0∞sin⁡xx ⁣dx=π2\int_0^\infty\frac{1 - \cos x}{x^2}\dd x = \int_0^\infty\frac{\sin x}x\dd x = \frac\pi2. Consistency: 1−cos⁡x=2sin⁡2x21 - \cos x = 2\sin^2\frac x2, and the substitution x=2ux = 2u turns ∫2sin⁡2(x/2)x2 ⁣dx\int\frac{2\sin^2(x/2)}{x^2}\dd x into ∫sin⁡2uu2 ⁣du\int\frac{\sin^2u}{u^2}\dd u: the two computations agree.

11. ∫0∞sin⁡(ax)x ⁣dx=π2\int_0^\infty\frac{\sin(ax)}x\dd x = \frac\pi2 for every a>0a > 0 (substitute u=axu = ax: the integral is scale-invariant); ∫Re−ax2 ⁣dx=π/a\int_\R\eu^{-ax^2}\dd x = \sqrt{\pi/a} (substitute u=a xu = \sqrt a\,x). Scaling in the argument leaves the Dirichlet integral fixed and divides the Gaussian by a\sqrt a.

12. A dominator valid for all t∈[0,1]t \in [0,1] must dominate sup⁡t∈[0,1]∣e−txsin⁡xx∣=∣sin⁡xx∣\sup_{t\in[0,1]}\abs{\eu^{-tx}\frac{\sin x}x} = \abs{\frac{\sin x}x}, which is not integrable (Exercise 10.6): DCT cannot cross t=0t = 0. The substitute of question 7 — tails uniformly small in the parameter, compact part handled by DCT — is the standard “uniform integrability” pattern, and reappears whenever conditional convergence meets limit interchange.

13. If g=log⁡fg = \log f is convex then f=exp⁡∘gf = \exp\circ g is convex (exp is convex increasing: f(λx+(1−λ)y)≤eλg(x)+(1−λ)g(y)≤λf(x)+(1−λ)f(y)f(\lambda x + (1-\lambda) y) \leq \eu^{\lambda g(x) + (1-\lambda)g(y)} \leq \lambda f(x) + (1-\lambda)f(y), the last step by convexity of exp between the points g(x),g(y)g(x), g(y)). Products and affine substitutions: logarithms turn them into sums and affine substitutions of convex functions. Hölder (1p+1q=1\frac1p + \frac1q = 1): for ∫∣u∣p=∫∣v∣q=1\int\abs u^p = \int\abs v^q = 1, Young gives ∣uv∣≤∣u∣pp+∣v∣qq\abs{uv} \leq \frac{\abs u^p}p + \frac{\abs v^q}q, integrate: ∫∣uv∣≤1\int\abs{uv} \leq 1; the general case by homogeneity. Then, for λ∈(0, 1)\lambda \in \intoo01, apply it with p=1λp = \frac1\lambda to the factorization

tλx+(1−λ)y−1e−t=(tx−1e−t)λ(ty−1e−t)1−λ:Γ(λx+(1−λ)y)≤Γ(x)λ Γ(y)1−λ.t^{\lambda x + (1-\lambda)y - 1}\eu^{-t} = \bigl(t^{x-1}\eu^{-t}\bigr)^{\lambda} \bigl(t^{y-1}\eu^{-t}\bigr)^{1-\lambda} : \qquad \Gamma(\lambda x + (1{-}\lambda)y) \leq \Gamma(x)^\lambda\,\Gamma(y)^{1-\lambda} .

14. For a convex gg, the slope of a chord increases with its endpoints (three-chord inequality). Comparing the chords over [n−1,n][n-1, n], [n,n+x][n, n+x], [n,n+1][n, n+1]:

log⁡(n−1)=g(n)−g(n−1)1≤g(n+x)−g(n)x≤g(n+1)−g(n)1=log⁡n,\log(n-1) = \frac{g(n) - g(n-1)}1 \leq \frac{g(n+x) - g(n)}x \leq \frac{g(n+1) - g(n)}1 = \log n,

and multiplying by x>0x > 0 gives the claim.

15. With f(n)=(n−1)!f(n) = (n-1)! (recursion from f(1)=1f(1) = 1) and f(n+x)=x(x+1)⋯(x+n−1) f(x)f(n + x) = x(x+1)\cdots(x+n-1)\,f(x), question 14 reads

(n−1)x (n−1)!  ≤  x(x+1)⋯(x+n−1) f(x)  ≤  nx (n−1)!.(n-1)^x\,(n-1)! \;\leq\; x(x+1)\cdots(x+n-1)\,f(x) \;\leq\; n^x\,(n-1)! .

The upper bound rewrites as f(x)≤n! nxx(x+1)⋯(x+n)⋅x+nnf(x) \leq \frac{n!\,n^x}{x(x+1)\cdots(x+n)}\cdot\frac{x+n}n, and the lower bound at rank n+1n+1 as f(x)≥n! nxx(x+1)⋯(x+n)f(x) \geq \frac{n!\,n^x}{x(x+1)\cdots(x+n)}. The correction factor x+nn→1\frac{x+n}n \to 1: the sandwich forces

f(x)=lim⁡nn! nxx(x+1)⋯(x+n)(x∈(0, 1]),f(x) = \lim_n\frac{n!\,n^x}{x(x+1)\cdots(x+n)} \qquad (x \in \intoc01),

an expression independent of ff: uniqueness on (0, 1]\intoc01, hence everywhere by the recursion. Since Γ\Gamma satisfies all three hypotheses (question 13), f=Γf = \Gamma and Gauss’s formula holds — for all x>0x > 0, as both sides obey the same recursion.

16. Near 00, the integrand is ∼tx−1\sim t^{x-1}, integrable iff x>0x > 0; near 11, symmetric with yy. Integration by parts on [ε, 1−ε]\intcc\varepsilon{1-\varepsilon}, letting ε→0\varepsilon \to 0 (boundary terms vanish for x,y>0x, y > 0):

B(x+1,y)=[−tx(1−t)yy]01+xy∫01tx−1(1−t)y  ⁣dt=xy(B(x,y)−B(x+1,y)),B(x{+}1, y) = \Bigl[-t^x\frac{(1-t)^y}y\Bigr]_0^1 + \frac xy\int_0^1t^{x-1}(1-t)^y\,\dd t = \frac xy\bigl(B(x, y) - B(x{+}1, y)\bigr),

using (1−t)y=(1−t)y−1(1−t)(1-t)^y = (1-t)^{y-1}(1 - t); solving, B(x+1,y)=xx+yB(x,y)B(x+1, y) = \frac{x}{x+y}B(x, y). And B(1,y)=∫01(1−t)y−1 ⁣dt=1yB(1, y) = \int_0^1(1-t)^{y-1}\dd t = \frac1y.

17. f(1)=1y⋅Γ(1+y)Γ(y)=1f(1) = \frac1y\cdot\frac{\Gamma(1+y)}{\Gamma(y)} = 1; f(x+1)=xx+yB(x,y)⋅(x+y)Γ(x+y)Γ(y)=x f(x)f(x+1) = \frac{x}{x+y}B(x,y)\cdot\frac{(x+y)\Gamma(x+y)} {\Gamma(y)} = x\,f(x); and ff is log-convex in xx as a product of the log-convex B(⋅,y)B(\cdot, y) (Hölder on the factorization t(λx1+(1−λ)x2)−1(1−t)y−1=(⋯ )λ(⋯ )1−λt^{(\lambda x_1 + (1-\lambda)x_2)-1}(1-t)^{y-1} = (\cdots)^\lambda(\cdots)^{1-\lambda}, as in question 13) and Γ(⋅+y)\Gamma(\cdot + y) (affine shift). Bohr–Mollerup: f=Γf = \Gamma, i.e. B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x, y) = \frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}.

18. With t=sin⁡2θt = \sin^2\theta,  ⁣dt=2sin⁡θcos⁡θ  ⁣dθ\dd t = 2\sin\theta\cos\theta\,\dd\theta and t−1/2(1−t)−1/2=1sin⁡θcos⁡θt^{-1/2}(1-t)^{-1/2} = \frac1{\sin\theta\cos\theta}:

B(12,12)=∫0π/22  ⁣dθ=π=Γ(12)2Γ(1)⟹Γ(12)=π.B\Bigl(\frac12, \frac12\Bigr) = \int_0^{\pi/2}2\,\dd\theta = \pi = \frac{\Gamma(\frac12)^2}{\Gamma(1)} \quad\Longrightarrow\quad \Gamma\Bigl(\frac12\Bigr) = \sqrt\pi .

Part I reached the same constant by differentiating a parameter and racing two functions to their limits; here convexity alone rigidified the problem until only one value survived. Analysis by motion versus analysis by shape.

19. g(1)=20πΓ(12)Γ(1)=1g(1) = \frac{2^0}{\sqrt\pi}\Gamma(\frac12) \Gamma(1) = 1. Recursion:

g(x+1)=2xπ Γ(x+12)Γ(x2+1)=2xπ⋅x2 Γ(x2)Γ(x+12)=x g(x).g(x+1) = \frac{2^{x}}{\sqrt\pi}\, \Gamma\Bigl(\frac{x+1}2\Bigr)\Gamma\Bigl(\frac x2 + 1\Bigr) = \frac{2^{x}}{\sqrt\pi}\cdot\frac x2\, \Gamma\Bigl(\frac x2\Bigr)\Gamma\Bigl(\frac{x+1}2\Bigr) = x\,g(x) .

Log-convexity: product of e(x−1)log⁡2\eu^{(x-1)\log2} (log-affine) and two affine reparametrizations of the log-convex Γ\Gamma. Bohr–Mollerup gives g=Γg = \Gamma; setting x=2zx = 2z: Γ(2z)=22z−1πΓ(z)Γ(z+12)\Gamma(2z) = \frac{2^{2z-1}}{\sqrt\pi}\Gamma(z)\Gamma(z + \frac12).

20. From Γ(12)=π\Gamma(\frac12) = \sqrt\pi and the recursion, Γ(n+12)=(n−12)(n−32)⋯12 π=(2n−1)(2n−3)⋯12nπ=(2n)!4nn!π\Gamma(n + \frac12) = (n - \frac12)(n - \frac32)\cdots\frac12\,\sqrt\pi = \frac{(2n-1)(2n-3)\cdots1}{2^n}\sqrt\pi = \frac{(2n)!}{4^nn!}\sqrt\pi (complete the odd product with the evens). Asymptotics: question 14 with g=log⁡Γg = \log\Gamma and x=12x = \frac12 gives n−1≤Γ(n+12)Γ(n)≤n\sqrt{n-1} \leq \frac{\Gamma(n+\frac12)}{\Gamma(n)} \leq \sqrt n, so the ratio to n\sqrt n is squeezed between 1−1n\sqrt{1 - \frac1n} and 11.

21. From question 20, (2n)!=4nn!πΓ(n+12)(2n)! = \frac{4^nn!}{\sqrt\pi}\Gamma(n + \tfrac12), so

(2nn)=(2n)!(n!)2=4n Γ(n+12)π  n!=4nπ⋅Γ(n+12)n Γ(n)∼4nπ⋅nn=4nπn.\binom{2n}n = \frac{(2n)!}{(n!)^2} = \frac{4^n\,\Gamma(n+\frac12)}{\sqrt\pi\;n!} = \frac{4^n}{\sqrt\pi}\cdot \frac{\Gamma(n+\frac12)}{n\,\Gamma(n)} \sim \frac{4^n}{\sqrt\pi}\cdot\frac{\sqrt n}{n} = \frac{4^n}{\sqrt{\pi n}} .

Numerically, 410/10π=1048576/5.6050≈187078.64^{10}/\sqrt{10\pi} = 1048576/5.6050 \approx 187078.6, against (2010)=184756\binom{20}{10} = 184756: ratio 0.98760.9876 — the error is O(1/n)O(1/n), visible at n=10n = 10.

22. Equivalences: Γ(12)=π⇔G=π\Gamma(\frac12) = \sqrt\pi \Leftrightarrow G = \sqrt\pi (the substitution x=u2x = u^2 of question 3) ⇔B(12,12)=π\Leftrightarrow B(\frac12, \frac12) = \pi (Euler’s formula); duplication at z=nz = n is the closed form of Γ(n+12)\Gamma(n + \frac12), which is the central binomial estimate up to the slope lemma. The parameter technique (Part I–II) computes limits of moving quantities and is indispensable when a genuine deformation is present (Dirichlet’s integral has no convexity proof); the convexity technique computes nothing but forbids everything — it excels at uniqueness and functional equations (Gauss, Euler, Legendre in three strokes), where differentiation would drown in computation. A complete analyst carries both.

23. With t=sin⁡2θt = \sin^2\theta,  ⁣dt=2sin⁡θcos⁡θ  ⁣dθ=2 t1/2(1−t)1/2  ⁣dθ\dd t = 2\sin\theta\cos\theta\,\dd\theta = 2\,t^{1/2}(1-t)^{1/2}\,\dd\theta, so

Wp=∫01tp/2  ⁣dt2 t1/2(1−t)1/2=12∫01tp+12−1(1−t)12−1 ⁣dt=12 B(p+12,12).W_p = \int_0^1 t^{p/2}\, \frac{\dd t}{2\,t^{1/2}(1-t)^{1/2}} = \frac12\int_0^1 t^{\frac{p+1}2 - 1}(1-t)^{\frac12 - 1}\dd t = \frac12\,B\Bigl(\frac{p+1}2, \frac12\Bigr).

The Beta recursion with x=n+12x = \frac{n+1}2, y=12y = \frac12 gives

Wn+2=12 (n+1)/2(n+2)/2 B(n+12,12)=n+1n+2 Wn.W_{n+2} = \frac12\, \frac{(n+1)/2}{(n+2)/2}\,B\Bigl(\frac{n+1}2, \frac12\Bigr) = \frac{n+1}{n+2}\,W_n .

Starting from W0=π2W_0 = \frac\pi2, W1=1W_1 = 1:

W2n=π2∏k=1n2k−12k,W2n+1=∏k=1n2k2k+1.W_{2n} = \frac\pi2\prod_{k=1}^n\frac{2k-1}{2k}, \qquad W_{2n+1} = \prod_{k=1}^n\frac{2k}{2k+1} .

Since sin⁡n+1≤sin⁡n\sin^{n+1} \leq \sin^n on [0, π/2]\intcc0{\pi/2}, the sequence (Wn)(W_n) is nonincreasing, so

1≥W2n+1W2n≥W2n+1W2n−1=2n2n+1⟶1.1 \geq \frac{W_{2n+1}}{W_{2n}} \geq \frac{W_{2n+1}}{W_{2n-1}} = \frac{2n}{2n+1} \longrightarrow 1 .

But the closed forms give

W2n+1W2n=2π∏k=1n(2k)2(2k−1)(2k+1)=2π∏k=1n4k24k2−1,\frac{W_{2n+1}}{W_{2n}} = \frac2\pi \prod_{k=1}^n\frac{(2k)^2}{(2k-1)(2k+1)} = \frac2\pi\prod_{k=1}^n\frac{4k^2}{4k^2-1},

and letting n→∞n \to \infty yields Wallis’s product. (Via Euler’s formula, Wp=π2 Γ(p+12)/Γ(p2+1)W_p = \frac{\sqrt\pi}2\, \Gamma(\frac{p+1}2)/\Gamma(\frac p2 + 1): Wallis is the Gaussian integral in yet another costume.)

24. The bb-derivative of the integrand is −2x e−x2sin⁡(2bx)-2x\,\eu^{-x^2}\sin(2bx), dominated by 2∣x∣ e−x2∈L1(R)2\abs x\,\eu^{-x^2} \in L^1(\R) uniformly in bb: Φ\Phi is C1\mathcal C^1 with

Φ′(b)=−∫−∞+∞2x e−x2sin⁡(2bx)  ⁣dx.\Phi'(b) = -\int_{-\infty}^{+\infty} 2x\,\eu^{-x^2}\sin(2bx)\,\dd x .

Integrating by parts with u=sin⁡(2bx)u = \sin(2bx),  ⁣dv=−2x e−x2 ⁣dx\dd v = -2x\,\eu^{-x^2}\dd x (so v=e−x2v = \eu^{-x^2}), the boundary terms vanish and

Φ′(b)=−2b∫−∞+∞e−x2cos⁡(2bx)  ⁣dx=−2b Φ(b).\Phi'(b) = -2b\int_{-\infty}^{+\infty} \eu^{-x^2}\cos(2bx)\,\dd x = -2b\,\Phi(b).

Hence (Φ(b) eb2)′=0\bigl(\Phi(b)\,\eu^{b^2}\bigr)' = 0 and Φ(b)=Φ(0) e−b2=π e−b2\Phi(b) = \Phi(0)\,\eu^{-b^2} = \sqrt\pi\,\eu^{-b^2} by Part I. Up to normalization this says the Fourier transform of e−x2\eu^{-x^2} is again a Gaussian — the fixed point on which the inversion theory of Chapter 14 pivots.

25. For 0<a<b0 < a < b and x>0x > 0,

0≤e−ax−e−bxx=∫abe−sx  ⁣ds≤(b−a) e−ax,0 \leq \frac{\eu^{-ax} - \eu^{-bx}}{x} = \int_a^b \eu^{-sx}\,\dd s \leq (b - a)\,\eu^{-ax},

so the integral I(a)I(a) converges (Lebesgue); the pointwise bound also shows the integrand extends continuously by b−ab - a at x=0x = 0. Fix bb; on [a0, +∞)\intco{a_0}{+\infty} the aa-derivative of the integrand is −e−ax-\eu^{-ax}, dominated by e−a0x∈L1((0, +∞))\eu^{-a_0x} \in L^1(\intoo0{+\infty}), so II is C1\mathcal C^1 on (0, b)\intoo0b with

I′(a)=−∫0+∞e−ax ⁣dx=−1a,henceI(a)=−log⁡a+c.I'(a) = -\int_0^{+\infty}\eu^{-ax}\dd x = -\frac1a, \qquad\text{hence}\qquad I(a) = -\log a + c .

The two-sided bound gives 0≤I(a)≤(b−a)/a→00 \leq I(a) \leq (b-a)/a \to 0 as a→b−a \to b^-, so c=log⁡bc = \log b and I(a)=log⁡baI(a) = \log\frac ba. The removable singularity is needed at 00: each term e−ax/x\eu^{-ax}/x separately has a divergent (logarithmic) integral near 00, and only the first-order cancellation e−ax−e−bx=O(x)\eu^{-ax} - \eu^{-bx} = O(x) makes the difference integrable there; at infinity each term is harmless on its own.

Terms defined in this chapter

See all 395 terms in the glossary