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 (xi1,xi)\intoo{x_{i-1}}{x_i} (values at the nodes are unconstrained). Its integral is

abφ=i=1nci(xixi1),\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(xi1,xi)t \in \intoo{x_{i-1}}{x_i} into a subdivision replaces the term ci(xixi1)c_i(x_i - x_{i-1}) by ci(txi1)+ci(xit)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: cidic_i \leq d_i on each piece gives ciΔidiΔi\sum c_i \Delta_i \leq \sum d_i \Delta_i (lengths Δi0\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>baδn > \frac{b - a}{\delta} pieces). On each closed piece [xi1,xi]\intcc{x_{i-1}}{x_i}, ff attains a minimum mim_i and a maximum MiM_i (Theorem 13.13), and MimiεM_i - m_i \leq \varepsilon (the two extremal points are within δ\delta). Define φ=mi\varphi = m_i and ψ=Mi\psi = M_i on (xi1,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 ψφε(ba)\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,

03t ⁣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 φ=k1n\varphi = \frac{k-1}{n} and ψ=kn\psi = \frac kn on the kk-th piece, with

01φ=k=1nk1n1n=n12n,01ψ=k=1nkn1n=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 n12nI(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: fg    abfabgf \leq g \implies \int_a^b f \leq \int_a^b g; and abfabf(ba)supf\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, f0f \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 φffψf\varphi_f \leq f \leq \psi_f and φggψg\varphi_g \leq g \leq \psi_g with gaps ε\leq \varepsilon (Theorem 15.3). For λ0\lambda \geq 0, φf+λφgf+λ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 ff\abs{\int f} \leq \int \abs f comes from fff-\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 fm2f \geq \frac m2; the step function worth m2\frac m2 there and 00 elsewhere is f\leq f, so fmη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.

02x1 ⁣dx=01(1x) ⁣dx+12(x1) ⁣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πsint ⁣dt=0πsint ⁣dtπ2πsint ⁣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πsint ⁣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 f0f \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 ff\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 aIa \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 x0Ix_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 tx0δ\abs{t - x_0} \leq \delta. For 0<hδ0 < \abs{h} \leq \delta (and x0+hIx_0 + h \in I), Chasles gives

F(x0+h)F(x0)hf(x0)=1hx0x0+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 1hhε=ε\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 (GF)=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 ttt \mapsto -t sends a0f\int_{-a}^{0} f to 0af-\int_0^a f, so

aaf(t) ⁣dt=0;if f is even,aaf=20af.\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 11t3cost1+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 ππt2cost ⁣dt=20πt2cost ⁣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+cosx\displaystyle\int_0^{\pi/2} \frac{\dd x}{1 + \cos x}. The half-angle identity 1+cosx=2cos2x21 + \cos x = 2\cos^2\frac x2 turns the integrand into 12(1+tan2x2)\frac{1}{2}\bigl(1 + \tan^2\frac x2\bigr), which is exactly the derivative of tanx2\tan\frac x2:

0π/2 ⁣dx1+cosx=[tanx2]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=tanx2t = \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)=0xet2 ⁣dt.F(x) = \int_0^x \eu^{-t^2}\,\dd t .

No combination of classical functions has derivative et2\eu^{-t^2} (a theorem of Liouville, admitted); yet FF exists, is C1C^1 with F(x)=ex2>0F'(x) = \eu^{-x^2} > 0, strictly increasing, odd (substitute ttt \mapsto -t), and bounded: for x1x \geq 1,

F(x)F(1)=1xet2 ⁣dt1xet ⁣dte1,F(x) - F(1) = \int_1^x \eu^{-t^2}\dd t \leq \int_1^x \eu^{-t}\dd t \leq \eu^{-1} ,

so FF(1)+e11+e1F \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 ⁣dxxx2et2 ⁣dt=2xex4ex2\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}, 1211+t1\frac12 \leq \frac{1}{1+t} \leq 1, so

12(n+1)=1201tn ⁣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 (RnR_n \sim a multiple of 1n\frac1n, in fact Rn12nR_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 1baabf\frac{1}{b-a}\int_a^b f. For the arch of the sine:

1π0πsint ⁣dt=1π[cost]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: sinc=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}:

    abuv=[uv]ababuv.\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)=uv+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

