Mathematics · Book 3 · Bachelor Year 1

University Mathematics — Year 1

University Mathematics — Year 1 · Bachelor Year 1

15Integration on a Segment

The integral of the High School volume was founded on areas taken intuitively. This chapter constructs it: first for step functions, where the integral is a finite sum, then for continuous (and piecewise continuous) functions by uniform approximation — the place where Heine’s theorem (Theorem 13.22) earns its keep. The fundamental theorem of calculus then connects the construction to primitives, and Riemann sums connect it to discrete averages.

Throughout, a<ba < b are reals.

15.1 Step functions

Definition 15.1

φ ⁣:[a, b]→R\varphi \colon \intcc{a}{b} \to \R is a step function when there is a subdivision a=x0<x1<⋯<xn=ba = x_0 < x_1 < \dots < x_n = b such that φ\varphi is constant, equal to cic_i, on each open interval (xi−1, xi)\intoo{x_{i-1}}{x_i} (values at the nodes are unconstrained). Its integral is

∫abφ=∑i=1nci (xi−xi−1),\int_a^b \varphi = \sum_{i=1}^{n} c_i\,(x_i - x_{i-1}),

independent of the chosen subdivision (refine two subdivisions by their common one: each side is unchanged under refinement).

Proposition 15.2

On step functions, the integral is linear, increasing (φ≤ψ  ⟹  ∫φ≤∫ψ\varphi \leq \psi \implies \int\varphi \leq \int\psi), and satisfies the Chasles relation ∫ab=∫ac+∫cb\int_a^b = \int_a^c + \int_c^b for a<c<ba < c < b.

Proof. The engine is refinement invariance, stated in the definition: inserting one extra node t∈(xi−1, xi)t \in \intoo{x_{i-1}}{x_i} into a subdivision replaces the term ci(xi−xi−1)c_i(x_i - x_{i-1}) by ci(t−xi−1)+ci(xi−t)c_i(t - x_{i-1}) + c_i(x_i - t) — the same number — so the integral is unchanged by any finite refinement. Now take φ\varphi with subdivision σ\sigma and ψ\psi with subdivision σ′\sigma': on the common refinement σ∪σ′\sigma \cup \sigma' both are step functions with the same nodes, and on each piece φ+λψ\varphi + \lambda\psi is constant equal to ci+λdic_i + \lambda d_i: linearity reduces to linearity of finite sums. Increase: ci≤dic_i \leq d_i on each piece gives ∑ciΔi≤∑diΔi\sum c_i \Delta_i \leq \sum d_i \Delta_i (lengths Δi≥0\Delta_i \geq 0). Chasles: insert the node cc and split the sum at it. ∎

15.2 Integral of a continuous function

Theorem 15.3 (Uniform approximation)

Let ff be continuous on [a, b]\intcc{a}{b}. For every ε>0\varepsilon > 0 there are step functions φ,ψ\varphi, \psi with

φ≤f≤ψandψ−φ≤ε on [a, b].\varphi \leq f \leq \psi \qquad\text{and}\qquad \psi - \varphi \leq \varepsilon \text{ on } \intcc{a}{b}.

Proof. By Heine’s theorem (Theorem 13.22), ff is uniformly continuous: pick δ\delta for ε\varepsilon, and a subdivision of mesh <δ< \delta (equally spaced, say, with n>b−aδn > \frac{b - a}{\delta} pieces). On each closed piece [xi−1, xi]\intcc{x_{i-1}}{x_i}, ff attains a minimum mim_i and a maximum MiM_i (Theorem 13.13), and Mi−mi≤εM_i - m_i \leq \varepsilon (the two extremal points are within δ\delta). Define φ=mi\varphi = m_i and ψ=Mi\psi = M_i on (xi−1, xi)\intoo{x_{i-1}}{x_i} (and φ=ψ=f\varphi = \psi = f at the nodes). ∎

Theorem 15.4 (Definition of the integral)

Let ff be continuous on [a, b]\intcc{a}{b}. The two numbers

I−(f)=sup⁡{∫abφ:φ step, φ≤f},I+(f)=inf⁡{∫abψ:ψ step, ψ≥f}I_-(f) = \sup\Bigl\{\int_a^b \varphi : \varphi \text{ step},\ \varphi \leq f\Bigr\}, \qquad I_+(f) = \inf\Bigl\{\int_a^b \psi : \psi \text{ step},\ \psi \geq f\Bigr\}

are equal; their common value is the integral ∫abf\int_a^b f (also written ∫abf(t)  ⁣dt\int_a^b f(t)\,\dd t). It coincides with the previous notion on step functions, and extends to piecewise continuous functions by splitting [a, b]\intcc{a}{b} at the discontinuities (Chasles as a definition there).

Proof. Both sets are nonempty (ff is bounded) and every lower step integral is ≤\leq every upper one (monotonicity on step functions): so I−(f)≤I+(f)I_-(f) \leq I_+(f). By Theorem 15.3, for every ε\varepsilon there is a pair with ∫ψ−∫φ≤ε(b−a)\int\psi - \int\varphi \leq \varepsilon(b - a): the sup and inf are squeezed together, I−=I+I_- = I_+. ∎

Example 15.5 (Piecewise continuous, no drama)

The floor function on [0, 3]\intcc{0}{3} is a step function in disguise: splitting at its jumps,

∫03⌊t⌋  ⁣dt=∫010+∫121+∫232=0+1+2=3,\int_0^3 \lfloor t \rfloor\,\dd t = \int_0^1 0 + \int_1^2 1 + \int_2^3 2 = 0 + 1 + 2 = 3 ,

and the values at the jump points 1,21, 2 are irrelevant: changing a function at finitely many points changes no integral (the framing step functions are unaffected). This is the whole content of the “piecewise continuous” extension: cut at the finitely many discontinuities, integrate each continuous piece, add — Chasles as a definition.

Example 15.6 (The definition computes, once)

Let f(x)=xf(x) = x on [0, 1]\intcc{0}{1} and cut into nn equal pieces. The best step functions constant on the pieces are φ=k−1n\varphi = \frac{k-1}{n} and ψ=kn\psi = \frac kn on the kk-th piece, with

∫01φ=∑k=1nk−1n⋅1n=n−12n,∫01ψ=∑k=1nkn⋅1n=n+12n.\int_0^1 \varphi = \sum_{k=1}^{n} \frac{k-1}{n}\cdot\frac1n = \frac{n-1}{2n}, \qquad \int_0^1 \psi = \sum_{k=1}^{n} \frac{k}{n}\cdot\frac1n = \frac{n+1}{2n} .

Every lower integral is ≤I−(f)≤I+(f)≤\leq I_-(f) \leq I_+(f) \leq every upper one, so n−12n≤I−(f)≤I+(f)≤n+12n\frac{n-1}{2n} \leq I_-(f) \leq I_+(f) \leq \frac{n+1}{2n} for all nn: both squeeze onto 12\frac12, and ∫01x  ⁣dx=12\int_0^1 x\,\dd x = \frac12 straight from the definition. The closing insight: this is the first and last time we integrate from the definition — the fundamental theorem below replaces all such computations by one antiderivative lookup, which is the entire economic point of this chapter.

Theorem 15.7 (Properties)

For f,gf, g continuous (or piecewise continuous) on [a, b]\intcc{a}{b} and λ∈R\lambda \in \R:

  1. linearity: ∫(f+λg)=∫f+λ∫g\int (f + \lambda g) = \int f + \lambda \int g;
  2. monotonicity: f≤g  ⟹  ∫abf≤∫abgf \leq g \implies \int_a^b f \leq \int_a^b g; and ∣∫abf∣≤∫ab∣f∣≤(b−a) sup⁡∣f∣\bigl|\int_a^b f\bigr| \leq \int_a^b \abs f \leq (b - a)\, \sup\abs f;
  3. Chasles: ∫ab=∫ac+∫cb\int_a^b = \int_a^c + \int_c^b (with the convention ∫ba=−∫ab\int_b^a = -\int_a^b, valid for any order of the bounds);
  4. strict positivity: if ff is continuous, f≥0f \geq 0 and ∫abf=0\int_a^b f = 0, then f=0f = 0 everywhere on [a, b]\intcc{a}{b}.

Proof. (1)–(3) pass from step functions to the limit through the sup/inf definition. Linearity deserves the details once: given ε>0\varepsilon > 0, frame φf≤f≤ψf\varphi_f \leq f \leq \psi_f and φg≤g≤ψg\varphi_g \leq g \leq \psi_g with gaps ≤ε\leq \varepsilon (Theorem 15.3). For λ≥0\lambda \geq 0, φf+λφg≤f+λg≤ψf+λψg\varphi_f + \lambda\varphi_g \leq f + \lambda g \leq \psi_f + \lambda\psi_g is a framing by step functions with gap ≤(1+λ)ε\leq (1 + \lambda)\varepsilon, and its step integrals equal ∫φf+λ∫φg\int \varphi_f + \lambda\int\varphi_g etc. (Proposition 15.2): letting ε→0\varepsilon \to 0 squeezes ∫(f+λg)\int(f + \lambda g) onto ∫f+λ∫g\int f + \lambda\int g. For λ<0\lambda < 0, multiplying by λ\lambda reverses the framing of gg — the lower step function of λg\lambda g is λψg\lambda \psi_g — and the same squeeze runs with the roles swapped. The bound ∣∫f∣≤∫∣f∣\abs{\int f} \leq \int \abs f comes from −∣f∣≤f≤∣f∣-\abs f \leq f \leq \abs f and monotonicity.

(4) Contrapositive: if f(x0)=m>0f(x_0) = m > 0, continuity provides a subinterval of length η>0\eta > 0 on which f≥m2f \geq \frac m2; the step function worth m2\frac m2 there and 00 elsewhere is ≤f\leq f, so ∫f≥mη2>0\int f \geq \frac{m\eta}{2} > 0. ∎

Example 15.8 (Chasles at work: integrals with absolute values)

To integrate an absolute value, cut where the sign changes.

∫02∣x−1∣  ⁣dx=∫01(1−x)  ⁣dx+∫12(x−1)  ⁣dx=12+12=1,\int_0^2 \abs{x - 1}\,\dd x = \int_0^1 (1 - x)\,\dd x + \int_1^2 (x - 1)\,\dd x = \frac12 + \frac12 = 1 ,

and, cutting [0, 2π]\intcc{0}{2\pi} at π\pi:

∫02π∣sin⁡t∣  ⁣dt=∫0πsin⁡t  ⁣dt−∫π2πsin⁡t  ⁣dt=2+2=4,\int_0^{2\pi} \abs{\sin t}\,\dd t = \int_0^{\pi} \sin t\,\dd t - \int_{\pi}^{2\pi} \sin t\,\dd t = 2 + 2 = 4 ,

