Mathematics · Book 3 · Bachelor Year 1

Universitaire wiskunde — Bachelor jaar 1

Universitaire wiskunde — Bachelor jaar 1 · Bachelor Year 1

15Integratie op een segment

De integraal uit het middelbareschoolvolume was gebaseerd op oppervlakten die intuïtief werden opgevat. Dit hoofdstuk construeert haar: eerst voor trapfuncties, waar de integraal een eindige som is, daarna voor continue (en stuksgewijs continue) functies door uniforme benadering — de plek waar de stelling van Heine (Stelling 13.22) haar nut bewijst. De hoofdstelling van de integraalrekening verbindt vervolgens de constructie met primitieven, en Riemannsommen verbinden haar met discrete gemiddelden.

Overal zijn a<ba < b reëel.

15.1 Trapfuncties

Definitie 15.1

φ ⁣:[a,b]R\varphi \colon \intcc{a}{b} \to \R is een trapfunctie wanneer er een onderverdeling a=x0<x1<<xn=ba = x_0 < x_1 < \dots < x_n = b bestaat zodat φ\varphi constant is, gelijk aan cic_i, op elk open interval (xi1,xi)\intoo{x_{i-1}}{x_i} (de waarden in de knooppunten zijn onbepaald). Haar integraal is

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

onafhankelijk van de gekozen onderverdeling (verfijn twee onderverdelingen tot hun gemeenschappelijke: elk lid blijft ongewijzigd onder verfijning).

Propositie 15.2

Op trapfuncties is de integraal lineair, stijgend (φψ    φψ\varphi \leq \psi \implies \int\varphi \leq \int\psi), en voldoet aan de Chasles-relatie ab=ac+cb\int_a^b = \int_a^c + \int_c^b voor a<c<ba < c < b.

Bewijs. De motor is de verfijningsinvariantie, gesteld in de definitie: het invoegen van één extra knooppunt t(xi1,xi)t \in \intoo{x_{i-1}}{x_i} in een onderverdeling vervangt de term ci(xixi1)c_i(x_i - x_{i-1}) door ci(txi1)+ci(xit)c_i(t - x_{i-1}) + c_i(x_i - t) — hetzelfde getal — zodat de integraal ongewijzigd blijft onder elke eindige verfijning. Neem nu φ\varphi met onderverdeling σ\sigma en ψ\psi met onderverdeling σ\sigma': op de gemeenschappelijke verfijning σσ\sigma \cup \sigma' zijn beide trapfuncties met dezelfde knooppunten, en op elk stuk is φ+λψ\varphi + \lambda\psi constant gelijk aan ci+λdic_i + \lambda d_i: lineariteit reduceert tot lineariteit van eindige sommen. Stijgend: cidic_i \leq d_i op elk stuk geeft ciΔidiΔi\sum c_i \Delta_i \leq \sum d_i \Delta_i (de lengten Δi0\Delta_i \geq 0). Chasles: voeg het knooppunt cc in en splits de som daar.

15.2 Integraal van een continue functie

Stelling 15.3 (Uniforme benadering)

Zij ff continu op [a,b]\intcc{a}{b}. Voor elke ε>0\varepsilon > 0 bestaan er trapfuncties φ,ψ\varphi, \psi met

φfψenψφε op [a,b].\varphi \leq f \leq \psi \qquad\text{en}\qquad \psi - \varphi \leq \varepsilon \text{ op } \intcc{a}{b}.

Bewijs. Volgens de stelling van Heine (Stelling 13.22) is ff uniform continu: kies δ\delta voor ε\varepsilon, en een onderverdeling met maaswijdte <δ< \delta (bijvoorbeeld gelijkmatig verdeeld, met n>baδn > \frac{b - a}{\delta} stukken). Op elk gesloten stuk [xi1,xi]\intcc{x_{i-1}}{x_i} bereikt ff een minimum mim_i en een maximum MiM_i (Stelling 13.13), en MimiεM_i - m_i \leq \varepsilon (de twee extreme punten liggen binnen δ\delta van elkaar). Definieer φ=mi\varphi = m_i en ψ=Mi\psi = M_i op (xi1,xi)\intoo{x_{i-1}}{x_i} (en φ=ψ=f\varphi = \psi = f in de knooppunten).

Stelling 15.4 (Definitie van de integraal)

Zij ff continu op [a,b]\intcc{a}{b}. De twee getallen

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

zijn gelijk; hun gemeenschappelijke waarde is de integraal abf\int_a^b f (ook geschreven als abf(t) ⁣dt\int_a^b f(t)\,\dd t). Zij valt samen met het vorige begrip op trapfuncties, en strekt zich uit tot stuksgewijs continue functies door [a,b]\intcc{a}{b} te splitsen bij de discontinuïteiten (Chasles als definitie daar).

Bewijs. Beide verzamelingen zijn niet-leeg (ff is begrensd) en elke onderste trapintegraal is \leq elke bovenste (monotonie op trapfuncties): dus I(f)I+(f)I_-(f) \leq I_+(f). Volgens Stelling 15.3 bestaat er voor elke ε\varepsilon een paar met ψφε(ba)\int\psi - \int\varphi \leq \varepsilon(b - a): het supremum en infimum worden samengeknepen, I=I+I_- = I_+.

Voorbeeld 15.5 (Stuksgewijs continu, zonder drama)

De vloerfunctie op [0,3]\intcc{0}{3} is een vermomde trapfunctie: door te splitsen bij haar sprongen,

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 ,

en de waarden in de sprongpunten 1,21, 2 zijn irrelevant: een functie in eindig veel punten wijzigen verandert geen enkele integraal (de omkaderende trapfuncties blijven onaangetast). Dit is de volledige inhoud van de uitbreiding “stuksgewijs continu”: snijd bij de eindig veel discontinuïteiten, integreer elk continu stuk, tel op — Chasles als definitie.

Voorbeeld 15.6 (De definitie rekent, eenmalig)

Zij f(x)=xf(x) = x op [0,1]\intcc{0}{1} en verdeel in nn gelijke stukken. De beste trapfuncties die constant zijn op de stukken zijn φ=k1n\varphi = \frac{k-1}{n} en ψ=kn\psi = \frac kn op het kk-de stuk, met

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} .

Elke onderste integraal is I(f)I+(f)\leq I_-(f) \leq I_+(f) \leq elke bovenste, dus n12nI(f)I+(f)n+12n\frac{n-1}{2n} \leq I_-(f) \leq I_+(f) \leq \frac{n+1}{2n} voor alle nn: beide knijpen naar 12\frac12, en 01x ⁣dx=12\int_0^1 x\,\dd x = \frac12 rechtstreeks uit de definitie. Het afsluitende inzicht: dit is de eerste en laatste keer dat we vanuit de definitie integreren — de hoofdstelling hieronder vervangt al zulke berekeningen door één opzoeking van een primitieve, wat het hele economische punt van dit hoofdstuk is.

Stelling 15.7 (Eigenschappen)

Voor f,gf, g continu (of stuksgewijs continu) op [a,b]\intcc{a}{b} en λR\lambda \in \R:

  1. lineariteit: (f+λg)=f+λg\int (f + \lambda g) = \int f + \lambda \int g;
  2. monotonie: fg    abfabgf \leq g \implies \int_a^b f \leq \int_a^b g; en 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 (met de conventie ba=ab\int_b^a = -\int_a^b, geldig voor elke volgorde van de grenzen);
  4. strikte positiviteit: als ff continu is, f0f \geq 0 en abf=0\int_a^b f = 0, dan is f=0f = 0 overal op [a,b]\intcc{a}{b}.

