Mathematics · Boek 3 · Bachelor Year 1

Universitaire wiskunde — Bachelor jaar 1

Universitaire wiskunde — Bachelor jaar 1 · Bachelor Year 1

15Integratie op een segment

De integraal van het bovenbouwvolume steunde op intuïtief opgevatte oppervlakten. Dit hoofdstuk construeert hem: eerst voor trapfuncties, waar de integraal een eindige som is, en daarna voor continue (en stuksgewijs continue) functies door uniforme benadering — de plek waar de stelling van Heine (Stelling 13.22) haar kost verdient. De hoofdstelling van de integraalrekening verbindt de constructie vervolgens met primitieven, en de Riemannsommen verbinden hem met discrete gemiddelden.

Overal zijn a<ba < b reële getallen.

15.1 Trapfuncties

Definitie 15.1

φ ⁣:[a,b]R\varphi \colon \intcc{a}{b} \to \R heet een trapfunctie wanneer er een verdeling a=x0<x1<<xn=ba = x_0 < x_1 < \dots < x_n = b bestaat waarvoor φ\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 vrij). 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 verdeling (verfijn twee verdelingen tot hun gemeenschappelijke verdeling: elk van beide leden blijft bij verfijning ongewijzigd).

Propositie 15.2

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

Bewijs. De motor is de invariantie onder verfijning, vermeld in de definitie: het invoegen van één extra knooppunt t(xi1,xi)t \in \intoo{x_{i-1}}{x_i} in een verdeling 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 onder elke eindige verfijning onveranderd blijft. Neem nu φ\varphi met verdeling σ\sigma en ψ\psi met verdeling σ\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: de lineariteit herleidt zich tot die van eindige sommen. Stijging: cidic_i \leq d_i op elk stuk geeft ciΔidiΔi\sum c_i \Delta_i \leq \sum d_i \Delta_i (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 bij ε\varepsilon, en een verdeling met maaswijdte <δ< \delta (bijvoorbeeld gelijkmatig, 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 is MimiεM_i - m_i \leq \varepsilon (de twee extremale 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φ:φ trapfunctie, φf},I+(f)=inf{abψ:ψ trapfunctie, ψf}I_-(f) = \sup\Bigl\{\int_a^b \varphi : \varphi \text{ trapfunctie},\ \varphi \leq f\Bigr\}, \qquad I_+(f) = \inf\Bigl\{\int_a^b \psi : \psi \text{ trapfunctie},\ \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). Hij valt op trapfuncties samen met het vorige begrip, en breidt uit tot stuksgewijs continue functies door [a,b]\intcc{a}{b} in de discontinuïteitspunten te splitsen (daar is Chasles een definitie).

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 is er voor elke ε\varepsilon een paar met ψφε(ba)\int\psi - \int\varphi \leq \varepsilon(b - a): het supremum en het infimum worden samengeknepen, I=I+I_- = I_+.

Voorbeeld 15.5 (Stuksgewijs continu, geen drama)

Het gehele deel op [0,3]\intcc{0}{3} is een vermomde trapfunctie: door in zijn sprongen te splitsen krijgen we

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 doen er niet toe: een functie in eindig veel punten wijzigen verandert geen enkele integraal (de omsluitende trapfuncties blijven onaangetast). Dat is de hele inhoud van de uitbreiding “stuksgewijs continu”: snijd in 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 snijd 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 worden op 12\frac12 samengeknepen, en 01x ⁣dx=12\int_0^1 x\,\dd x = \frac12 rechtstreeks uit de definitie. Het afsluitende inzicht: dit is de eerste en de laatste keer dat we uit de definitie integreren — de hoofdstelling hieronder vervangt al zulke berekeningen door het opzoeken van één primitieve, en dat is het hele economische punt van dit hoofdstuk.

Stelling 15.7 (Eigenschappen)

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

  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 afspraak ba=ab\int_b^a = -\int_a^b, geldig voor elke volgorde van de grenzen);
  4. strikte positiviteit: is ff continu, 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 over op de limiet via de definitie met supremum en infimum. De lineariteit verdient het detail eenmaal: zij ε>0\varepsilon > 0 gegeven; omsluit φffψf\varphi_f \leq f \leq \psi_f en φggψg\varphi_g \leq g \leq \psi_g met openingen ε\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 omsluiting door trapfuncties met opening (1+λ)ε\leq (1 + \lambda)\varepsilon, en haar trapintegralen zijn gelijk aan φf+λφg\int \varphi_f + \lambda\int\varphi_g enzovoort (Propositie 15.2): door ε0\varepsilon \to 0 te laten gaan wordt (f+λg)\int(f + \lambda g) op f+λg\int f + \lambda\int g samengeknepen. Voor λ<0\lambda < 0 keert vermenigvuldiging met λ\lambda de omsluiting van gg om — de onderste trapfunctie van λg\lambda g is λψg\lambda \psi_g — en dezelfde insluiting loopt met omgewisselde rollen. De grens ff\abs{\int f} \leq \int \abs f komt uit fff-\abs f \leq f \leq \abs f en de monotonie.

(4) Contrapositie: is f(x0)=m>0f(x_0) = m > 0, dan levert de continuïteit een deelinterval van lengte η>0\eta > 0 waarop fm2f \geq \frac m2; de trapfunctie die daar m2\frac m2 waard is en elders 00, 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 je 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} in π\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: de opheffing is echt, en daarom draagt de uitspraak over strikte positiviteit (Stelling 15.7 (4)) de hypothese f0f \geq 0 — zonder haar bewijst een verdwijnende integraal niets over ff. Het afsluitende inzicht: f\int \abs f meet oppervlakte, f\int f meet tekenbalans; de ongelijkheid ff\abs{\int f} \leq \int\abs f is het exacte verslag van wat de 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. Leg x0Ix_0 \in I en ε>0\varepsilon > 0 vast; de continuïteit in x0x_0 levert een δ\delta 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 beide volgordes van de grenzen). Dus F(x0)=f(x0)F'(x_0) = f(x_0); F=fF' = f is continu, dus FF is C1C^1. Is 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 de pariteit het werk: is ff oneven, dan stuurt de substitutie ttt \mapsto -t de integraal a0f\int_{-a}^{0} f naar 0af-\int_0^a f, dus