while ∫02πsin⁡t  ⁣dt=0\int_0^{2\pi} \sin t\,\dd t = 0: cancellation is real, and this is why the strict positivity statement (Theorem 15.7 (4)) carries the hypothesis f≥0f \geq 0 — without it, a vanishing integral proves nothing about ff. The closing insight: ∫∣f∣\int \abs f measures area, ∫f\int f measures signed balance; the inequality ∣∫f∣≤∫∣f∣\abs{\int f} \leq \int\abs f is the exact record of what cancellation can destroy.

15.3 The fundamental theorem of calculus

Theorem 15.9 (Fundamental theorem of calculus)

Let ff be continuous on an interval II and a∈Ia \in I. The function

F(x)=∫axf(t)  ⁣dtF(x) = \int_a^x f(t)\, \dd t

is of class C1C^1 on II, with F′=fF' = f: every continuous function on an interval has primitives. Consequently, for any primitive GG of ff:

∫abf(t)  ⁣dt=G(b)−G(a).\int_a^b f(t)\,\dd t = G(b) - G(a) .

Proof. Fix x0∈Ix_0 \in I and ε>0\varepsilon > 0; continuity at x0x_0 provides δ\delta with ∣f(t)−f(x0)∣≤ε\abs{f(t) - f(x_0)} \leq \varepsilon for ∣t−x0∣≤δ\abs{t - x_0} \leq \delta. For 0<∣h∣≤δ0 < \abs{h} \leq \delta (and x0+h∈Ix_0 + h \in I), Chasles gives

F(x0+h)−F(x0)h−f(x0)=1h∫x0x0+h(f(t)−f(x0)) ⁣dt,\frac{F(x_0 + h) - F(x_0)}{h} - f(x_0) = \frac 1h \int_{x_0}^{x_0+h} \bigl(f(t) - f(x_0)\bigr)\dd t ,

whose absolute value is ≤1∣h∣⋅∣h∣ ε=ε\leq \frac{1}{\abs h}\cdot \abs h\, \varepsilon = \varepsilon (bound (2), valid for either order of the bounds). So F′(x0)=f(x0)F'(x_0) = f(x_0); F′=fF' = f is continuous: FF is C1C^1. If G′=fG' = f too, then (G−F)′=0(G - F)' = 0 on the interval, so G=F+cG = F + c (Corollary 14.12), and G(b)−G(a)=F(b)−F(a)=∫abfG(b) - G(a) = F(b) - F(a) = \int_a^b f. ∎

Example 15.10 (Symmetry before computation)

On a symmetric interval, parity does the work: if ff is odd, the substitution t↦−tt \mapsto -t sends ∫−a0f\int_{-a}^{0} f to −∫0af-\int_0^a f, so

∫−aaf(t)  ⁣dt=0;if f is even,∫−aaf=2∫0af.\int_{-a}^{a} f(t)\,\dd t = 0 ; \qquad\text{if $f$ is even,}\quad \int_{-a}^{a} f = 2\int_0^a f .

Thus ∫−11t3cos⁡t1+t4  ⁣dt=0\int_{-1}^{1} \frac{t^3\cos t}{1 + t^4}\,\dd t = 0 with no primitive in sight (the integrand is odd), and ∫−ππt2cos⁡t  ⁣dt=2∫0πt2cos⁡t  ⁣dt\int_{-\pi}^{\pi} t^2\cos t\,\dd t = 2\int_0^\pi t^2\cos t\, \dd t. Check symmetry before reaching for techniques: the fastest integral is the one never computed.

Example 15.11 (Recognizing a derivative on sight)

Compute ∫0π/2 ⁣dx1+cos⁡x\displaystyle\int_0^{\pi/2} \frac{\dd x}{1 + \cos x}. The half-angle identity 1+cos⁡x=2cos⁡2x21 + \cos x = 2\cos^2\frac x2 turns the integrand into 12(1+tan⁡2x2)\frac{1}{2}\bigl(1 + \tan^2\frac x2\bigr), which is exactly the derivative of tan⁡x2\tan\frac x2:

∫0π/2 ⁣dx1+cos⁡x=[tan⁡x2]0π/2=tan⁡π4=1.\int_0^{\pi/2} \frac{\dd x}{1 + \cos x} = \Bigl[\tan\frac x2\Bigr]_0^{\pi/2} = \tan\frac\pi4 = 1 .

No substitution machinery was needed — only the reflex of reading an integrand as somebody’s derivative, the fundamental theorem doing the rest. (The systematic tool behind such trigonometric integrals, the substitution t=tan⁡x2t = \tan\frac x2, belongs to the standard toolkit built from Theorem 15.15 (2).)

Example 15.12 (Functions defined by integrals)

The fundamental theorem manufactures functions. Let

F(x)=∫0xe−t2  ⁣dt.F(x) = \int_0^x \eu^{-t^2}\,\dd t .

No combination of classical functions has derivative e−t2\eu^{-t^2} (a theorem of Liouville, admitted); yet FF exists, is C1C^1 with F′(x)=e−x2>0F'(x) = \eu^{-x^2} > 0, strictly increasing, odd (substitute t↦−tt \mapsto -t), and bounded: for x≥1x \geq 1,

F(x)−F(1)=∫1xe−t2 ⁣dt≤∫1xe−t ⁣dt≤e−1,F(x) - F(1) = \int_1^x \eu^{-t^2}\dd t \leq \int_1^x \eu^{-t}\dd t \leq \eu^{-1} ,

so F≤F(1)+e−1≤1+e−1F \leq F(1) + \eu^{-1} \leq 1 + \eu^{-1}. (The exact limit, π2\frac{\sqrt\pi}{2}, is computed with double integrals in the Year 3 volume.) Chain rule for moving bounds:  ⁣d ⁣dx∫xx2e−t2 ⁣dt=2x e−x4−e−x2\frac{\dd}{\dd x}\int_x^{x^2} \eu^{-t^2}\dd t = 2x\,\eu^{-x^4} - \eu^{-x^2}. The closing insight: integration creates new functions from old ones, with all their properties readable from the integrand — the primitive you cannot write down is still a function you fully control.

Example 15.13 (Estimating without evaluating)

The integrals Rn=∫01tn1+t  ⁣dtR_n = \int_0^1 \frac{t^n}{1 + t}\,\dd t have no pleasant closed form, yet monotonicity pins them precisely: on [0, 1]\intcc{0}{1}, 12≤11+t≤1\frac12 \leq \frac{1}{1+t} \leq 1, so

12(n+1)=12∫01tn  ⁣dt  ≤  Rn  ≤  ∫01tn  ⁣dt=1n+1:\frac{1}{2(n+1)} = \frac12\int_0^1 t^n\,\dd t \;\leq\; R_n \;\leq\; \int_0^1 t^n \,\dd t = \frac{1}{n+1} :

the exact order of decay (Rn∼R_n \sim a multiple of 1n\frac1n, in fact Rn∼12nR_n \sim \frac{1}{2n}) with two lines and no antiderivative. The weekend problems of this chapter and the next run on exactly such framings — the analyst’s first instinct before an integral should be bound it, and only then, if needed, compute it.

Example 15.14 (Average values)

The mean of a continuous ff over [a, b]\intcc{a}{b} is 1b−a∫abf\frac{1}{b-a}\int_a^b f. For the arch of the sine:

1π∫0πsin⁡t  ⁣dt=1π[−cos⁡t]0π=2π≈0.637:\frac{1}{\pi}\int_0^\pi \sin t\,\dd t = \frac{1}{\pi}\bigl[-\cos t\bigr]_0^\pi = \frac{2}{\pi} \approx 0.637 :

a full positive arch averages not to 12\frac12 but to 2π\frac2\pi — the curve spends more time high than a triangle would. By Exercise 15.11 (mean value theorem for integrals, g=1g = 1), the mean is a value: sin⁡c=2π\sin c = \frac2\pi for some c∈(0, π)c \in \intoo{0}{\pi}. And by the Riemann sums of this chapter, the mean is the limit of ordinary averages of nn samples — the bridge between the discrete mean of data and the continuous mean of a signal, which is how the integral enters physics.

Theorem 15.15 (Integration by parts; substitution)

  1. If u,vu, v are C1C^1 on [a, b]\intcc{a}{b}:

    ∫abu′v=[uv]ab−∫abuv′.\int_a^b u'v = \bigl[uv\bigr]_a^b - \int_a^b uv' .
  2. If φ\varphi is C1C^1 on [α, β]\intcc{\alpha}{\beta} and ff continuous on φ([α, β])\varphi(\intcc{\alpha}{\beta}):

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

Proof. (1) (uv)′=u′v+uv′(uv)' = u'v + uv'; integrate over [a, b]\intcc{a}{b} and apply the fundamental theorem to the C1C^1 function uvuv.

(2) Let FF be a primitive of ff on the image interval (Theorem 15.9). Then (F∘φ)′=(f∘φ) φ′(F \circ \varphi)' = (f \circ \varphi)\,\varphi' (chain rule), so both sides equal F(φ(β))−F(φ(α))F(\varphi(\beta)) - F(\varphi(\alpha)). ∎

Example 15.16

∫01t et ⁣dt=[t et]01−∫01et ⁣dt=e−(e−1)=1\int_0^1 t\,\eu^t \dd t = \bigl[t\,\eu^t\bigr]_0^1 - \int_0^1 \eu^t\dd t = \eu - (\eu - 1) = 1. And with the substitution x=sin⁡tx = \sin t (t∈[0, π2]t \in \intcc{0}{\frac\pi2}):

∫011−x2  ⁣dx=∫0π/2cos⁡t⋅cos⁡t  ⁣dt=∫0π/21+cos⁡2t2  ⁣dt=π4,\int_0^1 \sqrt{1 - x^2}\, \dd x = \int_0^{\pi/2} \cos t \cdot \cos t \, \dd t = \int_0^{\pi/2} \frac{1 + \cos 2t}{2}\, \dd t = \frac\pi4 ,

— a quarter of the unit disk, as geometry demands. (Linearization from Method 3.11 at work.)

Remark 15.17 (Common pitfalls in integral calculus)