Bewijs. (1)–(3) gaan van trapfuncties naar de limiet via de supremum/infimum-definitie. Lineariteit verdient de details één keer: gegeven ε>0\varepsilon > 0, kader φffψf\varphi_f \leq f \leq \psi_f en φggψg\varphi_g \leq g \leq \psi_g met marges ε\leq \varepsilon (Stelling 15.3). Voor λ0\lambda \geq 0 is φf+λφgf+λgψf+λψg\varphi_f + \lambda\varphi_g \leq f + \lambda g \leq \psi_f + \lambda\psi_g een omkadering door trapfuncties met marge (1+λ)ε\leq (1 + \lambda)\varepsilon, en haar trapintegralen zijn gelijk aan φf+λφg\int \varphi_f + \lambda\int\varphi_g enz. (Propositie 15.2): ε0\varepsilon \to 0 laten gaan knijpt (f+λg)\int(f + \lambda g) naar f+λg\int f + \lambda\int g. Voor λ<0\lambda < 0 keert vermenigvuldiging met λ\lambda de omkadering van gg om — de onderste trapfunctie van λg\lambda g is λψg\lambda \psi_g — en dezelfde knijp verloopt met verwisselde rollen. De grens ff\abs{\int f} \leq \int \abs f komt uit fff-\abs f \leq f \leq \abs f en monotonie.

(4) Contrapositie: als f(x0)=m>0f(x_0) = m > 0, dan levert continuïteit een deelinterval van lengte η>0\eta > 0 op waar fm2f \geq \frac m2; de trapfunctie die daar m2\frac m2 waard is en 00 elders is f\leq f, dus fmη2>0\int f \geq \frac{m\eta}{2} > 0.

Voorbeeld 15.8 (Chasles aan het werk: integralen met absolute waarden)

Om een absolute waarde te integreren, snijd waar het teken verandert.

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 ,

en, door [0,2π]\intcc{0}{2\pi} bij π\pi te snijden:

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 ,

terwijl 02πsint ⁣dt=0\int_0^{2\pi} \sin t\,\dd t = 0: opheffing is reëel, en dit is waarom de uitspraak over strikte positiviteit (Stelling 15.7 (4)) de hypothese f0f \geq 0 met zich meedraagt — zonder haar bewijst een verdwijnende integraal niets over ff. Het afsluitende inzicht: f\int \abs f meet oppervlakte, f\int f meet getekende balans; de ongelijkheid ff\abs{\int f} \leq \int\abs f is het exacte verslag van wat opheffing kan vernietigen.

15.3 De hoofdstelling van de integraalrekening

Stelling 15.9 (Hoofdstelling van de integraalrekening)

Zij ff continu op een interval II en aIa \in I. De functie

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

is van klasse C1C^1 op II, met F=fF' = f: elke continue functie op een interval heeft primitieven. Bijgevolg geldt voor elke primitieve GG van ff:

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

Bewijs. Fixeer x0Ix_0 \in I en ε>0\varepsilon > 0; continuïteit in x0x_0 levert δ\delta op met f(t)f(x0)ε\abs{f(t) - f(x_0)} \leq \varepsilon voor tx0δ\abs{t - x_0} \leq \delta. Voor 0<hδ0 < \abs{h} \leq \delta (en x0+hIx_0 + h \in I) geeft Chasles

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 ,

waarvan de absolute waarde 1hhε=ε\leq \frac{1}{\abs h}\cdot \abs h\, \varepsilon = \varepsilon is (grens (2), geldig voor elke volgorde van de grenzen). Dus F(x0)=f(x0)F'(x_0) = f(x_0); F=fF' = f is continu: FF is C1C^1. Als ook G=fG' = f, dan is (GF)=0(G - F)' = 0 op het interval, dus G=F+cG = F + c (Gevolg 14.12), en G(b)G(a)=F(b)F(a)=abfG(b) - G(a) = F(b) - F(a) = \int_a^b f.

Voorbeeld 15.10 (Symmetrie vóór berekening)

Op een symmetrisch interval doet pariteit het werk: als ff oneven is, stuurt de substitutie ttt \mapsto -t de integraal a0f\int_{-a}^{0} f naar 0af-\int_0^a f, dus

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

Zo is 11t3cost1+t4 ⁣dt=0\int_{-1}^{1} \frac{t^3\cos t}{1 + t^4}\,\dd t = 0 zonder een primitieve in zicht (de integrand is oneven), en ππt2cost ⁣dt=20πt2cost ⁣dt\int_{-\pi}^{\pi} t^2\cos t\,\dd t = 2\int_0^\pi t^2\cos t\, \dd t. Controleer symmetrie voordat je naar technieken grijpt: de snelste integraal is die welke nooit berekend wordt.

Voorbeeld 15.11 (Een afgeleide op het zicht herkennen)

Bereken 0π/2 ⁣dx1+cosx\displaystyle\int_0^{\pi/2} \frac{\dd x}{1 + \cos x}. De halvehoekidentiteit 1+cosx=2cos2x21 + \cos x = 2\cos^2\frac x2 maakt van de integrand 12(1+tan2x2)\frac{1}{2}\bigl(1 + \tan^2\frac x2\bigr), wat precies de afgeleide van tanx2\tan\frac x2 is:

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 .

Geen substitutiemachinerie was nodig — enkel de reflex om een integrand als iemands afgeleide te lezen, waarbij de hoofdstelling de rest doet. (Het systematische instrument achter zulke goniometrische integralen, de substitutie t=tanx2t = \tan\frac x2, hoort tot de standaardgereedschapskist opgebouwd vanuit Stelling 15.15 (2).)

Voorbeeld 15.12 (Functies gedefinieerd door integralen)

De hoofdstelling vervaardigt functies. Zij

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

Geen enkele combinatie van klassieke functies heeft afgeleide et2\eu^{-t^2} (een stelling van Liouville, toegegeven); toch bestaat FF, is C1C^1 met F(x)=ex2>0F'(x) = \eu^{-x^2} > 0, strikt stijgend, oneven (substitueer ttt \mapsto -t), en begrensd: voor 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} ,

dus FF(1)+e11+e1F \leq F(1) + \eu^{-1} \leq 1 + \eu^{-1}. (De exacte limiet, π2\frac{\sqrt\pi}{2}, wordt berekend met dubbelintegralen in het volume van Jaar 3.) Kettingregel voor bewegende grenzen:  ⁣d ⁣dxxx2et2 ⁣dt=2xex4ex2\frac{\dd}{\dd x}\int_x^{x^2} \eu^{-t^2}\dd t = 2x\,\eu^{-x^4} - \eu^{-x^2}. Het afsluitende inzicht: integratie schept nieuwe functies uit oude, met al hun eigenschappen afleesbaar uit de integrand — de primitieve die je niet kunt opschrijven is nog steeds een functie die je volledig beheerst.

Voorbeeld 15.13 (Schatten zonder evalueren)

De integralen Rn=01tn1+t ⁣dtR_n = \int_0^1 \frac{t^n}{1 + t}\,\dd t hebben geen prettige gesloten vorm, toch pinnen monotonie ze precies vast: op [0,1]\intcc{0}{1} is 1211+t1\frac12 \leq \frac{1}{1+t} \leq 1, dus

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} :

de exacte orde van verval (RnR_n \sim een veelvoud van 1n\frac1n, in feite Rn12nR_n \sim \frac{1}{2n}) met twee regels en zonder primitieve. De weekendproblemen van dit hoofdstuk en het volgende draaien op precies zulke omkaderingen — het eerste instinct van de analyst vóór een integraal moet zijn: begrens haar, en pas dan, indien nodig, bereken haar.

Voorbeeld 15.14 (Gemiddelde waarden)

Het gemiddelde van een continue ff over [a,b]\intcc{a}{b} is 1baabf\frac{1}{b-a}\int_a^b f. Voor de boog van de sinus:

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 :

een volledige positieve boog heeft niet als gemiddelde 12\frac12 maar 2π\frac2\pi — de kromme brengt meer tijd hoog door dan een driehoek zou doen. Volgens Oefening 15.11 (middelwaardestelling voor integralen, g=1g = 1) is het gemiddelde een waarde: sinc=2π\sin c = \frac2\pi voor een zekere c(0,π)c \in \intoo{0}{\pi}. En volgens de Riemannsommen van dit hoofdstuk is het gemiddelde de limiet van gewone gemiddelden van nn steekproeven — de brug tussen het discrete gemiddelde van gegevens en het continue gemiddelde van een signaal, wat de manier is waarop de integraal de fysica binnenkomt.