01tet ⁣dt=[tet]0101et ⁣dt=e(e1)=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=sintx = \sin t (t[0,π2]t \in \intcc{0}{\frac\pi2}):

011x2 ⁣dx=0π/2costcost ⁣dt=0π/21+cos2t2 ⁣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:  ⁣dxx2=lnx2+C\int \frac{\dd x}{x - 2} = \ln\abs{x - 2} + C on each side of 22 separately — writing ln(x2)\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 banf(a+kban)\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(lnt)2 ⁣dt\int_1^x (\ln t)^2\,\dd t? Guess a primitive of the form tP(lnt)t\,P(\ln t) with PP polynomial and differentiate:

(tP(lnt))=P(lnt)+P(lnt).\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)=u22u+2P(u) = u^2 - 2u + 2 (matching coefficients downward from u2u^2). Hence

1x(lnt)2 ⁣dt=[t((lnt)22lnt+2)]1x=x(lnx)22xlnx+2x2,\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 lnt\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π/2excosx ⁣dxI = \int_0^{\pi/2} \eu^x \cos x\,\dd x. Integrate by parts twice, differentiating the trigonometric factor each time:

I=[exsinx]0π/20π/2exsinx ⁣dx=eπ/2J,J=[excosx]0π/2+0π/2excosx ⁣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π/21II = \eu^{\pi/2} - 1 - I, whence

I=eπ/212.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}, cosbx\cos bx, sinbx\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=bank=0n1f(a+kban)nabf(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+kban)f(a + k\frac{b-a}{n}) on the kk-th subinterval. Given ε>0\varepsilon > 0, uniform continuity (Heine) provides δ\delta; for n>baδ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

abfSn=ab(fφn)(ba)ε.\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=1nk=1n11+k/nn01 ⁣dx1+x=ln2\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 limnk=1nn(n+k)2\lim_{n\to\infty} \sum_{k=1}^{n} \dfrac{n}{(n+k)^2}. Calibrate:

k=1nn(n+k)2=1nk=1nn2(n+k)2=1nk=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 ⁣dxx24\displaystyle\int_0^1 \frac{\dd x}{x^2 - 4} (partial fractions, Chapter 9); 1elnt ⁣dt\displaystyle\int_1^{\eu} \ln t \,\dd t; 0πtsint ⁣dt\displaystyle\int_0^{\pi} t \sin t\, \dd t.

Solution

Solution of Exercise 15.1.

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

01 ⁣dxx24=14[lnx2lnx+2]01=14(ln13ln1)=ln34.\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=lntv = \ln t): 1elnt ⁣dt=[tlnt]1e1e ⁣dt=e(e1)=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=sintu' = \sin t, v=tv = t): 0πtsint ⁣dt=[tcost]0π+0πcost ⁣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π/2cos3t ⁣dt\displaystyle\int_0^{\pi/2} \cos^3 t\, \dd t (write cos3=cos(1sin2)\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=1212 ⁣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=sintu = \sin t: 0π/2cos3t ⁣dt=0π/2(1sin2t)cost ⁣dt=[sintsin3t3]0π/2=113=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=1nk=1n11+(k/n)2u_n = \dfrac1n \sum_{k=1}^{n} \dfrac{1}{1 + (k/n)^2}: a Riemann sum of x11+x2x \mapsto \frac{1}{1+x^2} on [0,1]\intcc{0}{1}, so un01 ⁣dx1+x2=arctan1=π4u_n \to \int_0^1 \frac{\dd x}{1+x^2} = \arctan 1 = \dfrac\pi4.

