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<b reële getallen.
15.1 Trapfuncties
Definitie 15.1
φ:[a,b]→R heet een trapfunctie wanneer er een verdeling a=x0<x1<⋯<xn=b bestaat waarvoor φ constant is, gelijk aan ci, op elk openinterval(xi−1,xi) (de waarden in de knooppunten zijn vrij). Haar integraal is
∫abφ=i=1∑nci(xi−xi−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 (φ≤ψ⟹∫φ≤∫ψ), en voldoet hij aan de relatie van Chasles ∫ab=∫ac+∫cb voor a<c<b.
Bewijs. De motor is de invariantie onder verfijning, vermeld in de definitie: het invoegen van één extra knooppunt t∈(xi−1,xi) in een verdeling vervangt de term ci(xi−xi−1) door ci(t−xi−1)+ci(xi−t) — hetzelfde getal — zodat de integraal onder elke eindige verfijning onveranderd blijft. Neem nu φ met verdeling σ en ψ met verdeling σ′: op de gemeenschappelijke verfijning σ∪σ′ zijn beide trapfuncties met dezelfde knooppunten, en op elk stuk is φ+λψ constant gelijk aan ci+λdi: de lineariteit herleidt zich tot die van eindige sommen. Stijging: ci≤di op elk stuk geeft ∑ciΔi≤∑diΔi (lengten Δi≥0). Chasles: voeg het knooppunt c in en splits de som daar. ∎
Bewijs. Volgens de stelling van Heine (Stelling 13.22) is funiform continu: kies δ bij ε, en een verdeling met maaswijdte <δ (bijvoorbeeld gelijkmatig, met n>δb−a stukken). Op elk gesloten stuk [xi−1,xi] bereikt f een minimum mi en een maximum Mi (Stelling 13.13), en is Mi−mi≤ε (de twee extremale punten liggen binnen δ van elkaar). Definieer φ=mi en ψ=Mi op (xi−1,xi) (en φ=ψ=f in de knooppunten). ∎
zijn gelijk; hun gemeenschappelijke waarde is de integraal∫abf (ook geschreven als ∫abf(t)dt). Hij valt op trapfuncties samen met het vorige begrip, en breidt uit tot stuksgewijs continue functies door [a,b] in de discontinuïteitspunten te splitsen (daar is Chasles een definitie).
Bewijs. Beide verzamelingen zijn niet leeg (f is begrensd) en elke onderste trapintegraal is ≤ elke bovenste (monotonie op trapfuncties): dus I−(f)≤I+(f). Volgens Stelling 15.3 is er voor elke ε een paar met ∫ψ−∫φ≤ε(b−a): het supremum en het infimum worden samengeknepen, I−=I+. ∎
Voorbeeld 15.5(Stuksgewijs continu, geen drama)
Het gehele deel op [0,3] is een vermomde trapfunctie: door in zijn sprongen te splitsen krijgen we
∫03⌊t⌋dt=∫010+∫121+∫232=0+1+2=3,
en de waarden in de sprongpunten 1,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)=x op [0,1] en snijd in n gelijke stukken. De beste trapfuncties die constant zijn op de stukken, zijn φ=nk−1 en ψ=nk op het k-de stuk, met
Elke onderste integraal is ≤I−(f)≤I+(f)≤ elke bovenste, dus 2nn−1≤I−(f)≤I+(f)≤2nn+1 voor alle n: beide worden op 21 samengeknepen, en ∫01xdx=21 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,gcontinu (of stuksgewijs continu) op [a,b] en λ∈R geldt:
lineariteit: ∫(f+λg)=∫f+λ∫g;
monotonie: f≤g⟹∫abf≤∫abg; en ∫abf≤∫ab∣f∣≤(b−a)sup∣f∣;
Chasles: ∫ab=∫ac+∫cb (met de afspraak ∫ba=−∫ab, geldig voor elke volgorde van de grenzen);
strikte positiviteit: is fcontinu, f≥0 en ∫abf=0, dan is f=0 overal op [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 gegeven; omsluit φf≤f≤ψf en φg≤g≤ψg met openingen ≤ε (Stelling 15.3). Voor λ≥0 is φf+λφg≤f+λg≤ψf+λψg een omsluiting door trapfuncties met opening ≤(1+λ)ε, en haar trapintegralen zijn gelijk aan ∫φf+λ∫φg enzovoort (Propositie 15.2): door ε→0 te laten gaan wordt ∫(f+λg) op ∫f+λ∫g samengeknepen. Voor λ<0keert vermenigvuldiging met λ de omsluiting van g om — de onderste trapfunctie van λg is λψg — en dezelfde insluiting loopt met omgewisselde rollen. De grens ∫f≤∫∣f∣ komt uit −∣f∣≤f≤∣f∣ en de monotonie.
(4) Contrapositie: is f(x0)=m>0, dan levert de continuïteit een deelinterval van lengte η>0 waarop f≥2m; de trapfunctie die daar 2m waard is en elders 0, is ≤f, dus ∫f≥2mη>0. ∎
Voorbeeld 15.8(Chasles aan het werk: integralen met absolute waarden)
Om een absolute waarde te integreren, snijd je waar het teken verandert.
∫02∣x−1∣dx=∫01(1−x)dx+∫12(x−1)dx=21+21=1,
en door [0,2π] in π te snijden:
∫02π∣sint∣dt=∫0πsintdt−∫π2πsintdt=2+2=4,
terwijl ∫02πsintdt=0: de opheffing is echt, en daarom draagt de uitspraak over strikte positiviteit (Stelling 15.7 (4)) de hypothese f≥0 — zonder haar bewijst een verdwijnende integraal niets over f. Het afsluitende inzicht: ∫∣f∣ meet oppervlakte, ∫f meet tekenbalans; de ongelijkheid ∫f≤∫∣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)
waarvan de absolute waarde ≤∣h∣1⋅∣h∣ε=ε is (grens (2), geldig voor beide volgordes van de grenzen). Dus F′(x0)=f(x0); F′=f is continu, dus F is C1. Is ook G′=f, dan is (G−F)′=0 op het interval, dus G=F+c (Gevolg 14.12), en G(b)−G(a)=F(b)−F(a)=∫abf. ∎
Voorbeeld 15.10(Symmetrie vóór berekening)
Op een symmetrisch interval doet de pariteit het werk: is f oneven, dan stuurt de substitutie t↦−t de integraal∫−a0f naar −∫0af, dus
∫−aaf(t)dt=0;is f even, dan∫−aaf=2∫0af.
Zo is ∫−111+t4t3costdt=0 zonder dat er een primitieve in zicht is (de integrand is oneven), en ∫−ππt2costdt=2∫0πt2costdt. 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π/21+cosxdx. De halveringsformule 1+cosx=2cos22x maakt van de integrand 21(1+tan22x), en dat is precies de afgeleide van tan2x:
∫0π/21+cosxdx=[tan2x]0π/2=tan4π=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=tan2x, 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)=∫0xe−t2dt.
Geen enkele combinatie van klassieke functies heeft afgeleidee−t2 (een stelling van Liouville, aangenomen); toch bestaat F, is zij C1 met F′(x)=e−x2>0, strikt stijgend, oneven (substitueer t↦−t) en begrensd: voor x≥1 is
F(x)−F(1)=∫1xe−t2dt≤∫1xe−tdt≤e−1,
dus F≤F(1)+e−1≤1+e−1. (De exacte limiet, 2π, wordt met dubbele integralen berekend in het volume van bachelorjaar 3.) Kettingregel voor bewegende grenzen: dxd∫xx2e−t2dt=2xe−x4−e−x2. 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 integralenRn=∫011+ttndt hebben geen aangename gesloten vorm, en toch pint de monotonie ze nauwkeurig vast: op [0,1] is 21≤1+t1≤1, dus
2(n+1)1=21∫01tndt≤Rn≤∫01tndt=n+11:
de exacte orde van afname (Rn∼ een veelvoud van n1, in feite Rn∼2n1) 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 continuef over [a,b] is b−a1∫abf. Voor de boog van de sinus:
π1∫0πsintdt=π1[−cost]0π=π2≈0.637:
een volledige positieve boog heeft niet 21 maar π2 als gemiddelde — de kromme brengt meer tijd hoog door dan een driehoek zou doen. Volgens Oefening 15.11 (middelwaardestelling voor integralen, g=1) is het gemiddelde een waarde: sinc=π2 voor zekere c∈(0,π). En volgens de Riemannsommen van dit hoofdstuk is het gemiddelde de limiet van gewone gemiddelden van n 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)
Zijn u,v van klasse C1 op [a,b], dan geldt
∫abu′v=[uv]ab−∫abuv′.
Is φ van klasse C1 op [α,β] en fcontinu op φ([α,β]), dan geldt
∫αβf(φ(t))φ′(t)dt=∫φ(α)φ(β)f(x)dx.
Bewijs. (1) (uv)′=u′v+uv′; integreer over [a,b] en pas de hoofdstelling toe op de C1-functie uv.
(2) Zij F een primitieve van f op het beeldinterval (Stelling 15.9). Dan is (F∘φ)′=(f∘φ)φ′ (kettingregel), dus zijn beide leden gelijk aan F(φ(β))−F(φ(α)). ∎
Voorbeeld 15.16
∫01tetdt=[tet]01−∫01etdt=e−(e−1)=1. En met de substitutie x=sint (t∈[0,2π]):
— 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 C1 zijn op het hele interval: de verandering x=t1 is ongeoorloofd over 0 heen; blind toegepast op ∫−111+x2dx “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: ∫x−2dx=ln∣x−2∣+C, aan elke kant van 2 apart — ln(x−2) schrijven op (0,1) is de logaritme van een negatief getal schrijven; en de constante C mag aan weerszijden van de singulariteit verschillen. (iii) Een verdwijnende integraal doodt de functie niet: ∫02πsin=0; positiviteit van de integrand is vereist voordat men f=0 mag besluiten (Voorbeeld 15.8). (iv) Riemannsommen moeten worden geijkt: in nb−a∑f(a+knb−a) moeten de stap buiten en de punten binnen bij dezelfde verdeling horen — de veelgemaakte fout is een som ∑k=1nf(nk)zonder de factor n1, die divergeert in plaats van naar ∫01f te convergeren. Controlelijst voordat je Stelling 15.20 inroept: zonder n1 af, herschrijf de sommand als f van nk, benoem f en controleer haar continuïteit.
Voorbeeld 15.18(Raden, afleiden, bijstellen)
Wat is ∫1x(lnt)2dt? Raad een primitieve van de vorm tP(lnt) met P een veelterm en leid af:
(tP(lnt))′=P(lnt)+P′(lnt).
We hebben P(u)+P′(u)=u2 nodig: neem P(u)=u2−2u+2 (coëfficiënten van u2 neerwaarts passend maken). Bijgevolg
een resultaat dat men anders met twee partiële integraties bereikt. Het afsluitende inzicht: voor integranden van de vorm (veelterm in lnt) of (veelterm maal eλ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π/2excosxdx. Integreer twee keer partieel, waarbij je telkens de goniometrische factor afleidt:
De integraal is bij zichzelf teruggekeerd: I=eπ/2−1−I, waaruit
I=2eπ/2−1.
Het afsluitende inzicht: is de integrand een product van twee functies die zich onder afleiden reproduceren (eax, cosbx, sinbx), 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 (Hoofdstuk 3) en neem het reële deel. Beide wegen geven hetzelfde antwoord, en nagaan dat dit zo is, is een gratis controle.
en evenzo met willekeurige evaluatiepunten binnen de deelintervallen.
Bewijs.Sn is de integraal van de trapfunctieφn die gelijk is aan f(a+knb−a) op het k-de deelinterval. Zij ε>0 gegeven; de uniforme continuïteit (Heine) levert een δ; voor n>δb−a ligt elk punt van een deelinterval binnen δ van zijn evaluatiepunt, dus is ∣f−φn∣≤ε op [a,b], waaruit
∫abf−Sn=∫ab(f−φn)≤(b−a)ε.
∎
Een linkse Riemannsom met n=8 rechthoeken: naarmate de maaswijdte krimpt, dwingt de uniforme continuïteit de oppervlakte van de trap naar ∫abf.
Voorbeeld 15.21
k=1∑nn+k1=n1k=1∑n1+k/n1n→∞∫011+xdx=ln2: een limiet die onzichtbaar is voor elementaire afschattingen en doorzichtig als Riemannsom.
Voorbeeld 15.22(Een tweede Riemannsom, met ijking)
Het afsluitende inzicht: de hele kunst zit in de middelste regel — dwing de sommand in de gedaante f(nk) ten koste van het afzonderen van precies één factor n1; 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 ∑nα1 met ∫tαdt); 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.
15.5 Oefeningen
Oefening 15.1★
Bereken: ∫01x2−4dx(splitsing in partieelbreuken, Hoofdstuk 9); ∫1elntdt; ∫0πtsintdt.
un=n1∑k=1n1+(k/n)21: een Riemannsom van x↦1+x21 op [0,1], dus un→∫011+x2dx=arctan1=4π.
lnvn=n1∑k=1nlnnn+k=n1∑k=1nln(1+nk)→∫01ln(1+x)dx=[(1+x)ln(1+x)−x]01=2ln2−1. Bijgevolg vn→e2ln2−1=e4. (Controle van de identificatie: n!nn(2n)!=∏k=1nnn+k.)
Oefening 15.4★
Zij fcontinu op [0,1]. Bereken limn→∞∫01xnf(x)dx. (Snijd [0,1] in 1−δ.)
Oplossing
Oplossing van Oefening 15.4.
De limiet is 0. Zij M=sup∣f∣ en ε∈(0,1). Snijd in 1−ε:
∫01xnf≤∫01−εxn∣f∣+∫1−ε1xn∣f∣≤M(1−ε)n+Mε.
Omdat (1−ε)n→0 (Oefening 11.3), is de limsup van het linkerlid ≤Mε voor elke ε: de integraal nadert tot 0.
Oefening 15.5★★
(Cauchy–Schwarz) Bewijs voor f,gcontinu op [a,b] dat
(∫abfg)2≤∫abf2⋅∫abg2,
door ∫ab(f+λg)2≥0 te ontwikkelen als een tweedegraadsuitdrukking in λ. Wanneer geldt de gelijkheid?
Oplossing
Oplossing van Oefening 15.5.
Q(λ)=∫ab(f+λg)2=∫f2+2λ∫fg+λ2∫g2≥0 voor alle λ. Is ∫g2=0, dan is g=0 (strikte positiviteit, Stelling 15.7 (4)) en luidt de ongelijkheid 0≤0. Anders is Q een echte tweedegraadsuitdrukking, overal ≥0: haar discriminant is ≤0, dat wil zeggen (∫fg)2≤∫f2∫g2.
Gelijkheid geldt dan en slechts dan als de discriminant nul is, dus als Q(λ0)=0 voor zekere λ0, dat wil zeggen ∫(f+λ0g)2=0, dus (opnieuw strikte positiviteit) f=−λ0g: de gelijkheid geldt precies wanneer f en g evenredig zijn.
Oefening 15.6★★
Zij fcontinu op R en T-periodiek. Bewijs dat ∫aa+Tf niet van a afhangt, en dat x1∫0xf(t)dt→T1∫0Tf als x→+∞.
Oplossing
Oplossing van Oefening 15.6.
Zij Φ(a)=∫aa+Tf. Volgens de hoofdstelling (Stelling 15.9) is Φafleidbaar met Φ′(a)=f(a+T)−f(a)=0: constant.
Voor x>0 schrijf x=nT+r, 0≤r<T (n=⌊x/T⌋). Chasles:
Dan is x1∫0xf=xnT⋅T1∫0Tf+O(x1), en xnT→1: de limiet is T1∫0Tf.
Oefening 15.7★★
Bewijs dat een continuef op [0,1] met ∫01f=21 een vast punt heeft in [0,1]. (Integreer f(x)−x en gebruik de strikte positiviteit, Stelling 15.7 (4), via haar contrapositie gecombineerd met de tussenwaardestelling.)
Werd g nooit nul, dan zou de tussenwaardestelling een vast teken afdwingen (een continue functie op een interval die beide tekens aanneemt, wordt nul); zeg g>0. Dan geeft de strikte positiviteit (Stelling 15.7 (4), toegepast op g>0) dat ∫g>0: in tegenspraak met ∫g=0. Dus is g(c)=0 voor zekere c: f(c)=c.
Bewijs met partiële integratie de recursie nWn=(n−1)Wn−2 (n≥2), en bereken W0,W1, en daarna W2p en W2p+1 in gesloten vorm.
Bewijs dat (Wn) dalend is met WnWn+1→1; bewijs dat de grootheid (n+1)Wn+1Wn constant is, gelijk aan 2π; en leid daaruit de asymptotische gelijkwaardigheid Wn∼2nπ af.
Op (0,2π) is 0<sint<1, dus sinn+1<sinn en is (Wn) (strikt) dalend en positief. Insluiting met de recursie:
n+1n=Wn−1Wn+1≤WnWn+1≤1⟹WnWn+1→1.
Invariant: an=(n+1)Wn+1Wn voldoet aan an=an−1 wegens de recursie (n+1)Wn+1=nWn−1, dus an=a0=1⋅W1W0=2π. Dan
nWn2∼(n+1)Wn+1Wn=2π⟹Wn∼2nπ.
Oefening 15.9★★★
(Niven: π is irrationaal) Stel dat π=ba met a,b∈N∗, en stel, voor een nog te kiezen n,
P(x)=n!xn(a−bx)n,In=∫0πP(x)sinxdx.
Bewijs dat 0<In≤πn!(πa)n, wat <1 is voor grote n.
Bewijs dat P en al haar afgeleidengehele waarden aannemen in 0 en in π=ba. (Binomiaalontwikkeling: de coëfficiënten van P maal k! zijn geheel voor k≥n; en P(π−x)=P(x).)
Stel Q=P−P′′+P(4)−… (een eindige som). Ga na dat (Q′sinx−Qcosx)′=Psinx, en leid af dat In=Q(π)+Q(0) een geheel getal is.
Besluit.
Oplossing
Oplossing van Oefening 15.9.
Op (0,π) is x>0, a−bx=b(ba−x)=b(π−x)>0 en sinx>0, dus is de integrand >0 en In>0 (strikte positiviteit). Grens: op [0,π] is x≤π en a−bx≤a, dus P≤n!πnan en In≤πn!(πa)n, wat naar 0 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<1 voor grote n.
Ontwikkel xn(a−bx)n=∑j=0n(jn)an−j(−b)jxn+j: dus P=n!1∑jcjxn+j met gehele cj. Dan is P(k)(0)=0 voor k<n (valuatie) en, voor n≤k≤2n, P(k)(0)=n!k!ck−n, een geheel getal omdat n!∣k!. Bovendien is P(π−x)=P(x) (substitueer: π−x verwisselt de factoren, met a−b(π−x)=bx), dus P(k)(π)=±P(k)(0): eveneens gehele getallen.
Met Q=P−P′′+P(4)−… (eindig, want P heeft graad 2n) is Q+Q′′=P, en
(Q′sinx−Qcosx)′=(Q+Q′′)sinx=Psinx.
Bijgevolg is In=[Q′sinx−Qcosx]0π=Q(π)+Q(0), een som van waarden P(2k) in 0 en π: een geheel getal volgens (2).
Voor grote n is In een geheel getal met 0<In<1: onmogelijk. De aanname π=ba faalt: π is irrationaal.
Oefening 15.10★★★
Zij f van klasse C1 op [a,b]. Bewijs met partiële integratie de limiet van het type Riemann–Lebesgue
en elke ∫sinλt=λcosλxi−1−cosλxi≤λ2: de tweede term nadert tot 0. Bijgevolg is de limsup ≤(b−a)ε voor elke ε: de limiet is 0.
Oefening 15.11★★
(Middelwaardestelling voor integralen) Zij f,gcontinu op [a,b] met g≥0. Bewijs dat er een c∈[a,b] bestaat met
∫abf(t)g(t)dt=f(c)∫abg(t)dt,
en toon met een voorbeeld aan dat de hypothese g≥0 niet kan worden weggelaten.
Oplossing
Oplossing van Oefening 15.11.
Zij m=minf en M=maxf, bereikt volgens de extremumstelling. Omdat g≥0, is mg≤fg≤Mg, dus volgens de monotonie
m∫abg≤∫abfg≤M∫abg.
Is ∫abg=0, dan dwingt de strikte positiviteit (Stelling 15.7 (4)) af dat g≡0, zijn beide leden nul, en voldoet elke c. Anders ligt t=∫g∫fg in [m,M]=f([a,b]) (Stellingen 13.13 en 13.10), dus is t=f(c) voor zekere c.
Het teken doet ertoe: op [−1,1] met f(t)=g(t)=t is ∫fg=∫−11t2=32, terwijl f(c)∫−11tdt=0 voor elke c.
Oefening 15.12★★★
(Momenten dwingen nulpunten af) Zij fcontinu op [a,b] met
∫abf(t)tkdt=0voor k=0,1,…,n.
Bewijs dat f nul wordt in n+1 verschillende punten van (a,b). (Wisselt f alleen van teken in z1<⋯<zm met m≤n, integreer f dan tegen P(t)=(t−z1)⋯(t−zm) en gebruik de strikte positiviteit.)
Oplossing
Oplossing van Oefening 15.12.
Is f≡0, dan is de bewering leeg (elk punt is een nulpunt). Neem dus aan dat f≡0 en stel dat zij hoogstens n verschillende nulpunten in (a,b) heeft; zij z1<⋯<zm (m≤n) die nulpunten waar fvan teken wisselt (eventueel geen enkel). Stel P(t)=∏i=1m(t−zi) (het lege product is 1), van graad m≤n. Op elk deelinterval dat door de zi wordt afgesneden hebben f en P een vast teken, en beide klappen om bij het passeren van een zi: het product fP heeft één vast teken op heel (a,b). Omdat het continu is, niet identiek nul en van vast teken, is ∫abfP>0 (strikte positiviteit toegepast op ∣fP∣). Maar ∫fP is een lineaire combinatie van de momenten ∫ftk, k≤n, die alle nul zijn: tegenspraak. Bijgevolg heeft f minstens n+1 verschillende nulpunten in (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 n keer herhaalde partiële integratie. Hoofdstuk 17 zet de omsluiting van sommen door integralen om in het beslissende criterium voor ∑n−α, 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)2dt, 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≤(b−a)sup∣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 en π2 zijn irrationaal, e tot op zes decimalen, en 722>π met bewijs
Eén mechanisme drijft deze hele opgave aan: een uitdrukking die een positief geheel getal moet zijn maar aantoonbaar kleiner is dan 1, kan niet bestaan. Oefening 15.9 (Niven) liet het één keer draaien om π∈/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 n!1 die de integraal verplettert) en een brug (partiële integratie) die ze verbindt. We bewijzen dat e irrationaal is en berekenen het met gecertificeerde fout, bewijzen de scherpere stelling van Legendre dat π2 irrationaal is, en eindigen met de charmantste integraal uit de analyse: ∫011+x2x4(1−x)4dx=722−π, die π met de hand insluit.
Deel I — Brandstof.
Bewijs dat n!cn→0 voor elke vaste c>0(voorbij n≥2c halveert elke stap de term minstens).
Leid af dat ∫01(x(1−x))ndx=(2n+1)(n2n)1, en — door te vergelijken met de grens x(1−x)≤41 — de schatting (n2n)≥2n+14n, in overeenstemming met (n2n)1/n→4 uit Probleem 11.1.
Bewijs het kleinheidslemma dat hieronder twee keer wordt gebruikt: voor elke continueg>0 op (0,1) geldt
0<∫01(x(1−x))ng(x)dx≤4nsup[0,1]∣g∣.
Deel II — e: eerst irrationaliteit, dan zes decimalen. Stel An=∫01xnexdx.
Bereken A0 en A1, bewijs de recursie An=e−nAn−1, en de grenzen 0<An≤n+1e.
Toon met inductie aan dat An=αn+βne met αn,βn∈Z.
Leid af dat e irrationaal is (is e=qp, dan is qAn voor grote n een geheel getal gevangen in (0,1)). Vergelijk met het bewijs van Oefening 11.9: dezelfde clou, andere brandstof.
Bewijs, met inductie en partiële integratie, de exacte formule voor de restterm
Neem n=9: begrens R9 met e<2.75 (uit b2=2.75 in Voorbeeld 11.12), evalueer de som, en besluit tot de gecertificeerde insluiting 2.7182818≤e≤2.7182823 — zes decimalen, e≈2.718282, met bewijs.
Deel III — De stelling van Legendre: π2 is irrationaal. Zij f(x)=n!xn(1−x)n, en stel dat π2=ba met a,b∈N∗.
Toon aan dat f(1−x)=f(x) en 0<f≤4nn!1 op (0,1).
Toon aan dat f(k)(0) en f(k)(1) gehele getallen zijn voor elke k≥0(ontwikkel xn(1−x)n met gehele coëfficiënten; n!k!∈Z voor k≥n; gebruik daarna de symmetrie).
Definieer
G=bnk=0∑n(−1)kπ2n−2kf(2k).
Toon aan dat G(0) en G(1) gehele getallen zijn (elke bnπ2n−2k=an−kbk).
Ga de telescopering G′′+π2G=bnπ2n+2f=π2anf na, en vervolgens
dxd(G′(x)sinπx−πG(x)cosπx)=π2anf(x)sinπx.
Integreer over [0,1] en besluit dat
πan∫01f(x)sin(πx)dx=G(0)+G(1)∈Z,
een positief geheel getal begrensd door 4nn!πan.
Besluit met vraag 1 dat π2 irrationaal is (Legendre, 1794), en dat dit Oefening 15.9 versterkt: waarom impliceert de irrationaliteit van π2 die van π, en niet omgekeerd?
Deel IV — De machine begrijpen.
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 n!1 in f de kern is: als hij wordt weggehaald, overleeft de geheeltalligheid, maar welke ongelijkheid sterft, en voor welke beweerde breuken ba faalt het bewijs dan?
Effectiviteit: stel dat iemand beweert dat π2=ba met a≤10. Toon aan dat de tegenspraak al bij n=7 landt: bereken π(10/4)7/7!≈0.38<1. De machine weerlegt niet alleen; zij weerlegt in een vast, berekenbaar stadium.
en breng dit in overeenstemming met vraag 14 bij n=1 (houd π2 symbolisch: de getelescopeerde identiteit luidt π3∫01f1sinπx=−(f1′′(0)+f1′′(1))=4).
Wat maakt ex en sinπ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 en πtranscendent zijn — voorbij dit volume.
Deel V — 722 tegenover π, en de moraal.
Voer de veeltermdeling
1+x2x4(1−x)4=x6−4x5+5x4−4x2+4−1+x24
uit, en leid de beroemde identiteit
∫011+x2x4(1−x)4dx=722−π
af.
De integrand is positief: besluit dat π<722. Bewijs vervolgens, door 1+x21 tussen 21 en 1 te begrenzen en ∫01(x(1−x))4=6301 (vraag 3) te gebruiken, dat
722−6301≤π≤722−12601,
dat wil zeggen 3.14126≤π≤3.14207: twee juiste decimalen, met de hand.
Veralgemeen: toon aan, door x4m(1−x)4m door 1+x2 te delen, dat de rest de constante (−4)m is (werk modulo x2+1: (1−x)2≡−2x), leid rationale getallen rm af met
∣π−rm∣≤41−5m,
en ga na dat m=1 de vragen 20–21 reproduceert.
Confronteer deze rationale getallen met de benaderingstheorie van Probleem 14.1: bereken π−722≈1.26⋅10−3 tegenover de garantie van Dirichlet 491, en vermeld π−113355≈2.7⋅10−7 tegenover 11321≈7.8⋅10−5: er bestaan uitzonderlijk goede rationale benaderingen van π — consistent, want van π is niet bekend dat het slecht benaderbaar is.
(De gehele-getallenval, geabstraheerd) Bewijs het lemma dat alles verenigt: bestaan er voor x∈R gehele getallen an,bn met 0<∣an+bnx∣→0, dan is x irrationaal. Som de gevallen ervan op in deze opgave, in Oefening 15.9, in Oefening 11.9 en in Probleem 14.1.
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 π, één binomiale grens); (iv) de grens (Hermite, Lindemann; en dezelfde val, losgelaten op ζ(2) en ζ(3), in de rekenkunde van de twintigste eeuw).
Oplossing
Oplossing van Probleem 15.1.
1. Zij N=⌈2c⌉. Voor n≥N is cn/n!cn+1/(n+1)!=n+1c≤21, dus 0<n!cn≤N!cN2−(n−N)→0: insluiting.
2. Leg k vast; inductie op m. Voor m=0: ∫01xk=k+11=(k+1)!k!0!. Stap, partieel (u=(1−x)m, v′=xk):
3.k=m=n: ∫01(x(1−x))n=(2n+1)!(n!)2=(2n+1)(n2n)1. Omdat x(1−x)≤41 op [0,1], is de integraal≤4−n, waaruit (n2n)≥2n+14n — in overeenstemming met (n2n)1/n→4 (Probleem 11.1).
4. De integrand is continu, ≥0 en positief op (0,1), dus niet identiek nul: zijn integraal is >0 (Stelling 15.7 (4)). Bovengrens: (x(1−x))n≤4−n en g≤sup∣g∣, daarna de monotonie.
5.A0=e−1; A1=[xex]01−∫01ex=e−(e−1)=1. Partieel: An=[xnex]01−n∫01xn−1ex=e−nAn−1. Grenzen: de integrand is positief, dus An>0; en ex≤e geeft An≤e∫01xn=n+1e.
6.A0=−1+1⋅e. Is An−1=αn−1+βn−1e met gehele ingangen, dan is
7. Is e=qp, dan is qAn=qαn+pβn∈Z, en 0<qAn≤n+1qe<1 voor grote n: een geheel getal strikt tussen 0 en 1 — onmogelijk. Dus e∈/Q. In Oefening 11.9 was het gevangen gehele getal q!qp−q!aq; hier is het qAn: dezelfde val, integrale brandstof.
8.n=0: R0=∫01et=e−1, dus e=1+R0. Partieel (u=et, v=−n+1(1−t)n+1):
zodat de formule zich van n naar n+1 voortplant. Grenzen: 1≤et≤e op [0,1] en ∫01(1−t)n=n+11 geven (n+1)!1≤Rn≤(n+1)!e.
9.∑k=09k!1=362880986410=2.71828152…, en
10!1=2.76⋅10−7≤R9≤10!2.75=7.58⋅10−7,
dus 2.7182818≤e≤2.7182823: met bewijs is e=2.718282 tot op zes decimalen (werkelijke waarde 2.7182818…).
10.f(1−x)=n!(1−x)nxn=f(x). Op (0,1) is 0<x(1−x)≤41, dus 0<f≤4nn!1.
11.xn(1−x)n=∑j=0n(−1)j(jn)xn+j, dus f=n!1∑jcjxn+j met cj∈Z. Bijgevolg is f(k)(0)=0 voor k<n of k>2n, en voor n≤k≤2n is f(k)(0)=n!k!ck−n, een geheel getal omdat n!∣k!. De symmetrie geeft f(k)(1)=(−1)kf(k)(0)∈Z.
12.bnπ2n−2k=bn(ba)n−k=an−kbk∈Z, dus zijn G(0)=∑k(−1)kan−kbkf(2k)(0) en evenzo G(1) gehele getallen volgens vraag 11.
13. In π2G+G′′ draagt de term k van π2G de bijdrage π2n−2k+2f(2k) en de term j=k−1 van G′′ de bijdrage (−1)k−1π2n−2k+2f(2k): alles heft elkaar op behalve k=0 in de eerste som en j=n in de tweede, dat wil zeggen
G′′+π2G=bn(π2n+2f+(−1)nf(2n+2))=bnπ2n+2f=π2anf
(f heeft graad 2n, dus f(2n+2)=0; en bnπ2n=an). Dan is
dus πan∫01fsinπx=G(0)+G(1)∈Z. Op (0,1) is f>0 en sinπx>0: het linkerlid is positief, dus G(0)+G(1)≥1; en sin≤1 samen met vraag 10 begrenst het door 4nn!πan.
15. Volgens vraag 1 (met c=4a) nadert 4nn!πan tot 0: voor grote n is het <1, in tegenspraak met G(0)+G(1)≥1. Geen enkele breuk ba is dus gelijk aan π2: de stelling van Legendre. Was π rationaal, dan zou π2 dat ook zijn: dus π∈/Q — en de implicatie loopt alleen deze kant op (2 is irrationaal met een rationaal kwadraat), en daarom is π2∈/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 n!1 blijven de randgegevens geheel (nog gemakkelijker zelfs), maar wordt de grens 4nπan, die alleen naar 0 nadert wanneer a<4 — en elke kandidaat heeft a=bπ2>9. De faculteit is precies wat de meetkundige groei an voorbijstreeft: geen faculteit, geen stelling.
17. Voor a≤10 is het gehele getal G(0)+G(1) positief en hoogstens π(10/4)n/n!. Bij n=7: 2.57=610.35…, dus is de grens 5040π×610.35≈0.38<1 (bij n=6 is zij nog 1.07): de tegenspraak landt expliciet in het zevende stadium.
(de haaktermen verdwijnen: x(1−x) in 0,1, en sinπx in 0,1). Symbolisch luidt de telescopering bij n=1 (zonder enige aanname over π) π3∫01f1sinπx=−(f1′′(0)+f1′′(1)) met f1=x(1−x), f1′′=−2: rechterlid 4 — de twee berekeningen stemmen overeen.
19.ex lost y′=y op en sinπx lost y′′=−π2y op: lineaire vergelijkingen met constante coëfficiënten, zodat partiële integratie de kern naar zichzelf terugvoert en alle randgegevens binnen Z+Ze houdt (respectievelijk binnen de gehele veeltermen in π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 (etranscendent, 1873) en die van Lindemann (πtranscendent, 1882) — voorbij dit volume.
20. Veeltermdeling (of vermenigvuldig terug en controleer):
x4(1−x)4=(x6−4x5+5x4−4x2+4)(1+x2)−4.
Integratie van de getoonde identiteit gedeeld door 1+x2 geeft
21. De integrand is continu en positief op (0,1): de integraal is >0, dus π<722. Bovendien is 21≤1+x21≤1 op [0,1] en ∫01(x(1−x))4=9!(4!)2=6301 (vraag 3):
22. Modulo x2+1: x2≡−1, dus x4m=(x2)2m≡1 en (1−x)2=1−2x+x2≡−2x, zodat (1−x)4m≡(−2x)2m=4m(x2)m≡(−4)m: de rest is de constante (−4)m, en het quotiënt Qm heeft gehele coëfficiënten (deling door een monische gehele veelterm). Deling van de identiteit door 1+x2, gevolgd door integratie, geeft
Oplossen naar π: met rm=(−1)m+141−msm∈Q is ∣π−rm∣=41−mJm≤41−m⋅4−4m=41−5m. Voor m=1: s1=722, r1=722, grens 4−4=2561 — opnieuw de vragen 20–21.
23.π−722=722−π≈1.26⋅10−3, zestien keer beter dan de ijkwaarde van orde 2, 721≈2.0⋅10−2, die Dirichlet garandeert (Probleem 14.1, vraag 4); en π−113355≈2.7⋅10−7 verslaat 11321≈7.8⋅10−5 met een factor ≈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 π is op dit niveau niet beschikbaar — π mag spectaculair goed benaderd worden.
24. Lemma: stel x=qp en 0<∣an+bnx∣→0. Dan is ∣an+bnx∣=q∣qan+pbn∣, waarbij qan+pbn een niet-nul geheel getal is (niet nul omdat de absolute waarde >0 is): dus is ∣an+bnx∣≥q1 voor elke n, in tegenspraak met de convergentie naar 0. Gevallen: vraag 7 (x=e, an=αn, bn=βn); Oefening 11.9 (opnieuw x=e, met aq=−q!∑k≤qk!1, bq=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) 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 c), en daarom reikt de machine tot π2 terwijl Probleem 14.1 slechts benaderingsexponenten bereikte. (iii) Geoogst: e∈/Q, π2∈/Q (en dus π∈/Q), e=2.718282 gecertificeerd, 722−6301≤π≤722−12601, en (n2n)≥2n+14n. (iv) De grens: Hermite en Lindemann duwen dezelfde machine door tot transcendentie; en Apéry (1979) liet de gehele-getallenval draaien op ζ(3) — de machine produceert nog steeds wiskunde van de twintigste eeuw.