Stelling 15.15 (Partiële integratie; substitutie)

  1. Als u,vu, v van klasse C1C^1 zijn op [a,b]\intcc{a}{b}:

    abuv=[uv]ababuv.\int_a^b u'v = \bigl[uv\bigr]_a^b - \int_a^b uv' .
  2. Als φ\varphi van klasse C1C^1 is op [α,β]\intcc{\alpha}{\beta} en ff continu op φ([α,β])\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 .

Bewijs. (1) (uv)=uv+uv(uv)' = u'v + uv'; integreer over [a,b]\intcc{a}{b} en pas de hoofdstelling toe op de C1C^1-functie uvuv.

(2) Zij FF een primitieve van ff op het beeldinterval (Stelling 15.9). Dan is (Fφ)=(fφ)φ(F \circ \varphi)' = (f \circ \varphi)\,\varphi' (kettingregel), dus beide leden zijn gelijk aan F(φ(β))F(φ(α))F(\varphi(\beta)) - F(\varphi(\alpha)).

Voorbeeld 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. En met de substitutie 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 ,

— een kwart van de eenheidsschijf, zoals de meetkunde eist. (Linearisatie uit Methode 3.11 aan het werk.)

Opmerking 15.17 (Veelvoorkomende valkuilen in de integraalrekening)

(i) Substituties moeten C1C^1 zijn op het hele interval: de verandering x=1tx = \frac1t is ongeoorloofd over 00 heen; blindelings toegepast op 11 ⁣dx1+x2\int_{-1}^{1}\frac{\dd x}{1 + x^2} “bewijst” zij dat de integraal gelijk is aan haar eigen tegengestelde. Wanneer een substitutie een singulariteit heeft, snijd het interval eerst (Chasles), substitueer op elk stuk, en herenig pas daarna. (ii) Logaritmische primitieven hebben absolute waarden nodig:  ⁣dxx2=lnx2+C\int \frac{\dd x}{x - 2} = \ln\abs{x - 2} + C op elke zijde van 22 afzonderlijk — ln(x2)\ln(x - 2) schrijven op (0,1)\intoo{0}{1} is de logaritme van een negatief getal schrijven; en de constante CC mag verschillen op de twee zijden van de singulariteit. (iii) Een verdwijnende integraal doodt de functie niet: 02πsin=0\int_0^{2\pi}\sin = 0; positiviteit van de integrand is vereist voordat men f=0f = 0 besluit (Voorbeeld 15.8). (iv) Riemannsommen moeten gekalibreerd zijn: in banf(a+kban)\frac{b-a}{n}\sum f\bigl(a + k\frac{b-a}{n}\bigr) moeten de stap buiten en de punten binnen bij dezelfde onderverdeling passen — de frequente fout is een som k=1nf(kn)\sum_{k=1}^{n} f\bigl(\frac kn\bigr) zonder de factor 1n\frac1n, die divergeert in plaats van naar 01f\int_0^1 f te convergeren. Checklist vóór het inroepen van Stelling 15.20: haal 1n\frac1n buiten, herschrijf de somterm als ff van kn\frac kn, benoem ff en controleer haar continuïteit.

Voorbeeld 15.18 (Raad, differentieer, pas aan)

Wat is 1x(lnt)2 ⁣dt\int_1^x (\ln t)^2\,\dd t? Raad een primitieve van de vorm tP(lnt)t\,P(\ln t) met PP een veelterm en differentieer:

(tP(lnt))=P(lnt)+P(lnt).\bigl(t\,P(\ln t)\bigr)' = P(\ln t) + P'(\ln t) .

We hebben P(u)+P(u)=u2P(u) + P'(u) = u^2 nodig: neem P(u)=u22u+2P(u) = u^2 - 2u + 2 (door de coëfficiënten neerwaarts vanaf u2u^2 te matchen). Bijgevolg

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 ,

een resultaat dat anders bereikt wordt met twee partiële integraties. Het afsluitende inzicht: voor integranden van de vorm (veelterm in lnt\ln t) of (veelterm maal eλt\eu^{\lambda t}) heeft de primitieve dezelfde gedaante — het differentiëren van een gevormde gok zet integratie om in lineaire algebra op coëfficiënten, sneller en minder foutgevoelig dan herhaalde partiële integratie.

Voorbeeld 15.19 (De boemerangintegraal)

Bereken I=0π/2excosx ⁣dxI = \int_0^{\pi/2} \eu^x \cos x\,\dd x. Integreer tweemaal partieel, waarbij je telkens de goniometrische factor differentieert:

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 .

De integraal is naar zichzelf teruggekeerd: I=eπ/21II = \eu^{\pi/2} - 1 - I, waaruit

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

Het afsluitende inzicht: wanneer de integrand een product is van twee functies die zichzelf reproduceren onder differentiatie (eax\eu^{ax}, cosbx\cos bx, sinbx\sin bx), leveren twee partiële integraties een lineaire vergelijking voor de onbekende integraal — los haar op in plaats van te integreren; equivalent, ga langs e(a+ib)x\eu^{(a + \iu b)x} (Hoofdstuk 3) en neem de reële delen. Beide wegen geven hetzelfde antwoord, en controleren dat ze dat doen is een gratis gezondheidstest.

15.4 Riemannsommen

Stelling 15.20 (Riemannsommen)

Zij ff continu op [a,b]\intcc{a}{b}. Dan geldt

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 ,

en evenzo met willekeurige evaluatiepunten binnen de deelintervallen.

Bewijs. SnS_n is de integraal van de trapfunctie φn\varphi_n die gelijk is aan f(a+kban)f(a + k\frac{b-a}{n}) op het kk-de deelinterval. Gegeven ε>0\varepsilon > 0 levert uniforme continuïteit (Heine) δ\delta op; voor n>baδn > \frac{b-a}{\delta} ligt elk punt van een deelinterval binnen δ\delta van zijn evaluatiepunt, dus fφnε\abs{f - \varphi_n} \leq \varepsilon op [a,b]\intcc{a}{b}, waaruit

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

Een linkse Riemannsom met n = 8 rechthoeken: naarmate de maaswijdte krimpt, dwingt uniforme continuïteit de trapoppervlakte naar ∈t_ab f.
Een linkse Riemannsom met n=8n = 8 rechthoeken: naarmate de maaswijdte krimpt, dwingt uniforme continuïteit de trapoppervlakte naar abf\int_a^b f.

Voorbeeld 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: een limiet onzichtbaar voor elementaire grenzen, doorzichtig als Riemannsom.

Voorbeeld 15.22 (Een tweede Riemannsom, met kalibratie)

Bepaal limnk=1nn(n+k)2\lim_{n\to\infty} \sum_{k=1}^{n} \dfrac{n}{(n+k)^2}. Kalibreer:

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} ,

een Riemannsom van de continue f(x)=1(1+x)2f(x) = \frac{1}{(1+x)^2} op [0,1]\intcc{0}{1}: de limiet 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 .

Het afsluitende inzicht: de hele kunst zit in de middelste regel — dwing de somterm in de gedaante f(kn)f(\frac kn) ten koste van het uithalen van precies één factor 1n\frac1n; zodra de gedaante klopt, doet de stelling de analyse en doet de hoofdstelling het rekenwerk.

Opmerking 15.23 (Waar de integraal hierna werkt)

De constructies van dit hoofdstuk hebben elk een vervolg. Riemannsommen keren terug in Hoofdstuk 17 als de brug tussen reeksen en integralen (vergelijking van 1nα\sum \frac{1}{n^\alpha} met  ⁣dttα\int \frac{\dd t}{t^\alpha}); de integraalrest is de scherpste vorm van de formule van Taylor in Hoofdstuk 16; de op sup\sup gebaseerde definitie is het prototype voor de Lebesgue-integraal van het volume van Jaar 3, waar dezelfde drie eigenschappen (lineariteit, monotonie, een convergentiestelling) opnieuw worden opgebouwd op een veel grotere klasse van functies. En het weekendprobleem hieronder maakt van partiële integratie rekenkunde: de irrationaliteit van π2\pi^2.