(i) Substitutions must be C1C^1 on the whole interval: the change x=1tx = \frac1t is illegal across 00; applied blindly to ∫−11 ⁣dx1+x2\int_{-1}^{1}\frac{\dd x}{1 + x^2} it “proves” that the integral equals its own negative. When a substitution has a singularity, cut the interval first (Chasles), substitute on each piece, and only then recombine. (ii) Logarithmic primitives need absolute values: ∫ ⁣dxx−2=ln⁡∣x−2∣+C\int \frac{\dd x}{x - 2} = \ln\abs{x - 2} + C on each side of 22 separately — writing ln⁡(x−2)\ln(x - 2) on (0, 1)\intoo{0}{1} is writing the logarithm of a negative number; and the constant CC may differ on the two sides of the singularity. (iii) A vanishing integral does not kill the function: ∫02πsin⁡=0\int_0^{2\pi}\sin = 0; positivity of the integrand is required before concluding f=0f = 0 (Example 15.8). (iv) Riemann sums must be calibrated: in b−an∑f(a+kb−an)\frac{b-a}{n}\sum f\bigl(a + k\frac{b-a}{n}\bigr), the step outside and the points inside must match the same subdivision — the frequent error is a sum ∑k=1nf(kn)\sum_{k=1}^{n} f\bigl(\frac kn\bigr) without the factor 1n\frac1n, which diverges instead of converging to ∫01f\int_0^1 f. Checklist before invoking Theorem 15.20: factor out 1n\frac1n, rewrite the summand as ff of kn\frac kn, name ff and check its continuity.

Example 15.18 (Guess, differentiate, adjust)

What is ∫1x(ln⁡t)2  ⁣dt\int_1^x (\ln t)^2\,\dd t? Guess a primitive of the form t P(ln⁡t)t\,P(\ln t) with PP polynomial and differentiate:

(t P(ln⁡t))′=P(ln⁡t)+P′(ln⁡t).\bigl(t\,P(\ln t)\bigr)' = P(\ln t) + P'(\ln t) .

We need P(u)+P′(u)=u2P(u) + P'(u) = u^2: take P(u)=u2−2u+2P(u) = u^2 - 2u + 2 (matching coefficients downward from u2u^2). Hence

∫1x(ln⁡t)2  ⁣dt=[t((ln⁡t)2−2ln⁡t+2)]1x=x(ln⁡x)2−2xln⁡x+2x−2,\int_1^x (\ln t)^2\,\dd t = \bigl[t\bigl((\ln t)^2 - 2\ln t + 2\bigr)\bigr]_1^x = x(\ln x)^2 - 2x\ln x + 2x - 2 ,

a result otherwise reached by two integrations by parts. The closing insight: for integrands of the form (polynomial in ln⁡t\ln t) or (polynomial times eλt\eu^{\lambda t}), the primitive has the same shape — differentiating a shaped guess converts integration into linear algebra on coefficients, faster and less error-prone than iterated parts.

Example 15.19 (The boomerang integral)

Compute I=∫0π/2excos⁡x  ⁣dxI = \int_0^{\pi/2} \eu^x \cos x\,\dd x. Integrate by parts twice, differentiating the trigonometric factor each time:

I=[exsin⁡x]0π/2−∫0π/2exsin⁡x  ⁣dx=eπ/2−J,J=[−excos⁡x]0π/2+∫0π/2excos⁡x  ⁣dx=1+I.I = \bigl[\eu^x \sin x\bigr]_0^{\pi/2} - \int_0^{\pi/2} \eu^x \sin x\,\dd x = \eu^{\pi/2} - J, \qquad J = \bigl[-\eu^x\cos x\bigr]_0^{\pi/2} + \int_0^{\pi/2} \eu^x\cos x\,\dd x = 1 + I .

The integral has returned to itself: I=eπ/2−1−II = \eu^{\pi/2} - 1 - I, whence

I=eπ/2−12.I = \frac{\eu^{\pi/2} - 1}{2} .

The closing insight: when the integrand is a product of two functions that reproduce themselves under differentiation (eax\eu^{ax}, cos⁡bx\cos bx, sin⁡bx\sin bx), two integrations by parts produce a linear equation for the unknown integral — solve it instead of integrating; equivalently, pass through e(a+ib)x\eu^{(a + \iu b)x} (Chapter 3) and take real parts. Both roads give the same answer, and checking that they do is a free sanity test.

15.4 Riemann sums

Theorem 15.20 (Riemann sums)

Let ff be continuous on [a, b]\intcc{a}{b}. Then

Sn=b−an∑k=0n−1f(a+k b−an)→n→∞∫abf(t)  ⁣dt,S_n = \frac{b - a}{n} \sum_{k=0}^{n-1} f\Bigl(a + k\,\frac{b-a}{n}\Bigr) \xrightarrow[n \to \infty]{} \int_a^b f(t)\, \dd t ,

and likewise with any evaluation points inside the subintervals.

Proof. SnS_n is the integral of the step function φn\varphi_n equal to f(a+kb−an)f(a + k\frac{b-a}{n}) on the kk-th subinterval. Given ε>0\varepsilon > 0, uniform continuity (Heine) provides δ\delta; for n>b−aδn > \frac{b-a}{\delta}, every point of a subinterval is within δ\delta of its evaluation point, so ∣f−φn∣≤ε\abs{f - \varphi_n} \leq \varepsilon on [a, b]\intcc{a}{b}, whence

∣∫abf−Sn∣=∣∫ab(f−φn)∣≤(b−a) ε.\Bigl| \int_a^b f - S_n \Bigr| = \Bigl| \int_a^b (f - \varphi_n) \Bigr| \leq (b-a)\,\varepsilon . \qedhere

∎

A left Riemann sum with n = 8 rectangles: as the mesh shrinks, uniform continuity forces the staircase area toward ∈t_ab f.
A left Riemann sum with n=8n = 8 rectangles: as the mesh shrinks, uniform continuity forces the staircase area toward ∫abf\int_a^b f.

Example 15.21

∑k=1n1n+k=1n∑k=1n11+k/n→n→∞∫01 ⁣dx1+x=ln⁡2\displaystyle\sum_{k=1}^{n} \frac{1}{n + k} = \frac 1n \sum_{k=1}^{n} \frac{1}{1 + k/n} \xrightarrow[n\to\infty]{} \int_0^1 \frac{\dd x}{1 + x} = \ln 2: a limit invisible to elementary bounds, transparent as a Riemann sum.

Example 15.22 (A second Riemann sum, with calibration)

Find lim⁡n→∞∑k=1nn(n+k)2\lim_{n\to\infty} \sum_{k=1}^{n} \dfrac{n}{(n+k)^2}. Calibrate:

∑k=1nn(n+k)2=1n∑k=1nn2(n+k)2=1n∑k=1n1(1+kn)2,\sum_{k=1}^{n} \frac{n}{(n+k)^2} = \frac{1}{n}\sum_{k=1}^{n} \frac{n^2}{(n+k)^2} = \frac1n \sum_{k=1}^{n} \frac{1}{\bigl(1 + \frac kn\bigr)^2} ,

a Riemann sum of the continuous f(x)=1(1+x)2f(x) = \frac{1}{(1+x)^2} on [0, 1]\intcc{0}{1}: the limit is

∫01 ⁣dx(1+x)2=[−11+x]01=12.\int_0^1 \frac{\dd x}{(1 + x)^2} = \Bigl[-\frac{1}{1+x}\Bigr]_0^1 = \frac12 .

The closing insight: the whole art is the middle line — force the summand into the shape f(kn)f(\frac kn) at the cost of extracting exactly one factor 1n\frac1n; once the shape is right, the theorem does the analysis and the fundamental theorem does the arithmetic.

Remark 15.23 (Where the integral works next)

The chapter’s constructions each have a sequel. Riemann sums return in Chapter 17 as the bridge between series and integrals (comparison of ∑1nα\sum \frac{1}{n^\alpha} with ∫ ⁣dttα\int \frac{\dd t}{t^\alpha}); the integral remainder is the sharpest form of Taylor’s formula in Chapter 16; the sup⁡\sup-based definition is the prototype for the Lebesgue integral of the Year 3 volume, where the same three properties (linearity, monotonicity, a convergence theorem) are rebuilt on a far larger class of functions. And the weekend problem below turns integration by parts into arithmetic: the irrationality of π2\pi^2.

15.5 Exercises

Exercise 15.1 ★

Compute: ∫01 ⁣dxx2−4\displaystyle\int_0^1 \frac{\dd x}{x^2 - 4} (partial fractions, Chapter 9); ∫1eln⁡t  ⁣dt\displaystyle\int_1^{\eu} \ln t \,\dd t; ∫0πtsin⁡t  ⁣dt\displaystyle\int_0^{\pi} t \sin t\, \dd t.

Solution

Solution of Exercise 15.1.

1x2−4=1/4x−2−1/4x+2\dfrac{1}{x^2-4} = \dfrac{1/4}{x - 2} - \dfrac{1/4}{x+2} (cover-up), so

∫01 ⁣dxx2−4=14[ln⁡∣x−2∣−ln⁡∣x+2∣]01=14(ln⁡13−ln⁡1)=−ln⁡34.\int_0^1 \frac{\dd x}{x^2 - 4} = \frac14\Bigl[\ln\abs{x-2} - \ln\abs{x+2}\Bigr]_0^1 = \frac14\Bigl(\ln\frac{1}{3} - \ln 1\Bigr) = -\frac{\ln 3}{4}.