lnvn=1nk=1nlnn+kn=1nk=1nln(1+kn)01ln(1+x) ⁣dx=[(1+x)ln(1+x)x]01=2ln21\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 vne2ln21=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 limn01xnf(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=supfM = \sup \abs f and ε(0,1)\varepsilon \in \intoo{0}{1}. Cut at 1ε1 - \varepsilon:

01xnf01εxnf+1ε1xnfM(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ε)n0(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) ⁣2abf2abg2,\Bigl(\int_a^b fg\Bigr)^{\!2} \leq \int_a^b f^2 \cdot \int_a^b g^2 ,

by expanding ab(f+λg)20\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+λ2g20Q(\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 000 \leq 0. Otherwise QQ is a genuine quadratic, everywhere 0\geq 0: its discriminant is 0\leq 0, i.e. (fg)2f2g2\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 1x0xf(t) ⁣dt1T0Tf\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, 0r<T0 \leq r < T (n=x/Tn = \lfloor x/T \rfloor). Chasles:

0xf=n0Tf+nTnT+rf,nTnT+rfTsup[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 1x0xf=nTx1T0Tf+O(1x)\frac 1x \int_0^x f = \frac{nT}{x}\cdot\frac 1T \int_0^T f + O\bigl(\frac 1x\bigr), and nTx1\frac{nT}{x} \to 1: the limit is 1T0Tf\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=01f12=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π/2sinnt ⁣dtW_n = \int_0^{\pi/2} \sin^n t\,\dd t.

  1. Prove the recurrence nWn=(n1)Wn2n W_n = (n-1) W_{n-2} (n2n \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+1Wn1\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=sintu' = \sin t, v=sinn1tv = \sin^{n-1} t:

    Wn=[costsinn1t]0π/2+(n1)0π/2cos2tsinn2t ⁣dt=(n1)(Wn2Wn),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=(n1)Wn2nW_n = (n-1)W_{n-2}. From W0=π2W_0 = \frac\pi2, W1=1W_1 = 1:

    W2p=(2p1)(2p3)1(2p)(2p2)2π2=(2p)!4p(p!)2π2,W2p+1=(2p)(2p2)2(2p+1)(2p1)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<sint<10 < \sin t < 1, so sinn+1<sinn\sin^{n+1} < \sin^n and (Wn)(W_n) is (strictly) decreasing, positive. Sandwiching with the recurrence:

    nn+1=Wn+1Wn1Wn+1Wn1    Wn+1Wn1.\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=an1a_n = a_{n-1} by the recurrence (n+1)Wn+1=nWn1(n+1)W_{n+1} = nW_{n-1}, so an=a0=1W1W0=π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,bNa, b \in \N^*, and set, for nn to be chosen,

P(x)=xn(abx)nn!,In=0πP(x)sinx ⁣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 knk \geq n; and P(πx)=P(x)P(\pi - x) = P(x).)
  3. Set Q=PP+P(4)Q = P - P'' + P^{(4)} - \dots (a finite sum). Check that (QsinxQcosx)=Psinx\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, abx=b(abx)=b(πx)>0a - bx = b(\frac ab - x) = b(\pi - x) > 0 and sinx>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 abxaa - 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(abx)n=j=0n(nj)anj(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 nk2nn \leq k \leq 2n, P(k)(0)=k!n!cknP^{(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 ab(π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=PP+P(4)Q = P - P'' + P^{(4)} - \dots (finite: PP has degree 2n2n): Q+Q=PQ + Q'' = P, and

    (QsinxQcosx)=(Q+Q)sinx=Psinx.\bigl(Q'\sin x - Q\cos x\bigr)' = (Q + Q'')\sin x = P\sin x .

    Hence In=[QsinxQcosx]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 2supf+(ba)supfλ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λtabfφ+abφsinλt(ba)ε+icixi1xisinλ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λxi1cosλ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 (ba)ε\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 g0g \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 g0g \geq 0 cannot be dropped.

Solution

Solution of Exercise 15.11.

Let m=minfm = \min f and M=maxfM = \max f, attained by the extreme value theorem. Since g0g \geq 0: mgfgMgm\,g \leq fg \leq M\,g, so by monotonicity

mabg    abfg    Mabg.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 g0g \equiv 0, both sides vanish, and any cc works. Otherwise t=fggt = \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 mnm \leq n, integrate ff against P(t)=(tz1)(tzm)P(t) = (t - z_1) \cdots (t - z_m) and use strict positivity.)

Solution

Solution of Exercise 15.12.

If f0f \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 (mnm \leq n) be those zeros where ff changes sign (possibly none). Set P(t)=i=1m(tzi)P(t) = \prod_{i=1}^{m}(t - z_i) (empty product =1= 1), of degree mnm \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 ftk\int f\,t^k, knk \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(ba)supf\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(1x)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 cnn!0\dfrac{c^{\,n}}{n!} \to 0 for every fixed c>0c > 0 (beyond n2cn \geq 2c, each step at least halves the term).
  2. (Beta integrals) Prove, by induction on mm with integration by parts:

    01xk(1x)m ⁣dx=k!m!(k+m+1)!(k,mN).\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(1x))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(1x)14x(1 - x) \leq \frac14 — the estimate (2nn)4n2n+1\binom{2n}{n} \geq \dfrac{4^n}{2n+1}, matching (2nn)1/n4\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(1x))ng(x) ⁣dxsup[0,1]g4n.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=enAn1A_n = \eu - n A_{n-1}, and the bounds 0<Anen+10 < A_n \leq \dfrac{\eu}{n+1}.
  2. Show by induction that An=αn+βneA_n = \alpha_n + \beta_n \eu with αn,βnZ\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(1t)net ⁣dt,1(n+1)!Rne(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.7182818e2.71828232.7182818 \leq \eu \leq 2.7182823 — six decimals, e2.718282\eu \approx 2.718282, with proof.

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

  1. Show f(1x)=f(x)f(1 - x) = f(x) and 0<f14nn!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 k0k \geq 0 (expand xn(1x)nx^n(1-x)^n with integer coefficients; k!n!Z\frac{k!}{n!} \in \Z for knk \geq n; then use the symmetry).
  3. Define

    G=bnk=0n(1)kπ2n2kf(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π2n2k=ankbkb^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

    πan01f(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 πan4nn!\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 a10a \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(1x)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 π301f1sinπ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(1x)41+x2=x64x5+5x44x2+441+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(1x)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(1x))4=1630\int_0^1 (x(1-x))^4 = \frac{1}{630} (question 3), prove

    2271630    π    22711260,\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(1x)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: (1x)22x(1-x)^2 \equiv -2x), deduce rationals rmr_m with

    πrm415m,\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 π2271.26103\abs{\pi - \frac{22}{7}} \approx 1.26\cdot10^{-3} against the Dirichlet guarantee 149\frac{1}{49}, and quote π3551132.7107\abs{\pi - \frac{355}{113}} \approx 2.7\cdot10^{-7} against 111327.8105\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 xRx \in \R and there exist integers an,bna_n, b_n with 0<an+bnx00 < \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=2cN = \lceil 2c \rceil. For nNn \geq N: cn+1/(n+1)!cn/n!=cn+112\frac{c^{n+1}/(n+1)!}{c^n/n!} = \frac{c}{n+1} \leq \frac12, so 0<cnn!cNN!2(nN)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=(1x)mu = (1-x)^m, v=xkv' = x^k):

01xk(1x)m ⁣dx=[xk+1k+1(1x)m]01+mk+101xk+1(1x)m1 ⁣dx=mk+1(k+1)!(m1)!(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(1x))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(1x)14x(1-x) \leq \frac14 on [0,1]\intcc{0}{1}, the integral is 4n\leq 4^{-n}, whence (2nn)4n2n+1\binom{2n}{n} \geq \frac{4^n}{2n+1} — consistent with (2nn)1/n4\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(1x))n4n(x(1-x))^n \leq 4^{-n} and gsupgg \leq \sup\abs g, then monotonicity.

5. A0=e1A_0 = \eu - 1; A1=[xex]0101ex=e(e1)=1A_1 = [x\eu^x]_0^1 - \int_0^1 \eu^x = \eu - (\eu - 1) = 1. By parts: An=[xnex]01n01xn1ex=enAn1A_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 exe\eu^x \leq \eu gives Ane01xn=en+1A_n \leq \eu\int_0^1 x^n = \frac{\eu}{n+1}.

6. A0=1+1eA_0 = -1 + 1\cdot\eu. If An1=αn1+βn1eA_{n-1} = \alpha_{n-1} + \beta_{n-1}\eu with integer entries, then

An=enαn1nβn1e=(nαn1)αn+(1nβn1)βne,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βnZq A_n = q\alpha_n + p\beta_n \in \Z, and 0<qAnqen+1<10 < qA_n \leq \frac{q\eu}{n+1} < 1 for nn large: an integer strictly between 00 and 11 — impossible. So eQ\eu \notin \Q. In Exercise 11.9 the trapped integer was q!pqq!aqq!\,\frac pq - q!\,a_q; here it is qAnqA_n: same trap, integral fuel.

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

Rn=1n!([(1t)n+1n+1et]01+1n+101(1t)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: 1ete1 \leq \eu^t \leq \eu on [0,1]\intcc{0}{1} and 01(1t)n=1n+1\int_0^1 (1-t)^n = \frac{1}{n+1} give 1(n+1)!Rne(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.76107R92.7510!=7.58107,\frac{1}{10!} = 2.76\cdot10^{-7} \leq R_9 \leq \frac{2.75}{10!} = 7.58\cdot10^{-7} ,

so 2.7182818e2.71828232.7182818 \leq \eu \leq 2.7182823: with proof, e=2.718282\eu = 2.718282 to six decimals (true value 2.71828182.7182818\dots).

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

11. xn(1x)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!jcjxn+jf = \frac{1}{n!}\sum_j c_j\,x^{n+j} with cjZc_j \in \Z. Hence f(k)(0)=0f^{(k)}(0) = 0 for k<nk < n or k>2nk > 2n, and for nk2nn \leq k \leq 2n: f(k)(0)=k!n!cknf^{(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π2n2k=bn(ab)nk=ankbkZb^n \pi^{2n-2k} = b^n\bigl(\frac ab\bigr)^{\,n-k} = a^{\,n-k}\,b^{\,k} \in \Z, so G(0)=k(1)kankbkf(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 π2n2k+2f(2k)\pi^{2n-2k+2}f^{(2k)} and the term j=k1j = k - 1 of GG'' carries (1)k1π2n2k+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

(Gsinπ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}:

π2an01fsin(πx) ⁣dx=[Gsinπ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 πan01fsinπ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 sin1\sin \leq 1 with question 10 bounds it by πan4nn!\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 π2Q\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 a10a \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.352.5^7 = 610.35\dots, so the bound is π×610.3550400.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(1x)sinπx ⁣dx=1π01(12x)cosπx ⁣dx=2π201sinπx ⁣dx=2π22π=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(1x)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 π301f1sinπx=(f1(0)+f1(1))\pi^3\int_0^1 f_1\sin\pi x = -(f_1''(0) + f_1''(1)) with f1=x(1x)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(1x)4=(x64x5+5x44x2+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(1x)41+x2 ⁣dx=(1746+143+4)4arctan1=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 1723+143+4=3+17\frac17 - \frac23 + 1 - \frac43 + 4 = 3 + \frac17 and arctan1=π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 1211+x21\frac12 \leq \frac{1}{1+x^2} \leq 1 on [0,1]\intcc{0}{1} and 01(x(1x))4=(4!)29!=1630\int_0^1 (x(1-x))^4 = \frac{(4!)^2}{9!} = \frac{1}{630} (question 3):

11260227π16303.14126<2271630π22711260<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: x21x^2 \equiv -1, so x4m=(x2)2m1x^{4m} = (x^2)^{2m} \equiv 1 and (1x)2=12x+x22x(1 - x)^2 = 1 - 2x + x^2 \equiv -2x, hence (1x)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(1x))4m1+x2 ⁣dx=sm+(4)mπ4,sm=01QmQ.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+141msmQr_m = (-1)^{m+1}\,4^{\,1-m} s_m \in \Q,   πrm=41mJm41m44m=415m\;\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 44=12564^{-4} = \frac{1}{256} — questions 20–21 again.

23. π227=227π1.26103\bigl|\pi - \frac{22}{7}\bigr| = \frac{22}{7} - \pi \approx 1.26\cdot10^{-3}, sixteen times better than the order-22 benchmark 1722.0102\frac{1}{7^2} \approx 2.0\cdot10^{-2} guaranteed by Dirichlet (Problem 14.1, question 4); and π3551132.7107\bigl|\pi - \frac{355}{113}\bigr| \approx 2.7\cdot10^{-7} beats 111327.8105\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+bnx00 < \abs{a_n + b_n x} \to 0. Then an+bnx=qan+pbnq\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+bnx1q\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!kq1k!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: eQ\eu \notin \Q, π2Q\pi^2 \notin \Q (hence πQ\pi \notin \Q), e=2.718282\eu = 2.718282 certified, 2271630π22711260\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.