15.5 Oefeningen

Oefening 15.1

Bereken: 01 ⁣dxx24\displaystyle\int_0^1 \frac{\dd x}{x^2 - 4} (partiële breuken, Hoofdstuk 9); 1elnt ⁣dt\displaystyle\int_1^{\eu} \ln t \,\dd t; 0πtsint ⁣dt\displaystyle\int_0^{\pi} t \sin t\, \dd t.

Oplossing

Oplossing van Oefening 15.1.

1x24=1/4x21/4x+2\dfrac{1}{x^2-4} = \dfrac{1/4}{x - 2} - \dfrac{1/4}{x+2} (afdekmethode), dus

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}.

Met partiële integratie (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.

Met partiële integratie (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.

Oefening 15.2

Bereken 01t(t2+1)2 ⁣dt\displaystyle\int_0^{1} \frac{t}{(t^2+1)^2}\,\dd t (substitutie) en 0π/2cos3t ⁣dt\displaystyle\int_0^{\pi/2} \cos^3 t\, \dd t (schrijf cos3=cos(1sin2)\cos^3 = \cos(1 - \sin^2)).

Oplossing

Oplossing van Oefening 15.2.

Substitutie 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 .

Met 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.

Oefening 15.3

Bepaal de limieten, als Riemannsommen:

un=k=1nnn2+k2,vn=1n(2n)!n!nnn (neem logaritmen).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{(neem logaritmen)}.
Oplossing

Oplossing van Oefening 15.3.

un=1nk=1n11+(k/n)2u_n = \dfrac1n \sum_{k=1}^{n} \dfrac{1}{1 + (k/n)^2}: een Riemannsom van x11+x2x \mapsto \frac{1}{1+x^2} op [0,1]\intcc{0}{1}, dus 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. Bijgevolg vne2ln21=4ev_n \to \eu^{2\ln 2 - 1} = \dfrac 4\eu. (Controle van de identificatie: (2n)!n!nn=k=1nn+kn\frac{(2n)!}{n!\,n^n} = \prod_{k=1}^{n} \frac{n+k}{n}.)

Oefening 15.4

Zij ff continu op [0,1]\intcc{0}{1}. Bereken limn01xnf(x) ⁣dx\lim_{n\to\infty} \int_0^1 x^n f(x)\,\dd x. (Snijd [0,1]\intcc{0}{1} bij 1δ1 - \delta.)

Oplossing

Oplossing van Oefening 15.4.

De limiet is 00. Zij M=supfM = \sup \abs f en ε(0,1)\varepsilon \in \intoo{0}{1}. Snijd bij 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 .

Aangezien (1ε)n0(1 - \varepsilon)^n \to 0 (Oefening 11.3), is de limsup van het linkerlid Mε\leq M\varepsilon voor elke ε\varepsilon: de integraal streeft naar 00.

Oefening 15.5 ★★

(Cauchy–Schwarz) Voor f,gf, g continu op [a,b]\intcc{a}{b}, bewijs

(abfg) ⁣2abf2abg2,\Bigl(\int_a^b fg\Bigr)^{\!2} \leq \int_a^b f^2 \cdot \int_a^b g^2 ,

door ab(f+λg)20\int_a^b (f + \lambda g)^2 \geq 0 als een kwadratische vergelijking in λ\lambda te ontwikkelen. Wanneer geldt gelijkheid?

Oplossing

Oplossing van Oefening 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 voor alle λ\lambda. Als g2=0\int g^2 = 0, dan is g=0g = 0 (strikte positiviteit, Stelling 15.7 (4)) en de ongelijkheid is 000 \leq 0. Anders is QQ een echte kwadratische vergelijking, overal 0\geq 0: haar discriminant is 0\leq 0, d.w.z. (fg)2f2g2\bigl(\int fg\bigr)^2 \leq \int f^2 \int g^2.

Gelijkheid dan en slechts dan als de discriminant verdwijnt dan en slechts dan als Q(λ0)=0Q(\lambda_0) = 0 voor een zekere λ0\lambda_0, d.w.z. (f+λ0g)2=0\int (f + \lambda_0 g)^2 = 0, d.w.z. (opnieuw strikte positiviteit) f=λ0gf = -\lambda_0 g: gelijkheid geldt precies wanneer ff en gg evenredig zijn.

Oefening 15.6 ★★

Zij ff continu op R\R, TT-periodiek. Bewijs dat aa+Tf\int_a^{a+T} f niet afhangt van aa, en dat 1x0xf(t) ⁣dt1T0Tf\frac1x \int_0^x f(t)\,\dd t \to \frac 1T \int_0^T f als x+x \to +\infty.

Oplossing

Oplossing van Oefening 15.6.

Zij Φ(a)=aa+Tf\Phi(a) = \int_a^{a+T} f. Volgens de hoofdstelling (Stelling 15.9) is Φ\Phi afleidbaar met Φ(a)=f(a+T)f(a)=0\Phi'(a) = f(a + T) - f(a) = 0: constant.

Voor x>0x > 0, schrijf 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 .

Dan is 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), en nTx1\frac{nT}{x} \to 1: de limiet is 1T0Tf\frac1T \int_0^T f.

Oefening 15.7 ★★

Voor ff continu op [0,1]\intcc{0}{1} met 01f=12\int_0^1 f = \frac12, bewijs dat ff een vast punt heeft in [0,1]\intcc{0}{1}. (Integreer f(x)xf(x) - x en gebruik strikte positiviteit, Stelling 15.7 (4), via haar contrapositie gecombineerd met de tussenwaardestelling.)

Oplossing

Oplossing van Oefening 15.7.

Zij g(x)=f(x)xg(x) = f(x) - x: continu, met

01g=01f12=0.\int_0^1 g = \int_0^1 f - \frac12 = 0 .

Als gg nooit verdween, zou de tussenwaardestelling een constant teken afdwingen (een continue functie op een interval die beide tekens aanneemt, verdwijnt); zeg g>0g > 0. Dan, door strikte positiviteit (Stelling 15.7 (4) toegepast op g>0g > 0, wat g>0\int g > 0 geeft): tegenspraak met g=0\int g = 0. Dus g(c)=0g(c) = 0 voor een zekere cc: f(c)=cf(c) = c.

Oefening 15.8 ★★★

(Wallis-integralen) Zij Wn=0π/2sinnt ⁣dtW_n = \int_0^{\pi/2} \sin^n t\,\dd t.

  1. Bewijs de recurrentie nWn=(n1)Wn2n W_n = (n-1) W_{n-2} (n2n \geq 2) met partiële integratie, en bereken W0,W1W_0, W_1, daarna W2pW_{2p} en W2p+1W_{2p+1} in gesloten vorm.
  2. Bewijs dat (Wn)(W_n) dalend is met Wn+1Wn1\frac{W_{n+1}}{W_n} \to 1; bewijs dat de grootheid (n+1)Wn+1Wn(n+1)\,W_{n+1} W_n constant is, gelijk aan π2\frac\pi2; en leid het equivalent Wnπ2nW_n \sim \sqrt{\dfrac{\pi}{2n}} af.
Oplossing

Oplossing van Oefening 15.8.

  1. Met partiële integratie met 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),

    dus nWn=(n1)Wn2nW_n = (n-1)W_{n-2}. Uit 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. Op (0,π2)\intoo{0}{\frac\pi2} is 0<sint<10 < \sin t < 1, dus sinn+1<sinn\sin^{n+1} < \sin^n en (Wn)(W_n) is (strikt) dalend, positief. Insluiten met de recurrentie:

    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 voldoet aan an=an1a_n = a_{n-1} door de recurrentie (n+1)Wn+1=nWn1(n+1)W_{n+1} = nW_{n-1}, dus an=a0=1W1W0=π2a_n = a_0 = 1 \cdot W_1 W_0 = \frac\pi2. Dan

    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}} .

