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 0sinxx ⁣dx\int_0^\infty\frac{\sin x}x\,\dd x and Rex2 ⁣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 ⁣:XYf \colon X \to Y is measurable if f1(B)Af^{-1}(B) \in \mathcal A for every BBB \in \mathcal B. For real (or [,+][-\infty,+\infty]-valued) functions, Y=RY = \R carries its Borel σ\sigma-algebra, and it suffices to check f1((t,+))={f>t}Af^{-1}(\intoo t{+\infty}) = \{f > t\} \in \mathcal A for all tRt \in \R: the good sets {B:f1(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 ⁣:XRf, 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 supnfn\sup_nf_n, infnfn\inf_nf_n, lim supfn\limsup f_n, lim inffn\liminf f_n are measurable; if fnff_n \to f pointwise, ff is measurable.

Proof. (a) (gf)1(B)=f1(g1(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) ⁣:XR2(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)=f1(U)g1(V)(f,g)^{-1}(U\times V) = f^{-1}(U)\cap g^{-1}(V) — and +,×,max+, \times, \max are continuous R2R\R^2 \to \R: compose. (c) {supfn>t}=n{fn>t}\{\sup f_n > t\} = \bigcup_n\{f_n > t\}; inf=sup()\inf = -\sup(-); lim sup=infNsupnN\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=1nci1Ais = \sum_{i=1}^n c_i\,\mathbf 1_{A_i}, AiAA_i \in \mathcal A disjoint, ci0c_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=1n2nk12n1{k12nf<k2n}+n1{fn}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 k12nf(x)<k2n\frac{k-1}{2^n} \leq f(x) < \frac k{2^n} gets either 2k22n+1\frac{2k-2}{2^{n+1}} or 2k12n+1\frac{2k-1}{2^{n+1}}, both k12n\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)2nf(x) - s_n(x) \leq 2^{-n}; if f(x)=f(x) = \infty, sn(x)=ns_n(x) = n \to \infty.

10.2 The integral and the convergence theorems

Definition 10.5

For measurable f0f \geq 0:

f ⁣dμ=sup{s ⁣dμ:s simple, 0sf}[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 0fnf0 \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 LfL \leq \int f (monotonicity). Conversely, fix a simple s=ci1Aifs = \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

fnEnθs ⁣dμ=θiciμ(AiEn)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 sfs \leq f: LfL \geq \int f.

Corollary 10.7

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

Proof. For simple functions, additivity is a computation on a common refinement. In general take snfs_n \nearrow f, tngt_n \nearrow g (Theorem 10.4): sn+tnf+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 fn0f_n \geq 0:

lim infnfn ⁣dμ    lim infnfn ⁣dμ.\int \liminf_n f_n \,\dd\mu \;\leq\; \liminf_n \int f_n\,\dd\mu .

Proof. Let gN=infnNfng_N = \inf_{n\geq N}f_n: measurable, 0gNlim inffn0 \leq g_N \nearrow \liminf f_n, and gNfng_N \leq f_n for every nNn \geq N, so gNinfnNfn\int g_N \leq \inf_{n \geq N}\int f_n. Apply MCT to the left side: lim inffn=limNgNlimNinfnNfn=lim inffn\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 ⁣:XRf \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 ff\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 fnff_n \to f a.e., with fng\abs{f_n} \leq g a.e. for a fixed integrable gg. Then ff is integrable and

fn ⁣dμf ⁣dμ,indeedfnf ⁣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. fg\abs f \leq g: ff is integrable. The functions hn=2gfnf0h_n = 2g - \abs{f_n - f} \geq 0 satisfy lim infhn=2g\liminf h_n = 2g; Fatou gives

2glim inf(2gfnf)=2glim supfnf,\int 2g \leq \liminf\int\bigl(2g - \abs{f_n - f}\bigr) = \int 2g - \limsup\int\abs{f_n - f},

so lim supfnf0\limsup\int\abs{f_n - f} \leq 0 (the subtraction is legal: 2g<\int 2g < \infty). Finally fnffnf0\abs{\int f_n - \int f} \leq \int\abs{f_n - f} \to 0.

Method 10.11

Faced with limnfn\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 gfng \geq \abs{f_n}, found by crude bounds (“supn\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 fn0f_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 (xi1,xi)\intoo{x_{i-1}}{x_i}, to mi=inf[xi1,xi]fm_i = \inf_{[x_{i-1}, x_i]}f and Mi=supM_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 xDx \notin D, writing In(x)I_n(x) for the open σn\sigma_n-interval containing xx: (x)=supninfIn(x)f\ell(x) = \sup_n\inf_{I_n(x)}f and u(x)=infnsupIn(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 xxu(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λ=limnLσ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. (u0u - \ell \geq 0; Exercise 10.5)     \iff ff continuous a.e. In that case fu\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: 0sinxx ⁣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 sinxxL1((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×XCf \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: tf(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 tnt0t_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, tf(t,x)t \mapsto f(t,x) is differentiable on TT, with

ft(t,x)g(x)for all tT, 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)=Xft(t,x) ⁣dμ(x)F'(t) = \int_X \frac{\partial f}{\partial t}(t, x)\,\dd\mu(x).

Proof. Fix tt and hn0h_n \to 0: the difference quotients

φn(x)=f(t+hn,x)f(t,x)hnft(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)supstf(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=φntf(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+xt1ex ⁣dx.\Gamma(t) = \int_0^{+\infty} x^{t-1}\eu^{-x}\,\dd x .

The integral converges: near 00, xt1x^{t-1} is integrable (t>0t > 0); at infinity, xt1exCex/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(xt1ex)=lnxxt1ex\partial_t\bigl(x^{t-1}\eu^{-x}\bigr) = \ln x\cdot x^{t-1}\eu^{-x} is dominated by lnx(xa1+xb1)ex\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(lnx)kxt1ex ⁣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 sinxx\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: sinxxL1\frac{\sin x}x \notin L^1. Lebesgue integrability is absolute integrability.

10.5 Exercises

Exercise 10.1

(a) Show that a monotone function RR\R \to \R is Borel-measurable, and that a derivative (of an everywhere differentiable function) is Borel-measurable. (b) Show that f ⁣:XRf \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)=limnn(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}=qQ,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:

limn0+cosx(1+x/n)n ⁣dx,limn01nxn11+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)nex(1 + x/n)^n \nearrow \eu^x for x0x \geq 0, so the integrand tends pointwise to excosx\eu^{-x}\cos x; for n2n \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:

limn0cosx(1+x/n)n ⁣dx=0excosx ⁣dx=Re0e(1i)x ⁣dx=Re11i=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):

01nxn11+x ⁣dx=01 ⁣du1+u1/n01 ⁣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/n1u^{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 fn0f_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=n1(0,1/n)f_n = n\,\mathbf 1_{\intoo0{1/n}}: lim inffn=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: supnfn\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 supfnlim supfn\int\limsup f_n \geq \limsup\int f_n” fails for the translating bump: left side 00, right side 11. “lim suplim 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+xex1 ⁣dx=n11n2=π26\displaystyle\int_0^{+\infty}\frac{x}{\eu^x - 1}\,\dd x = \sum_{n\geq1}\frac1{n^2} = \frac{\pi^2}6 (expand 1ex1\frac1{\eu^x - 1} in a geometric series and integrate term by term — which theorem permits it?). (b) (Sophomore’s dream) Show 01xx ⁣dx=n1nn\displaystyle\int_0^1 x^{-x}\,\dd x = \sum_{n\geq1}n^{-n}. (Write xx=exlnx=k(xlnx)kk!x^{-x} = \eu^{-x\ln x} = \sum_k\frac{(-x\ln x)^k}{k!} and compute 01(xlnx)k ⁣dx\int_0^1(-x\ln x)^k\dd x by substituting x=eu/(k+1)x = \eu^{-u/(k+1)}, recognizing Γ\Gamma.)

Solution

Solution of Exercise 10.4.

(a) For x>0x > 0: 1ex1=ex1ex=n1enx\frac1{\eu^x - 1} = \frac{\eu^{-x}}{1 - \eu^{-x}} = \sum_{n\geq1}\eu^{-nx}, so xex1=n1xenx\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:

0x ⁣dxex1=n10xenx ⁣dx=n11n2=π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

(0xenx ⁣dx=n2\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, xlnx0-x\ln x \geq 0, so xx=exlnx=k(xlnx)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=eu/(k+1)x = \eu^{-u/(k+1)}:

01(xlnx)k ⁣dx=0(uk+1)kekuk+1  euk+1k+1 ⁣du=1(k+1)k+10ukeu ⁣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 01xx ⁣dx=k01(k+1)k+1=n1nn\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 f0f \geq 0 measurable with f ⁣dμ=0\int f\,\dd\mu = 0 satisfies f=0f = 0 a.e. (Consider {f1/n}\{f \geq 1/n\} and Markov’s inequality: μ({fa})1af\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: a1{fa}fa\,\mathbf 1_{\{f \geq a\}} \leq f, integrate: μ({fa})1af\mu(\{f \geq a\}) \leq \frac1a\int f. If f=0\int f = 0: μ({f1n})=0\mu(\{f \geq \frac1n\}) = 0 for every nn, and {f>0}=n{f1n}\{f > 0\} = \bigcup_n\{f \geq \frac1n\} is null.

(b) μ({f=})μ({fn})1nf0\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+sinxx ⁣dx=+\int_1^{+\infty}\abs{\frac{\sin x}x}\,\dd x = +\infty, while limA1Asinxx ⁣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)πsinxx ⁣dx1(k+1)πkπ(k+1)πsinx ⁣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,

πAsinxx ⁣dx=[cosxx]πAπAcosxx2 ⁣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 AA \to \infty (1x2\frac1{x^2} is integrable): conditional convergence without absolute integrability.

Exercise 10.7 ★★

Justify that F(t)=0+ex2cos(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)et2/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(ex2cos(tx))=xex2sin(tx)xex2\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)=0xex2sin(tx) ⁣dx=[12ex2sin(tx)]0t20ex2cos(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=xex2 ⁣dx\dd v = x\eu^{-x^2}\dd x). The linear ODE gives F(t)=F(0)et2/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+eaxebxx ⁣dx=lnba,\int_0^{+\infty}\frac{\eu^{-ax} - \eu^{-bx}}{x}\,\dd x = \ln\frac ba,

by writing the integrand as abext ⁣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 bab - a (bounded), and it decays like eax\eu^{-ax} at infinity. Fix aa and view I(b)=0eaxebxx ⁣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 bab \geq a: b(integrand)=ebxeax\abs{\partial_b(\text{integrand})} = \eu^{-bx} \leq \eu^{-ax}, and 0eax ⁣dx=1a<\int_0^\infty\eu^{-ax}\dd x = \frac1a < \infty: an integrable dominator. So Theorem 10.15 gives I(b)=0ebx ⁣dx=1bI'(b) = \int_0^\infty\eu^{-bx}\dd x = \frac1b, and I(a)=0I(a) = 0:

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

(Equivalently, the hint’s route: the integrand is abext ⁣dt0\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 f0f \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 g0g \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), f1An=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 g0g \geq 0 measurable, take simple sngs_n \nearrow g (Theorem 10.4): snfgfs_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)=0cos(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)=π2etf'(t) = \frac\pi2\eu^{-t} for t>0t > 0 (proved in Chapter 17), what is 0xsin(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)tx\abs{\sin(tx)} \leq \abs tx gives sin(tx)x(1+x2)t1+x2\abs{\frac{\sin(tx)}{x(1+x^2)}} \leq \frac{\abs t}{1+x^2}: the integral converges, and on tT\abs t \leq T the dominator T1+x2\frac{T}{1+x^2} yields continuity (Theorem 10.14). Differentiation: t=cos(tx)1+x211+x2\abs{\partial_t} = \abs{\frac{\cos(tx)}{1+x^2}} \leq \frac1{1+x^2}, integrable: f(t)=0cos(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)=π2etf'(t) = \frac\pi2\eu^{-t} for t>0t > 0: by oddness of ff, ff' is even, so f(t)=π2etf'(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 0xsin(tx)1+x2 ⁣dx\int_0^\infty\frac{x\sin(tx)}{1+x^2}\dd x equals f(t)=π2et-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,f0f_n, f \geq 0 be integrable with fnff_n \to f a.e. and fnf\int f_n \to \int f. (a) Show that fnf0\int\abs{f_n - f} \to 0. (Apply dominated convergence to gn=(ffn)+fg_n = (f - f_n)^+ \leq f, and write fnf=2gn(ffn)\int\abs{f_n - f} = 2\int g_n - \int(f - f_n).) (b) Show by example that the hypothesis fnf\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 fn0f_n \to 0 a.e., fn=00\int f_n = 0 \to 0, yet fn=2\int\abs{f_n} = 2. (c) Application (densities): if probability densities pnpp_n \to p a.e., then automatically pnp0\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=(ffn)+g_n = (f - f_n)^+: then 0gnf0 \leq g_n \leq f (positivity of fnf_n), gn0g_n \to 0 a.e., and ff is an integrable dominator: gn0\int g_n \to 0 (DCT). Since fnf=2(ffn)+(ffn)\abs{f_n - f} = 2(f - f_n)^+ - (f - f_n),

fnf=2gn(ffn)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 fn0f_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: fn0f_n \to 0 at every x0x \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 gnfg_n \leq f.

(c) Densities satisfy pn=1=p\int p_n = 1 = \int p: the hypothesis of (a) is automatic, so pnpp_n \to p a.e. forces pnpL10\norm{p_n - p}_{L^1} \to 0 — and hence convergence of the probabilities ApnAp\int_Ap_n \to \int_Ap uniformly over all measurable AA (A(pnp)pnp1\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) limn0n(1xn)nex/2 ⁣dx,(b) limn01nxn11+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) limn0 ⁣dx(1+x/n)nx1/n.\text{(c)}\ \lim_{n\to\infty}\int_0^\infty \frac{\dd x}{(1 + x/n)^n\,x^{1/n}} .

(For (a): (1x/n)nex(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 n2n \geq 2 by splitting at x=1x = 1.)

Solution

Solution of Exercise 10.12.

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

limn0n(1xn)nex/2 ⁣dx=0ex/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}, nxn1 ⁣dx= ⁣dun x^{n-1}\dd x = \dd u):

01nxn11+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/n1u^{1/n} \to 1, so the integrand tends to 12\frac12, dominated by 11: the limit is 12\frac12 (DCT). (The mass of nxn1nx^{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)nex(1 + x/n)^n \nearrow \eu^x and x1/n1x^{1/n} \to 1 (x>0x > 0): the integrand tends to ex\eu^{-x}. Dominator for n2n \geq 2: on (0,1]\intoc01, x1/nx1/2x^{-1/n} \leq x^{-1/2} and (1+x/n)n1(1 + x/n)^{-n} \leq 1: bound x1/2x^{-1/2}, integrable; on (1,)\intoo1\infty, x1/n1x^{-1/n} \leq 1 and (1+x/n)n1+(n2)x2n21+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:

limn0 ⁣dx(1+x/n)nx1/n=0ex ⁣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=+ex2 ⁣dx=π,D=0+sinxx ⁣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 t0t \geq 0 set

A(t)=(0tex2 ⁣dx)2,B(t)=01et2(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 AA' and BB'; show A(t)+B(t)=0A'(t) + B'(t) = 0. (In BB', substitute u=txu = tx.)
  2. Compute A(0)+B(0)A(0) + B(0) and limt+(A+B)(t)\lim_{t\to+\infty}(A + B)(t) — justify the limit under the integral in BB.
  3. Conclude 0ex2 ⁣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 t0t \geq 0 set

F(t)=0+etxsinxx ⁣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=limA0Asinxx ⁣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+etxsinx ⁣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)=π2arctantF(t) = \frac\pi2 - \arctan t on (0,+)\intoo0{+\infty}.
  4. The delicate point: D=limt0+F(t)D = \lim_{t\to0^+}F(t). Prove it by uniform control of the tail: for 0t10 \leq t \leq 1 and AπA \geq \pi, integrate by parts to show

    A+etxsinxx ⁣dxCA\Bigl|\int_A^{+\infty}\eu^{-tx}\frac{\sin x}x\,\dd x\Bigr| \leq \frac{C}{A}

    with CC independent of tt (differentiate etxx\frac{\eu^{-tx}}x and bound cos\abs{\cos} by 11; note tetx1/x(txetx)t\eu^{-tx} \leq 1/x\cdot(tx\eu^{-tx}) with supu0ueu<1\sup_{u\geq0}u\eu^{-u} < 1); then split F(t)DF(t) - D into [0,A][0, A] (where DCT applies as t0t \to 0) and [A,)[A, \infty).

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

Part III — Dividends.

  1. Compute 0+sin2xx2 ⁣dx\int_0^{+\infty}\frac{\sin^2x}{x^2}\,\dd x (integrate by parts and reduce to DD via sin2x=2sinxcosx\sin 2x = 2\sin x\cos x).
  2. Compute 0+1cosxx2 ⁣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 +eax2 ⁣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 logf\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(up)1/p(vq)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)=logng(n+1) - g(n) = \log n for every integer n1n \geq 1. For x(0,1]x \in \intoc01 and n2n \geq 2, compare the slopes of gg over [n1,n][n-1, n], [n,n+x][n, n+x] and [n,n+1][n, n+1], and deduce

    xlog(n1)    g(n+x)g(n)    xlogn.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 logf\log f convex. Unwinding the recursion into f(n+x)=x(x+1)(x+n1)f(x)f(n + x) = x(x+1)\cdots(x + n - 1)\,f(x) and f(n)=(n1)!f(n) = (n-1)!, deduce from question 14 that for x(0,1]x \in \intoc01

    f(x)=limnn!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)=01tx1(1t)y1 ⁣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+yB(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 xB(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=sin2θ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)=2x1πΓ(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)=22z1πΓ(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)!4nn!π\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)n1.\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=sin2θt = \sin^2\theta to show

    Wp=0π/2sinpθ ⁣dθ=12B(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+2WnW_{n+2} = \frac{n+1}{n+2}\,W_n for integers n0n \geq 0. Compute W2nW_{2n} and W2n+1W_{2n+1} in closed form, show W2n+1/W2n1W_{2n+1}/W_{2n} \to 1 by squeezing, and conclude with Wallis’s product

    π2=limnk=1n4k24k21.\frac\pi2 = \lim_{n\to\infty} \prod_{k=1}^{n}\frac{4k^2}{4k^2 - 1} .
  2. (The Gaussian meets a frequency) For bRb \in \R set

    Φ(b)=+ex2cos(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)=πeb2:\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+eaxebxx ⁣dx=logba,\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 aba \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)=2et20tex2 ⁣dxA'(t) = 2\eu^{-t^2}\int_0^t\eu^{-x^2}\dd x. For BB: t(et2(1+x2)1+x2)=2tet2(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)=2t01et2(1+x2) ⁣dx=2et201te(tx)2 ⁣dx=2et20teu2 ⁣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: 0B(t)et200 \leq B(t) \leq \eu^{-t^2} \to 0, and A(t)(0ex2 ⁣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 (0ex2 ⁣dx)2=π4\bigl(\int_0^\infty\eu^{-x^2}\dd x\bigr)^2 = \frac\pi4: 0ex2 ⁣dx=π2\int_0^\infty\eu^{-x^2}\dd x = \frac{\sqrt\pi}2, and by evenness G=πG = \sqrt\pi. Also Γ(12)=0x1/2ex ⁣dx=x=u220eu2 ⁣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: etxsinxxetx\abs{\eu^{-tx}\frac{\sin x}x} \leq \eu^{-tx}, integrable. For t=0t = 0, improper convergence: on [π,A][\pi, A],

πAsinxx ⁣dx=[cosxx]πAπAcosxx2 ⁣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 AA \to \infty; near 00 the integrand extends continuously by 11.

5. On [t0,+)[t_0, +\infty) (t0>0t_0 > 0): t(etxsinxx)=etxsinxet0x\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)=0etxsinx ⁣dx=Im0e(ti)x ⁣dx=Im1ti=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)0etx ⁣dx=1t0\abs{F(t)} \leq \int_0^\infty\eu^{-tx}\dd x = \frac1t \to 0. Integrating F=11+t2F' = -\frac1{1+t^2}: F(t)=CarctantF(t) = C - \arctan t, and tt \to \infty forces C=π2C = \frac\pi2: F(t)=π2arctantF(t) = \frac\pi2 - \arctan t on (0,+)\intoo0{+\infty}.

7. Integrate by parts on [A,R][A, R] with sinx=(cosx)\sin x = (-\cos x)' and let RR \to \infty:

Aetxsinxx ⁣dx=cosA  etAAAcosx(tx+1x2)etx ⁣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 cos1\abs{\cos} \leq 1: the first term is 1A\leq \frac1A; the integral is at most Atetx ⁣dxx+A ⁣dxx21A0tetx ⁣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)D0A(etx1)sinxx ⁣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]: etx1txtA\abs{\eu^{-tx} - 1} \leq tx \leq tA and sinxx1\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=limt0+F(t)=limt0+(π2arctant)=π2D = \lim_{t\to0^+}F(t) = \lim_{t\to0^+}\bigl(\frac\pi2 - \arctan t\bigr) = \frac\pi2.

9. By parts (u=sin2xu = \sin^2x, v=x2v' = x^{-2}):

0sin2xx2 ⁣dx=[sin2xx]0+02sinxcosxx ⁣dx=0sin2xx ⁣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: sin2x/x0\sin^2 x/x \to 0 at both ends).

10. By parts (u=1cosxu = 1 - \cos x, v=x2v' = x^{-2}): 01cosxx2 ⁣dx=0sinxx ⁣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: 1cosx=2sin2x21 - \cos x = 2\sin^2\frac x2, and the substitution x=2ux = 2u turns 2sin2(x/2)x2 ⁣dx\int\frac{2\sin^2(x/2)}{x^2}\dd x into sin2uu2 ⁣du\int\frac{\sin^2u}{u^2}\dd u: the two computations agree.

11. 0sin(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); Reax2 ⁣dx=π/a\int_\R\eu^{-ax^2}\dd x = \sqrt{\pi/a} (substitute u=axu = \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 supt[0,1]etxsinxx=sinxx\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=logfg = \log f is convex then f=expgf = \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 up=vq=1\int\abs u^p = \int\abs v^q = 1, Young gives uvupp+vqq\abs{uv} \leq \frac{\abs u^p}p + \frac{\abs v^q}q, integrate: uv1\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λ)y1et=(tx1et)λ(ty1et)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 [n1,n][n-1, n], [n,n+x][n, n+x], [n,n+1][n, n+1]:

log(n1)=g(n)g(n1)1g(n+x)g(n)xg(n+1)g(n)1=logn,\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)=(n1)!f(n) = (n-1)! (recursion from f(1)=1f(1) = 1) and f(n+x)=x(x+1)(x+n1)f(x)f(n + x) = x(x+1)\cdots(x+n-1)\,f(x), question 14 reads

(n1)x(n1)!    x(x+1)(x+n1)f(x)    nx(n1)!.(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+nn1\frac{x+n}n \to 1: the sandwich forces

f(x)=limnn!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 tx1\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(1t)yy]01+xy01tx1(1t)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 (1t)y=(1t)y1(1t)(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(1t)y1 ⁣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)=xf(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(1t)y1=()λ()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=sin2θt = \sin^2\theta,  ⁣dt=2sinθcosθ ⁣dθ\dd t = 2\sin\theta\cos\theta\,\dd\theta and t1/2(1t)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)=xg(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(x1)log2\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)=22z1πΓ(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)=(n12)(n32)12π=(2n1)(2n3)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 n1Γ(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 11n\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.6050187078.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=sin2θt = \sin^2\theta,  ⁣dt=2sinθcosθ ⁣dθ=2t1/2(1t)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 ⁣dt2t1/2(1t)1/2=1201tp+121(1t)121 ⁣dt=12B(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)/2B(n+12,12)=n+1n+2Wn.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=π2k=1n2k12k,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 sinn+1sinn\sin^{n+1} \leq \sin^n on [0,π/2]\intcc0{\pi/2}, the sequence (Wn)(W_n) is nonincreasing, so

1W2n+1W2nW2n+1W2n1=2n2n+11.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(2k1)(2k+1)=2πk=1n4k24k21,\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 nn \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 2xex2sin(2bx)-2x\,\eu^{-x^2}\sin(2bx), dominated by 2xex2L1(R)2\abs x\,\eu^{-x^2} \in L^1(\R) uniformly in bb: Φ\Phi is C1\mathcal C^1 with

Φ(b)=+2xex2sin(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=2xex2 ⁣dx\dd v = -2x\,\eu^{-x^2}\dd x (so v=ex2v = \eu^{-x^2}), the boundary terms vanish and

Φ(b)=2b+ex2cos(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)eb2=πeb2\Phi(b) = \Phi(0)\,\eu^{-b^2} = \sqrt\pi\,\eu^{-b^2} by Part I. Up to normalization this says the Fourier transform of ex2\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,

0eaxebxx=abesx ⁣ds(ba)eax,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 bab - a at x=0x = 0. Fix bb; on [a0,+)\intco{a_0}{+\infty} the aa-derivative of the integrand is eax-\eu^{-ax}, dominated by ea0xL1((0,+))\eu^{-a_0x} \in L^1(\intoo0{+\infty}), so II is C1\mathcal C^1 on (0,b)\intoo0b with

I(a)=0+eax ⁣dx=1a,henceI(a)=loga+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 0I(a)(ba)/a00 \leq I(a) \leq (b-a)/a \to 0 as aba \to b^-, so c=logbc = \log b and I(a)=logbaI(a) = \log\frac ba. The removable singularity is needed at 00: each term eax/x\eu^{-ax}/x separately has a divergent (logarithmic) integral near 00, and only the first-order cancellation eaxebx=O(x)\eu^{-ax} - \eu^{-bx} = O(x) makes the difference integrable there; at infinity each term is harmless on its own.