aaf(t) ⁣dt=0;is f even, danaaf=20af.\int_{-a}^{a} f(t)\,\dd t = 0 ; \qquad\text{is $f$ even, dan}\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 dat er een primitieve in zicht is (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 de symmetrie voordat je naar technieken grijpt: de snelste integraal is die welke nooit wordt berekend.

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 halveringsformule 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), en dat is precies de afgeleide van tanx2\tan\frac x2:

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

Er was geen substitutiemachinerie nodig — alleen de reflex om een integrand te lezen als iemands afgeleide, waarna de hoofdstelling de rest doet. (Het systematische gereedschap achter zulke goniometrische integralen, de substitutie t=tanx2t = \tan\frac x2, hoort bij de standaarduitrusting die op Stelling 15.15 (2) wordt gebouwd.)

Voorbeeld 15.12 (Functies gedefinieerd door integralen)

De hoofdstelling fabriceert 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, aangenomen); toch bestaat FF, is zij 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 is

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 met dubbele integralen berekend in het volume van bachelorjaar 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 te berekenen)

De integralen Rn=01tn1+t ⁣dtR_n = \int_0^1 \frac{t^n}{1 + t}\,\dd t hebben geen aangename gesloten vorm, en toch pint de monotonie ze nauwkeurig 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 afname (RnR_n \sim een veelvoud van 1n\frac1n, in feite Rn12nR_n \sim \frac{1}{2n}) met twee regels en zonder primitieve. De weekendopgaven van dit hoofdstuk en het volgende draaien op precies zulke omsluitingen — de eerste reflex van de analist tegenover een integraal moet begrens hem zijn, en pas daarna, indien nodig, hem berekenen.

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 12\frac12 maar 2π\frac2\pi als gemiddelde — 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 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, en zo doet de integraal zijn intrede in de fysica.

Stelling 15.15 (Partiële integratie; substitutie)

  1. Zijn u,vu, v van klasse C1C^1 op [a,b]\intcc{a}{b}, dan geldt

    abuv=[uv]ababuv.\int_a^b u'v = \bigl[uv\bigr]_a^b - \int_a^b uv' .
  2. Is φ\varphi van klasse C1C^1 op [α,β]\intcc{\alpha}{\beta} en ff continu op φ([α,β])\varphi(\intcc{\alpha}{\beta}), dan geldt

    αβ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 zijn beide leden 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. (De linearisatie uit Methode 3.11 aan het werk.)