Oefening 15.9 ★★★

(Niven: π\pi is irrationaal) Veronderstel π=ab\pi = \frac ab met a,bNa, b \in \N^*, en stel, voor een nader te kiezen nn,

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. Bewijs dat 0<Inπ(πa)nn!0 < I_n \leq \pi\,\frac{(\pi a)^n}{n!}, wat <1< 1 is voor grote nn.
  2. Bewijs dat PP en al haar afgeleiden gehele waarden aannemen in 00 en in π=ab\pi = \frac ab. (Binomiale ontwikkeling: de coëfficiënten van PP maal k!k! zijn geheel voor knk \geq n; en P(πx)=P(x)P(\pi - x) = P(x).)
  3. Stel Q=PP+P(4)Q = P - P'' + P^{(4)} - \dots (een eindige som). Controleer dat (QsinxQcosx)=Psinx\bigl(Q'\sin x - Q\cos x\bigr)' = P \sin x, en leid af dat In=Q(π)+Q(0)I_n = Q(\pi) + Q(0) een geheel getal is.
  4. Besluit.
Oplossing

Oplossing van Oefening 15.9.

  1. Op (0,π)\intoo{0}{\pi}: x>0x > 0, abx=b(abx)=b(πx)>0a - bx = b(\frac ab - x) = b(\pi - x) > 0 en sinx>0\sin x > 0, dus de integrand is >0> 0 en In>0I_n > 0 (strikte positiviteit). Grens: op [0,π]\intcc{0}{\pi} is xπx \leq \pi en abxaa - bx \leq a, dus Pπnann!P \leq \frac{\pi^n a^n}{n!} en Inπ(πa)nn!I_n \leq \pi\,\frac{(\pi a)^n}{n!}, wat naar 00 streeft (de faculteit verslaat de meetkundige term: het is de algemene term van de convergente exponentiële reeks, cf. Voorbeeld 11.12); in het bijzonder In<1I_n < 1 voor grote nn.
  2. Ontwikkel 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}: dus P=1n!jcjxn+jP = \frac{1}{n!}\sum_j c_j x^{n+j} met gehele cjc_j. Dan is P(k)(0)=0P^{(k)}(0) = 0 voor k<nk < n (valuatie) en, voor nk2nn \leq k \leq 2n, P(k)(0)=k!n!cknP^{(k)}(0) = \frac{k!}{n!} c_{k-n}, een geheel getal omdat n!k!n! \mid k!. Bovendien P(πx)=P(x)P(\pi - x) = P(x) (substitueer: πx\pi - x verwisselt de factoren, met ab(πx)=bxa - b(\pi - x) = bx), dus P(k)(π)=±P(k)(0)P^{(k)}(\pi) = \pm P^{(k)}(0): eveneens gehele getallen.
  3. Met Q=PP+P(4)Q = P - P'' + P^{(4)} - \dots (eindig: PP heeft graad 2n2n): Q+Q=PQ + Q'' = P, en

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

    Bijgevolg In=[QsinxQcosx]0π=Q(π)+Q(0)I_n = \bigl[Q'\sin x - Q\cos x\bigr]_0^{\pi} = Q(\pi) + Q(0), een som van waarden P(2k)P^{(2k)} in 00 en π\pi: een geheel getal volgens (2).

  4. Voor grote nn is InI_n een geheel getal met 0<In<10 < I_n < 1: onmogelijk. De veronderstelling π=ab\pi = \frac ab faalt: π\pi is irrationaal.

Oefening 15.10 ★★★

Zij ff van klasse C1C^1 op [a,b]\intcc{a}{b}. Bewijs de limiet van het Riemann–Lebesgue-type

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

door partiële integratie. Bewijs haar dan opnieuw voor ff enkel continu, door uniforme benadering met trapfuncties (Stelling 15.3).

Oplossing

Oplossing van Oefening 15.10.

C1C^1-geval: met partiële integratie,

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,

in absolute waarde begrensd door 2supf+(ba)supfλ0\frac{2\sup\abs f + (b - a)\sup\abs{f'}}{\lambda} \to 0.

Continu geval: zij ε>0\varepsilon > 0 en kies een trapfunctie φ\varphi met fφε\abs{f - \varphi} \leq \varepsilon (Stelling 15.3 levert φfψ\varphi \leq f \leq \psi met marge ε\leq\varepsilon; neem φ\varphi). Dan

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| ,

en elke 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}: de tweede term streeft naar 00. Bijgevolg is de limsup (ba)ε\leq (b-a)\varepsilon voor elke ε\varepsilon: de limiet is 00.

Oefening 15.11 ★★

(Middelwaardestelling voor integralen) Zij f,gf, g continu op [a,b]\intcc{a}{b} met g0g \geq 0. Bewijs dat er een c[a,b]c \in \intcc{a}{b} bestaat met

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 ,

en toon met een voorbeeld dat de hypothese g0g \geq 0 niet kan worden weggelaten.

Oplossing

Oplossing van Oefening 15.11.

Zij m=minfm = \min f en M=maxfM = \max f, bereikt volgens de extremumstelling. Aangezien g0g \geq 0: mgfgMgm\,g \leq fg \leq M\,g, dus door monotonie

mabg    abfg    Mabg.m \int_a^b g \;\leq\; \int_a^b fg \;\leq\; M \int_a^b g .

Als abg=0\int_a^b g = 0: strikte positiviteit (Stelling 15.7 (4)) dwingt g0g \equiv 0 af, beide leden verdwijnen, en elke cc voldoet. Anders ligt t=fggt = \frac{\int fg}{\int g} in [m,M]=f([a,b])\intcc{m}{M} = f(\intcc{a}{b}) (Stellingen 13.13 en 13.10), dus t=f(c)t = f(c) voor een zekere cc.

Het teken doet ertoe: op [1,1]\intcc{-1}{1} met f(t)=g(t)=tf(t) = g(t) = t: fg=11t2=23\int fg = \int_{-1}^1 t^2 = \frac23, terwijl f(c)11t ⁣dt=0f(c)\int_{-1}^1 t\,\dd t = 0 voor elke cc.

Oefening 15.12 ★★★

(Momenten forceren nulpunten) Zij ff continu op [a,b]\intcc{a}{b} met

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

Bewijs dat ff verdwijnt in n+1n + 1 verschillende punten van (a,b)\intoo{a}{b}. (Als ff enkel van teken wisselt in z1<<zmz_1 < \dots < z_m met mnm \leq n, integreer ff tegen P(t)=(tz1)(tzm)P(t) = (t - z_1) \cdots (t - z_m) en gebruik strikte positiviteit.)

Oplossing

Oplossing van Oefening 15.12.

Als f0f \equiv 0 is de bewering triviaal (elk punt is een nulpunt). Veronderstel dus f≢0f \not\equiv 0 en neem aan dat het hoogstens nn verschillende nulpunten heeft in (a,b)\intoo{a}{b}; zij z1<<zmz_1 < \dots < z_m (mnm \leq n) die nulpunten waar ff van teken wisselt (mogelijk geen). Stel P(t)=i=1m(tzi)P(t) = \prod_{i=1}^{m}(t - z_i) (leeg product =1= 1), van graad mnm \leq n. Op elk deelinterval gesneden door de ziz_i hebben zowel ff als PP een constant teken, en beide keren van teken bij het passeren van een zekere ziz_i: het product fPfP heeft één constant teken op heel (a,b)\intoo{a}{b}. Continu zijnde, niet identiek nul, van constant teken, heeft het abfP>0\bigl|\int_a^b fP\bigr| > 0 (strikte positiviteit toegepast op fP\abs{fP}). Maar fP\int f P is een lineaire combinatie van de momenten ftk\int f\,t^k, knk \leq n, alle nul: tegenspraak. Bijgevolg heeft ff minstens n+1n + 1 verschillende nulpunten in (a,b)\intoo{a}{b}.

Opmerking 15.24 (Perspectieven binnen dit volume)

Drie verdere hoofdstukken steunen rechtstreeks op dit hoofdstuk. Hoofdstuk 16 draagt de integraalrest — de scherpste van de drie Taylor-formules is een nn maal herhaalde partiële integratie. Hoofdstuk 17 zet de omkadering van sommen door integralen om in de beslissende test voor nα\sum n^{-\alpha}, en haar weekendprobleem verfijnt die omkadering tot de constante van Euler. Hoofdstuk 24 maakt de integraal meetkundig: de lengte van een geparametriseerde boog is x(t)2+y(t)2 ⁣dt\int \sqrt{x'(t)^2 + y'(t)^2}\,\dd t, een integraal van een continue functie op een segment — precies het object dat hier is opgebouwd, zonder oneigenlijke theorie nodig. Het meest herbruikte feit zal het nederigste zijn: f(ba)supf\bigl|\int f\bigr| \leq (b - a)\sup\abs{f}, de ongelijkheid die elke puntsgewijze schatting omzet in een integraalschatting.

15.6 Probleem: De integraal-irrationaliteitsmachine

Probleem 15.1

Weekendprobleem — e\eu en π2\pi^2 zijn irrationaal, e\eu tot op zes decimalen, en 227>π\frac{22}{7} > \pi met bewijs

Eén mechanisme drijft dit hele probleem aan: een uitdrukking die een positief geheel getal moet zijn, maar aantoonbaar kleiner is dan 11, kan niet bestaan. Oefening 15.9 (Niven) liet het één keer draaien om te bewijzen dat πQ\pi \notin \Q; hier industrialiseren we het. De machine heeft drie onderdelen nodig: een geheeltalligheids-invoer (eindpuntwaarden van welgekozen veeltermen), een kleinheids-invoer (een factor 1n!\frac{1}{n!} die de integraal verplettert), en een brug (partiële integratie) die ze verbindt. We bewijzen dat e\eu irrationaal is en berekenen haar met gecertificeerde fout, bewijzen Legendres scherpere stelling dat π2\pi^2 irrationaal is, en eindigen met de charmantste integraal in de analyse: 01x4(1x)41+x2 ⁣dx=227π\int_0^1 \frac{x^4(1 - x)^4}{1 + x^2}\dd x = \frac{22}{7} - \pi, die π\pi met de hand omsluit.

Deel I — Brandstof.

  1. Bewijs dat cnn!0\dfrac{c^{\,n}}{n!} \to 0 voor elke vaste c>0c > 0 (voorbij n2cn \geq 2c halveert elke stap de term minstens).
  2. (Beta-integralen) Bewijs, door inductie op mm met partiële integratie:

    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. Leid 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}} af, en — door te vergelijken met de grens x(1x)14x(1 - x) \leq \frac14 — de schatting (2nn)4n2n+1\binom{2n}{n} \geq \dfrac{4^n}{2n+1}, in overeenstemming met (2nn)1/n4\binom{2n}{n}^{1/n} \to 4 uit Probleem 11.1.
  4. Bewijs het kleinheidslemma dat hieronder tweemaal wordt gebruikt: voor elke continue g>0g > 0 op (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}} .