By parts (u′=1u' = 1, v=ln⁡tv = \ln t): ∫1eln⁡t  ⁣dt=[tln⁡t]1e−∫1e ⁣dt=e−(e−1)=1\int_1^{\eu} \ln t\,\dd t = \bigl[t\ln t\bigr]_1^{\eu} - \int_1^{\eu} \dd t = \eu - (\eu - 1) = 1.

By parts (u′=sin⁡tu' = \sin t, v=tv = t): ∫0πtsin⁡t  ⁣dt=[−tcos⁡t]0π+∫0πcos⁡t  ⁣dt=π+0=π\int_0^\pi t\sin t\,\dd t = \bigl[-t\cos t\bigr]_0^\pi + \int_0^\pi \cos t\,\dd t = \pi + 0 = \pi.

Exercise 15.2 ★

Compute ∫01t(t2+1)2  ⁣dt\displaystyle\int_0^{1} \frac{t}{(t^2+1)^2}\,\dd t (substitution) and ∫0π/2cos⁡3t  ⁣dt\displaystyle\int_0^{\pi/2} \cos^3 t\, \dd t (write cos⁡3=cos⁡(1−sin⁡2)\cos^3 = \cos(1 - \sin^2)).

Solution

Solution of Exercise 15.2.

Substitution u=t2+1u = t^2 + 1,  ⁣du=2t  ⁣dt\dd u = 2t\,\dd t:

∫01t  ⁣dt(t2+1)2=12∫12 ⁣duu2=12[−1u]12=14.\int_0^1 \frac{t\,\dd t}{(t^2+1)^2} = \frac12 \int_1^2 \frac{\dd u}{u^2} = \frac12\Bigl[-\frac1u\Bigr]_1^2 = \frac14 .

With u=sin⁡tu = \sin t: ∫0π/2cos⁡3t  ⁣dt=∫0π/2(1−sin⁡2t)cos⁡t  ⁣dt=[sin⁡t−sin⁡3t3]0π/2=1−13=23\int_0^{\pi/2} \cos^3 t\,\dd t = \int_0^{\pi/2} (1 - \sin^2 t)\cos t\,\dd t = \bigl[\sin t - \frac{\sin^3 t}{3}\bigr]_0^{\pi/2} = 1 - \frac13 = \frac23.

Exercise 15.3 ★

Find the limits, as Riemann sums:

un=∑k=1nnn2+k2,vn=1n(2n)!n! nnn (take logarithms).u_n = \sum_{k=1}^{n} \frac{n}{n^2 + k^2}, \qquad v_n = \frac{1}{n}\sqrt[n]{\frac{(2n)!}{n!\,n^n}} \ \emph{(take logarithms)}.
Solution

Solution of Exercise 15.3.

un=1n∑k=1n11+(k/n)2u_n = \dfrac1n \sum_{k=1}^{n} \dfrac{1}{1 + (k/n)^2}: a Riemann sum of x↦11+x2x \mapsto \frac{1}{1+x^2} on [0, 1]\intcc{0}{1}, so un→∫01 ⁣dx1+x2=arctan⁡1=π4u_n \to \int_0^1 \frac{\dd x}{1+x^2} = \arctan 1 = \dfrac\pi4.

ln⁡vn=1n∑k=1nln⁡n+kn=1n∑k=1nln⁡(1+kn)→∫01ln⁡(1+x)  ⁣dx=[(1+x)ln⁡(1+x)−x]01=2ln⁡2−1\ln v_n = \dfrac1n \sum_{k=1}^{n} \ln\dfrac{n+k}{n} = \dfrac 1n \sum_{k=1}^n \ln\Bigl(1 + \dfrac kn\Bigr) \to \int_0^1 \ln(1+x)\,\dd x = \bigl[(1+x)\ln(1+x) - x\bigr]_0^1 = 2\ln 2 - 1. Hence vn→e2ln⁡2−1=4ev_n \to \eu^{2\ln 2 - 1} = \dfrac 4\eu. (Check of the identification: (2n)!n! nn=∏k=1nn+kn\frac{(2n)!}{n!\,n^n} = \prod_{k=1}^{n} \frac{n+k}{n}.)

Exercise 15.4 ★

Let ff be continuous on [0, 1]\intcc{0}{1}. Compute lim⁡n→∞∫01xnf(x)  ⁣dx\lim_{n\to\infty} \int_0^1 x^n f(x)\,\dd x. (Cut [0, 1]\intcc{0}{1} at 1−δ1 - \delta.)

Solution

Solution of Exercise 15.4.

The limit is 00. Let M=sup⁡∣f∣M = \sup \abs f and ε∈(0, 1)\varepsilon \in \intoo{0}{1}. Cut at 1−ε1 - \varepsilon:

∣∫01xnf∣≤∫01−εxn∣f∣+∫1−ε1xn∣f∣≤M (1−ε)n+Mε.\Bigl| \int_0^1 x^n f \Bigr| \leq \int_0^{1 - \varepsilon} x^n \abs f + \int_{1-\varepsilon}^1 x^n \abs f \leq M\,(1-\varepsilon)^n + M\varepsilon .

Since (1−ε)n→0(1 - \varepsilon)^n \to 0 (Exercise 11.3), the limsup of the left side is ≤Mε\leq M\varepsilon for every ε\varepsilon: the integral tends to 00.

Exercise 15.5 ★★

(Cauchy–Schwarz) For f,gf, g continuous on [a, b]\intcc{a}{b}, prove

(∫abfg) ⁣2≤∫abf2⋅∫abg2,\Bigl(\int_a^b fg\Bigr)^{\!2} \leq \int_a^b f^2 \cdot \int_a^b g^2 ,

by expanding ∫ab(f+λg)2≥0\int_a^b (f + \lambda g)^2 \geq 0 as a quadratic in λ\lambda. When is it an equality?

Solution

Solution of Exercise 15.5.

Q(λ)=∫ab(f+λg)2=∫f2+2λ∫fg+λ2∫g2≥0Q(\lambda) = \int_a^b (f + \lambda g)^2 = \int f^2 + 2\lambda \int fg + \lambda^2 \int g^2 \geq 0 for all λ\lambda. If ∫g2=0\int g^2 = 0, then g=0g = 0 (strict positivity, Theorem 15.7 (4)) and the inequality is 0≤00 \leq 0. Otherwise QQ is a genuine quadratic, everywhere ≥0\geq 0: its discriminant is ≤0\leq 0, i.e. (∫fg)2≤∫f2∫g2\bigl(\int fg\bigr)^2 \leq \int f^2 \int g^2.

Equality iff the discriminant vanishes iff Q(λ0)=0Q(\lambda_0) = 0 for some λ0\lambda_0, i.e. ∫(f+λ0g)2=0\int (f + \lambda_0 g)^2 = 0, i.e. (strict positivity again) f=−λ0gf = -\lambda_0 g: equality holds exactly when ff and gg are proportional.

Exercise 15.6 ★★

Let ff be continuous on R\R, TT-periodic. Prove that ∫aa+Tf\int_a^{a+T} f does not depend on aa, and that 1x∫0xf(t)  ⁣dt→1T∫0Tf\frac1x \int_0^x f(t)\,\dd t \to \frac 1T \int_0^T f as x→+∞x \to +\infty.

Solution

Solution of Exercise 15.6.

Let Φ(a)=∫aa+Tf\Phi(a) = \int_a^{a+T} f. By the fundamental theorem (Theorem 15.9), Φ\Phi is differentiable with Φ′(a)=f(a+T)−f(a)=0\Phi'(a) = f(a + T) - f(a) = 0: constant.

For x>0x > 0, write x=nT+rx = nT + r, 0≤r<T0 \leq r < T (n=⌊x/T⌋n = \lfloor x/T \rfloor). Chasles:

∫0xf=n∫0Tf+∫nTnT+rf,∣∫nTnT+rf∣≤Tsup⁡[0, T]∣f∣=C.\int_0^x f = n \int_0^T f + \int_{nT}^{nT + r} f, \qquad \Bigl| \int_{nT}^{nT+r} f \Bigr| \leq T \sup_{\intcc{0}{T}} \abs f = C .

Then 1x∫0xf=nTx⋅1T∫0Tf+O(1x)\frac 1x \int_0^x f = \frac{nT}{x}\cdot\frac 1T \int_0^T f + O\bigl(\frac 1x\bigr), and nTx→1\frac{nT}{x} \to 1: the limit is 1T∫0Tf\frac1T \int_0^T f.

Exercise 15.7 ★★

For ff continuous on [0, 1]\intcc{0}{1} with ∫01f=12\int_0^1 f = \frac12, prove that ff has a fixed point in [0, 1]\intcc{0}{1}. (Integrate f(x)−xf(x) - x and use strict positivity, Theorem 15.7 (4), through its contrapositive combined with the intermediate value theorem.)

Solution

Solution of Exercise 15.7.

Let g(x)=f(x)−xg(x) = f(x) - x: continuous, with

∫01g=∫01f−12=0.\int_0^1 g = \int_0^1 f - \frac12 = 0 .

If gg never vanished, the intermediate value theorem would force a constant sign (a continuous function on an interval taking both signs vanishes); say g>0g > 0. Then, by strict positivity (Theorem 15.7 (4) applied to g>0g > 0, giving ∫g>0\int g > 0): contradiction with ∫g=0\int g = 0. So g(c)=0g(c) = 0 for some cc: f(c)=cf(c) = c.

Exercise 15.8 ★★★

(Wallis integrals) Let Wn=∫0π/2sin⁡nt  ⁣dtW_n = \int_0^{\pi/2} \sin^n t\,\dd t.

  1. Prove the recurrence nWn=(n−1)Wn−2n W_n = (n-1) W_{n-2} (n≥2n \geq 2) by parts, and compute W0,W1W_0, W_1, then W2pW_{2p} and W2p+1W_{2p+1} in closed form.
  2. Prove that (Wn)(W_n) is decreasing with Wn+1Wn→1\frac{W_{n+1}}{W_n} \to 1; prove that the quantity (n+1) Wn+1Wn(n+1)\,W_{n+1} W_n is constant, equal to π2\frac\pi2; and deduce the equivalent Wn∼π2nW_n \sim \sqrt{\dfrac{\pi}{2n}}.
Solution

Solution of Exercise 15.8.

  1. By parts with u′=sin⁡tu' = \sin t, v=sin⁡n−1tv = \sin^{n-1} t:

    Wn=[−cos⁡tsin⁡n−1t]0π/2+(n−1)∫0π/2cos⁡2t sin⁡n−2t  ⁣dt=(n−1)(Wn−2−Wn),W_n = \bigl[-\cos t\sin^{n-1}t\bigr]_0^{\pi/2} + (n-1)\int_0^{\pi/2} \cos^2 t\,\sin^{n-2} t\,\dd t = (n-1)(W_{n-2} - W_n),

    so nWn=(n−1)Wn−2nW_n = (n-1)W_{n-2}. From W0=π2W_0 = \frac\pi2, W1=1W_1 = 1:

    W2p=(2p−1)(2p−3)⋯1(2p)(2p−2)⋯2 π2=(2p)!4p(p!)2 π2,W2p+1=(2p)(2p−2)⋯2(2p+1)(2p−1)⋯3=4p(p!)2(2p+1)!.W_{2p} = \frac{(2p-1)(2p-3)\cdots 1}{(2p)(2p-2)\cdots 2}\, \frac{\pi}{2} = \frac{(2p)!}{4^p (p!)^2}\,\frac\pi2, \qquad W_{2p+1} = \frac{(2p)(2p-2)\cdots 2}{(2p+1)(2p-1)\cdots 3} = \frac{4^p (p!)^2}{(2p+1)!} .
  2. On (0, π2)\intoo{0}{\frac\pi2}, 0<sin⁡t<10 < \sin t < 1, so sin⁡n+1<sin⁡n\sin^{n+1} < \sin^n and (Wn)(W_n) is (strictly) decreasing, positive. Sandwiching with the recurrence:

    nn+1=Wn+1Wn−1≤Wn+1Wn≤1  ⟹  Wn+1Wn→1.\frac{n}{n+1} = \frac{W_{n+1}}{W_{n-1}} \leq \frac{W_{n+1}}{W_n} \leq 1 \quad\implies\quad \frac{W_{n+1}}{W_n} \to 1 .

    Invariant: an=(n+1)Wn+1Wna_n = (n+1)W_{n+1}W_n satisfies an=an−1a_n = a_{n-1} by the recurrence (n+1)Wn+1=nWn−1(n+1)W_{n+1} = nW_{n-1}, so an=a0=1⋅W1W0=π2a_n = a_0 = 1 \cdot W_1 W_0 = \frac\pi2. Then

    nWn2∼(n+1)Wn+1Wn=π2  ⟹  Wn∼π2n.n W_n^2 \sim (n+1) W_{n+1} W_n = \frac\pi2 \quad\implies\quad W_n \sim \sqrt{\frac{\pi}{2n}} .

Exercise 15.9 ★★★

(Niven: π\pi is irrational) Suppose π=ab\pi = \frac ab with a,b∈N∗a, b \in \N^*, and set, for nn to be chosen,

P(x)=xn(a−bx)nn!,In=∫0πP(x)sin⁡x  ⁣dx.P(x) = \frac{x^n (a - bx)^n}{n!}, \qquad I_n = \int_0^{\pi} P(x)\sin x\, \dd x .
  1. Prove that 0<In≤π (πa)nn!0 < I_n \leq \pi\,\frac{(\pi a)^n}{n!}, which is <1< 1 for nn large.
  2. Prove that PP and all its derivatives take integer values at 00 and at π=ab\pi = \frac ab. (Binomial expansion: the coefficients of PP times k!k! are integers for k≥nk \geq n; and P(π−x)=P(x)P(\pi - x) = P(x).)
  3. Set Q=P−P′′+P(4)−…Q = P - P'' + P^{(4)} - \dots (a finite sum). Check that (Q′sin⁡x−Qcos⁡x)′=Psin⁡x\bigl(Q'\sin x - Q\cos x\bigr)' = P \sin x, and deduce that In=Q(π)+Q(0)I_n = Q(\pi) + Q(0) is an integer.
  4. Conclude.
Solution

Solution of Exercise 15.9.

  1. On (0, π)\intoo{0}{\pi}: x>0x > 0, a−bx=b(ab−x)=b(π−x)>0a - bx = b(\frac ab - x) = b(\pi - x) > 0 and sin⁡x>0\sin x > 0, so the integrand is >0> 0 and In>0I_n > 0 (strict positivity). Bound: on [0, π]\intcc{0}{\pi}, x≤πx \leq \pi and a−bx≤aa - bx \leq a, so P≤πnann!P \leq \frac{\pi^n a^n}{n!} and In≤π (πa)nn!I_n \leq \pi\,\frac{(\pi a)^n}{n!}, which tends to 00 (the factorial beats the geometric term: it is the general term of the convergent exponential series, cf. Example 11.12); in particular In<1I_n < 1 for large nn.
  2. Expand xn(a−bx)n=∑j=0n(nj)an−j(−b)jxn+jx^n(a - bx)^n = \sum_{j=0}^{n} \binom nj a^{n-j} (-b)^j x^{n+j}: so P=1n!∑jcjxn+jP = \frac{1}{n!}\sum_j c_j x^{n+j} with integer cjc_j. Then P(k)(0)=0P^{(k)}(0) = 0 for k<nk < n (valuation) and, for n≤k≤2nn \leq k \leq 2n, P(k)(0)=k!n!ck−nP^{(k)}(0) = \frac{k!}{n!} c_{k-n}, an integer since n!∣k!n! \mid k!. Moreover P(π−x)=P(x)P(\pi - x) = P(x) (substitute: π−x\pi - x swaps the factors, using a−b(π−x)=bxa - b(\pi - x) = bx), so P(k)(π)=±P(k)(0)P^{(k)}(\pi) = \pm P^{(k)}(0): integers as well.
  3. With Q=P−P′′+P(4)−…Q = P - P'' + P^{(4)} - \dots (finite: PP has degree 2n2n): Q+Q′′=PQ + Q'' = P, and

    (Q′sin⁡x−Qcos⁡x)′=(Q+Q′′)sin⁡x=Psin⁡x.\bigl(Q'\sin x - Q\cos x\bigr)' = (Q + Q'')\sin x = P\sin x .

    Hence In=[Q′sin⁡x−Qcos⁡x]0π=Q(π)+Q(0)I_n = \bigl[Q'\sin x - Q\cos x\bigr]_0^{\pi} = Q(\pi) + Q(0), a sum of values P(2k)P^{(2k)} at 00 and π\pi: an integer by (2).

  4. For nn large, InI_n is an integer with 0<In<10 < I_n < 1: impossible. The assumption π=ab\pi = \frac ab fails: π\pi is irrational.

Exercise 15.10 ★★★

Let ff be C1C^1 on [a, b]\intcc{a}{b}. Prove the Riemann–Lebesgue-type limit

∫abf(t)sin⁡(λt)  ⁣dt→λ→+∞0\int_a^b f(t)\sin(\lambda t)\,\dd t \xrightarrow[\lambda \to +\infty]{} 0

by integrating by parts. Then prove it again for ff merely continuous, by uniform approximation with step functions (Theorem 15.3).

Solution

Solution of Exercise 15.10.

C1C^1 case: by parts,

∫abf(t)sin⁡λt  ⁣dt=[−f(t)cos⁡λtλ]ab+1λ∫abf′(t)cos⁡λt  ⁣dt,\int_a^b f(t)\sin\lambda t\,\dd t = \Bigl[-f(t)\frac{\cos\lambda t}{\lambda}\Bigr]_a^b + \frac{1}{\lambda}\int_a^b f'(t)\cos\lambda t\,\dd t,

bounded in absolute value by 2sup⁡∣f∣+(b−a)sup⁡∣f′∣λ→0\frac{2\sup\abs f + (b - a)\sup\abs{f'}} {\lambda} \to 0.

Continuous case: let ε>0\varepsilon > 0 and pick a step function φ\varphi with ∣f−φ∣≤ε\abs{f - \varphi} \leq \varepsilon (Theorem 15.3 provides φ≤f≤ψ\varphi \leq f \leq \psi with gap ≤ε\leq\varepsilon; take φ\varphi). Then

∣∫abfsin⁡λt∣≤∫ab∣f−φ∣+∣∫abφsin⁡λt∣≤(b−a)ε+∑i∣ci∣ ∣∫xi−1xisin⁡λt  ⁣dt∣,\Bigl|\int_a^b f\sin\lambda t\Bigr| \leq \int_a^b \abs{f - \varphi} + \Bigl|\int_a^b \varphi \sin\lambda t\Bigr| \leq (b-a)\varepsilon + \sum_i \abs{c_i}\,\Bigl|\int_{x_{i-1}}^{x_i} \sin\lambda t\,\dd t\Bigr| ,

and each ∣∫sin⁡λt∣=∣cos⁡λxi−1−cos⁡λxiλ∣≤2λ\bigl|\int \sin \lambda t\bigr| = \bigl|\frac{\cos\lambda x_{i-1} - \cos\lambda x_i}{\lambda}\bigr| \leq \frac{2}{\lambda}: the second term tends to 00. Hence the limsup is ≤(b−a)ε\leq (b-a)\varepsilon for every ε\varepsilon: the limit is 00.

Exercise 15.11 ★★

(Mean value theorem for integrals) Let f,gf, g be continuous on [a, b]\intcc{a}{b} with g≥0g \geq 0. Prove that there exists c∈[a, b]c \in \intcc{a}{b} with

∫abf(t) g(t)  ⁣dt=f(c)∫abg(t)  ⁣dt,\int_a^b f(t)\,g(t)\,\dd t = f(c)\int_a^b g(t)\,\dd t ,

and show by an example that the hypothesis g≥0g \geq 0 cannot be dropped.

Solution

Solution of Exercise 15.11.

Let m=min⁡fm = \min f and M=max⁡fM = \max f, attained by the extreme value theorem. Since g≥0g \geq 0: m g≤fg≤M gm\,g \leq fg \leq M\,g, so by monotonicity

m∫abg  ≤  ∫abfg  ≤  M∫abg.m \int_a^b g \;\leq\; \int_a^b fg \;\leq\; M \int_a^b g .

If ∫abg=0\int_a^b g = 0: strict positivity (Theorem 15.7 (4)) forces g≡0g \equiv 0, both sides vanish, and any cc works. Otherwise t=∫fg∫gt = \frac{\int fg}{\int g} lies in [m, M]=f([a, b])\intcc{m}{M} = f(\intcc{a}{b}) (Theorems 13.13 and 13.10), so t=f(c)t = f(c) for some cc.

Sign matters: on [−1, 1]\intcc{-1}{1} with f(t)=g(t)=tf(t) = g(t) = t: ∫fg=∫−11t2=23\int fg = \int_{-1}^1 t^2 = \frac23, while f(c)∫−11t  ⁣dt=0f(c)\int_{-1}^1 t\,\dd t = 0 for every cc.

Exercise 15.12 ★★★

(Moments force zeros) Let ff be continuous on [a, b]\intcc{a}{b} with

∫abf(t) tk  ⁣dt=0for k=0,1,…,n.\int_a^b f(t)\,t^k\,\dd t = 0 \qquad \text{for } k = 0, 1, \dots, n .

Prove that ff vanishes at n+1n + 1 distinct points of (a, b)\intoo{a}{b}. (If ff changes sign only at z1<⋯<zmz_1 < \dots < z_m with m≤nm \leq n, integrate ff against P(t)=(t−z1)⋯(t−zm)P(t) = (t - z_1) \cdots (t - z_m) and use strict positivity.)

Solution

Solution of Exercise 15.12.

If f≡0f \equiv 0 the claim is vacuous (every point is a zero). So assume f≢0f \not\equiv 0 and suppose it has at most nn distinct zeros in (a, b)\intoo{a}{b}; let z1<⋯<zmz_1 < \dots < z_m (m≤nm \leq n) be those zeros where ff changes sign (possibly none). Set P(t)=∏i=1m(t−zi)P(t) = \prod_{i=1}^{m}(t - z_i) (empty product =1= 1), of degree m≤nm \leq n. On each subinterval cut by the ziz_i, both ff and PP have constant sign, and both flip sign when crossing some ziz_i: the product fPfP has one constant sign on all of (a, b)\intoo{a}{b}. Being continuous, not identically zero, of constant sign, it has ∣∫abfP∣>0\bigl|\int_a^b fP\bigr| > 0 (strict positivity applied to ∣fP∣\abs{fP}). But ∫fP\int f P is a linear combination of the moments ∫f tk\int f\,t^k, k≤nk \leq n, all zero: contradiction. Hence ff has at least n+1n + 1 distinct zeros in (a, b)\intoo{a}{b}.

Remark 15.24 (Perspectives inside this volume)

Three chapters ahead lean directly on this one. Chapter 16 carries the integral remainder — the sharpest of the three Taylor formulas is an integration by parts iterated nn times. Chapter 17 converts the framing of sums by integrals into the decisive test for ∑n−α\sum n^{-\alpha}, and its weekend problem refines that framing into Euler’s constant. Chapter 24 makes the integral geometric: the length of a parametrized arc is ∫x′(t)2+y′(t)2  ⁣dt\int \sqrt{x'(t)^2 + y'(t)^2}\,\dd t, an integral of a continuous function on a segment — precisely the object built here, no improper theory needed. The single most reused fact will be the humblest: ∣∫f∣≤(b−a)sup⁡∣f∣\bigl|\int f\bigr| \leq (b - a)\sup\abs{f}, the inequality that turns every pointwise estimate into an integral estimate.

15.6 Problem: The integral irrationality machine

Problem 15.1

Weekend problem — e\eu and π2\pi^2 are irrational, e\eu to six decimals, and 227>π\frac{22}{7} > \pi with proof

One mechanism powers this whole problem: an expression that must be a positive integer, yet is provably smaller than 11, cannot exist. Exercise 15.9 (Niven) ran it once to prove π∉Q\pi \notin \Q; here we industrialize it. The machine needs three parts: an integrality input (endpoint values of well-chosen polynomials), a smallness input (a factor 1n!\frac{1}{n!} crushing the integral), and a bridge (integration by parts) connecting them. We prove that e\eu is irrational and compute it with certified error, prove Legendre’s sharper theorem that π2\pi^2 is irrational, and end with the most charming integral in analysis: ∫01x4(1−x)41+x2 ⁣dx=227−π\int_0^1 \frac{x^4(1 - x)^4}{1 + x^2}\dd x = \frac{22}{7} - \pi, which brackets π\pi by hand.

Part I — Fuel.

  1. Prove that c nn!→0\dfrac{c^{\,n}}{n!} \to 0 for every fixed c>0c > 0 (beyond n≥2cn \geq 2c, each step at least halves the term).
  2. (Beta integrals) Prove, by induction on mm with integration by parts:

    ∫01x k (1−x) m  ⁣dx=k! m!(k+m+1)!(k,m∈N).\int_0^1 x^{\,k}\,(1 - x)^{\,m}\,\dd x = \frac{k!\,m!}{(k + m + 1)!} \qquad (k, m \in \N).
  3. Deduce ∫01(x(1−x))n ⁣dx=1(2n+1)(2nn)\displaystyle\int_0^1 \bigl(x(1-x)\bigr)^n \dd x = \frac{1}{(2n+1)\binom{2n}{n}}, and — comparing with the bound x(1−x)≤14x(1 - x) \leq \frac14 — the estimate (2nn)≥4n2n+1\binom{2n}{n} \geq \dfrac{4^n}{2n+1}, matching (2nn)1/n→4\binom{2n}{n}^{1/n} \to 4 from Problem 11.1.
  4. Prove the smallness lemma used twice below: for every continuous g>0g > 0 on (0, 1)\intoo{0}{1},

    0<∫01(x(1−x))ng(x)  ⁣dx≤sup⁡[0, 1]∣g∣4 n.0 < \int_0^1 \bigl(x(1-x)\bigr)^n g(x)\,\dd x \leq \frac{\sup_{\intcc{0}{1}}\abs g}{4^{\,n}} .

Part II — e\eu: irrationality, then six decimals. Set An=∫01xnex ⁣dxA_n = \displaystyle\int_0^1 x^n \eu^x \dd x.

  1. Compute A0A_0 and A1A_1, prove the recurrence An=e−nAn−1A_n = \eu - n A_{n-1}, and the bounds 0<An≤en+10 < A_n \leq \dfrac{\eu}{n+1}.
  2. Show by induction that An=αn+βneA_n = \alpha_n + \beta_n \eu with αn,βn∈Z\alpha_n, \beta_n \in \Z.
  3. Deduce that e\eu is irrational (if e=pq\eu = \frac pq, then qAnq A_n is an integer trapped in (0, 1)\intoo{0}{1} for nn large). Compare with the proof of Exercise 11.9: same punchline, different fuel.
  4. Prove, by induction and integration by parts, the exact remainder formula

    e=∑k=0n1k!+Rn,Rn=1n!∫01(1−t)n e t  ⁣dt,1(n+1)!≤Rn≤e(n+1)!.\eu = \sum_{k=0}^{n} \frac{1}{k!} + R_n, \qquad R_n = \frac{1}{n!}\int_0^1 (1 - t)^n\,\eu^{\,t}\,\dd t, \qquad \frac{1}{(n+1)!} \leq R_n \leq \frac{\eu}{(n+1)!} .
  5. Take n=9n = 9: bound R9R_9 using e<2.75\eu < 2.75 (from b2=2.75b_2 = 2.75 in Example 11.12), evaluate the sum, and conclude the certified bracketing 2.7182818≤e≤2.71828232.7182818 \leq \eu \leq 2.7182823 — six decimals, e≈2.718282\eu \approx 2.718282, with proof.

Part III — Legendre’s theorem: π2\pi^2 is irrational. Let f(x)=xn(1−x)nn!f(x) = \dfrac{x^n (1 - x)^n}{n!}, and suppose π2=ab\pi^2 = \frac ab with a,b∈N∗a, b \in \N^*.

  1. Show f(1−x)=f(x)f(1 - x) = f(x) and 0<f≤14n n!0 < f \leq \dfrac{1}{4^n\,n!} on (0, 1)\intoo{0}{1}.
  2. Show that f(k)(0)f^{(k)}(0) and f(k)(1)f^{(k)}(1) are integers for every k≥0k \geq 0 (expand xn(1−x)nx^n(1-x)^n with integer coefficients; k!n!∈Z\frac{k!}{n!} \in \Z for k≥nk \geq n; then use the symmetry).
  3. Define

    G=b n∑k=0n(−1)k π2n−2kf(2k).G = b^{\,n} \sum_{k=0}^{n} (-1)^k\, \pi^{2n - 2k} f^{(2k)} .

    Show that G(0)G(0) and G(1)G(1) are integers (each bnπ2n−2k=a n−k b kb^n \pi^{2n-2k} = a^{\,n-k}\,b^{\,k}).

  4. Verify the telescoping G′′+π2G=bnπ2n+2f=π2anfG'' + \pi^2 G = b^n \pi^{2n+2} f = \pi^2 a^n f, then

     ⁣d ⁣dx(G′(x)sin⁡πx−π G(x)cos⁡πx)=π2anf(x)sin⁡πx.\frac{\dd}{\dd x}\Bigl(G'(x)\sin \pi x - \pi\,G(x)\cos \pi x\Bigr) = \pi^2 a^n f(x)\sin \pi x .
  5. Integrate over [0, 1]\intcc{0}{1} and conclude

    π an∫01f(x)sin⁡(πx)  ⁣dx=G(0)+G(1)∈Z,\pi\,a^n \int_0^1 f(x)\sin(\pi x)\,\dd x = G(0) + G(1) \in \Z ,

    a positive integer bounded by πan4n n!\dfrac{\pi a^n}{4^n\,n!}.

  6. Conclude with question 1 that π2\pi^2 is irrational (Legendre, 1794), and that this strengthens Exercise 15.9: why does the irrationality of π2\pi^2 imply that of π\pi, and not conversely?

Part IV — Understanding the machine.

  1. Locate the two opposing forces (integrality of endpoint data; analytic smallness of the integral) and the bridge, in Parts II and III. Then explain why the factor 1n!\frac{1}{n!} in ff is the crux: if it is removed, integrality survives, but which inequality dies, and for which claimed fractions ab\frac ab does the proof then fail?
  2. Effectivity: suppose someone claims π2=ab\pi^2 = \frac ab with a≤10a \leq 10. Show that the contradiction already lands at n=7n = 7: compute π (10/4)7/7!≈0.38<1\pi\,(10/4)^7/7! \approx 0.38 < 1. The machine does not merely refute; it refutes by a fixed, computable stage.
  3. Sanity check the bridge unconditionally: prove by two integrations by parts that

    ∫01x(1−x)sin⁡(πx)  ⁣dx=4π3,\int_0^1 x(1 - x)\sin(\pi x)\,\dd x = \frac{4}{\pi^3},

    and reconcile with question 14 at n=1n = 1 (keep π2\pi^2 symbolic: the telescoped identity reads π3∫01f1sin⁡πx=−(f1′′(0)+f1′′(1))=4\pi^3 \int_0^1 f_1 \sin \pi x = -(f_1''(0) + f_1''(1)) = 4).

  4. What makes ex\eu^x and sin⁡πx\sin \pi x eligible as the machine’s kernels? Identify the property (each satisfies a linear differential equation with constant coefficients, so repeated integration by parts cycles back to the start), and name the frontier: the same machine, refined by Hermite and Lindemann, proves e\eu and π\pi transcendental — beyond this volume.

Part V — 227\frac{22}{7} versus π\pi, and the moral.

  1. Establish the polynomial division

    x4(1−x)41+x2=x6−4x5+5x4−4x2+4−41+x2,\frac{x^4(1-x)^4}{1 + x^2} = x^6 - 4x^5 + 5x^4 - 4x^2 + 4 - \frac{4}{1 + x^2},

    and deduce the celebrated identity

    ∫01x4(1−x)41+x2  ⁣dx=227−π.\int_0^1 \frac{x^4 (1-x)^4}{1 + x^2}\,\dd x = \frac{22}{7} - \pi .
  2. The integrand is positive: conclude π<227\pi < \frac{22}{7}. Then, bounding 11+x2\frac{1}{1+x^2} between 12\frac12 and 11 and using ∫01(x(1−x))4=1630\int_0^1 (x(1-x))^4 = \frac{1}{630} (question 3), prove

    227−1630  ≤  π  ≤  227−11260,\frac{22}{7} - \frac{1}{630} \;\leq\; \pi \;\leq\; \frac{22}{7} - \frac{1}{1260} ,

    i.e. 3.14126≤π≤3.142073.14126 \leq \pi \leq 3.14207: two correct decimals, by hand.

  3. Generalize: dividing x4m(1−x)4mx^{4m}(1-x)^{4m} by 1+x21 + x^2, show the remainder is the constant (−4)m(-4)^m (work modulo x2+1x^2 + 1: (1−x)2≡−2x(1-x)^2 \equiv -2x), deduce rationals rmr_m with

    ∣π−rm∣≤4 1−5m,\abs{\pi - r_m} \leq 4^{\,1 - 5m} ,

    and check that m=1m = 1 reproduces question 20–21.

  4. Confront these rationals with the approximation theory of Problem 14.1: compute ∣π−227∣≈1.26⋅10−3\abs{\pi - \frac{22}{7}} \approx 1.26\cdot10^{-3} against the Dirichlet guarantee 149\frac{1}{49}, and quote ∣π−355113∣≈2.7⋅10−7\abs{\pi - \frac{355}{113}} \approx 2.7\cdot10^{-7} against 11132≈7.8⋅10−5\frac{1}{113^2} \approx 7.8\cdot10^{-5}: exceptionally good rational approximations exist for π\pi — consistent, since π\pi is not known to be badly approximable.
  5. (The integer trap, abstracted) Prove the lemma that unifies everything: if x∈Rx \in \R and there exist integers an,bna_n, b_n with 0<∣an+bnx∣→00 < \abs{a_n + b_n x} \to 0, then xx is irrational. List its instances in this problem, in Exercise 15.9, in Exercise 11.9, and in Problem 14.1.
  6. Synthesis, one sentence each: (i) the machine’s three parts and where each lives in the toolbox of this chapter; (ii) what the integral contributes that the mean value theorem of Problem 14.1 could not; (iii) the inventory of results extracted (two irrationalities, one six-decimal constant, one bracketing of π\pi, one binomial bound); (iv) the frontier (Hermite, Lindemann; and the same trap, run on ζ(2)\zeta(2) and ζ(3)\zeta(3), in twentieth-century arithmetic).
Solution

Solution of Problem 15.1.

1. Let N=⌈2c⌉N = \lceil 2c \rceil. For n≥Nn \geq N: cn+1/(n+1)!cn/n!=cn+1≤12\frac{c^{n+1}/(n+1)!}{c^n/n!} = \frac{c}{n+1} \leq \frac12, so 0<cnn!≤cNN! 2−(n−N)→00 < \frac{c^n}{n!} \leq \frac{c^N}{N!}\,2^{-(n - N)} \to 0: squeeze.

2. Fix kk; induction on mm. For m=0m = 0: ∫01xk=1k+1=k! 0!(k+1)!\int_0^1 x^k = \frac{1}{k+1} = \frac{k!\,0!}{(k+1)!}. Step, by parts (u=(1−x)mu = (1-x)^m, v′=xkv' = x^k):

∫01xk(1−x)m ⁣dx=[xk+1k+1(1−x)m]01+mk+1∫01xk+1(1−x)m−1 ⁣dx=mk+1⋅(k+1)! (m−1)!(k+m+1)!,\int_0^1 x^k (1-x)^m \dd x = \Bigl[\frac{x^{k+1}}{k+1}(1-x)^m\Bigr]_0^1 + \frac{m}{k+1}\int_0^1 x^{k+1}(1-x)^{m-1}\dd x = \frac{m}{k+1}\cdot\frac{(k+1)!\,(m-1)!}{(k+m+1)!} ,

which is k! m!(k+m+1)!\frac{k!\,m!}{(k+m+1)!}.

3. k=m=nk = m = n: ∫01(x(1−x))n=(n!)2(2n+1)!=1(2n+1)(2nn)\int_0^1 (x(1-x))^n = \frac{(n!)^2}{(2n+1)!} = \frac{1}{(2n+1)\binom{2n}{n}}. Since x(1−x)≤14x(1-x) \leq \frac14 on [0, 1]\intcc{0}{1}, the integral is ≤4−n\leq 4^{-n}, whence (2nn)≥4n2n+1\binom{2n}{n} \geq \frac{4^n}{2n+1} — consistent with (2nn)1/n→4\binom{2n}{n}^{1/n} \to 4 (Problem 11.1).

4. The integrand is continuous, ≥0\geq 0, and positive on (0, 1)\intoo{0}{1}, hence not identically zero: its integral is >0> 0 (Theorem 15.7 (4)). Upper bound: (x(1−x))n≤4−n(x(1-x))^n \leq 4^{-n} and g≤sup⁡∣g∣g \leq \sup\abs g, then monotonicity.

5. A0=e−1A_0 = \eu - 1; A1=[xex]01−∫01ex=e−(e−1)=1A_1 = [x\eu^x]_0^1 - \int_0^1 \eu^x = \eu - (\eu - 1) = 1. By parts: An=[xnex]01−n∫01xn−1ex=e−nAn−1A_n = [x^n \eu^x]_0^1 - n\int_0^1 x^{n-1}\eu^x = \eu - n A_{n-1}. Bounds: the integrand is positive, so An>0A_n > 0; and ex≤e\eu^x \leq \eu gives An≤e∫01xn=en+1A_n \leq \eu\int_0^1 x^n = \frac{\eu}{n+1}.

6. A0=−1+1⋅eA_0 = -1 + 1\cdot\eu. If An−1=αn−1+βn−1eA_{n-1} = \alpha_{n-1} + \beta_{n-1}\eu with integer entries, then

An=e−nαn−1−nβn−1e=(−n αn−1)⏟αn+(1−n βn−1)⏟βn e,A_n = \eu - n\alpha_{n-1} - n\beta_{n-1}\eu = \underbrace{(-n\,\alpha_{n-1})}_{\alpha_n} + \underbrace{(1 - n\,\beta_{n-1})}_{\beta_n}\,\eu ,

both integers.

7. If e=pq\eu = \frac pq: qAn=qαn+pβn∈Zq A_n = q\alpha_n + p\beta_n \in \Z, and 0<qAn≤qen+1<10 < qA_n \leq \frac{q\eu}{n+1} < 1 for nn large: an integer strictly between 00 and 11 — impossible. So e∉Q\eu \notin \Q. In Exercise 11.9 the trapped integer was q! pq−q! aqq!\,\frac pq - q!\,a_q; here it is qAnqA_n: same trap, integral fuel.

8. n=0n = 0: R0=∫01et=e−1R_0 = \int_0^1 \eu^t = \eu - 1, so e=1+R0\eu = 1 + R_0. By parts (u=etu = \eu^t, v=−(1−t)n+1n+1v = -\frac{(1-t)^{n+1}}{n+1}):

Rn=1n!([−(1−t)n+1n+1et]01+1n+1∫01(1−t)n+1et)=1(n+1)!+Rn+1,R_n = \frac{1}{n!}\Bigl(\Bigl[-\frac{(1-t)^{n+1}}{n+1} \eu^t\Bigr]_0^1 + \frac{1}{n+1}\int_0^1 (1-t)^{n+1}\eu^t\Bigr) = \frac{1}{(n+1)!} + R_{n+1} ,

so the formula propagates from nn to n+1n + 1. Bounds: 1≤et≤e1 \leq \eu^t \leq \eu on [0, 1]\intcc{0}{1} and ∫01(1−t)n=1n+1\int_0^1 (1-t)^n = \frac{1}{n+1} give 1(n+1)!≤Rn≤e(n+1)!\frac{1}{(n+1)!} \leq R_n \leq \frac{\eu}{(n+1)!}.

9. ∑k=091k!=986410362880=2.71828152…\sum_{k=0}^{9} \frac{1}{k!} = \frac{986410}{362880} = 2.71828152\dots, and

110!=2.76⋅10−7≤R9≤2.7510!=7.58⋅10−7,\frac{1}{10!} = 2.76\cdot10^{-7} \leq R_9 \leq \frac{2.75}{10!} = 7.58\cdot10^{-7} ,

so 2.7182818≤e≤2.71828232.7182818 \leq \eu \leq 2.7182823: with proof, e=2.718282\eu = 2.718282 to six decimals (true value 2.7182818…2.7182818\dots).

10. f(1−x)=(1−x)nxnn!=f(x)f(1 - x) = \frac{(1-x)^n x^n}{n!} = f(x). On (0, 1)\intoo{0}{1}: 0<x(1−x)≤140 < x(1-x) \leq \frac14, so 0<f≤14n n!0 < f \leq \frac{1}{4^n\,n!}.

11. xn(1−x)n=∑j=0n(−1)j(nj) xn+jx^n(1-x)^n = \sum_{j=0}^{n} (-1)^j\binom nj\,x^{n+j}, so f=1n!∑jcj xn+jf = \frac{1}{n!}\sum_j c_j\,x^{n+j} with cj∈Zc_j \in \Z. Hence f(k)(0)=0f^{(k)}(0) = 0 for k<nk < n or k>2nk > 2n, and for n≤k≤2nn \leq k \leq 2n: f(k)(0)=k!n! ck−nf^{(k)}(0) = \frac{k!}{n!}\, c_{k-n}, an integer because n!∣k!n! \mid k!. The symmetry gives f(k)(1)=(−1)kf(k)(0)∈Zf^{(k)}(1) = (-1)^k f^{(k)}(0) \in \Z.

12. bnπ2n−2k=bn(ab) n−k=a n−k b k∈Zb^n \pi^{2n-2k} = b^n\bigl(\frac ab\bigr)^{\,n-k} = a^{\,n-k}\,b^{\,k} \in \Z, so G(0)=∑k(−1)kan−kbkf(2k)(0)G(0) = \sum_k (-1)^k a^{n-k} b^k f^{(2k)}(0) and likewise G(1)G(1) are integers by question 11.

13. In π2G+G′′\pi^2 G + G'', the term kk of π2G\pi^2 G carries π2n−2k+2f(2k)\pi^{2n-2k+2}f^{(2k)} and the term j=k−1j = k - 1 of G′′G'' carries (−1)k−1π2n−2k+2f(2k)(-1)^{k-1}\pi^{2n-2k+2}f^{(2k)}: everything cancels except k=0k = 0 in the first sum and j=nj = n in the second, i.e.

G′′+π2G=bn(π2n+2f+(−1)nf(2n+2))=bnπ2n+2f=π2anfG'' + \pi^2 G = b^n\bigl(\pi^{2n+2} f + (-1)^n f^{(2n+2)}\bigr) = b^n \pi^{2n+2} f = \pi^2 a^n f

(ff has degree 2n2n, so f(2n+2)=0f^{(2n+2)} = 0; and bnπ2n=anb^n\pi^{2n} = a^n). Then

(G′sin⁡πx−πGcos⁡πx)′=(G′′+π2G)sin⁡πx=π2anfsin⁡πx.\bigl(G'\sin\pi x - \pi G\cos\pi x\bigr)' = (G'' + \pi^2 G)\sin \pi x = \pi^2 a^n f\sin\pi x .

14. Integrating over [0, 1]\intcc{0}{1}:

π2an∫01fsin⁡(πx)  ⁣dx=[G′sin⁡πx−πGcos⁡πx]01=π(G(1)+G(0)),\pi^2 a^n \int_0^1 f\sin(\pi x)\,\dd x = \bigl[G'\sin\pi x - \pi G\cos\pi x\bigr]_0^1 = \pi\bigl(G(1) + G(0)\bigr) ,

so πan∫01fsin⁡πx=G(0)+G(1)∈Z\pi a^n \int_0^1 f\sin\pi x = G(0) + G(1) \in \Z. On (0, 1)\intoo{0}{1}, f>0f > 0 and sin⁡πx>0\sin\pi x > 0: the left side is positive, so G(0)+G(1)≥1G(0) + G(1) \geq 1; and sin⁡≤1\sin \leq 1 with question 10 bounds it by πan4n n!\frac{\pi a^n}{4^n\,n!}.

15. By question 1 (with c=a4c = \frac a4), πan4nn!→0\frac{\pi a^n}{4^n n!} \to 0: for large nn it is <1< 1, contradicting G(0)+G(1)≥1G(0) + G(1) \geq 1. So no fraction ab\frac ab equals π2\pi^2: Legendre’s theorem. If π\pi were rational, π2\pi^2 would be too: so π∉Q\pi \notin \Q — and the implication only runs this way (2\sqrt2 is irrational with rational square), which is why π2∉Q\pi^2 \notin \Q is strictly stronger than Exercise 15.9.

16. Integrality: questions 11–12 (endpoint derivatives); smallness: questions 10 and 1; bridge: questions 13–14 (the telescoped double integration by parts). Without 1n!\frac{1}{n!}, the endpoint data remain integers (even more easily), but the bound becomes πan4n\frac{\pi a^n}{4^n}, which tends to 00 only when a<4a < 4 — and every candidate has a=b π2>9a = b\,\pi^2 > 9. The factorial is exactly what outruns the geometric growth ana^n: no factorial, no theorem.

17. For a≤10a \leq 10, the integer G(0)+G(1)G(0) + G(1) is positive and at most π (10/4)n/n!\pi\,(10/4)^n/n!. At n=7n = 7: 2.57=610.35…2.5^7 = 610.35\dots, so the bound is π×610.355040≈0.38<1\frac{\pi \times 610.35}{5040} \approx 0.38 < 1 (at n=6n = 6 it is still 1.071.07): the contradiction lands at the seventh stage, explicitly.

18. Two integrations by parts:

∫01x(1−x)sin⁡πx  ⁣dx=1π∫01(1−2x)cos⁡πx  ⁣dx=2π2∫01sin⁡πx  ⁣dx=2π2⋅2π=4π3\int_0^1 x(1-x)\sin\pi x\,\dd x = \frac1\pi\int_0^1 (1 - 2x)\cos\pi x\,\dd x = \frac{2}{\pi^2}\int_0^1 \sin\pi x\,\dd x = \frac{2}{\pi^2}\cdot\frac{2}{\pi} = \frac{4}{\pi^3}

(the bracket terms vanish: x(1−x)x(1-x) at 0,10, 1, and sin⁡πx\sin\pi x at 0,10, 1). Symbolically, the n=1n = 1 telescope (no assumption on π\pi) reads π3∫01f1sin⁡πx=−(f1′′(0)+f1′′(1))\pi^3\int_0^1 f_1\sin\pi x = -(f_1''(0) + f_1''(1)) with f1=x(1−x)f_1 = x(1 - x), f1′′=−2f_1'' = -2: right side 44 — the two computations agree.

19. ex\eu^x solves y′=yy' = y and sin⁡πx\sin\pi x solves y′′=−π2yy'' = -\pi^2 y: linear equations with constant coefficients, so integration by parts cycles the kernel back to itself and keeps all boundary data inside Z+Ze\Z + \Z\eu (resp. integer polynomials in π2\pi^2). That closure property is what the machine needs. Refined with kernels adapted to several points at once, the same mechanism yields Hermite’s theorem (e\eu transcendental, 1873) and Lindemann’s (π\pi transcendental, 1882) — beyond this volume.

20. Polynomial division (or multiply back and check):

x4(1−x)4=(x6−4x5+5x4−4x2+4)(1+x2)−4.x^4(1-x)^4 = (x^6 - 4x^5 + 5x^4 - 4x^2 + 4)(1 + x^2) - 4 .

Integrating the displayed identity divided by 1+x21 + x^2:

∫01x4(1−x)41+x2 ⁣dx=(17−46+1−43+4)−4arctan⁡1=227−π,\int_0^1 \frac{x^4(1-x)^4}{1+x^2}\dd x = \Bigl(\frac17 - \frac46 + 1 - \frac43 + 4\Bigr) - 4\arctan 1 = \frac{22}{7} - \pi ,

using 17−23+1−43+4=3+17\frac17 - \frac23 + 1 - \frac43 + 4 = 3 + \frac17 and arctan⁡1=π4\arctan 1 = \frac\pi4.

21. The integrand is continuous, positive on (0, 1)\intoo{0}{1}: the integral is >0> 0, so π<227\pi < \frac{22}{7}. Moreover 12≤11+x2≤1\frac12 \leq \frac{1}{1+x^2} \leq 1 on [0, 1]\intcc{0}{1} and ∫01(x(1−x))4=(4!)29!=1630\int_0^1 (x(1-x))^4 = \frac{(4!)^2}{9!} = \frac{1}{630} (question 3):

11260≤227−π≤1630⟹3.14126<227−1630≤π≤227−11260<3.14207.\frac{1}{1260} \leq \frac{22}{7} - \pi \leq \frac{1}{630} \quad\Longrightarrow\quad 3.14126 < \frac{22}{7} - \frac{1}{630} \leq \pi \leq \frac{22}{7} - \frac{1}{1260} < 3.14207 .

22. Modulo x2+1x^2 + 1: x2≡−1x^2 \equiv -1, so x4m=(x2)2m≡1x^{4m} = (x^2)^{2m} \equiv 1 and (1−x)2=1−2x+x2≡−2x(1 - x)^2 = 1 - 2x + x^2 \equiv -2x, hence (1−x)4m≡(−2x)2m=4m(x2)m≡(−4)m(1-x)^{4m} \equiv (-2x)^{2m} = 4^m (x^2)^m \equiv (-4)^m: the remainder is the constant (−4)m(-4)^m, and the quotient QmQ_m has integer coefficients (division by a monic integer polynomial). Dividing the identity by 1+x21 + x^2 and integrating:

Jm:=∫01(x(1−x))4m1+x2 ⁣dx=sm+(−4)m π4,sm=∫01Qm∈Q.J_m := \int_0^1 \frac{(x(1-x))^{4m}}{1+x^2}\dd x = s_m + (-4)^m\,\frac{\pi}{4}, \qquad s_m = \int_0^1 Q_m \in \Q .

Solving for π\pi: with rm=(−1)m+1 4 1−msm∈Qr_m = (-1)^{m+1}\,4^{\,1-m} s_m \in \Q,   ∣π−rm∣=4 1−mJm≤4 1−m⋅4−4m=4 1−5m\;\abs{\pi - r_m} = 4^{\,1-m} J_m \leq 4^{\,1-m}\cdot4^{-4m} = 4^{\,1-5m}. For m=1m = 1: s1=227s_1 = \frac{22}{7}, r1=227r_1 = \frac{22}{7}, bound 4−4=12564^{-4} = \frac{1}{256} — questions 20–21 again.

23. ∣π−227∣=227−π≈1.26⋅10−3\bigl|\pi - \frac{22}{7}\bigr| = \frac{22}{7} - \pi \approx 1.26\cdot10^{-3}, sixteen times better than the order-22 benchmark 172≈2.0⋅10−2\frac{1}{7^2} \approx 2.0\cdot10^{-2} guaranteed by Dirichlet (Problem 14.1, question 4); and ∣π−355113∣≈2.7⋅10−7\bigl|\pi - \frac{355}{113}\bigr| \approx 2.7\cdot10^{-7} beats 11132≈7.8⋅10−5\frac{1}{113^2} \approx 7.8\cdot10^{-5} by a factor ≈300\approx 300. No contradiction with anything proved: Liouville inequalities lower-bound approximation errors only for algebraic numbers, and no such bound for π\pi is available at this level — π\pi is free to be approximated spectacularly well.

24. Lemma: suppose x=pqx = \frac pq and 0<∣an+bnx∣→00 < \abs{a_n + b_n x} \to 0. Then ∣an+bnx∣=∣qan+pbn∣q\abs{a_n + b_n x} = \frac{\abs{q a_n + p b_n}}{q}, with qan+pbnq a_n + p b_n a nonzero integer (nonzero because the absolute value is >0> 0): so ∣an+bnx∣≥1q\abs{a_n + b_n x} \geq \frac1q for every nn, contradicting the convergence to 00. Instances: question 7 (x=ex = \eu, an=αna_n = \alpha_n, bn=βnb_n = \beta_n); Exercise 11.9 (x=ex = \eu again, with aq=−q!∑k≤q1k!a_q = -q!\sum_{k \leq q}\frac{1}{k!}, bq=q!b_q = q!); and Problem 14.1, question 1, is its geometric form. In Part III and in Exercise 15.9 the trap runs inside the contradiction: assuming rationality converts an expression into an integer, which the analysis then squeezes into (0, 1)\intoo{0}{1} — the same principle, transposed.

25. (i) Integrality lives in the endpoint calculus of polynomials (questions 6, 11–12), smallness in the bounds sup⁡\sup-monotonicity gives (questions 4, 10), the bridge in integration by parts (questions 8, 13–14) — all three are theorems of this chapter. (ii) The integral supplies what the mean value theorem could not: an exact identity between the analytic object and the arithmetic data (equality, not just an inequality with an unknown cc), which is why the machine reaches π2\pi^2 while Problem 14.1 reached only approximation exponents. (iii) Extracted: e∉Q\eu \notin \Q, π2∉Q\pi^2 \notin \Q (hence π∉Q\pi \notin \Q), e=2.718282\eu = 2.718282 certified, 227−1630≤π≤227−11260\frac{22}{7} - \frac{1}{630} \leq \pi \leq \frac{22}{7} - \frac{1}{1260}, and (2nn)≥4n2n+1\binom{2n}{n} \geq \frac{4^n}{2n+1}. (iv) Frontier: Hermite and Lindemann push the same machine to transcendence; and Apéry (1979) ran the integer trap on ζ(3)\zeta(3) — the machine is still producing twentieth-century mathematics.

Terms defined in this chapter

See all 395 terms in the glossary