Opmerking 15.17 (Veelgemaakte fouten in de integraalrekening)

(i) Substituties moeten C1C^1 zijn op het hele interval: de verandering x=1tx = \frac1t is ongeoorloofd over 00 heen; blind toegepast op 11 ⁣dx1+x2\int_{-1}^{1}\frac{\dd x}{1 + x^2} “bewijst” zij dat de integraal gelijk is aan zijn eigen tegengestelde. Heeft een substitutie een singulariteit, snijd dan eerst het interval (Chasles), substitueer op elk stuk, en voeg pas daarna weer samen. (ii) Logaritmische primitieven vragen absolute waarden:  ⁣dxx2=lnx2+C\int \frac{\dd x}{x - 2} = \ln\abs{x - 2} + C, aan elke kant van 22 apart — 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 aan weerszijden van de singulariteit verschillen. (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 mag besluiten (Voorbeeld 15.8). (iv) Riemannsommen moeten worden geijkt: 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 verdeling horen — de veelgemaakte 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. Controlelijst voordat je Stelling 15.20 inroept: zonder 1n\frac1n af, herschrijf de sommand als ff van kn\frac kn, benoem ff en controleer haar continuïteit.

Voorbeeld 15.18 (Raden, afleiden, bijstellen)

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 leid af:

(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 (coëfficiënten van u2u^2 neerwaarts passend maken). 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 men anders met twee partiële integraties bereikt. 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 afleiden van een geraden gedaante zet integreren 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 twee keer partieel, waarbij je telkens de goniometrische factor afleidt:

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 bij zichzelf teruggekeerd: I=eπ/21II = \eu^{\pi/2} - 1 - I, waaruit

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

Het afsluitende inzicht: is de integrand een product van twee functies die zich onder afleiden reproduceren (eax\eu^{ax}, cosbx\cos bx, sinbx\sin bx), dan leveren twee partiële integraties een lineaire vergelijking voor de onbekende integraal — los die op in plaats van te integreren; equivalent: ga langs e(a+ib)x\eu^{(a + \iu b)x} (Hoofdstuk 3) en neem het reële deel. Beide wegen geven hetzelfde antwoord, en nagaan dat dit zo is, is een gratis controle.

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. Zij ε>0\varepsilon > 0 gegeven; de uniforme continuïteit (Heine) levert een δ\delta; voor n>baδn > \frac{b-a}{\delta} ligt elk punt van een deelinterval binnen δ\delta van zijn evaluatiepunt, dus is 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 de uniforme continuïteit de oppervlakte van de trap naar ∈t_ab f.
Een linkse Riemannsom met n=8n = 8 rechthoeken: naarmate de maaswijdte krimpt, dwingt de uniforme continuïteit de oppervlakte van de trap 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 die onzichtbaar is voor elementaire afschattingen en doorzichtig als Riemannsom.

Voorbeeld 15.22 (Een tweede Riemannsom, met ijking)

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

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 sommand in de gedaante f(kn)f(\frac kn) ten koste van het afzonderen van precies één factor 1n\frac1n; is de gedaante eenmaal goed, dan doet de stelling de analyse en de hoofdstelling het rekenwerk.

Opmerking 15.23 (Waar de integraal hierna werkt)

Elk van de constructies uit dit hoofdstuk heeft een vervolg. De 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 integraalrestterm is de scherpste vorm van de formule van Taylor in Hoofdstuk 16; de op suprema gebaseerde definitie is het prototype van de Lebesgue-integraal van het volume van bachelorjaar 3, waar dezelfde drie eigenschappen (lineariteit, monotonie, een convergentiestelling) worden herbouwd op een veel grotere klasse functies. En de weekendopgave hieronder maakt van de 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} (splitsing in partieelbreuken, 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}.

Partieel (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.

Partieel (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} in 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 in 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 .

Omdat (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 nadert tot 00.

Oefening 15.5 ★★

(Cauchy–Schwarz) Bewijs voor f,gf, g continu op [a,b]\intcc{a}{b} dat

(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 te ontwikkelen als een tweedegraadsuitdrukking in λ\lambda. Wanneer geldt de 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. Is g2=0\int g^2 = 0, dan is g=0g = 0 (strikte positiviteit, Stelling 15.7 (4)) en luidt de ongelijkheid 000 \leq 0. Anders is QQ een echte tweedegraadsuitdrukking, overal 0\geq 0: haar discriminant is 0\leq 0, dat wil zeggen (fg)2f2g2\bigl(\int fg\bigr)^2 \leq \int f^2 \int g^2.

Gelijkheid geldt dan en slechts dan als de discriminant nul is, dus als Q(λ0)=0Q(\lambda_0) = 0 voor zekere λ0\lambda_0, dat wil zeggen (f+λ0g)2=0\int (f + \lambda_0 g)^2 = 0, dus (opnieuw strikte positiviteit) f=λ0gf = -\lambda_0 g: de gelijkheid geldt precies wanneer ff en gg evenredig zijn.

Oefening 15.6 ★★

Zij ff continu op R\R en TT-periodiek. Bewijs dat aa+Tf\int_a^{a+T} f niet van aa afhangt, 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 ★★

Bewijs dat een continue ff op [0,1]\intcc{0}{1} met 01f=12\int_0^1 f = \frac12 een vast punt heeft in [0,1]\intcc{0}{1}. (Integreer f(x)xf(x) - x en gebruik de 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 .

Werd gg nooit nul, dan zou de tussenwaardestelling een vast teken afdwingen (een continue functie op een interval die beide tekens aanneemt, wordt nul); zeg g>0g > 0. Dan geeft de strikte positiviteit (Stelling 15.7 (4), toegepast op g>0g > 0) dat g>0\int g > 0: in tegenspraak met g=0\int g = 0. Dus is g(c)=0g(c) = 0 voor zekere cc: f(c)=cf(c) = c.

Oefening 15.8 ★★★

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

  1. Bewijs met partiële integratie de recursie nWn=(n1)Wn2n W_n = (n-1) W_{n-2} (n2n \geq 2), en bereken W0,W1W_0, W_1, en 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 daaruit de asymptotische gelijkwaardigheid Wnπ2nW_n \sim \sqrt{\dfrac{\pi}{2n}} af.
Oplossing

Oplossing van Oefening 15.8.

  1. Partieel 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 is (Wn)(W_n) (strikt) dalend en positief. Insluiting met de recursie:

    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} wegens de recursie (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) Stel dat π=ab\pi = \frac ab met a,bNa, b \in \N^*, en stel, voor een nog 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. (Binomiaalontwikkeling: 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). Ga na 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} is 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 is de integrand >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 nadert (de faculteit verslaat de meetkundige term: het is de algemene term van de convergente exponentiële reeks, vergelijk Voorbeeld 11.12); in het bijzonder is 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 is 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, want PP heeft graad 2n2n) is 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 is 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 aanname π=ab\pi = \frac ab faalt: π\pi is irrationaal.

Oefening 15.10 ★★★

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

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

Bewijs haar daarna opnieuw voor ff alleen continu, met uniforme benadering door trapfuncties (Stelling 15.3).

Oplossing

Oplossing van Oefening 15.10.

Het C1C^1-geval: partieel,

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.

Het continue 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 opening ε\leq\varepsilon; neem φ\varphi). Dan is

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 nadert tot 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 aan 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. Omdat g0g \geq 0, is mgfgMgm\,g \leq fg \leq M\,g, dus volgens de monotonie

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

Is abg=0\int_a^b g = 0, dan dwingt de strikte positiviteit (Stelling 15.7 (4)) af dat g0g \equiv 0, zijn beide leden nul, en voldoet elke cc. 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 is t=f(c)t = f(c) voor zekere cc.

Het teken doet ertoe: op [1,1]\intcc{-1}{1} met f(t)=g(t)=tf(t) = g(t) = t is 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 dwingen nulpunten af) 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 nul wordt in n+1n + 1 verschillende punten van (a,b)\intoo{a}{b}. (Wisselt ff alleen van teken in z1<<zmz_1 < \dots < z_m met mnm \leq n, integreer ff dan tegen P(t)=(tz1)(tzm)P(t) = (t - z_1) \cdots (t - z_m) en gebruik de strikte positiviteit.)