Deel II — e\eu: irrationaliteit, dan zes decimalen. Stel An=01xnex ⁣dxA_n = \displaystyle\int_0^1 x^n \eu^x \dd x.

  1. Bereken A0A_0 en A1A_1, bewijs de recurrentie An=enAn1A_n = \eu - n A_{n-1}, en de grenzen 0<Anen+10 < A_n \leq \dfrac{\eu}{n+1}.
  2. Toon door inductie dat An=αn+βneA_n = \alpha_n + \beta_n \eu met αn,βnZ\alpha_n, \beta_n \in \Z.
  3. Leid af dat e\eu irrationaal is (als e=pq\eu = \frac pq, dan is qAnq A_n een geheel getal gevangen in (0,1)\intoo{0}{1} voor grote nn). Vergelijk met het bewijs van Oefening 11.9: dezelfde clou, andere brandstof.
  4. Bewijs, door inductie en partiële integratie, de exacte restformule

    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. Neem n=9n = 9: begrens R9R_9 met behulp van e<2.75\eu < 2.75 (uit b2=2.75b_2 = 2.75 in Voorbeeld 11.12), evalueer de som, en besluit de gecertificeerde omsluiting 2.7182818e2.71828232.7182818 \leq \eu \leq 2.7182823 — zes decimalen, e2.718282\eu \approx 2.718282, met bewijs.

Deel III — Legendres stelling: π2\pi^2 is irrationaal. Zij f(x)=xn(1x)nn!f(x) = \dfrac{x^n (1 - x)^n}{n!}, en veronderstel π2=ab\pi^2 = \frac ab met a,bNa, b \in \N^*.

  1. Toon f(1x)=f(x)f(1 - x) = f(x) en 0<f14nn!0 < f \leq \dfrac{1}{4^n\,n!} op (0,1)\intoo{0}{1}.
  2. Toon dat f(k)(0)f^{(k)}(0) en f(k)(1)f^{(k)}(1) gehele getallen zijn voor elke k0k \geq 0 (ontwikkel xn(1x)nx^n(1-x)^n met gehele coëfficiënten; k!n!Z\frac{k!}{n!} \in \Z voor knk \geq n; gebruik dan de symmetrie).
  3. Definieer

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

    Toon dat G(0)G(0) en G(1)G(1) gehele getallen zijn (elke bnπ2n2k=ankbkb^n \pi^{2n-2k} = a^{\,n-k}\,b^{\,k}).

  4. Verifieer de telescopering G+π2G=bnπ2n+2f=π2anfG'' + \pi^2 G = b^n \pi^{2n+2} f = \pi^2 a^n f, daarna

     ⁣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. Integreer over [0,1]\intcc{0}{1} en besluit

    π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 ,

    een positief geheel getal begrensd door πan4nn!\dfrac{\pi a^n}{4^n\,n!}.

  6. Besluit met vraag 1 dat π2\pi^2 irrationaal is (Legendre, 1794), en dat dit Oefening 15.9 versterkt: waarom impliceert de irrationaliteit van π2\pi^2 die van π\pi, en niet omgekeerd?

