Mathematics · Book 4 · Bachelor Year 2

University Mathematics — Year 2

University Mathematics — Year 2 · Bachelor Year 2

10Sequences and Series of Functions

When functions converge to a function, which properties survive the passage to the limit? Pointwise convergence preserves almost nothing; uniform convergence — convergence in the sup norm — preserves continuity, integrals on segments and, with one twist, derivatives. This chapter proves the three transfer theorems, their series versions, and crowns them with the Weierstrass approximation theorem, proved by Bernstein’s beautiful probabilistic polynomials.

10.1 Pointwise and uniform convergence

Definition 10.1

Let fn,f ⁣:XRf_n, f \colon X \to \R (or C\C, or a normed space), XX any set. (fn)(f_n) converges to ff pointwise when fn(x)f(x)f_n(x) \to f(x) for every xx; uniformly when

fnf=supxXfn(x)f(x)n0.\norm{f_n - f}_\infty = \sup_{x \in X}\, \abs{f_n(x) - f(x)} \xrightarrow[n \to \infty]{} 0 .

Uniform implies pointwise; on C([a,b])C(\intcc{a}{b}), uniform convergence is exactly convergence in the Banach space (C([a,b]),)\bigl(C(\intcc{a}{b}), \norm\cdot_\infty\bigr) of Chapter 5.

Example 10.2

On [0,1]\intcc{0}{1}, fn(x)=xnf_n(x) = x^n converges pointwise to the discontinuous limit f=1{1}f = \mathbf{1}_{\{1\}}; the convergence is not uniform: fnffn(11n)=(11n)ne10\norm{f_n - f}_\infty \geq f_n\bigl(1 - \tfrac1n\bigr) = (1 - \tfrac1n)^n \to \eu^{-1} \neq 0. On [0,a]\intcc{0}{a} with a<1a < 1 it is uniform (sup=an0\sup = a^n \to 0): uniformity is a property of the domain as much as of the sequence.

The sequence xn on [0, 1]: the graphs sag toward 0 but all must climb to 1 at x = 1 — the sup distance to the discontinuous pointwise limit never shrinks below a constant.
The sequence xnx^n on [0,1]\intcc{0}{1}: the graphs sag toward 00 but all must climb to 11 at x=1x = 1 — the sup distance to the discontinuous pointwise limit never shrinks below a constant.

Example 10.3 (Two limits that refuse to commute)

The whole chapter is about interchanging limits, so here is the smallest possible failure. Let an,m=nn+ma_{n,m} = \dfrac{n}{n+m} for n,m1n, m \geq 1. Then

limm(limnan,m)=limm1=1,limn(limman,m)=limn0=0:\lim_{m\to\infty}\Bigl(\lim_{n\to\infty}a_{n,m}\Bigr) = \lim_{m\to\infty} 1 = 1, \qquad \lim_{n\to\infty}\Bigl(\lim_{m\to\infty}a_{n,m}\Bigr) = \lim_{n\to\infty} 0 = 0 :

both iterated limits exist and they differ. Every transfer theorem of this chapter is a licence to commute two limits — limn\lim_n with limxa\lim_{x\to a} (continuity), with \int (integration), with  ⁣d ⁣dx\frac{\dd}{\dd x} (differentiation) — and uniform convergence is precisely the fee that makes the commutation legal. Closing insight: whenever a “proof” silently swaps two limit operations, this two-line array is the counterexample to hold against it; the sliding bumps of Exercise 10.2 are the same phenomenon wearing an integral sign.

10.2 The three transfer theorems

Theorem 10.4 (Continuity)

If each fnf_n is continuous at aa and fnff_n \to f uniformly on a neighborhood of aa, then ff is continuous at aa. A uniform limit of continuous functions is continuous.

Proof. The 3ε3\varepsilon argument already used in Theorem 4.9: pick nn with fnfε\norm{f_n - f}_\infty \leq \varepsilon, then δ\delta from continuity of fnf_n at aa; for xaδ\abs{x - a} \leq \delta,

f(x)f(a)f(x)fn(x)+fn(x)fn(a)+fn(a)f(a)3ε.\abs{f(x) - f(a)} \leq \abs{f(x) - f_n(x)} + \abs{f_n(x) - f_n(a)} + \abs{f_n(a) - f(a)} \leq 3\varepsilon . \qedhere

Example 10.5 (Uniformity fails exactly where the limit breaks)

On [0,2]\intcc{0}{2}, let fn(x)=xn1+xnf_n(x) = \dfrac{x^n}{1 + x^n}. The pointwise limit is a three-piece function:

f(x)={00x<1,12x=1,11<x2,f(x) = \begin{cases} 0 & 0 \leq x < 1,\\[2pt] \tfrac12 & x = 1,\\[2pt] 1 & 1 < x \leq 2, \end{cases}

discontinuous at 11, so by Theorem 10.4 the convergence cannot be uniform on [0,2]\intcc{0}{2}. On the closed pieces avoiding the threshold it is: for 0xa<10 \leq x \leq a < 1,

sup[0,a]fn0=an1+anan0,\sup_{\intcc{0}{a}}\abs{f_n - 0} = \frac{a^n}{1 + a^n} \leq a^n \to 0 ,

and for 1<bx21 < b \leq x \leq 2,

sup[b,2]fn1=11+bnbn0,\sup_{\intcc{b}{2}}\abs{f_n - 1} = \frac{1}{1 + b^n} \leq b^{-n} \to 0 ,

both suprema computed by monotonicity of uu1+uu \mapsto \frac{u}{1+u} and of xxnx \mapsto x^n. Closing insight: the failure of uniformity is localized at the discontinuity of the limit — the same geometry as Example 10.2, and the reason the “uniform on every segment inside” discipline recurs all chapter long.

Theorem 10.6 (Integration on a segment)

If fnff_n \to f uniformly on [a,b]\intcc{a}{b}, with fnf_n piecewise continuous (ff likewise), then

abfnabf.\int_a^b f_n \longrightarrow \int_a^b f .

Proof. Linearity and the triangle inequality for integrals give

abfnabf=ab(fnf)abfnf(ba)fnf0.\Bigl|\int_a^b f_n - \int_a^b f\Bigr| = \Bigl|\int_a^b (f_n - f)\Bigr| \leq \int_a^b\abs{f_n - f} \leq (b - a)\,\norm{f_n - f}_\infty \longrightarrow 0 .

The length factor (ba)(b - a) is where compactness of the segment enters: on non-compact intervals the same estimate produces the useless bound 0\infty\cdot0, and the conclusion genuinely fails without domination — the flat bumps fn=1n1[0,n]f_n = \frac1n\mathbf 1_{\intcc{0}{n}} converge uniformly to 00 on [0,)\intco{0}{\infty} yet keep fn=1\int f_n = 1 (pitfalls remark below), and the sliding bumps of Chapter 9’s remark do the same with pointwise convergence; uniform convergence controls heights, never widths.

The bumps g_n(x) = nx\, -nx2 of : they converge to 0 at every point, but the peaks (height √n/(2 ), drifting toward 0) grow without bound — pointwise convergence with g_n_∈fty ∈fty, and ∈t_01 g_n 1/2 ≠ 0: the mass hides under the moving peak.
The bumps gn(x)=nxenx2g_n(x) = nx\,\eu^{-nx^2} of Exercise 10.1: they converge to 00 at every point, but the peaks (height n/(2e)\sim\sqrt{n/(2\eu)}, drifting toward 00) grow without bound — pointwise convergence with gn\norm{g_n}_\infty \to \infty, and 01gn120\int_0^1 g_n \to \frac12 \neq 0: the mass hides under the moving peak.

Theorem 10.7 (Differentiation)

Let fnf_n be C1C^1 on an interval II, with: (fn)(f_n') converging uniformly on II (or on every segment of II) to some gg, and (fn(x0))(f_n(x_0)) converging at one point x0x_0. Then (fn)(f_n) converges (uniformly on segments) to a C1C^1 function ff, and f=gf' = g: one may differentiate the limit.

Proof. Define f(x)=limfn(x0)+x0xgf(x) = \lim f_n(x_0) + \int_{x_0}^x g: legitimate, gg being continuous — indeed gg is the uniform limit on segments of the continuous fnf_n', so Theorem 10.4 applies, and the integral of a continuous function is well defined with, by the fundamental theorem of calculus,

f(x)=g(x)(xI):f'(x) = g(x) \qquad (x \in I) :

the candidate limit is C1C^1 with the right derivative by construction, before any convergence is proved. By the fundamental theorem again, fn(x)=fn(x0)+x0xfnf_n(x) = f_n(x_0) + \int_{x_0}^x f_n'; subtracting,

fn(x)f(x)fn(x0)limfn(x0)+xx0fng,\abs{f_n(x) - f(x)} \leq \abs{f_n(x_0) - \lim f_n(x_0)} + \abs{x - x_0}\,\norm{f_n' - g}_{\infty} ,

which tends to 00 uniformly on every segment. And ff is C1C^1 with f=gf' = g by construction.

Example 10.8 (Why the hypothesis sits on the derivatives)

Let Fn(x)=x2+1nF_n(x) = \sqrt{x^2 + \frac1n} on R\R. Each FnF_n is C1C^1 (indeed CC^\infty), and the convergence to x\abs x is uniform on all of R\R:

0Fn(x)x=(x2+1n)x2x2+1n+x=1/nx2+1n+x1/n1/n=1n.0 \leq F_n(x) - \abs x = \frac{(x^2 + \frac1n) - x^2}{\sqrt{x^2+\frac1n} + \abs x} = \frac{1/n}{\sqrt{x^2 + \frac1n} + \abs x} \leq \frac{1/n}{1/\sqrt n} = \frac{1}{\sqrt n} .

Yet the limit x\abs x is not differentiable at 00: uniform convergence of the functions, however fast, transfers no differentiability. The failure is visible on the derivatives:

Fn(x)=xx2+1n{1x>0,0x=0,1x<0,F_n'(x) = \frac{x}{\sqrt{x^2 + \frac1n}} \longrightarrow \begin{cases} 1 & x > 0,\\ 0 & x = 0,\\ -1 & x < 0, \end{cases}

a discontinuous pointwise limit, so (Fn)(F_n') cannot converge uniformly near 00 (Theorem 10.4 again). Closing insight: Theorem 10.7 deliberately assumes uniform convergence of the fnf_n', not of the fnf_n — this example is the reason.

10.3 Series of functions

Definition 10.9

A series of functions un\sum u_n converges pointwise/uniformly when its partial sums do. It converges normally (on XX) when un<\sum \norm{u_n}_\infty < \infty. Normal convergence implies uniform convergence (in the Banach space of bounded functions: Theorem 5.21), which implies pointwise; both implications are strict.

Example 10.10 (One series, three verdicts)

Take un(x)=xnnu_n(x) = \frac{x^n}{n} on [0,1)\intco{0}{1}. Pointwise: converges for every x[0,1)x \in \intco01 (comparison with the geometric series). Normal on [0,a]\intcc{0}{a}, a<1a < 1: un,[0,a]=ann\norm{u_n}_{\infty,\intcc0a} = \frac{a^n}{n}, summable. Not normal on [0,1)\intco{0}{1}: un,[0,1)=1n\norm{u_n}_{\infty,\intco01} = \frac1n, and 1n\sum\frac1n diverges. Not even uniform on [0,1)\intco{0}{1}: the remainder resists near 11,

RN(x)=n>Nxnnn=N+12NxnnNx2N2N=x2N2x112,R_N(x) = \sum_{n>N}\frac{x^n}{n} \geq \sum_{n=N+1}^{2N}\frac{x^n}{n} \geq \frac{N\,x^{2N}}{2N} = \frac{x^{2N}}{2} \xrightarrow[x\to1^-]{} \frac12 ,

so sup[0,1)RN12\sup_{\intco01}\abs{R_N} \geq \frac12 for every NN. Closing insight: all four verdicts coexist peacefully — the sum ln(1x)-\ln(1-x) is continuous on [0,1)\intco{0}{1} because continuity only needs uniformity near each point, i.e. on the segments [0,a]\intcc0a; blowing up at the edge is the sum’s right.

Theorem 10.11 (Transfer for series)

If un\sum u_n converges uniformly (e.g. normally) on the relevant set: continuity of all unu_n at aa passes to the sum; integration on a segment can be done term by term; and if un(x0)\sum u_n(x_0) converges while un\sum u_n' converges uniformly on segments, the sum is C1C^1 with derivative un\sum u_n'.

Proof. Everything is the corresponding theorem applied to the partial sums SN=nNunS_N = \sum_{n\leq N}u_n, which are finite sums of functions with the relevant regularity. Continuity: each SNS_N is continuous at aa and SNunS_N \to \sum u_n uniformly: Theorem 10.4. Integration: on the segment,

abn0un=limNabSN=limNn=0Nabun=n0abun\int_a^b \sum_{n\geq0} u_n = \lim_N \int_a^b S_N = \lim_N \sum_{n=0}^{N}\int_a^b u_n = \sum_{n\geq0}\int_a^b u_n

by Theorem 10.6 (first equality) and linearity of the integral (second). Differentiation: the SNS_N are C1C^1, SN(x0)S_N(x_0) converges, and SN=nNunS_N' = \sum_{n\leq N}u_n' converges uniformly on segments: Theorem 10.7 gives that the sum is C1C^1 with derivative limSN=un\lim S_N' = \sum u_n'.

Remark 10.12 (Common pitfalls)

Four traps, all seen in examination copies. (i) Half-checked suprema: evaluating fnf_n along a well-chosen sequence xnx_n only bounds fnf\norm{f_n - f}_\infty from below — enough to disprove uniformity (as in Example 10.2), never to prove it; to prove it, bound the sup by a computation valid for all xx. (ii) Uniformity on the wrong set: normal or uniform convergence often holds on every [a,a]\intcc{-a}{a} or [δ,)\intco\delta\infty but fails on the open union; that is no obstacle — continuity and differentiability are local, so the segment-by-segment discipline of Example 10.13 gives them on the whole open set. (iii) Integrating over non-segments: Theorem 10.6 is a statement about segments; on [0,)\intco0\infty, uniform convergence does not prevent mass from escaping to infinity (fn=1n1[0,n]f_n = \frac1n\mathbf 1_{\intcc{0}{n}} converges uniformly to 00, with fn=1\int f_n = 1) — use dominated convergence there. (iv) Differentiating the limit: Example 10.8; the derivative hypothesis is on (fn)(f_n'), and no rate of convergence of (fn)(f_n) can replace it.

Example 10.13 (The Riemann ζ\zeta function)

ζ(s)=n1ns\zeta(s) = \sum_{n\geq1} n^{-s} converges normally on every half-line [a,+)\intco{a}{+\infty}, a>1a > 1 (ns=na\norm{n^{-s}}_\infty = n^{-a}, summable): ζ\zeta is continuous on (1,+)\intoo{1}{+\infty}; differentiating termwise (the derived series lnn  ns\sum -\ln n\; n^{-s} also converges normally on [a,)\intco{a}{\infty}), ζ\zeta is C1C^1 — and, iterating, CC^\infty — with ζ(s)=lnnns\zeta'(s) = -\sum \frac{\ln n}{n^s}. Note the discipline: normal convergence is checked on sub-half-lines, never on the open (1,)\intoo{1}{\infty} itself, where it fails.

Example 10.14 (A logarithmic series, worked to the end)

Let F(x)=n1enxnF(x) = \sum_{n\geq1} \frac{\eu^{-nx}}{n} on (0,)\intoo{0}{\infty}. Each term is bounded on [δ,)\intco{\delta} \infty by enδnenδ\frac{\eu^{-n\delta}}{n} \leq \eu^{-n\delta}, a convergent geometric series: normal convergence on every [δ,)\intco\delta\infty, so FF is continuous on (0,)\intoo{0}{\infty}. The derived series enx\sum -\eu^{-nx} is likewise normally convergent on [δ,)\intco\delta\infty (enx,[δ,)=enδ\norm{\eu^{-nx}}_{\infty,\intco\delta\infty} = \eu^{-n\delta}), so FF is C1C^1 with a geometric derivative:

F(x)=n1enx=ex1ex=1ex1.F'(x) = -\sum_{n\geq1}\eu^{-nx} = \frac{-\eu^{-x}}{1 - \eu^{-x}} = \frac{-1}{\eu^{x} - 1} .

Iterating, FF is CC^\infty. Integrating FF' (both FF and xln(1ex)x \mapsto -\ln(1 - \eu^{-x}) vanish at ++\infty and have the same derivative on (0,)\intoo0\infty):

F(x)=ln(1ex),F(x) = -\ln\bigl(1 - \eu^{-x}\bigr),

the logarithmic series at t=ext = \eu^{-x}. Closing insight: as x0+x \to 0^+, F(x)=ln(x+O(x2))=ln1x+O(x)F(x) = -\ln(x + O(x^2)) = \ln\frac1x + O(x) — the series diverges logarithmically at the boundary, exactly like the harmonic series it becomes at x=0x = 0; normal convergence on [δ,)\intco\delta\infty but not on (0,)\intoo0\infty is the symptom.

Method 10.15 (Proving or disproving uniform convergence)

For fnff_n \to f pointwise on XX:

  1. Compute or bound fnf\norm{f_n - f}_\infty: study the function xfn(x)f(x)x \mapsto \abs{f_n(x) - f(x)} (derivative, monotonicity) to locate its maximum; a bound valid for all xx that tends to 00 proves uniformity.
  2. To disprove: exhibit points xnx_n with fn(xn)f(xn)↛0\abs{f_n(x_n) - f(x_n)} \not\to 0 (often xnx_n tracks the moving bump, as in Exercise 10.1); or invoke a transfer theorem in contrapositive — a discontinuous limit of continuous functions (Example 10.5), or fn↛f\int f_n \not\to \int f on a segment.
  3. For series, try normal convergence first (supun<\sum\sup\abs{u_n} < \infty); if it fails globally, test it on the sub-segments that matter (Example 10.10); if it fails everywhere, uniform convergence may still hold via the alternating remainder bound (Exercise 10.4) or summation by parts.

10.4 The Weierstrass approximation theorem

Theorem 10.16 (Weierstrass, via Bernstein)

Every continuous f ⁣:[0,1]Rf \colon \intcc{0}{1} \to \R is a uniform limit of polynomials — explicitly, of its Bernstein polynomials

Bn(f)(x)=k=0nf(kn)(nk)xk(1x)nk.B_n(f)(x) = \sum_{k=0}^{n} f\Bigl(\frac kn\Bigr)\binom nk x^k (1-x)^{n-k} .

Proof. Fix x[0,1]x \in \intcc{0}{1} and set pk(x)=(nk)xk(1x)nkp_k(x) = \binom nk x^k(1 - x)^{n-k}. Three binomial identities, obtained by evaluating (x+y)n(x + y)^n and its two xx-derivatives at y=1xy = 1 - x:

kpk=1,kkpk=nx,kk(k1)pk=n(n1)x2.\sum_k p_k = 1, \qquad \sum_k k\,p_k = nx, \qquad \sum_k k(k-1) p_k = n(n-1)x^2 .

In detail: (x+y)n=k(nk)xkynk(x+y)^n = \sum_k\binom nk x^ky^{n-k} at y=1xy = 1-x is the first; differentiating in xx,

n(x+y)n1=kk(nk)xk1ynk,n(x+y)^{n-1} = \sum_k k\binom nk x^{k-1}y^{n-k} ,

then multiplying by xx and setting y=1xy = 1 - x gives the second; differentiating twice and multiplying by x2x^2 gives the third. Expanding (knx)2=k(k1)+k(12nx)+n2x2(k - nx)^2 = k(k-1) + k(1 - 2nx) + n^2x^2 and combining the three:

k(knx)2pk=n(n1)x2+nx(12nx)+n2x2=nx(1x)n4,\sum_k (k - nx)^2 p_k = n(n-1)x^2 + nx(1 - 2nx) + n^2x^2 = nx(1 - x) \leq \frac n4 ,

the variance identity.

Now estimate, using pk=1\sum p_k = 1:

Bn(f)(x)f(x)kf(kn)f(x)pk(x)=Σnear+Σfar,\abs{B_n(f)(x) - f(x)} \leq \sum_{k} \Bigl| f\Bigl(\frac kn\Bigr) - f(x)\Bigr|\, p_k(x) = \Sigma_{\text{near}} + \Sigma_{\text{far}} ,

splitting according to knxδ\abs{\frac kn - x} \leq \delta or not. Given ε>0\varepsilon > 0, uniform continuity of ff (Heine) provides δ\delta with Σnearε\Sigma_{\text{near}} \leq \varepsilon. For the far sum, with M=fM = \norm f_\infty: by the variance identity and Chebyshev’s counting trick,

Σfar2Mknx>nδpk2Mk(knx)2pkn2δ22M4nδ21=M2nδ2n0,\Sigma_{\text{far}} \leq 2M \sum_{\abs{k - nx} > n\delta} p_k \leq 2M\,\frac{\sum_k (k - nx)^2 p_k}{n^2\delta^2} \leq \frac{2M}{4 n \delta^2} \cdot 1 = \frac{M}{2n\delta^2} \xrightarrow[n\to\infty]{} 0 ,

uniformly in xx. So Bn(f)fε+M2nδ22ε\norm{B_n(f) - f}_\infty \leq \varepsilon + \frac{M}{2n\delta^2} \leq 2\varepsilon for large nn.

Remark 10.17

By affine substitution the theorem holds on any segment [a,b]\intcc{a}{b}. It fails on R\R (a uniform limit of polynomials on R\R is a polynomial: Exercise 10.8). The probabilistic reading — Bn(f)(x)B_n(f)(x) is the expected value of ff at a binomial average, and the variance bound is Chebyshev’s inequality — is made honest in Chapter 23.

Bernstein approximation of f(x) = x2 (red), using the exact formula B_nf = x2 + x(1-x)/n of : B_1f is the chord, and each doubling of n halves the gap. Reliable but slow — the 1n saturation that Voronovskaya’s theorem (weekend problem) makes exact.
Bernstein approximation of f(x)=x2f(x) = x^2 (red), using the exact formula Bnf=x2+x(1x)nB_nf = x^2 + \frac{x(1-x)}{n} of Exercise 10.7: B1fB_1f is the chord, and each doubling of nn halves the gap. Reliable but slow — the 1n\frac1n saturation that Voronovskaya’s theorem (weekend problem) makes exact.

Remark 10.18 (Where this is used)

Weierstrass approximation is the density theorem of classical analysis: it makes C([a,b])C(\intcc ab) separable, lets one check integral identities on polynomials only (moment problems), and underlies the trigonometric version proved in the Fourier chapter via Fejér’s kernel. The weekend problem of this chapter extracts the quantitative content of Bernstein’s proof — convergence rates governed by the modulus of continuity — and then isolates what really made it work, in Korovkin’s theorem: positivity plus three test functions. The Year 3 volume generalizes the density statement to arbitrary subalgebras (Stone–Weierstrass) and to compact spaces.

Example 10.19 (Polygonal approximation, with a rate)

For LL-Lipschitz ff on [0,1]\intcc{0}{1}, let InfI_nf be the piecewise-affine interpolant at the nodes kn\frac kn. On a cell [kn,k+1n]\intcc{\frac kn}{\frac{k+1}n}, both f(x)f(x) and Inf(x)I_nf(x) lie between the extreme values that an LL-Lipschitz function can take given the two nodal values, so for xx in the cell, writing xk=knx_k = \frac kn:

Inf(x)f(x)Inf(x)f(xk)+f(xk)f(x)Lxxk+Lxxk2Ln\abs{I_nf(x) - f(x)} \leq \abs{I_nf(x) - f(x_k)} + \abs{f(x_k) - f(x)} \leq L\,\abs{x - x_k} + L\,\abs{x - x_k} \leq \frac{2L}{n}

(the interpolant is itself LL-Lipschitz on the cell: its slope is a difference quotient of ff). Hence Inff2Ln\norm{I_nf - f}_\infty \leq \frac{2L}{n}: polygonal approximation of Lipschitz functions converges at speed 1n\frac1nfaster than Bernstein’s 1n\frac{1}{\sqrt n} for the same class (weekend problem, Part II). Closing insight: the polygon interpolates but is not smooth, Bernstein is smooth but slow; there is no free lunch between regularity of the approximant and speed — a trade-off made precise by the saturation results of the weekend problem.

Remark 10.20 (Perspectives within this volume)

Uniform convergence is this book’s workhorse from here on. The power-series chapter runs entirely on normal convergence on compact subdisks — every term-by-term theorem there is a special case of this chapter’s transfer theorems. The Fourier chapter lives one floor up: its partial sums SNS_N fail exactly where this chapter warns they might (pointwise but not uniform at jumps), and its Fejér means succeed by the same 3ε3\varepsilon mechanics that proved Theorem 10.4. The differential-equations chapter defines etA\eu^{tA} by a normally convergent series and differentiates it termwise — literally Theorem 10.11 applied to matrix entries. When in doubt later in the book about “why can we do this”, the answer is usually a theorem of this chapter.

10.5 Exercises

Exercise 10.1

Study the pointwise and uniform convergence on [0,1]\intcc{0}{1}, then on [0,a]\intcc{0}{a} (a<1a < 1) or [δ,1)\intco{\delta}{1} as relevant, of:

fn(x)=x1+nx,gn(x)=nxenx2,hn(x)=xn(1xn).f_n(x) = \frac{x}{1 + nx}, \qquad g_n(x) = n x\,\eu^{-n x^2}, \qquad h_n(x) = x^n(1 - x^n).
Solution

Solution of Exercise 10.1.

fn(x)=x1+nxf_n(x) = \frac{x}{1 + nx}: pointwise limit 00 on [0,1]\intcc{0}{1}. Uniformly: fnf_n increases on [0,1]\intcc{0}{1} (derivative 1(1+nx)2>0\frac{1}{(1+nx)^2} > 0), so fn=fn(1)=11+n0\norm{f_n}_\infty = f_n(1) = \frac{1}{1+n} \to 0: uniform on [0,1]\intcc{0}{1}.

gn(x)=nxenx2g_n(x) = nx\,\eu^{-nx^2}: pointwise limit 00 (exponential beats). Sup: gn=nenx2(12nx2)g_n' = n\eu^{-nx^2}(1 - 2nx^2) vanishes at xn=12nx_n = \frac{1}{\sqrt{2n}}, where gn(xn)=n2e1/2g_n(x_n) = \sqrt{\frac n2}\,\eu^{-1/2} \to \infty: not uniform on [0,1]\intcc{0}{1} — but uniform on [δ,1)\intco{\delta}{1}, since there gn(x)nenδ20g_n(x) \leq n\,\eu^{-n\delta^2} \to 0.

hn(x)=xn(1xn)h_n(x) = x^n(1 - x^n): pointwise limit 00 on [0,1]\intcc{0}{1} (both factors; at x=1x = 1, hn=0h_n = 0). Sup: with u=xn[0,1]u = x^n \in \intcc{0}{1}, u(1u)14u(1-u) \leq \frac14 attained at u=12u = \frac12, i.e. x=21/n(0,1)x = 2^{-1/n} \in \intoo{0}{1}: hn=14↛0\norm{h_n}_\infty = \frac14 \not\to 0: not uniform on [0,1]\intcc{0}{1}; uniform on [0,a]\intcc{0}{a} (supan0\sup \leq a^n \to 0).

Exercise 10.2

Prove that 01gn↛01limgn\displaystyle\int_0^1 g_n \not\to \int_0^1 \lim g_n for the gng_n of Exercise 10.1, and reconcile with Theorem 10.6.

Solution

Solution of Exercise 10.2.

01nxenx2 ⁣dx=[12enx2]01=1en212\int_0^1 nx\,\eu^{-nx^2}\dd x = \bigl[-\tfrac12 \eu^{-nx^2}\bigr]_0^1 = \frac{1 - \eu^{-n}}{2} \to \frac12, while 01limgn=0\int_0^1 \lim g_n = 0. No contradiction: Theorem 10.6 requires uniform convergence on the segment, which fails here (the bump of height n\sim\sqrt n slides toward 00).

Exercise 10.3

Prove that S(x)=n1xnn2S(x) = \sum_{n\geq1} \dfrac{x^n}{n^2} is continuous on [1,1]\intcc{-1}{1}, and that SS is C1C^1 on (1,1)\intoo{-1}{1} with S(x)=ln(1x)xS'(x) = -\frac{\ln(1-x)}{x} for 0<x<10 < \abs x < 1.

Solution

Solution of Exercise 10.3.

Normal convergence on [1,1]\intcc{-1}{1}: xn/n2=1n2\norm{x^n/n^2}_\infty = \frac{1}{n^2}, summable: SS is continuous there (Theorem 10.11).

Derivative: the derived series xn1n\sum \frac{x^{n-1}}{n} converges normally on every [a,a]\intcc{-a}{a}, a<1a < 1 (sup=an1n\sup = \frac{a^{n-1}}{n}): SS is C1C^1 on (1,1)\intoo{-1}{1} with

S(x)=n1xn1n=1xn1xnn=ln(1x)x(0<x<1),S'(x) = \sum_{n\geq1} \frac{x^{n-1}}{n} = \frac1x \sum_{n\geq1} \frac{x^n}{n} = -\frac{\ln(1 - x)}{x} \qquad (0 < \abs x < 1),

the last identity being the Year 1 logarithmic series (re-derived honestly in Chapter 11).

Exercise 10.4 ★★

Let F(x)=n0(1)nn+xF(x) = \sum_{n \geq 0} \dfrac{(-1)^n}{n + x} on (0,+)\intoo{0}{+\infty}. Prove uniform (not normal) convergence on [δ,)\intco{\delta}{\infty} via the alternating series remainder bound, continuity, and the functional equation F(x)+F(x+1)=1xF(x) + F(x + 1) = \frac1x.

Solution

Solution of Exercise 10.4.

For fixed x>0x > 0 the series is alternating with 1n+x0\frac{1}{n + x} \downarrow 0: pointwise convergence, and the remainder bound RN(x)1N+1+x1N+1\abs{R_N(x)} \leq \frac{1}{N + 1 + x} \leq \frac{1}{N+1} is uniform on [δ,)\intco{\delta}{\infty} (indeed on (0,)\intoo{0}{\infty}): uniform convergence. (Not normal: (1)nn+x=1n+δ\norm{\frac{(-1)^n}{n+x}}_\infty = \frac{1}{n + \delta}, divergent.) Continuity follows from Theorem 10.11.

Functional equation: reindex F(x+1)F(x + 1) with m=n+1m = n + 1:

F(x+1)=n0(1)nn+1+x=m1(1)m1m+x,F(x+1) = \sum_{n\geq0} \frac{(-1)^n}{n + 1 + x} = \sum_{m\geq1}\frac{(-1)^{m-1}}{m+x} ,

so, isolating the m=0m = 0 term of F(x)F(x),

F(x)+F(x+1)=1x+m1(1)m+(1)m1m+x=1x.F(x) + F(x+1) = \frac{1}{x} + \sum_{m\geq1} \frac{(-1)^m + (-1)^{m-1}}{m+x} = \frac1x .

Exercise 10.5 ★★

(Dini) Let fn ⁣:KRf_n \colon K \to \R be continuous on a compact metric space, with fnff_n \to f pointwise, ff continuous, and (fn(x))(f_n(x)) decreasing in nn for each xx. Prove that the convergence is uniform. (Given ε\varepsilon, the open sets Un={x:fn(x)f(x)<ε}U_n = \{x : f_n(x) - f(x) < \varepsilon\} increase and cover KK; extract a finite subcover — Theorem 4.20.)

Solution

Solution of Exercise 10.5.

Set gn=fnf0g_n = f_n - f \geq 0 (decreasing in nn, by hypothesis; the limit is 00 pointwise); each gng_n is continuous. Fix ε>0\varepsilon > 0 and let Un={x:gn(x)<ε}U_n = \{x : g_n(x) < \varepsilon\}: open (preimage of an open set), increasing (gn+1gng_{n+1} \leq g_n), and covering KK (pointwise convergence). By Borel–Lebesgue (Theorem 4.20), finitely many Un1UnkU_{n_1} \subseteq \dots \subseteq U_{n_k} cover KK: so K=UnkK = U_{n_k}, i.e. gnkε\norm{g_{n_k}}_\infty \leq \varepsilon, and by monotonicity gnε\norm{g_n}_\infty \leq \varepsilon for all nnkn \geq n_k: uniform convergence. (Monotonicity is essential: the sliding bumps of Exercise 10.2 converge pointwise on a compact without uniformity.)

Exercise 10.6 ★★

Prove that limn01nf(x)1+n2x2 ⁣dx=π2f(0)\displaystyle\lim_{n\to\infty} \int_0^1 \frac{n\,f(x)}{1 + n^2x^2}\,\dd x = \frac{\pi}{2} f(0) for every continuous ff on [0,1]\intcc{0}{1}. (Substitute u=nxu = nx; split off f(0)f(0); dominate.)

Solution

Solution of Exercise 10.6.

Substitute u=nxu = nx:

01nf(x)1+n2x2 ⁣dx=0nf(u/n)1+u2 ⁣du.\int_0^1 \frac{n f(x)}{1 + n^2x^2}\dd x = \int_0^n \frac{f(u/n)}{1 + u^2}\,\dd u .

The integrands hn(u)=f(u/n)1+u21unh_n(u) = \frac{f(u/n)}{1+u^2}\mathbf{1}_{u \leq n} converge pointwise to f(0)1+u2\frac{f(0)}{1+u^2} (continuity of ff at 00) and are dominated by f1+u2\frac{\norm f_\infty}{1 + u^2}, integrable on [0,)\intco{0}{\infty}: dominated convergence (Theorem 9.6) gives the limit

0f(0)1+u2 ⁣du=π2f(0).\int_0^\infty \frac{f(0)}{1 + u^2}\dd u = \frac{\pi}{2} f(0) .

(The kernels concentrate at 00: an approximate identity.)

Exercise 10.7 ★★

Compute the Bernstein polynomials of f(x)=x2f(x) = x^2 explicitly and verify the uniform error Bnff=O(1n)\norm{B_n f - f}_\infty = O\bigl(\frac1n\bigr) predicted by the proof of Theorem 10.16 — here exactly x(1x)n\frac{x(1 - x)}{n} at each point.

Solution

Solution of Exercise 10.7.

For f(x)=x2f(x) = x^2, use the second binomial identity family from the proof: kk2pk=n(n1)x2+nx\sum_k k^2 p_k = n(n-1)x^2 + nx. Hence

Bn(f)(x)=kk2n2pk=n(n1)x2+nxn2=x2+x(1x)n:B_n(f)(x) = \sum_k \frac{k^2}{n^2}\,p_k = \frac{n(n-1)x^2 + nx}{n^2} = x^2 + \frac{x(1 - x)}{n} :

Bn(f)f=x(1x)nB_n(f) - f = \frac{x(1-x)}{n}, of sup norm 14n=O(1n)\frac{1}{4n} = O\bigl(\frac1n\bigr), as predicted.

Exercise 10.8 ★★

Prove that if polynomials PnP_n converge uniformly on all of R\R to ff, then ff is a polynomial. (For large m,nm, n, PnPmP_n - P_m is a bounded polynomial on R\R, hence constant; so the sequence stabilizes modulo constants.)

Solution

Solution of Exercise 10.8.

For ε=1\varepsilon = 1 there is NN with PnPm,R1\norm{P_n - P_m}_{\infty, \R} \leq 1 for m,nNm, n \geq N. A polynomial bounded on R\R is constant (a nonconstant one tends to ±\pm\infty): PnPm=cn,mP_n - P_m = c_{n,m}, constants. So for nNn \geq N: Pn=PN+cnP_n = P_N + c_n with cn=Pn(0)PN(0)c_n = P_n(0) - P_N(0) convergent (pointwise convergence at 00). Hence f=limPn=PN+limcnf = \lim P_n = P_N + \lim c_n: a polynomial.

Exercise 10.9 ★★★

(A continuous, nowhere differentiable function — guided) Let φ\varphi be the distance to the nearest integer (11-periodic, φ=12\norm{\varphi}_\infty = \frac12, 11-Lipschitz) and

W(x)=n=0(34) ⁣nφ(4nx).W(x) = \sum_{n=0}^{\infty} \Bigl(\frac{3}{4}\Bigr)^{\!n} \varphi(4^n x) .

Prove: (a) WW is continuous on R\R (normal convergence); (b) for every xx and every mm, choosing hm=±124mh_m = \pm\frac12\cdot 4^{-m} with the sign making φ\varphi affine on the segment from 4mx4^m x to 4m(x+hm)4^m(x + h_m), the difference quotient satisfies

W(x+hm)W(x)hm3mn<m3n3m+12m\Bigl|\frac{W(x + h_m) - W(x)}{h_m}\Bigr| \geq 3^m - \sum_{n<m} 3^n \geq \frac{3^m + 1}{2} \xrightarrow[m\to\infty]{} \infty

(terms n>mn > m vanish by periodicity; term n=mn = m contributes exactly 3m3^m; terms n<mn < m are bounded by the Lipschitz property). Conclude that WW is differentiable nowhere.

Solution

Solution of Exercise 10.9.

(a) (3/4)nφ(4n)=12(3/4)n\norm{(3/4)^n\varphi(4^n\cdot)}_\infty = \frac12 (3/4)^n: normal convergence, so WW is continuous (Theorem 10.11).

(b) Fix xx, mm; choose the sign of hm=±124mh_m = \pm\frac12 4^{-m} so that the segment [4mx,4m(x+hm)]\intcc{4^mx}{4^m(x + h_m)} (of length 12\frac12) contains no half-integer, making φ\varphi affine of slope ±1\pm1 on it (possible: an interval of length 12\frac12 meets at most one half-integer point; pick the side avoiding it).

For n>mn > m: 4nhm=±124nm4^n h_m = \pm\frac12 4^{n-m} is an integer, and φ\varphi is 11-periodic: the nn-th term of the difference vanishes.

For n=mn = m: φ(4mx+4mhm)φ(4mx)=4mhm=12\abs{\varphi(4^m x + 4^m h_m) - \varphi(4^m x)} = \abs{4^m h_m} = \frac12 (φ\varphi affine of slope ±1\pm 1 on the segment), so the term contributes exactly (3/4)m1/2hm=(3/4)m4m=3m(3/4)^m \cdot \frac{1/2}{\abs{h_m}} = (3/4)^m\,4^m = 3^m in the quotient.

For n<mn < m: the 11-Lipschitz φ\varphi gives (3/4)n(φ(4nx+4nhm)φ(4nx))(3/4)n4nhm=3nhm\bigl|(3/4)^n\bigl(\varphi(4^nx + 4^nh_m) - \varphi(4^nx)\bigr)\bigr| \leq (3/4)^n 4^n\abs{h_m} = 3^n\abs{h_m}: each contributes at most 3n3^n to the quotient.

Hence

W(x+hm)W(x)hm3mn=0m13n=3m3m12=3m+12.\Bigl|\frac{W(x + h_m) - W(x)}{h_m}\Bigr| \geq 3^m - \sum_{n=0}^{m-1} 3^n = 3^m - \frac{3^m - 1}{2} = \frac{3^m + 1}{2} \longrightarrow \infty .

If WW were differentiable at xx, every difference quotient along hm0h_m \to 0 would converge to W(x)W'(x): contradiction. WW is continuous everywhere, differentiable nowhere.

Exercise 10.10

Let un(x)=(1)nxn(1x)u_n(x) = (-1)^n x^n(1 - x) on [0,1]\intcc{0}{1}. Show that un\sum u_n converges pointwise on [0,1]\intcc{0}{1} and compute its sum; show that the convergence is uniform on [0,1]\intcc{0}{1} (bound the remainder RN(x)=n>Nun(x)R_N(x) = \sum_{n > N} u_n(x), a geometric tail, by its first term and maximize xN+1(1x)x^{N+1}(1-x)) but not normal (compute un\norm{u_n}_\infty): uniform convergence is strictly weaker than normal convergence. Contrast with xn(1x)\sum x^n(1-x), whose sum is discontinuous at 11: there even uniformity fails.

Solution

Solution of Exercise 10.10.

Pointwise: for x[0,1)x \in \intco{0}{1} the series is geometric of ratio x-x,

n0(1)nxn(1x)=1x1+x,\sum_{n\geq0}(-1)^n x^n(1-x) = \frac{1-x}{1+x},

and at x=1x = 1 every term vanishes: sum 0=1120 = \frac{1-1}{2}, consistent — the sum is continuous on [0,1]\intcc{0}{1}. Uniformity: the remainder is a geometric tail,

RN(x)=xN+1(1x)1+xxN+1(1x)max[0,1]tN+1(1t)=1N+2(N+1N+2) ⁣N+11N+20,\abs{R_N(x)} = \frac{x^{N+1}(1-x)}{1+x} \leq x^{N+1}(1 - x) \leq \max_{\intcc01} t^{N+1}(1-t) = \frac{1}{N+2}\Bigl(\frac{N+1}{N+2}\Bigr)^{\!N+1} \leq \frac{1}{N+2} \to 0 ,

uniformly in xx. Not normal: un=maxxn(1x)=1n+1(nn+1)n1en\norm{u_n}_\infty = \max x^n(1-x) = \frac{1}{n+1}\bigl(\frac{n}{n+1}\bigr)^n \sim \frac{1}{\eu\,n}, and 1en\sum \frac1{\eu n} diverges. Contrast: xn(1x)\sum x^n(1-x) has partial sums 1xN+11 - x^{N+1}, converging pointwise to the discontinuous 1[0,1)\mathbf 1_{\intco01}: by Theorem 10.4, that convergence cannot be uniform on [0,1]\intcc{0}{1}.

Exercise 10.11 ★★

Let fnff_n \to f uniformly on a metric space XX, each fnf_n continuous, and let xnxx_n \to x in XX. Prove fn(xn)f(x)f_n(x_n) \to f(x). Show by an example on X=[0,1]X = \intcc{0}{1} that pointwise convergence does not suffice, even with ff continuous (use the bumps gng_n of Exercise 10.1 and xn=12nx_n = \frac{1}{\sqrt{2n}}).

Solution

Solution of Exercise 10.11.

The limit ff is continuous (Theorem 10.4). Then

fn(xn)f(x)fn(xn)f(xn)+f(xn)f(x)fnf+f(xn)f(x),\abs{f_n(x_n) - f(x)} \leq \abs{f_n(x_n) - f(x_n)} + \abs{f(x_n) - f(x)} \leq \norm{f_n - f}_\infty + \abs{f(x_n) - f(x)} ,

and both terms tend to 00 (uniform convergence; continuity of ff at xx). Counterexample under mere pointwise convergence: gn(x)=nxenx20g_n(x) = nx\,\eu^{-nx^2} \to 0 pointwise on [0,1]\intcc{0}{1} with gng_n and the limit continuous, yet at xn=12n0x_n = \frac{1}{\sqrt{2n}} \to 0:

gn(xn)=n2e1/2+0=f(0).g_n(x_n) = \sqrt{\frac n2}\,\eu^{-1/2} \longrightarrow +\infty \neq 0 = f(0) .

Exercise 10.12 ★★★

(A Volterra integral equation by series) For fC([0,1])f \in C(\intcc{0}{1}) define Tf(x)=0xf(t) ⁣dtTf(x) = \int_0^x f(t)\,\dd t.

  1. Show by induction that for n1n \geq 1:

    Tnf(x)=0x(xt)n1(n1)!f(t) ⁣dt,Tnffn!.T^n f(x) = \int_0^x \frac{(x - t)^{n-1}}{(n-1)!}\,f(t)\,\dd t, \qquad \norm{T^n f}_\infty \leq \frac{\norm f_\infty}{n!} .
  2. Deduce that S=n0TnfS = \sum_{n\geq0} T^n f converges normally on [0,1]\intcc{0}{1} and solves the integral equation S=f+TSS = f + TS.
  3. Verify that S(x)=f(x)+0xextf(t) ⁣dtS(x) = f(x) + \int_0^x \eu^{x-t}f(t)\,\dd t solves the same equation, and prove uniqueness of continuous solutions (if S=TSS = TS then STnS0\norm{S}_\infty \leq \norm{T^nS}_\infty \to 0): conclude the closed form of the sum.
Solution

Solution of Exercise 10.12.

  1. Induction. n=1n = 1 is the definition. Assume the formula for nn and set g(x)=0x(xt)nn!f(t) ⁣dtg(x) = \int_0^x \frac{(x-t)^n}{n!} f(t)\dd t. For an integrand continuous in (x,t)(x,t) and C1C^1 in xx, the variable-limit parameter integral differentiates as

    g(x)=(xx)nn!f(x)+0x(xt)n1(n1)!f(t) ⁣dt=Tnf(x)g'(x) = \frac{(x-x)^n}{n!}f(x) + \int_0^x \frac{(x-t)^{n-1}}{(n-1)!}f(t)\dd t = T^nf(x)

    (split g(x+h)g(x)g(x+h) - g(x) into the strip xx+h\int_x^{x+h}, which is O(hsup)O(h\cdot\sup) with the integrand vanishing at t=xt = x like hnh^n, and the fixed integral of the xx-increment, handled by the mean value inequality and continuity). Also (Tn+1f)=Tnf(T^{n+1}f)' = T^nf (fundamental theorem of calculus) and g(0)=Tn+1f(0)=0g(0) = T^{n+1}f(0) = 0: two primitives of TnfT^nf vanishing at 00 coincide, so Tn+1f=gT^{n+1}f = g. The bound:

    Tnf(x)f0x(xt)n1(n1)! ⁣dt=fxnn!fn!.\abs{T^nf(x)} \leq \norm f_\infty \int_0^x \frac{(x-t)^{n-1}}{(n-1)!}\dd t = \norm f_\infty\,\frac{x^n}{n!} \leq \frac{\norm f_\infty}{n!} .
  2. nTnfef\sum_n \norm{T^nf}_\infty \leq \eu\,\norm f_\infty: normal, hence uniform, convergence; SS is continuous. The partial sums satisfy SN=f+TSN1S_N = f + T S_{N-1}, and TT is 11-Lipschitz for \norm\cdot_\infty (Tg(x)xg\abs{Tg(x)} \leq x\norm g_\infty): letting NN \to \infty in both sides gives S=f+TSS = f + TS.
  3. Set V(x)=f(x)+ex0xetf(t) ⁣dtV(x) = f(x) + \eu^x\int_0^x \eu^{-t}f(t)\dd t. Then VfV - f is C1C^1 with (Vf)(x)=ex0xetf+f(x)=V(x)(V-f)'(x) = \eu^x\int_0^x\eu^{-t}f + f(x) = V(x), and (TV)=V(TV)' = V with (Vf)(0)=TV(0)=0(V - f)(0) = TV(0) = 0: so Vf=TVV - f = TV, i.e. VV solves the equation. Uniqueness: if S1,S2S_1, S_2 are continuous solutions, D=S1S2D = S_1 - S_2 satisfies D=TDD = TD, hence D=TnDD = T^nD for all nn and DDn!0\norm D_\infty \leq \frac{\norm D_\infty}{n!} \to 0: D=0D = 0. Therefore

    n0Tnf(x)=f(x)+0xextf(t) ⁣dt.\sum_{n\geq0} T^nf(x) = f(x) + \int_0^x \eu^{x-t}f(t)\,\dd t .

    (The series Tn\sum T^n is a geometric series of operators: a first taste of the resolvent (IdT)1(\mathrm{Id} - T)^{-1}, developed in the Year 3 volume.)

10.6 Problem: Rates of approximation and Korovkin’s theorem

Problem 10.1

The Bernstein proof of Theorem 10.16 hides two treasures. First, it is quantitative: how fast BnffB_nf \to f is governed by the modulus of continuity of ff, with the sharp rate attained by x12\abs{x - \frac12}. Second, it is structural: all that mattered was that BnB_n is a positive linear operator behaving well on 11, xx, x2x^2 — that observation, isolated, is Korovkin’s theorem. This problem proves both, and closes with Voronovskaya’s exact asymptotic. Throughout, fC([0,1])f \in C(\intcc{0}{1}), M=fM = \norm f_\infty, pk(x)=(nk)xk(1x)nkp_k(x) = \binom nk x^k(1-x)^{n-k}, and eje_j denotes xxjx \mapsto x^j.

Part I — The Bernstein operator.

  1. Show that BnB_n is linear, positive (f0Bnf0f \geq 0 \Rightarrow B_nf \geq 0), hence monotone (fgBnfBngf \leq g \Rightarrow B_nf \leq B_ng), with Bnff\norm{B_nf}_\infty \leq \norm f_\infty, and that BnfB_nf interpolates ff at both endpoints.
  2. Rederive the identities Bne0=e0B_n e_0 = e_0, Bne1=e1B_n e_1 = e_1 and Bne2=e2+e1e2nB_n e_2 = e_2 + \frac{e_1 - e_2}{n} (differentiate (x+y)n(x + y)^n twice and set y=1xy = 1 - x).
  3. Deduce the variance identity k(knx)2pk(x)=x(1x)n\sum_k \bigl(\frac kn - x\bigr)^2 p_k(x) = \frac{x(1-x)}{n} and, by Cauchy–Schwarz, the first-moment bound

    k=0nknxpk(x)x(1x)n12n.\sum_{k=0}^{n}\Bigl|\frac kn - x\Bigr|\,p_k(x) \leq \sqrt{\frac{x(1-x)}{n}} \leq \frac{1}{2\sqrt n} .
  4. Show that if ff is convex then BnffB_nf \geq f on [0,1]\intcc{0}{1} (finite Jensen inequality for the weights pk(x)p_k(x)).
  5. (Chebyshev’s counting bound, restated) For δ>0\delta > 0 show

    k/nx>δpk(x)x(1x)nδ214nδ2,\sum_{\abs{k/n - x} > \delta} p_k(x) \leq \frac{x(1-x)}{n\delta^2} \leq \frac{1}{4n\delta^2} ,

    and give the probabilistic reading: Bnf(x)B_nf(x) averages ff over a binomial sample mean concentrating at xx.

Part II — Rates: the modulus of continuity. For δ>0\delta > 0 set ω(δ)=sup{f(s)f(t):s,t[0,1], stδ}\omega(\delta) = \sup\{\abs{f(s) - f(t)} : s, t \in \intcc{0}{1},\ \abs{s - t} \leq \delta\}.

  1. Show: ω\omega is finite, nondecreasing, ω(δ)0\omega(\delta) \to 0 as δ0+\delta \to 0^+ (Heine), subadditive (ω(δ1+δ2)ω(δ1)+ω(δ2)\omega(\delta_1 + \delta_2) \leq \omega(\delta_1) + \omega(\delta_2)), and ω(λδ)(1+λ)ω(δ)\omega(\lambda\delta) \leq (1 + \lambda)\,\omega(\delta) for all λ>0\lambda > 0.
  2. Prove the master estimate, for every δ>0\delta > 0:

    Bnf(x)f(x)kω(knx)pk(x)(1+1δkknxpk(x))ω(δ).\abs{B_nf(x) - f(x)} \leq \sum_k \omega\Bigl(\Bigl|\frac kn - x\Bigr|\Bigr)p_k(x) \leq \Bigl(1 + \frac1\delta\sum_k\Bigl|\frac kn - x\Bigr|p_k(x)\Bigr)\,\omega(\delta) .
  3. Choose δ=n1/2\delta = n^{-1/2} and conclude the quantitative Weierstrass theorem:

    Bnff32ω(1n)n0.\norm{B_nf - f}_\infty \leq \frac32\,\omega\Bigl(\frac{1}{\sqrt n}\Bigr) \xrightarrow[n\to\infty]{} 0 .
  4. Deduce the rates: Bnff3L2n\norm{B_nf - f}_\infty \leq \frac{3L}{2\sqrt n} for LL-Lipschitz ff, and 32Cnα/2\leq \frac32 C n^{-\alpha/2} for α\alpha-Hölder ff (f(s)f(t)Cstα\abs{f(s) - f(t)} \leq C\abs{s-t}^\alpha).
  5. (The sharp example — a binomial identity) For m1m \geq 1 prove

    k=m+12m(km)(2mk)=m2(2mm),hencek=02mkm(2mk)=m(2mm)\sum_{k=m+1}^{2m} (k - m)\binom{2m}{k} = \frac{m}{2}\binom{2m}{m}, \qquad\text{hence}\qquad \sum_{k=0}^{2m}\abs{k - m}\binom{2m}{k} = m\binom{2m}{m}

    (use k(2mk)=2m(2m1k1)k\binom{2m}k = 2m\binom{2m-1}{k-1} and the symmetry of the binomial line, which gives j=m2m1(2m1j)=22m2\sum_{j=m}^{2m-1}\binom{2m-1}{j} = 2^{2m-2}).

  6. For f(t)=t12f(t) = \abs{t - \frac12} deduce the exact value and its asymptotics (central binomial, Example 6.14):

    B2mf(12)f(12)=(2mm)24m    12πm:B_{2m}f\Bigl(\frac12\Bigr) - f\Bigl(\frac12\Bigr) = \frac{\binom{2m}{m}}{2\cdot4^{m}} \;\sim\; \frac{1}{2\sqrt{\pi m}} :

    the rate ω(n1/2)\omega(n^{-1/2}) of question 8 is attained (up to a constant) — for merely continuous ff, Bernstein’s n1/2n^{-1/2} is honest.

Part III — Korovkin’s theorem. Let (Ln)(L_n) be a sequence of positive linear operators from C([0,1])C(\intcc{0}{1}) to itself such that LnejejL_ne_j \to e_j uniformly for j=0,1,2j = 0, 1, 2.

  1. Show that a positive linear LL is monotone and satisfies LfLf\abs{Lf} \leq L\abs f pointwise.
  2. Show: for every ε>0\varepsilon > 0 there is δ>0\delta > 0 such that for all s,x[0,1]s, x \in \intcc{0}{1}:

    f(s)f(x)ε+2Mδ2(sx)2\abs{f(s) - f(x)} \leq \varepsilon + \frac{2M}{\delta^2}(s - x)^2

    (treat sxδ\abs{s - x} \leq \delta by Heine and sx>δ\abs{s-x} > \delta by the crude bound 2M2M).

  3. Fix xx, apply LnL_n to the inequality of question 13 in the variable ss, and derive

    Lnf(x)f(x)Lne0(x)εLne0(x)+2Mδ2(Lne2(x)2xLne1(x)+x2Lne0(x)).\abs{L_nf(x) - f(x)\,L_ne_0(x)} \leq \varepsilon\,L_ne_0(x) + \frac{2M}{\delta^2} \bigl(L_ne_2(x) - 2x\,L_ne_1(x) + x^2 L_ne_0(x)\bigr).
  4. Show that supx(Lne2(x)2xLne1(x)+x2Lne0(x))0\sup_x \bigl(L_ne_2(x) - 2x\,L_ne_1(x) + x^2 L_ne_0(x)\bigr) \to 0, then assemble Korovkin’s theorem: LnffL_nf \to f uniformly for every fC([0,1])f \in C(\intcc{0}{1}).
  5. Check that (Bn)(B_n) satisfies Korovkin’s hypotheses: Weierstrass a third time, from three monomials.
  6. Let InI_n be the piecewise-affine interpolation operator at the nodes kn\frac kn. Show InI_n is positive linear, Ine0=e0I_ne_0 = e_0, Ine1=e1I_ne_1 = e_1, and Ine2e2=14n2\norm{I_ne_2 - e_2}_\infty = \frac{1}{4n^2} (on each cell the error of affine interpolation of t2t^2 is (ta)(bt)(t - a)(b - t)). Conclude by Korovkin: polygonal interpolants converge uniformly for every continuous ff.

Part IV — Dividends: density, moments, derivatives.

  1. Show that polynomials with rational coefficients are dense in (C([0,1]),)\bigl(C(\intcc{0}{1}), \norm\cdot_\infty\bigr): this Banach space is separable.
  2. (Moments determine the function) Let fC([0,1])f \in C(\intcc{0}{1}) with 01f(t)tn ⁣dt=0\int_0^1 f(t)\,t^n \dd t = 0 for every nNn \in \N. Show 01fP=0\int_0^1 f P = 0 for every polynomial, then 01f2=0\int_0^1 f^2 = 0, then f=0f = 0.
  3. Prove the derivative identity

    (Bnf)(x)=nk=0n1(f(k+1n)f(kn))(n1k)xk(1x)n1k(B_nf)'(x) = n\sum_{k=0}^{n-1}\Bigl( f\Bigl(\frac{k+1}{n}\Bigr) - f\Bigl(\frac kn\Bigr)\Bigr)\, \binom{n-1}{k}x^k(1-x)^{n-1-k}

    (differentiate pkp_k and reindex — an Abel summation).

  4. Suppose ff is C1C^1. Using the mean value theorem in each increment and comparing with Bn1(f)B_{n-1}(f'), show (Bnf)f(B_nf)' \to f' uniformly on [0,1]\intcc{0}{1}. Deduce: for fC1f \in C^1 there are polynomials converging to ff together with their derivatives.
  5. Suppose ff is C2C^2. By Taylor–Lagrange at xx show

    Bnf(x)f(x)f2x(1x)nf8n:\abs{B_nf(x) - f(x)} \leq \frac{\norm{f''}_\infty}{2}\cdot\frac{x(1-x)}{n} \leq \frac{\norm{f''}_\infty}{8n} :

    smoothness upgrades the rate from n1/2n^{-1/2} to n1n^{-1}.

Part V — Saturation: Voronovskaya’s theorem.

  1. Prove the fourth-moment identity

    k(knx)4pk(x)=nx(1x)(1+3(n2)x(1x))n2(n1)\sum_k (k - nx)^4 p_k(x) = nx(1-x)\bigl(1 + 3(n-2)x(1-x)\bigr) \leq n^2 \quad (n \geq 1)

    (expand k4k^4 in falling factorials k(k1)k(k-1)\cdots and use the differentiation trick of question 2 twice more).

  2. (Voronovskaya) Let ff be C2C^2 and x[0,1]x \in \intcc{0}{1}. Writing f(t)=f(x)+f(x)(tx)+f(x)2(tx)2+η(t)(tx)2f(t) = f(x) + f'(x)(t-x) + \frac{f''(x)}2(t-x)^2 + \eta(t)(t-x)^2 with η\eta bounded and η(t)0\eta(t) \to 0 as txt \to x, prove

    n(Bnf(x)f(x))nx(1x)2f(x)n\bigl(B_nf(x) - f(x)\bigr) \xrightarrow[n\to\infty]{} \frac{x(1-x)}{2}\,f''(x)

    (split the η\eta sum at txδ\abs{t - x} \leq \delta; control the far part with question 23). So the error of question 22 is exact in order and constant: BnB_n saturates at 1n\frac1n, however smooth ff may be — compare Exercise 10.7.

  3. Synthesis. In one sentence each: (i) what positivity alone bought (Parts I and III); (ii) where compactness of [0,1]\intcc{0}{1} entered every part; (iii) why three test functions suffice in Korovkin’s theorem; (iv) the trade-off Bernstein makes (robust n1/2n^{-1/2} for rough ff, but a 1n\frac1n ceiling for smooth ff), and which chapter of this book will play the same game with trigonometric polynomials.
Solution

Solution of Problem 10.1.

1. Linearity is clear from the formula. Positivity: the weights pk(x)0p_k(x) \geq 0, so f0f \geq 0 forces Bnf0B_nf \geq 0; monotonicity follows applied to gfg - f. Bound: ±ff\pm f \leq \norm f_\infty gives ±BnffBne0=f\pm B_nf \leq \norm f_\infty B_ne_0 = \norm f_\infty. Endpoints: pk(0)=1k=0p_k(0) = \mathbf 1_{k=0} and pk(1)=1k=np_k(1) = \mathbf 1_{k=n}, so Bnf(0)=f(0)B_nf(0) = f(0), Bnf(1)=f(1)B_nf(1) = f(1).

2. Differentiate (x+y)n=k(nk)xkynk(x+y)^n = \sum_k\binom nk x^ky^{n-k} in xx, multiply by xx, and set y=1xy = 1 - x:

nx=kkpk(x);nx = \sum_k k\,p_k(x) ;

twice, multiplying by x2x^2: n(n1)x2=kk(k1)pk(x)n(n-1)x^2 = \sum_k k(k-1)p_k(x). Hence Bne0=1B_ne_0 = 1 (binomial theorem), Bne1(x)=nxn=xB_ne_1(x) = \frac{nx}{n} = x, and

Bne2(x)=kk2pkn2=n(n1)x2+nxn2=x2+x(1x)n.B_ne_2(x) = \frac{\sum_k k^2p_k}{n^2} = \frac{n(n-1)x^2 + nx}{n^2} = x^2 + \frac{x(1-x)}{n} .

3. Expand:

k(knx) ⁣2pk=Bne2(x)2xBne1(x)+x2=x(1x)n.\sum_k\Bigl(\frac kn - x\Bigr)^{\!2} p_k = B_ne_2(x) - 2x\,B_ne_1(x) + x^2 = \frac{x(1-x)}{n} .

Cauchy–Schwarz with the splitting k/nxpkpk\abs{k/n - x}\sqrt{p_k} \cdot \sqrt{p_k}:

kknxpk(k(knx)2pk) ⁣1/2=x(1x)n12n,\sum_k\Bigl|\frac kn - x\Bigr| p_k \leq \Bigl(\sum_k\Bigl(\frac kn - x\Bigr)^2 p_k\Bigr)^{\!1/2} = \sqrt{\frac{x(1-x)}{n}} \leq \frac{1}{2\sqrt n},

using x(1x)14x(1-x) \leq \frac14.

4. The weights pk(x)p_k(x) are nonnegative with sum 11 and barycenter kknpk(x)=x\sum_k \frac kn p_k(x) = x (question 2). The finite Jensen inequality for the convex ff (induction from the two-point definition, Year 1 volume) gives

f(x)=f(kknpk)kf(kn)pk=Bnf(x).f(x) = f\Bigl(\sum_k \frac kn\,p_k\Bigr) \leq \sum_k f\Bigl(\frac kn\Bigr)p_k = B_nf(x) .

5. On {k:k/nx>δ}\{k : \abs{k/n - x} > \delta\} one has (k/nxδ)2>1\bigl(\frac{k/n - x}{\delta}\bigr)^2 > 1, so

k/nx>δpk1δ2k(knx)2pk=x(1x)nδ214nδ2.\sum_{\abs{k/n-x}>\delta} p_k \leq \frac{1}{\delta^2}\sum_k\Bigl(\frac kn - x\Bigr)^2p_k = \frac{x(1-x)}{n\delta^2} \leq \frac{1}{4n\delta^2} .

Reading: pk(x)p_k(x) is the law of a sample frequency Sn/nS_n/n of nn coin flips of bias xx; its mean is xx, its variance x(1x)n0\frac{x(1-x)}n \to 0, and the display is Chebyshev’s inequality: the mass concentrates at xx, so averaging ff against it reproduces f(x)f(x) in the limit (Chapter 23 makes the vocabulary official).

6. ω2M<\omega \leq 2M < \infty; monotonicity is clear (sup over a larger set). Heine: ff continuous on a compact is uniformly continuous, which says exactly ω(δ)0\omega(\delta) \to 0 as δ0+\delta \to 0^+. Subadditivity: if stδ1+δ2\abs{s - t} \leq \delta_1 + \delta_2, the point uu on the segment [s,t]\intcc st at distance min(δ1,st)\min(\delta_1, \abs{s-t}) from ss satisfies suδ1\abs{s-u} \leq \delta_1, utδ2\abs{u-t} \leq \delta_2, and f(s)f(t)f(s)f(u)+f(u)f(t)\abs{f(s)-f(t)} \leq \abs{f(s)-f(u)} + \abs{f(u)-f(t)}. Iterating, ω(pδ)pω(δ)\omega(p\delta) \leq p\,\omega(\delta) for pNp \in \N^*; for λ>0\lambda > 0, with p=λ1+λp = \lceil\lambda\rceil \leq 1 + \lambda: ω(λδ)ω(pδ)pω(δ)(1+λ)ω(δ)\omega(\lambda\delta) \leq \omega(p\delta) \leq p\,\omega(\delta) \leq (1+\lambda)\omega(\delta).

7. Since pk=1\sum p_k = 1:

Bnf(x)f(x)=k(f(k/n)f(x))pkkω(k/nx)pk.\abs{B_nf(x) - f(x)} = \Bigl|\sum_k\bigl(f(k/n) - f(x)\bigr)p_k\Bigr| \leq \sum_k\omega\bigl(\abs{k/n - x}\bigr)p_k .

For each kk, question 6 with λ=k/nx/δ\lambda = \abs{k/n - x}/\delta gives ω(k/nx)(1+k/nxδ)ω(δ)\omega(\abs{k/n-x}) \leq \bigl(1 + \frac{\abs{k/n-x}}\delta\bigr)\omega(\delta); summing against the pkp_k yields the master estimate.

8. Insert the bound of question 3:

Bnf(x)f(x)(1+12δn)ω(δ),\abs{B_nf(x) - f(x)} \leq \Bigl(1 + \frac{1}{2\delta\sqrt n}\Bigr)\omega(\delta),

uniformly in xx; with δ=n1/2\delta = n^{-1/2} the parenthesis is 32\frac32: Bnff32ω(n1/2)0\norm{B_nf - f}_\infty \leq \frac32\omega(n^{-1/2}) \to 0 by question 6 (Heine). This reproves Theorem 10.16 with a rate.

9. LL-Lipschitz means ω(δ)Lδ\omega(\delta) \leq L\delta: rate 3L2n\frac{3L}{2\sqrt n}. α\alpha-Hölder means ω(δ)Cδα\omega(\delta) \leq C\delta^\alpha: rate 3C2nα/2\frac{3C}{2} n^{-\alpha/2}.

10. Using k(2mk)=2m(2m1k1)k\binom{2m}k = 2m\binom{2m-1}{k-1}:

k=m+12mk(2mk)=2mj=m2m1(2m1j)=2m22m2,\sum_{k=m+1}^{2m}k\binom{2m}k = 2m\sum_{j=m}^{2m-1}\binom{2m-1}{j} = 2m\cdot 2^{2m-2},

because j2m1jj \mapsto 2m-1-j bijects {m,,2m1}\{m,\dots,2m-1\} onto {0,,m1}\{0,\dots,m-1\}, so the sum is half of 22m12^{2m-1}. Also k=m+12m(2mk)=22m(2mm)2\sum_{k=m+1}^{2m}\binom{2m}k = \frac{2^{2m} - \binom{2m}m}{2} (same symmetry). Hence

k=m+12m(km)(2mk)=m22m1m22m(2mm)2=m2(2mm).\sum_{k=m+1}^{2m}(k-m)\binom{2m}k = m\,2^{2m-1} - m\,\frac{2^{2m} - \binom{2m}m}{2} = \frac m2\binom{2m}m .

The substitution k2mkk \mapsto 2m-k maps the terms with k<mk < m onto those with k>mk > m (equal binomials, equal km\abs{k-m}): the absolute sum is twice the one-sided sum, m(2mm)m\binom{2m}m.

11. At x=12x = \frac12, pk(12)=(2mk)22mp_k(\tfrac12) = \binom{2m}k2^{-2m} and f(12)=0f(\tfrac12) = 0:

B2mf(12)=kk2m12(2mk)22m=22m2mm(2mm)=(2mm)24m12πmB_{2m}f\Bigl(\frac12\Bigr) = \sum_k\Bigl|\frac{k}{2m} - \frac12\Bigr| \binom{2m}k 2^{-2m} = \frac{2^{-2m}}{2m}\,m\binom{2m}m = \frac{\binom{2m}m}{2\cdot4^m} \sim \frac{1}{2\sqrt{\pi m}}

by Example 6.14. Since ωf(δ)=δ\omega_f (\delta) = \delta here (the function is 11-Lipschitz and the bound is attained), question 8 predicts at most 32(2m)1/2\frac32(2m)^{-1/2}: the true error 12πm\frac{1}{2\sqrt{\pi m}} has exactly the order n1/2n^{-1/2} — the rate is sharp up to the constant.

12. fgf \leq g gives gf0g - f \geq 0, so L(gf)0L(g-f) \geq 0, i.e. LfLgLf \leq Lg. From fff-\abs f \leq f \leq \abs f: LfLfLf-L\abs f \leq Lf \leq L\abs f, i.e. LfLf\abs{Lf} \leq L\abs f.

13. By Heine choose δ\delta with f(s)f(x)ε\abs{f(s)-f(x)} \leq \varepsilon whenever sxδ\abs{s-x} \leq \delta. If sx>δ\abs{s - x} > \delta, then (sx)2δ2>1\frac{(s-x)^2}{\delta^2} > 1 and f(s)f(x)2M2Mδ2(sx)2\abs{f(s)-f(x)} \leq 2M \leq \frac{2M}{\delta^2}(s-x)^2. In both cases the claimed bound holds.

14. Fix xx; question 13 says, as functions of ss:

εe02Mδ2qx    ff(x)e0    εe0+2Mδ2qx,qx=e22xe1+x2e0.-\varepsilon e_0 - \frac{2M}{\delta^2}q_x \;\leq\; f - f(x)e_0 \;\leq\; \varepsilon e_0 + \frac{2M}{\delta^2}q_x, \qquad q_x = e_2 - 2x\,e_1 + x^2e_0 .

Apply the monotone linear LnL_n (question 12) and evaluate at xx:

Lnf(x)f(x)Lne0(x)εLne0(x)+2Mδ2(Lne2(x)2xLne1(x)+x2Lne0(x)).\abs{L_nf(x) - f(x)L_ne_0(x)} \leq \varepsilon L_ne_0(x) + \frac{2M}{\delta^2}\bigl(L_ne_2(x) - 2xL_ne_1(x) + x^2L_ne_0(x)\bigr) .

15. Write αj=Lnejej\alpha_j = L_ne_j - e_j, so αj0\norm{\alpha_j}_\infty \to 0. Since e2(x)2xe1(x)+x2e0(x)=0e_2(x) - 2xe_1(x) + x^2e_0(x) = 0:

Lne2(x)2xLne1(x)+x2Lne0(x)=α2(x)2xα1(x)+x2α0(x),L_ne_2(x) - 2xL_ne_1(x) + x^2L_ne_0(x) = \alpha_2(x) - 2x\,\alpha_1(x) + x^2\alpha_0(x),

of sup norm at most α2+2α1+α00\norm{\alpha_2} + 2\norm{\alpha_1} + \norm{\alpha_0} \to 0. Also Lne0e0L_ne_0 \to e_0 uniformly, so Lne02L_ne_0 \leq 2 for large nn, and f(x)Lne0(x)1Mα00\abs{f(x)}\abs{L_ne_0(x) - 1} \leq M\norm{\alpha_0} \to 0. Assembling with question 14: for nn large, uniformly in xx,

Lnf(x)f(x)2ε+2Mδ2o(1)+Mo(1)3ε:\abs{L_nf(x) - f(x)} \leq 2\varepsilon + \frac{2M}{\delta^2}\,o(1) + M\,o(1) \leq 3\varepsilon :

LnffL_nf \to f uniformly — Korovkin’s theorem.

16. Bne0=e0B_ne_0 = e_0 and Bne1=e1B_ne_1 = e_1 exactly, and Bne2e2=maxxx(1x)n=14n0\norm{B_ne_2 - e_2}_\infty = \max_x\frac{x(1-x)}{n} = \frac{1}{4n} \to 0 (question 2): Korovkin applies, and Weierstrass follows for the third time.

17. InfI_nf is linear in ff (nodal values are), and on each cell the affine interpolant of nonnegative nodal values is nonnegative: positive. Ine0=e0I_ne_0 = e_0 and Ine1=e1I_ne_1 = e_1 because an affine function equals its own interpolant. On a cell [a,b]\intcc ab (ba=1nb - a = \frac1n), the affine interpolant of e2e_2 is L(t)=(a+b)tabL(t) = (a+b)t - ab, and

L(t)t2=(ta)(bt)[0,(ba)24],L(t) - t^2 = (t-a)(b-t) \in \intcc{0}{\tfrac{(b-a)^2}{4}} ,

with the maximum at the midpoint: Ine2e2=14n20\norm{I_ne_2 - e_2}_\infty = \frac{1}{4n^2} \to 0. Korovkin: InffI_nf \to f uniformly for every continuous ff — polygonal approximation, with no further estimate needed.

18. Given ff and ε\varepsilon: Weierstrass provides a polynomial P=j=0dajxjP = \sum_{j=0}^d a_jx^j with fPε2\norm{f - P}_\infty \leq \frac\varepsilon2; replacing each aja_j by a rational bjb_j with ajbjε2(d+1)\abs{a_j - b_j} \leq \frac{\varepsilon}{2(d+1)} moves the sup norm on [0,1]\intcc{0}{1} by at most ε2\frac\varepsilon2. The set of rational-coefficient polynomials is a countable union (over dd) of countable sets, hence countable, and dense: C([0,1])C(\intcc{0}{1}) is separable.

19. By linearity 01fP=0\int_0^1 fP = 0 for every polynomial PP. Choose polynomials PnfP_n \to f uniformly (Weierstrass):

01f2=01f(fPn)ffPn0,\Bigl|\int_0^1 f^2\Bigr| = \Bigl|\int_0^1 f\,(f - P_n)\Bigr| \leq \norm f_\infty\,\norm{f - P_n}_\infty \longrightarrow 0 ,

so 01f2=0\int_0^1 f^2 = 0. If f(x0)0f(x_0) \neq 0, continuity gives f2c>0f^2 \geq c > 0 on a subinterval, contradicting the vanishing integral: f=0f = 0. Consequently two continuous functions with the same moments ftn\int f t^n coincide.

20. With pn,k(x)=(nk)xk(1x)nkp_{n,k}(x) = \binom nk x^k(1-x)^{n-k} and the conventions pn1,1=pn1,n=0p_{n-1,-1} = p_{n-1,n} = 0, the product rule and k(nk)=n(n1k1)k\binom nk = n\binom{n-1}{k-1}, (nk)(nk)=n(n1k)(n-k)\binom nk = n\binom{n-1}{k} give

pn,k(x)=n(pn1,k1(x)pn1,k(x)).p_{n,k}'(x) = n\bigl(p_{n-1,k-1}(x) - p_{n-1,k}(x)\bigr) .

Summing against f(k/n)f(k/n) and shifting the index in the first sum (Abel summation):

(Bnf)(x)=nj=0n1(f(j+1n)f(jn))pn1,j(x).(B_nf)'(x) = n\sum_{j=0}^{n-1}\Bigl(f\Bigl(\frac{j+1}n\Bigr) - f\Bigl(\frac jn\Bigr)\Bigr)p_{n-1,j}(x) .

21. By the mean value theorem, f(j+1n)f(jn)=1nf(ξj)f(\frac{j+1}n) - f(\frac jn) = \frac1n f'(\xi_j) with ξj(j/n,(j+1)/n)\xi_j \in \intoo{j/n}{(j+1)/n}, so (Bnf)(x)=jf(ξj)pn1,j(x)(B_nf)'(x) = \sum_j f'(\xi_j)\,p_{n-1,j}(x). The node jn1\frac{j}{n-1} also lies in [j/n,(j+1)/n]\intcc{j/n}{(j+1)/n} (both inequalities reduce to jn1j \leq n-1), hence ξjjn11n\abs{\xi_j - \frac j{n-1}} \leq \frac1n and

(Bnf)(x)Bn1(f)(x)jf(ξj)f(jn1)pn1,j(x)ωf(1n)0\bigl|(B_nf)'(x) - B_{n-1}(f')(x)\bigr| \leq \sum_j\Bigl|f'(\xi_j) - f'\Bigl(\frac{j}{n-1}\Bigr)\Bigr| p_{n-1,j}(x) \leq \omega_{f'}\Bigl(\frac1n\Bigr) \longrightarrow 0

uniformly. Since Bn1(f)fB_{n-1}(f') \to f' uniformly (Theorem 10.16 applied to the continuous ff'), the triangle inequality gives (Bnf)f(B_nf)' \to f' uniformly. The polynomials Pn=BnfP_n = B_nf then converge to ff in the C1C^1 sense.

22. Taylor–Lagrange at xx: f(kn)f(x)=f(x)(knx)+f(ξk)2(knx)2f(\frac kn) - f(x) = f'(x)(\frac kn - x) + \frac{f''(\xi_k)}2(\frac kn - x)^2. Summing against pkp_k, the linear term dies (question 2):

Bnf(x)f(x)f2k(knx)2pk=f2x(1x)nf8n.\abs{B_nf(x) - f(x)} \leq \frac{\norm{f''}_\infty}{2}\sum_k\Bigl(\frac kn - x\Bigr)^2p_k = \frac{\norm{f''}_\infty}{2}\cdot\frac{x(1-x)}{n} \leq \frac{\norm{f''}_\infty}{8n} .

23. Two more differentiations of (x+y)n(x+y)^n give the factorial moments, with n(j)=n(n1)(nj+1)n_{(j)} = n(n-1)\cdots(n-j+1):

kk(j)pk=n(j)xj(j=3,4),\sum_k k_{(j)}\,p_k = n_{(j)}\,x^j \qquad (j = 3, 4),

and k3=k(3)+3k(2)+kk^3 = k_{(3)} + 3k_{(2)} + k, k4=k(4)+6k(3)+7k(2)+kk^4 = k_{(4)} + 6k_{(3)} + 7k_{(2)} + k convert them to power moments:

kk3pk=n(3)x3+3n(2)x2+nx,kk4pk=n(4)x4+6n(3)x3+7n(2)x2+nx.\sum_k k^3p_k = n_{(3)}x^3 + 3n_{(2)}x^2 + nx, \qquad \sum_k k^4p_k = n_{(4)}x^4 + 6n_{(3)}x^3 + 7n_{(2)}x^2 + nx .

Expanding (knx)4(k - nx)^4 and collecting (a patient but purely mechanical computation with the four power moments):

k(knx)4pk=nx(1x)(1+3(n2)x(1x)).\sum_k(k-nx)^4p_k = nx(1-x)\bigl(1 + 3(n-2)x(1-x)\bigr) .

With x(1x)14x(1-x) \leq \frac14: the right side is at most n4(1+3n4)=3n216+n4n2\frac n4\bigl(1 + \frac{3n}4\bigr) = \frac{3n^2}{16} + \frac n4 \leq n^2 for n1n \geq 1.

24. The Peano form of Taylor at xx defines η(t)=f(t)f(x)f(x)(tx)12f(x)(tx)2(tx)2\eta(t) = \frac{f(t) - f(x) - f'(x)(t-x) - \frac12f''(x)(t-x)^2} {(t-x)^2} for txt \neq x, η(x)=0\eta(x) = 0: by Taylor–Lagrange η(t)=12(f(ξ)f(x))\eta(t) = \frac12\bigl(f''(\xi) - f''(x)\bigr) for some ξ\xi between tt and xx, so ηf\abs\eta \leq \norm{f''}_\infty and η(t)0\eta(t) \to 0 as txt \to x (continuity of ff''). Summing the expansion against pkp_k and using questions 2–3:

n(Bnf(x)f(x))=x(1x)2f(x)+nkη(kn)(knx)2pk.n\bigl(B_nf(x) - f(x)\bigr) = \frac{x(1-x)}{2}f''(x) + n\sum_k\eta\Bigl(\frac kn\Bigr)\Bigl(\frac kn - x\Bigr)^2p_k .

Given ε\varepsilon, choose δ\delta with ηε\abs\eta \leq \varepsilon on txδ\abs{t - x}\leq\delta. Near part: at most εnx(1x)nε\varepsilon\,n\cdot\frac{x(1-x)}n \leq \varepsilon. Far part: with C=fC = \norm{f''}_\infty and question 23,

nCk/nx>δ(knx)2pknCδ2k(knx)4pk=nCδ2n4k(knx)4pkCδ2n0.n\,C\sum_{\abs{k/n-x}>\delta}\Bigl(\frac kn - x\Bigr)^2p_k \leq \frac{nC}{\delta^2}\sum_k\Bigl(\frac kn - x\Bigr)^4p_k = \frac{nC}{\delta^2 n^4}\sum_k(k-nx)^4p_k \leq \frac{C}{\delta^2 n} \longrightarrow 0 .

Hence n(Bnf(x)f(x))x(1x)2f(x)n(B_nf(x) - f(x)) \to \frac{x(1-x)}2f''(x) — Voronovskaya’s theorem. For f=e2f = e_2 this is exact at every nn (Exercise 10.7): the 1n\frac1n ceiling is real.

25. (i) Positivity turned pointwise inequalities into operator inequalities: it gave the norm bound, Jensen, Chebyshev, and the whole of Korovkin — linearity alone proves nothing here. (ii) Compactness entered through Heine (questions 6, 13), through boundedness of ff, and through the very norm \norm\cdot_\infty being finite. (iii) Three test functions suffice because positivity reduces everything to controlling LnL_n on the single family (sx)2=e22xe1+x2e0(s-x)^2 = e_2 - 2xe_1 + x^2e_0, whose span is that of e0,e1,e2e_0, e_1, e_2. (iv) Bernstein converges at the honest rate ω(n1/2)\omega(n^{-1/2}) for every continuous ff (sharp, question 11) but saturates at 1n\frac1n for smooth ff (question 24); the Fourier chapter runs the same program for periodic functions with Fejér’s kernel — another positive operator with the same virtues and the same modesty.