Oplossing

Oplossing van Oefening 15.12.

Is f0f \equiv 0, dan is de bewering leeg (elk punt is een nulpunt). Neem dus aan dat f≢0f \not\equiv 0 en stel dat zij hoogstens nn verschillende nulpunten in (a,b)\intoo{a}{b} heeft; zij z1<<zmz_1 < \dots < z_m (mnm \leq n) die nulpunten waar ff van teken wisselt (eventueel geen enkel). Stel P(t)=i=1m(tzi)P(t) = \prod_{i=1}^{m}(t - z_i) (het lege product is 11), van graad mnm \leq n. Op elk deelinterval dat door de ziz_i wordt afgesneden hebben ff en PP een vast teken, en beide klappen om bij het passeren van een ziz_i: het product fPfP heeft één vast teken op heel (a,b)\intoo{a}{b}. Omdat het continu is, niet identiek nul en van vast teken, is 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, die alle nul zijn: tegenspraak. Bijgevolg heeft ff minstens n+1n + 1 verschillende nulpunten in (a,b)\intoo{a}{b}.

Opmerking 15.24 (Vooruitzichten binnen dit volume)

Drie komende hoofdstukken steunen rechtstreeks op dit hoofdstuk. Hoofdstuk 16 draagt de integraalrestterm — de scherpste van de drie formules van Taylor is een nn keer herhaalde partiële integratie. Hoofdstuk 17 zet de omsluiting van sommen door integralen om in het beslissende criterium voor nα\sum n^{-\alpha}, en haar weekendopgave verfijnt die omsluiting 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 gebouwd, zonder dat er oneigenlijke theorie nodig is. Het meest hergebruikte 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 schatting van de integraal.