Deel IV — De machine begrijpen.

  1. Lokaliseer de twee tegengestelde krachten (geheeltalligheid van de eindpuntgegevens; analytische kleinheid van de integraal) en de brug, in Deel II en III. Verklaar dan waarom de factor 1n!\frac{1}{n!} in ff de kern is: als hij wordt verwijderd, overleeft de geheeltalligheid, maar welke ongelijkheid sterft, en voor welke beweerde breuken ab\frac ab faalt het bewijs dan?
  2. Effectiviteit: veronderstel dat iemand beweert π2=ab\pi^2 = \frac ab met a10a \leq 10. Toon dat de tegenspraak reeds landt bij n=7n = 7: bereken π(10/4)7/7!0.38<1\pi\,(10/4)^7/7! \approx 0.38 < 1. De machine weerlegt niet enkel; ze weerlegt bij een vast, berekenbaar stadium.
  3. Controleer de brug onvoorwaardelijk op consistentie: bewijs door twee partiële integraties dat

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

    en verzoen met vraag 14 bij n=1n = 1 (houd π2\pi^2 symbolisch: de getelescopeerde identiteit luidt π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. Wat maakt ex\eu^x en sinπx\sin \pi x geschikt als de kernen van de machine? Identificeer de eigenschap (elk voldoet aan een lineaire differentiaalvergelijking met constante coëfficiënten, zodat herhaalde partiële integratie terugkeert naar het begin), en benoem de grens: dezelfde machine, verfijnd door Hermite en Lindemann, bewijst dat e\eu en π\pi transcendent zijn — voorbij dit volume.

Deel V — 227\frac{22}{7} tegenover π\pi, en de moraal.

  1. Stel de veeltermdeling vast

    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},

    en leid de befaamde identiteit af

    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. De integrand is positief: besluit π<227\pi < \frac{22}{7}. Begrens dan 11+x2\frac{1}{1+x^2} tussen 12\frac12 en 11 en gebruik 01(x(1x))4=1630\int_0^1 (x(1-x))^4 = \frac{1}{630} (vraag 3) om te bewijzen dat

    2271630    π    22711260,\frac{22}{7} - \frac{1}{630} \;\leq\; \pi \;\leq\; \frac{22}{7} - \frac{1}{1260} ,

    d.w.z. 3.14126π3.142073.14126 \leq \pi \leq 3.14207: twee correcte decimalen, met de hand.

  3. Veralgemeen: door x4m(1x)4mx^{4m}(1-x)^{4m} te delen door 1+x21 + x^2, toon dat de rest de constante (4)m(-4)^m is (werk modulo x2+1x^2 + 1: (1x)22x(1-x)^2 \equiv -2x), leid rationale getallen rmr_m af met

    πrm415m,\abs{\pi - r_m} \leq 4^{\,1 - 5m} ,

    en controleer dat m=1m = 1 vraag 20–21 reproduceert.

  4. Confronteer deze rationale getallen met de benaderingstheorie van Probleem 14.1: bereken π2271.26103\abs{\pi - \frac{22}{7}} \approx 1.26\cdot10^{-3} tegenover de Dirichlet-garantie 149\frac{1}{49}, en citeer π3551132.7107\abs{\pi - \frac{355}{113}} \approx 2.7\cdot10^{-7} tegenover 111327.8105\frac{1}{113^2} \approx 7.8\cdot10^{-5}: uitzonderlijk goede rationale benaderingen bestaan voor π\pi — consistent, aangezien π\pi niet bekendstaat als slecht benaderbaar.
  5. (De geheeltalligheidsval, geabstraheerd) Bewijs het lemma dat alles verenigt: als xRx \in \R en er gehele getallen an,bna_n, b_n bestaan met 0<an+bnx00 < \abs{a_n + b_n x} \to 0, dan is xx irrationaal. Somm haar gevallen op in dit probleem, in Oefening 15.9, in Oefening 11.9, en in Probleem 14.1.
  6. Synthese, één zin per stuk: (i) de drie onderdelen van de machine en waar elk woont in de gereedschapskist van dit hoofdstuk; (ii) wat de integraal bijdraagt dat de middelwaardestelling van Probleem 14.1 niet kon; (iii) de inventaris van geëxtraheerde resultaten (twee irrationaliteiten, één zes-decimalen-constante, één omsluiting van π\pi, één binomiale grens); (iv) de grens (Hermite, Lindemann; en dezelfde val, gedraaid op ζ(2)\zeta(2) en ζ(3)\zeta(3), in de twintigste-eeuwse rekenkunde).
Oplossing

Oplossing van Probleem 15.1.

1. Zij N=2cN = \lceil 2c \rceil. Voor 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, dus 0<cnn!cNN!2(nN)00 < \frac{c^n}{n!} \leq \frac{c^N}{N!}\,2^{-(n - N)} \to 0: knijp.

2. Fixeer kk; inductie op mm. Voor m=0m = 0: 01xk=1k+1=k!0!(k+1)!\int_0^1 x^k = \frac{1}{k+1} = \frac{k!\,0!}{(k+1)!}. Stap, met partiële integratie (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)!} ,

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

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}}. Aangezien x(1x)14x(1-x) \leq \frac14 op [0,1]\intcc{0}{1}, is de integraal 4n\leq 4^{-n}, waaruit (2nn)4n2n+1\binom{2n}{n} \geq \frac{4^n}{2n+1} — consistent met (2nn)1/n4\binom{2n}{n}^{1/n} \to 4 (Probleem 11.1).

4. De integrand is continu, 0\geq 0, en positief op (0,1)\intoo{0}{1}, dus niet identiek nul: zijn integraal is >0> 0 (Stelling 15.7 (4)). Bovengrens: (x(1x))n4n(x(1-x))^n \leq 4^{-n} en gsupgg \leq \sup\abs g, daarna monotonie.

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. Met partiële integratie: 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}. Grenzen: de integrand is positief, dus An>0A_n > 0; en exe\eu^x \leq \eu geeft Ane01xn=en+1A_n \leq \eu\int_0^1 x^n = \frac{\eu}{n+1}.

6. A0=1+1eA_0 = -1 + 1\cdot\eu. Als An1=αn1+βn1eA_{n-1} = \alpha_{n-1} + \beta_{n-1}\eu met gehele termen, dan

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 ,

beide geheel.

7. Als e=pq\eu = \frac pq: qAn=qαn+pβnZq A_n = q\alpha_n + p\beta_n \in \Z, en 0<qAnqen+1<10 < qA_n \leq \frac{q\eu}{n+1} < 1 voor grote nn: een geheel getal strikt tussen 00 en 11 — onmogelijk. Dus eQ\eu \notin \Q. In Oefening 11.9 was het gevangen geheel getal q!pqq!aqq!\,\frac pq - q!\,a_q; hier is het qAnqA_n: dezelfde val, integraalbrandstof.

8. n=0n = 0: R0=01et=e1R_0 = \int_0^1 \eu^t = \eu - 1, dus e=1+R0\eu = 1 + R_0. Met partiële integratie (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} ,

dus de formule plant zich voort van nn naar n+1n + 1. Grenzen: 1ete1 \leq \eu^t \leq \eu op [0,1]\intcc{0}{1} en 01(1t)n=1n+1\int_0^1 (1-t)^n = \frac{1}{n+1} geven 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, en

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} ,

dus 2.7182818e2.71828232.7182818 \leq \eu \leq 2.7182823: met bewijs, e=2.718282\eu = 2.718282 tot op zes decimalen (werkelijke waarde 2.71828182.7182818\dots).

10. f(1x)=(1x)nxnn!=f(x)f(1 - x) = \frac{(1-x)^n x^n}{n!} = f(x). Op (0,1)\intoo{0}{1}: 0<x(1x)140 < x(1-x) \leq \frac14, dus 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}, dus f=1n!jcjxn+jf = \frac{1}{n!}\sum_j c_j\,x^{n+j} met cjZc_j \in \Z. Bijgevolg f(k)(0)=0f^{(k)}(0) = 0 voor k<nk < n of k>2nk > 2n, en voor nk2nn \leq k \leq 2n: f(k)(0)=k!n!cknf^{(k)}(0) = \frac{k!}{n!}\, c_{k-n}, een geheel getal omdat n!k!n! \mid k!. De symmetrie geeft 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, dus G(0)=k(1)kankbkf(2k)(0)G(0) = \sum_k (-1)^k a^{n-k} b^k f^{(2k)}(0) en evenzo G(1)G(1) zijn gehele getallen volgens vraag 11.

13. In π2G+G\pi^2 G + G'' draagt de term kk van π2G\pi^2 G π2n2k+2f(2k)\pi^{2n-2k+2}f^{(2k)} en de term j=k1j = k - 1 van GG'' draagt (1)k1π2n2k+2f(2k)(-1)^{k-1}\pi^{2n-2k+2}f^{(2k)}: alles heft zich op behalve k=0k = 0 in de eerste som en j=nj = n in de tweede, d.w.z.

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 heeft graad 2n2n, dus f(2n+2)=0f^{(2n+2)} = 0; en bnπ2n=anb^n\pi^{2n} = a^n). Dan

(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. Door te integreren 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) ,

dus πan01fsinπx=G(0)+G(1)Z\pi a^n \int_0^1 f\sin\pi x = G(0) + G(1) \in \Z. Op (0,1)\intoo{0}{1} is f>0f > 0 en sinπx>0\sin\pi x > 0: het linkerlid is positief, dus G(0)+G(1)1G(0) + G(1) \geq 1; en sin1\sin \leq 1 met vraag 10 begrenst het door πan4nn!\frac{\pi a^n}{4^n\,n!}.