15.6 Opgave: de integrale irrationaliteitsmachine

Probleem 15.1

Weekendopgave — 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 deze hele opgave 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 πQ\pi \notin \Q te bewijzen; hier industrialiseren we het. De machine heeft drie onderdelen nodig: een invoer van geheeltalligheid (waarden in de randpunten van goedgekozen veeltermen), een invoer van kleinheid (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 het met gecertificeerde fout, bewijzen de scherpere stelling van Legendre dat π2\pi^2 irrationaal is, en eindigen met de charmantste integraal uit 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 insluit.

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. (Bèta-integralen) Bewijs, met inductie op mm en 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 af dat 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}}, 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 twee keer wordt gebruikt: voor elke continue g>0g > 0 op (0,1)\intoo{0}{1} geldt

    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: eerst 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 recursie An=enAn1A_n = \eu - n A_{n-1}, en de grenzen 0<Anen+10 < A_n \leq \dfrac{\eu}{n+1}.
  2. Toon met inductie aan 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 (is e=pq\eu = \frac pq, dan is qAnq A_n voor grote nn een geheel getal gevangen in (0,1)\intoo{0}{1}). Vergelijk met het bewijs van Oefening 11.9: dezelfde clou, andere brandstof.
  4. Bewijs, met inductie en partiële integratie, de exacte formule voor de restterm

    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 e<2.75\eu < 2.75 (uit b2=2.75b_2 = 2.75 in Voorbeeld 11.12), evalueer de som, en besluit tot de gecertificeerde insluiting 2.7182818e2.71828232.7182818 \leq \eu \leq 2.7182823 — zes decimalen, e2.718282\eu \approx 2.718282, met bewijs.

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

  1. Toon aan dat 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 aan 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 daarna 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 aan 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. Ga de telescopering G+π2G=bnπ2n+2f=π2anfG'' + \pi^2 G = b^n \pi^{2n+2} f = \pi^2 a^n f na, en vervolgens

     ⁣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 dat

    π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 randgegevens; analytische kleinheid van de integraal) en de brug, in de delen II en III. Leg daarna uit waarom de factor 1n!\frac{1}{n!} in ff de kern is: als hij wordt weggehaald, overleeft de geheeltalligheid, maar welke ongelijkheid sterft, en voor welke beweerde breuken ab\frac ab faalt het bewijs dan?
  2. Effectiviteit: stel dat iemand beweert dat π2=ab\pi^2 = \frac ab met a10a \leq 10. Toon aan dat de tegenspraak al bij n=7n = 7 landt: bereken π(10/4)7/7!0.38<1\pi\,(10/4)^7/7! \approx 0.38 < 1. De machine weerlegt niet alleen; zij weerlegt in een vast, berekenbaar stadium.
  3. Controleer de brug onvoorwaardelijk: bewijs met 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 breng dit in overeenstemming 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 kernen van de machine? Benoem de eigenschap (elk voldoet aan een lineaire differentiaalvergelijking met constante coëfficiënten, zodat herhaalde partiële integratie naar het begin terugkeert), 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. Voer de veeltermdeling

    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}

    uit, en leid de beroemde identiteit

    01x4(1x)41+x2 ⁣dx=227π\int_0^1 \frac{x^4 (1-x)^4}{1 + x^2}\,\dd x = \frac{22}{7} - \pi

    af.

  2. De integrand is positief: besluit dat π<227\pi < \frac{22}{7}. Bewijs vervolgens, door 11+x2\frac{1}{1+x^2} tussen 12\frac12 en 11 te begrenzen en 01(x(1x))4=1630\int_0^1 (x(1-x))^4 = \frac{1}{630} (vraag 3) te gebruiken, dat

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

    dat wil zeggen 3.14126π3.142073.14126 \leq \pi \leq 3.14207: twee juiste decimalen, met de hand.

  3. Veralgemeen: toon aan, door x4m(1x)4mx^{4m}(1-x)^{4m} door 1+x21 + x^2 te delen, 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 ga na dat m=1m = 1 de vragen 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 garantie van Dirichlet 149\frac{1}{49}, en vermeld π3551132.7107\abs{\pi - \frac{355}{113}} \approx 2.7\cdot10^{-7} tegenover 111327.8105\frac{1}{113^2} \approx 7.8\cdot10^{-5}: er bestaan uitzonderlijk goede rationale benaderingen van π\pi — consistent, want van π\pi is niet bekend dat het slecht benaderbaar is.
  5. (De gehele-getallenval, geabstraheerd) Bewijs het lemma dat alles verenigt: bestaan er voor xRx \in \R gehele getallen an,bna_n, b_n met 0<an+bnx00 < \abs{a_n + b_n x} \to 0, dan is xx irrationaal. Som de gevallen ervan op in deze opgave, in Oefening 15.9, in Oefening 11.9 en in Probleem 14.1.
  6. Synthese, telkens één zin: (i) de drie onderdelen van de machine en waar elk in het gereedschap van dit hoofdstuk huist; (ii) wat de integraal bijdraagt wat de middelwaardestelling van Probleem 14.1 niet kon; (iii) de inventaris van de geoogste resultaten (twee irrationaliteiten, één constante tot op zes decimalen, één insluiting van π\pi, één binomiale grens); (iv) de grens (Hermite, Lindemann; en dezelfde val, losgelaten op ζ(2)\zeta(2) en ζ(3)\zeta(3), in de rekenkunde van de twintigste eeuw).
Oplossing

Oplossing van Probleem 15.1.

1. Zij N=2cN = \lceil 2c \rceil. Voor nNn \geq N is 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: insluiting.

2. Leg kk vast; 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, partieel (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)!} ,

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

3. k=m=nk = m = n: 01(x(1x))n=(n!)2(2n+1)!=1(2n+1)(2nn)\int_0^1 (x(1-x))^n = \frac{(n!)^2}{(2n+1)!} = \frac{1}{(2n+1)\binom{2n}{n}}. Omdat 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} — in overeenstemming 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 de 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. Partieel: 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. Is An1=αn1+βn1eA_{n-1} = \alpha_{n-1} + \beta_{n-1}\eu met gehele ingangen, dan is

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. Is e=pq\eu = \frac pq, dan is 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 gehele getal q!pqq!aqq!\,\frac pq - q!\,a_q; hier is het qAnqA_n: dezelfde val, integrale brandstof.

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

zodat de formule zich van nn naar n+1n + 1 voortplant. 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 is 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} is 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 is f(k)(0)=0f^{(k)}(0) = 0 voor k<nk < n of k>2nk > 2n, en voor nk2nn \leq k \leq 2n is 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 zijn 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) gehele getallen volgens vraag 11.

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

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 is

(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. Integratie over [0,1]\intcc{0}{1} geeft

π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 samen met vraag 10 begrenst het door πan4nn!\frac{\pi a^n}{4^n\,n!}.

15. Volgens vraag 1 (met c=a4c = \frac a4) nadert πan4nn!\frac{\pi a^n}{4^n n!} tot 00: voor grote nn is het <1< 1, in tegenspraak met G(0)+G(1)1G(0) + G(1) \geq 1. Geen enkele breuk ab\frac ab is dus gelijk aan π2\pi^2: de stelling van Legendre. Was π\pi rationaal, dan zou π2\pi^2 dat ook zijn: dus πQ\pi \notin \Q — en de implicatie loopt alleen deze kant op (2\sqrt2 is irrationaal met een rationaal kwadraat), en daarom is π2Q\pi^2 \notin \Q strikt sterker dan Oefening 15.9.

16. Geheeltalligheid: de vragen 11–12 (afgeleiden in de randpunten); kleinheid: de vragen 10 en 1; brug: de vragen 13–14 (de getelescopeerde dubbele partiële integratie). Zonder 1n!\frac{1}{n!} blijven de randgegevens geheel (nog gemakkelijker zelfs), maar wordt de grens πan4n\frac{\pi a^n}{4^n}, die alleen naar 00 nadert 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 is de grens π×610.3550400.38<1\frac{\pi \times 610.35}{5040} \approx 0.38 < 1 (bij n=6n = 6 is zij nog 1.071.07): de tegenspraak landt expliciet in het zevende stadium.

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 telescopering bij n=1n = 1 (zonder enige 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 stemmen 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 (respectievelijk binnen de gehele veeltermen in π2\pi^2). Die geslotenheid is wat de machine nodig heeft. Verfijnd met kernen die op meerdere punten tegelijk zijn afgestemd, levert hetzelfde mechanisme de stelling van Hermite (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 .

Integratie van de getoonde identiteit gedeeld door 1+x21 + x^2 geeft

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 en positief op (0,1)\intoo{0}{1}: de integraal is >0> 0, dus π<227\pi < \frac{22}{7}. Bovendien is 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, zodat (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). Deling van de identiteit door 1+x21 + x^2, gevolgd door integratie, geeft

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 is   π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 de 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 ijkwaarde van orde 22, 1722.0102\frac{1}{7^2} \approx 2.0\cdot10^{-2}, die Dirichlet garandeert (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. Dat spreekt niets van het bewezene tegen: ongelijkheden van Liouville begrenzen benaderingsfouten alleen van onderen voor algebraïsche getallen, en zo’n grens voor π\pi is op dit niveau niet beschikbaar — π\pi mag spectaculair goed benaderd worden.

24. Lemma: stel 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}, waarbij qan+pbnq a_n + p b_n een niet-nul geheel getal is (niet nul omdat de absolute waarde >0> 0 is): dus is 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 (opnieuw x=ex = \eu, 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 de meetkundige vorm ervan. In deel III en in Oefening 15.9 loopt de val binnen de tegenspraak: de aanname van rationaliteit maakt van een uitdrukking een geheel getal, dat de analyse vervolgens in (0,1)\intoo{0}{1} perst — hetzelfde principe, getransponeerd.

25. (i) De geheeltalligheid huist in de randberekening met veeltermen (de vragen 6, 11–12), de kleinheid in de grenzen die de monotonie met suprema geeft (de vragen 4, 10), de brug in de partiële integratie (de 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 slechts een ongelijkheid met een onbekende cc), en daarom reikt de machine tot π2\pi^2 terwijl Probleem 14.1 slechts benaderingsexponenten bereikte. (iii) Geoogst: eQ\eu \notin \Q, π2Q\pi^2 \notin \Q (en 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) De grens: Hermite en Lindemann duwen dezelfde machine door tot transcendentie; en Apéry (1979) liet de gehele-getallenval draaien op ζ(3)\zeta(3) — de machine produceert nog steeds wiskunde van de twintigste eeuw.

Begrippen gedefinieerd in dit hoofdstuk

Bekijk alle 395 begrippen in de begrippenlijst