15. Volgens vraag 1 (met c=a4c = \frac a4), πan4nn!0\frac{\pi a^n}{4^n n!} \to 0: voor grote nn is het <1< 1, in tegenspraak met G(0)+G(1)1G(0) + G(1) \geq 1. Dus geen enkele breuk ab\frac ab is gelijk aan π2\pi^2: Legendres stelling. Als π\pi rationaal was, zou π2\pi^2 dat ook zijn: dus πQ\pi \notin \Q — en de implicatie loopt enkel deze kant op (2\sqrt2 is irrationaal met rationaal kwadraat), wat de reden is waarom π2Q\pi^2 \notin \Q strikt sterker is dan Oefening 15.9.

16. Geheeltalligheid: vragen 11–12 (eindpuntafgeleiden); kleinheid: vragen 10 en 1; brug: vragen 13–14 (de getelescopeerde dubbele partiële integratie). Zonder 1n!\frac{1}{n!} blijven de eindpuntgegevens gehele getallen (zelfs gemakkelijker), maar de grens wordt πan4n\frac{\pi a^n}{4^n}, die enkel naar 00 streeft wanneer a<4a < 4 — en elke kandidaat heeft a=bπ2>9a = b\,\pi^2 > 9. De faculteit is precies wat de meetkundige groei ana^n voorbijstreeft: geen faculteit, geen stelling.

17. Voor a10a \leq 10 is het gehele getal G(0)+G(1)G(0) + G(1) positief en hoogstens π(10/4)n/n!\pi\,(10/4)^n/n!. Bij n=7n = 7: 2.57=610.352.5^7 = 610.35\dots, dus de grens is π×610.3550400.38<1\frac{\pi \times 610.35}{5040} \approx 0.38 < 1 (bij n=6n = 6 is het nog steeds 1.071.07): de tegenspraak landt bij het zevende stadium, expliciet.

18. Twee partiële integraties:

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}

(de haaktermen verdwijnen: x(1x)x(1-x) in 0,10, 1, en sinπx\sin\pi x in 0,10, 1). Symbolisch luidt de n=1n = 1-telescoop (geen aanname over π\pi) π301f1sinπx=(f1(0)+f1(1))\pi^3\int_0^1 f_1\sin\pi x = -(f_1''(0) + f_1''(1)) met f1=x(1x)f_1 = x(1 - x), f1=2f_1'' = -2: rechterlid 44 — de twee berekeningen komen overeen.

19. ex\eu^x lost y=yy' = y op en sinπx\sin\pi x lost y=π2yy'' = -\pi^2 y op: lineaire vergelijkingen met constante coëfficiënten, zodat partiële integratie de kern naar zichzelf terugvoert en alle randgegevens binnen Z+Ze\Z + \Z\eu houdt (resp. gehele veeltermen in π2\pi^2). Die geslotenheidseigenschap is wat de machine nodig heeft. Verfijnd met kernen aangepast aan meerdere punten tegelijk levert hetzelfde mechanisme Hermites stelling (e\eu transcendent, 1873) en die van Lindemann (π\pi transcendent, 1882) — voorbij dit volume.

20. Veeltermdeling (of vermenigvuldig terug en controleer):

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 .

Door de weergegeven identiteit gedeeld door 1+x21 + x^2 te integreren:

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 ,

met 1723+143+4=3+17\frac17 - \frac23 + 1 - \frac43 + 4 = 3 + \frac17 en arctan1=π4\arctan 1 = \frac\pi4.

21. De integrand is continu, positief op (0,1)\intoo{0}{1}: de integraal is >0> 0, dus π<227\pi < \frac{22}{7}. Bovendien 1211+x21\frac12 \leq \frac{1}{1+x^2} \leq 1 op [0,1]\intcc{0}{1} en 01(x(1x))4=(4!)29!=1630\int_0^1 (x(1-x))^4 = \frac{(4!)^2}{9!} = \frac{1}{630} (vraag 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, dus x4m=(x2)2m1x^{4m} = (x^2)^{2m} \equiv 1 en (1x)2=12x+x22x(1 - x)^2 = 1 - 2x + x^2 \equiv -2x, bijgevolg (1x)4m(2x)2m=4m(x2)m(4)m(1-x)^{4m} \equiv (-2x)^{2m} = 4^m (x^2)^m \equiv (-4)^m: de rest is de constante (4)m(-4)^m, en het quotiënt QmQ_m heeft gehele coëfficiënten (deling door een monische gehele veelterm). Door de identiteit door 1+x21 + x^2 te delen en te integreren:

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 .

Oplossen naar π\pi: met 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}. Voor m=1m = 1: s1=227s_1 = \frac{22}{7}, r1=227r_1 = \frac{22}{7}, grens 44=12564^{-4} = \frac{1}{256} — opnieuw vragen 20–21.

23. π227=227π1.26103\bigl|\pi - \frac{22}{7}\bigr| = \frac{22}{7} - \pi \approx 1.26\cdot10^{-3}, zestien keer beter dan de orde-22-benchmark 1722.0102\frac{1}{7^2} \approx 2.0\cdot10^{-2} gegarandeerd door Dirichlet (Probleem 14.1, vraag 4); en π3551132.7107\bigl|\pi - \frac{355}{113}\bigr| \approx 2.7\cdot10^{-7} verslaat 111327.8105\frac{1}{113^2} \approx 7.8\cdot10^{-5} met een factor 300\approx 300. Geen tegenspraak met iets bewezen: Liouville- ongelijkheden onderbegrenzen benaderingsfouten enkel voor algebraïsche getallen, en geen zulke grens voor π\pi is beschikbaar op dit niveau — π\pi is vrij om spectaculair goed benaderd te worden.

24. Lemma: veronderstel x=pqx = \frac pq en 0<an+bnx00 < \abs{a_n + b_n x} \to 0. Dan is an+bnx=qan+pbnq\abs{a_n + b_n x} = \frac{\abs{q a_n + p b_n}}{q}, met qan+pbnq a_n + p b_n een niet-nul geheel getal (niet-nul omdat de absolute waarde >0> 0 is): dus an+bnx1q\abs{a_n + b_n x} \geq \frac1q voor elke nn, in tegenspraak met de convergentie naar 00. Gevallen: vraag 7 (x=ex = \eu, an=αna_n = \alpha_n, bn=βnb_n = \beta_n); Oefening 11.9 (x=ex = \eu opnieuw, met aq=q!kq1k!a_q = -q!\sum_{k \leq q}\frac{1}{k!}, bq=q!b_q = q!); en Probleem 14.1, vraag 1, is haar meetkundige vorm. In Deel III en in Oefening 15.9 draait de val binnenin de tegenspraak: rationaliteit aannemen zet een uitdrukking om in een geheel getal, dat de analyse vervolgens inknijpt in (0,1)\intoo{0}{1} — hetzelfde principe, getransponeerd.

25. (i) Geheeltalligheid woont in de eindpuntrekenkunde van veeltermen (vragen 6, 11–12), kleinheid in de grenzen die sup\sup-monotonie geeft (vragen 4, 10), de brug in partiële integratie (vragen 8, 13–14) — alle drie zijn stellingen van dit hoofdstuk. (ii) De integraal levert wat de middelwaardestelling niet kon: een exacte identiteit tussen het analytische object en de rekenkundige gegevens (gelijkheid, niet enkel een ongelijkheid met een onbekende cc), wat de reden is waarom de machine π2\pi^2 bereikt terwijl Probleem 14.1 enkel benaderingsexponenten bereikte. (iii) Geëxtraheerd: eQ\eu \notin \Q, π2Q\pi^2 \notin \Q (dus πQ\pi \notin \Q), e=2.718282\eu = 2.718282 gecertificeerd, 2271630π22711260\frac{22}{7} - \frac{1}{630} \leq \pi \leq \frac{22}{7} - \frac{1}{1260}, en (2nn)4n2n+1\binom{2n}{n} \geq \frac{4^n}{2n+1}. (iv) Grens: Hermite en Lindemann duwen dezelfde machine naar transcendentie; en Apéry (1979) liet de geheeltalligheidsval draaien op ζ(3)\zeta(3) — de machine produceert nog steeds twintigste-eeuwse wiskunde.