Mathematics · Book 4 · Bachelor Year 2

Universitaire wiskunde — Bachelor jaar 2

Universitaire wiskunde — Bachelor jaar 2 · Bachelor Year 2

6Vergelijking van functies

Asymptotische analyse — de kunst een gecompliceerde grootheid te vervangen door een eenvoudige plus een gecontroleerde fout — begon in het volume van Jaar 1 met Taylorontwikkelingen. Dit hoofdstuk maakt er een discipline op zich van: ontwikkelingen langs algemene schalen, de reeks–integraalvergelijking met haar volle asymptotische kracht, de formule van Stirling (volledig bewezen), en de systematische studie van impliciet gedefinieerde rijen. Deze technieken zijn het dagelijks brood van asymptotische analyse, en elk later hoofdstuk dat iets schat — reeksen, integralen, kansen — eet van deze tafel.

6.1 Vergelijkingsrelaties en schalen

Definitie 6.1

Nabij een punt aa (aRa \in \R of ±\pm\infty), voor functies (of rijen, met nn \to \infty): f=o(g)f = o(g), f=O(g)f = O(g), fgf \sim g zoals in het volume van Jaar 1. Een vergelijkingsschaal in aa is een familie van positieve functies, paarsgewijs vergelijkbaar, totaal geordend door o()o(\cdot) — de standaardschaal in ++\infty is

xα(lnx)β(α,βR),x^{\alpha} (\ln x)^{\beta} \qquad (\alpha, \beta \in \R),

lexicografisch geordend in (α,β)(\alpha, \beta), verfijnd wanneer nodig door exponentialen eγx\eu^{\gamma x}.

Definitie 6.2 (Asymptotische ontwikkeling)

ff laat de asymptotische ontwikkeling

f=c1φ1+c2φ2++ckφk+o(φk)(φi+1=o(φi) in de schaal)f = c_1 \varphi_1 + c_2\varphi_2 + \dots + c_k \varphi_k + o(\varphi_k) \qquad (\varphi_{i+1} = o(\varphi_i) \text{ in de schaal})

toe wanneer de opeenvolgende resten aan de weergegeven schattingen voldoen. De coëfficiënten zijn dan uniek: c1=limf/φ1c_1 = \lim f/\varphi_1, en inductief ci+1=lim(fjicjφj)/φi+1c_{i+1} = \lim\,(f - \sum_{j \leq i} c_j\varphi_j)/\varphi_{i+1}.

Voorbeeld 6.3

Taylorontwikkelingen zijn asymptotische ontwikkelingen langs de schaal (xa)k(x - a)^k in aa. Maar het begrip is strikt ruimer: in ++\infty,

1xlnx=1x11lnxx=1x+lnxx2+o(lnxx2),\frac{1}{x - \ln x} = \frac1x \cdot \frac{1}{1 - \frac{\ln x}{x}} = \frac1x + \frac{\ln x}{x^2} + o\Bigl(\frac{\ln x}{x^2}\Bigr),

een ontwikkeling langs de gemengde schaal — geen Taylorstelling is van toepassing, alleen de meetkundige ontwikkeling en de calculus van oo’s.

Voorbeeld 6.4 (De standaardschaal is echt geordend)

De lexicografische claim van Definitie 6.1 heeft één regel bewijs per geval nodig. Vergelijk xα(lnx)βx^{\alpha}(\ln x)^{\beta} en xα(lnx)βx^{\alpha'}(\ln x)^{\beta'} in ++\infty. Als α<α\alpha < \alpha': de ratio is xαα(lnx)ββ0x^{\alpha - \alpha'}(\ln x)^{\beta - \beta'} \to 0, omdat een negatieve macht van xx elke macht van lnx\ln x verplettert (stel x=etx = \eu^t: e(αα)ttββ0\eu^{(\alpha - \alpha')t}\,t^{\beta - \beta'} \to 0 door de exponentieel-verslaat-polynoom-limiet van het volume van Jaar 1). Als α=α\alpha = \alpha' en β<β\beta < \beta': de ratio is (lnx)ββ0(\ln x)^{\beta - \beta'} \to 0 direct. Dus ordenen de paren (α,β)(\alpha, \beta), lexicografisch geordend, de schaal door o()o(\cdot) — en de substitutie x=etx = \eu^t is de universele truc voor gemengde macht-log-vergelijkingen.

Voorbeeld 6.5 (Een menagerie rangschikken)

Schalen moeten geordend zijn; hier is de standaardoefening. In ++\infty, vergelijk n10n^{10}, elnnn\eu^{\sqrt{\ln n}\,\cdot\,\sqrt n}, 2n2^n en nlnnn^{\ln n} door logaritmen te nemen:

10lnn    (lnn)2    nlnn    nln2,10\ln n \;\ll\; (\ln n)^2 \;\ll\; \sqrt{n\ln n} \;\ll\; n\ln 2 ,

waar anbna_n \ll b_n betekent an=o(bn)a_n = o(b_n); de tweede entry is ln(nlnn)\ln(n^{\ln n}). Exponentialen bewaren deze strikte kloven (als lnunlnvn\ln u_n - \ln v_n \to -\infty dan un/vn0u_n/v_n \to 0), dus

n10=o(nlnn),nlnn=o(enlnn),enlnn=o(2n).n^{10} = o\bigl(n^{\ln n}\bigr), \qquad n^{\ln n} = o\bigl(\eu^{\sqrt{n\ln n}}\bigr), \qquad \eu^{\sqrt{n\ln n}} = o(2^n) .

De moraal, tweemaal: vergelijk altijd via logaritmen (verschillen van logs, niet ratio’s van logs), en concludeer nooit unvnu_n \sim v_n uit lnunlnvn\ln u_n \sim \ln v_n — het paar n10n^{10} en nlnnn^{\ln n} heeft ln\ln-ratio naar \infty, maar 2n2^n en 4n4^n hebben ln\ln-ratio precies 22 en zijn wild inequivalent.

6.2 Reeks–integraalvergelijking, asymptotisch

Stelling 6.6

Zij ff continu, positief, dalend op [1,+)\intco{1}{+\infty}.

  1. Als 1f\int_1^{\infty} f convergeert, voldoen de resten aan

    n+1f    k>nf(k)    nf.\int_{n+1}^{\infty} f \;\leq\; \sum_{k > n} f(k) \;\leq\; \int_{n}^{\infty} f .
  2. Als 1f\int_1^\infty f divergeert, voldoen de partiële sommen aan k=1nf(k)=1nf+C+o(1)\sum_{k=1}^{n} f(k) = \int_1^n f + C + o(1) voor zekere constante CC: het verschil knf(k)1nf\sum_{k \leq n} f(k) - \int_1^n f convergeert.

Bewijs. De inkleming f(k+1)kk+1ff(k)f(k+1) \leq \int_k^{k+1} f \leq f(k) (daling) was het apparaat van Jaar 1; sommeren over kn+1k \geq n+1, resp. knk \geq n, geeft (1). Voor (2), stel uk=f(k)kk+1fu_k = f(k) - \int_k^{k+1} f: door de inkleming, 0ukf(k)f(k+1)0 \leq u_k \leq f(k) - f(k+1), dus zijn de partiële sommen van uk\sum u_k begrensd door de telescoperende f(1)f(n+1)f(1)f(1) - f(n+1) \leq f(1): de reeks convergeert. Bovendien is de rij (nn+1f)n\bigl(\int_n^{n+1} f\bigr)_n niet-stijgend (ff daalt) en niet-negatief, dus convergent. Schrijvend

k=1nf(k)1nf=k=1nuk+nn+1f,\sum_{k=1}^{n} f(k) - \int_1^n f = \sum_{k=1}^{n} u_k + \int_n^{n+1} f ,

convergeert de rechterkant als nn \to \infty: het verschil convergeert naar een constante CC, wat uitspraak (2) is.

Voorbeeld 6.7 (De harmonische ontwikkeling)

Voor f(t)=1tf(t) = \frac1t: Hn=lnn+γ+o(1)H_n = \ln n + \gamma + o(1), herwint de constante van Euler (volume van Jaar 1) met een schoner bewijs. Eén orde verder duwen (Oefening 6.3):

Hn=lnn+γ+12n+o(1n).H_n = \ln n + \gamma + \frac{1}{2n} + o\Bigl(\frac1n\Bigr).

Getallen maken de winst zichtbaar bij n=10n = 10: H10=2.928968H_{10} = 2.928968\dots en ln10=2.302585\ln 10 = 2.302585\dots, dus de ruwe schatting van γ\gamma is H10ln10=0.626383H_{10} - \ln 10 = 0.626383, 0.0490.049 ernaast; aftrekken van de correctie 120\frac1{20} geeft 0.5763830.576383, naast γ=0.577216\gamma = 0.577216 met slechts 8.31048.3\cdot10^{-4} — wat zelf de volgende term 112100\frac{1}{12\cdot100} van de ontwikkeling is, zoals het weekendvraagstuk bewijst (vraag 8).

Voorbeeld 6.8 (Een ruwe ln(n!)\ln(n!) zonder Stirling)

Het inklemingsapparaat alleen al lokaliseert ln(n!)\ln(n!). Aangezien ln\ln stijgt,

k1klnt ⁣dt    lnk    kk+1lnt ⁣dt,\int_{k-1}^{k}\ln t\,\dd t \;\leq\; \ln k \;\leq\; \int_{k}^{k+1}\ln t\,\dd t ,

en sommeren over k=2,,nk = 2, \dots, n (met 1nln=nlnnn+1\int_1^n\ln = n\ln n - n + 1):

nlnnn+1    ln(n!)    (n+1)ln(n+1)n.n\ln n - n + 1 \;\leq\; \ln(n!) \;\leq\; (n+1)\ln(n+1) - n .

Beide hekken zijn nlnnn+O(lnn)n\ln n - n + O(\ln n): dus ln(n!)=nlnnn+O(lnn)\ln(n!) = n\ln n - n + O(\ln n), en in het bijzonder ln(n!)nlnn\ln(n!) \sim n\ln n. Wat Stirling toevoegt zijn de volgende twee treden — de 12lnn\frac12\ln n en de constante ln2π\ln\sqrt{2\pi} — die de fijnere telescopatie van Stelling 6.13 kosten. Weten welke precisie elk instrument koopt is de helft van asymptotisch vakmanschap.

Voorbeeld 6.9 (Cancellatie eist ontwikkelingen)

Bereken de limiet van n2+nn\sqrt{n^2 + n} - n. Beide termen zijn n\sim n, en “nn\sim n - n” is zinloos: equivalenten kunnen niet worden afgetrokken. Ontwikkel in plaats daarvan:

n2+nn=n(1+1n1)=n(12n18n2+O(1n3))=1218n+O(1n2):\sqrt{n^2 + n} - n = n\Bigl(\sqrt{1 + \tfrac1n} - 1\Bigr) = n\Bigl(\frac{1}{2n} - \frac{1}{8n^2} + O\Bigl(\frac1{n^3}\Bigr)\Bigr) = \frac12 - \frac{1}{8n} + O\Bigl(\frac1{n^2}\Bigr) :

limiet 12\frac12, met de naderingssnelheid 18n\frac1{8n} als bonus. Het mechanisme verdient een naam: een verschil van twee grote equivalente grootheden leeft volledig in hun volgende termen, dus moet men ontwikkelen tot de eerste orde waarop de twee kanten verschillen — en de rest meedragen om te certificeren dat niets anders op die orde overleeft.

Voorbeeld 6.10 (Een divergente vergelijking, uitgewerkt)

Voor f(t)=1tlntf(t) = \frac{1}{t\ln t} op [2,+)\intco{2}{+\infty} (continu, positief, dalend): 2xf=lnlnxlnln2\int_2^x f = \ln\ln x - \ln\ln 2 \to \infty, dus door Stelling 6.6 (2),

k=2n1klnk=lnlnn+C+o(1)\sum_{k=2}^{n}\frac{1}{k\ln k} = \ln\ln n + C + o(1)

voor zekere constante CC. Twee lessen. Eerst, de divergentie is echt maar gletsjerachtig: de partiële som overschrijdt 44 eerst rond nee4Cn \approx \eu^{\eu^{4 - C}}, astronomisch groot. Ten tweede werd de vorm lnlnn\ln\ln n geleverd door een antiderivaat, niet geraden: voor monotone termen is de integraal het canonieke sommeerapparaat, en de constante CC — zoals Eulers γ\gamma — is het geheugen van de begintermen.

6.3 De formule van Stirling

Lemma 6.11 (Wallisintegralen, herzien)

Zij Wn=0π/2sinnt ⁣dtW_n = \int_0^{\pi/2} \sin^n t\,\dd t. Dan nWnWn1=π2nW_nW_{n-1} = \frac\pi2 voor n1n \geq 1, (Wn)(W_n) daalt, en Wnπ2nW_n \sim \sqrt{\dfrac{\pi}{2n}}.

Bewijs. Partiële integratie geeft nWn=(n1)Wn2nW_n = (n-1)W_{n-2} (n2n \geq 2), dus is nWnWn1nW_nW_{n-1} constant in nn, gelijk aan 1W1W0=π21 \cdot W_1 W_0 = \frac\pi2. Daling: sinn+1sinn\sin^{n+1} \leq \sin^n op [0,π2]\intcc{0}{\frac\pi2}. De squeeze, in detail: monotoniciteit geeft Wn+1WnWn1W_{n+1} \leq W_n \leq W_{n-1}, en delen door Wn1>0W_{n-1} > 0,

nn+1=Wn+1Wn1WnWn11,\frac{n}{n+1} = \frac{W_{n+1}}{W_{n-1}} \leq \frac{W_n}{W_{n-1}} \leq 1 ,

de linkeridentiteit uit de recurrentie op index n+1n + 1. Beide grenzen gaan naar 11: WnWn1W_n \sim W_{n-1}, vandaar

nWn2nWnWn1=π2Wnπ2n.nW_n^2 \sim nW_nW_{n-1} = \frac\pi2 \qquad\Longrightarrow\qquad W_n \sim \sqrt{\frac{\pi}{2n}} .

Voorbeeld 6.12 (De eerste Wallisintegralen)

Uit W0=π2W_0 = \frac\pi2, W1=1W_1 = 1 en de recurrentie nWn=(n1)Wn2nW_n = (n-1)W_{n-2}:

W2=π4,W3=23,W4=3π16,W5=815,W6=5π32.W_2 = \frac\pi4, \qquad W_3 = \frac23, \qquad W_4 = \frac{3\pi}{16}, \qquad W_5 = \frac{8}{15}, \qquad W_6 = \frac{5\pi}{32}.

Even indices dragen een π\pi, oneven zijn rationaal — de twee verweven producten van de gesloten vormen. Numeriek W60.4909W_6 \approx 0.4909 tegen de asymptotiek π/120.5116\sqrt{\pi/12} \approx 0.5116: bij n=6n = 6 is het equivalent al binnen 5%5\%, en de productidentiteit is exact bij elke nn: 6W6W5=65π32815=π26\,W_6W_5 = 6\cdot\frac{5\pi}{32}\cdot\frac8{15} = \frac\pi2. Kleine tabellen als deze zijn de goedkoopste manier om een algebrafout te vangen voordat zij een asymptotisch argument infecteert.

Stelling 6.13 (Stirling)

n!    2πn(ne) ⁣n.n! \;\sim\; \sqrt{2\pi n}\, \Bigl(\frac{n}{\eu}\Bigr)^{\!n} .

Bewijs. Stap 1: n!Cn(n/e)nn! \sim C \sqrt n\, (n/\eu)^n voor zekere constante C>0C > 0. Stel

dn=ln(n!)(n+12)lnn+n.d_n = \ln(n!) - \Bigl(n + \frac12\Bigr)\ln n + n .

Dan

dndn+1=(n+12)lnn+1n1=(n+12)(1n12n2+13n3+o(n3))1=112n2+o(1n2),d_n - d_{n+1} = \Bigl(n + \frac12\Bigr) \ln\frac{n+1}{n} - 1 = \Bigl(n + \frac12\Bigr)\Bigl(\frac1n - \frac{1}{2n^2} + \frac{1}{3n^3} + o\bigl(n^{-3}\bigr)\Bigr) - 1 = \frac{1}{12n^2} + o\Bigl(\frac{1}{n^2}\Bigr),

door de Taylorontwikkeling van ln(1+1n)\ln(1 + \frac1n). De reeks (dndn+1)\sum (d_n - d_{n+1}) convergeert dus absoluut (vergelijking met n2\sum n^{-2}), dus convergeert (dn)(d_n), zeg naar dd; exponentiëren, n!Cn(n/e)nn! \sim C\sqrt n\,(n/\eu)^n met C=edC = \eu^{d}.

Stap 2: C=2πC = \sqrt{2\pi} via Wallis. De gesloten vorm W2p=(2p)!4p(p!)2π2W_{2p} = \frac{(2p)!}{4^p (p!)^2}\cdot\frac\pi2 (uit de recurrentie, berekening van Jaar 1 herdaan in de setting van Lemma 6.11) combineert met Stap 1:

W2pC2p(2p/e)2p4p(Cp(p/e)p)2π2=2pCpπ2=πC12p.W_{2p} \sim \frac{C\sqrt{2p}\,(2p/\eu)^{2p}} {4^p\,\bigl(C\sqrt p\,(p/\eu)^p\bigr)^2}\cdot\frac{\pi}{2} = \frac{\sqrt{2p}}{C\,p}\cdot\frac{\pi}{2} = \frac{\pi}{C}\cdot\frac{1}{\sqrt{2p}} .

Vergelijken met W2pπ4pW_{2p} \sim \sqrt{\frac{\pi}{4p}} (Lemma 6.11): πC2p=π4p(1+o(1))\frac{\pi}{C\sqrt{2p}} = \sqrt{\frac{\pi}{4p}}\,(1 + o(1)) forceert C=π4p2pπ=2πC = \pi \sqrt{\frac{4p}{2p\,\pi}} = \sqrt{2\pi}.

Voorbeeld 6.14 (Centraal binomiaalcoëfficiënt)

(2nn)=(2n)!(n!)24πn(2n/e)2n2πn(n/e)2n=4nπn:\binom{2n}{n} = \frac{(2n)!}{(n!)^2} \sim \frac{\sqrt{4\pi n}\,(2n/\eu)^{2n}}{2\pi n\,(n/\eu)^{2n}} = \frac{4^n}{\sqrt{\pi n}} :

de kans dat een symmetrische random walk op tijd 2n2n naar 00 terugkeert is 1πn\sim \frac{1}{\sqrt{\pi n}} — een aankondiging van Hoofdstuk 22.

Opmerking 6.15 (Perspectieven binnen dit volume)

Elk kwantitatief hoofdstuk vooruit spreekt de taal van dit hoofdstuk. Hoofdstuk 7 classificeert reeksen door termen te vergelijken met de schaal nα(lnn)βn^{-\alpha}(\ln n)^{-\beta} — haar weekendvraagstuk brengt die grens volledig in kaart. Hoofdstuk 9 doet hetzelfde voor oneigenlijke integralen, met de identieke schaal in de continue variabele. Hoofdstuk 11 berekent convergentiestralen uit lim supan1/n\limsup\abs{a_n}^{1/n}, een equivalent-van-nn-de-wortels-oefening waar Stirling de standaardsleutel is (n!nne\sqrt[n]{n!} \sim \frac n\eu, Oefening 6.4). En de kanshoofdstukken cashen Stirling direct: de lokale schattingen van Hoofdstuk 22 voor binomiaalcoëfficiënten zijn Voorbeeld 6.14 en Voorbeeld 6.21 letterlijk. Asymptotiek is hier geen hoofdstuk; het is het accent van het volume.

Methode 6.16 (De bootstrapchecklist)

Voor je een gebootstrapte ontwikkeling vertrouwt, audit vier punten. (1) Eerst bestaan: de wortel of rij moet vastgepind worden (monotoniciteit, tussenwaarden) vóór enige ontwikkeling — symbolen zonder referenten ontwikkelen prachtig en betekenen niets. (2) Eén orde per doorgang: elke substitutie mag alleen vertrouwd worden tot de orde van de ingevoerde schatting; twee nieuwe termen uit één doorgang extraheren is de klassieke bron van foute coëfficiënten. (3) Resten rijden mee: draag de o()o(\cdot) door elke algebraïsche stap en laat absorptie (kleinere termen opgeslokt door grotere resten) aan het eind expliciet gebeuren. (4) Numerieke audit: evalueer bij één eerlijke waarde van nn; een coëfficiëntfout overleeft algebraïsche herleiding verrassend vaak, en bijna nooit rekenkunde.

Opmerking 6.17 (Veelgemaakte valkuilen)

(i) Equivalenten optellen slecht: uit unn+lnnu_n \sim n + \ln n en vnnv_n \sim -n mag men niet concluderen un+vnlnnu_n + v_n \sim \ln n; cancellaties eisen ontwikkelingen met expliciete resten, nooit kale equivalenten. (ii) Exponentieer nooit een equivalentie: n+1nn + 1 \sim n maar en+1≁en\eu^{n+1} \not\sim \eu^n; de veilige richting is logaritmen nemen van equivalenten die naar ++\infty gaan (weekendvraagstuk van dit hoofdstuk, vraag 24). (iii) Een asymptotische ontwikkeling is gehecht aan een schaal: schrijven f=1x+o(1x2)f = \frac1x + o\bigl(\frac1{x^2}\bigr) claimt meer dan f=1x+o(1x)f = \frac1x + o\bigl(\frac1x\bigr), en de twee mengen invalideert latere algebra. (iv) In bootstraps, substitueer de hele huidige ontwikkeling, rest inbegrepen — een o()o(\cdot) mid-pass droppen produceert aannemelijke maar foute coëfficiënten. (v) De reeks–integraalvergelijking eist monotoniciteit: voor oscillerende termen faalt zij botweg (vergelijk sinkk\sum\frac{\sin k}k, Hoofdstuk 7).

Voorbeeld 6.18 (Stirling in getallen)

Bij n=10n = 10: de formule geeft 20π(10/e)103598696\sqrt{20\pi}\,(10/\eu)^{10} \approx 3\,598\,696 tegen 10!=362880010! = 3\,628\,800: relatieve fout 8.31038.3\cdot10^{-3}, opmerkelijk voor een “asymptotische” uitspraak bij n=10n = 10. De fout heeft structuur — de exacte verfijning n!=2πn(n/e)n(1+112n+O(n2))n! = \sqrt{2\pi n}\,(n/\eu)^n\bigl(1 + \frac1{12n} + O(n^{-2})\bigr) — waarvan de eerste correctie 11208.3103\frac1{120} \approx 8.3\cdot10^{-3} de waargenomen kloof bijna exact verklaart. De Euler–Maclaurin-machinerie van het weekendvraagstuk is precies de systematische bron van zulke correctietermen.

Opmerking 6.19 (Waar dit hoofdstuk gebruikt wordt)

Asymptotische vergelijking is de grammatica van alles kwantitatiefs stroomafwaarts: de convergentietests en het Bertrandpanorama van Hoofdstuk 7, de integreerbaarheidscriteria van Hoofdstuk 9, de convergentiestraalberekeningen van Hoofdstuk 11, en de limietstellingen van Hoofdstuk 22 (waar Stirling de de Moivre–Laplace-schattingen draait). Het volume van Jaar 3 industrialiseert het ene idee dat we hier met de hand bewijzen — extraheer de hoofdterm, begrens de rest — tot de Laplace-methode en gedomineerde convergentie.

Voorbeeld 6.20 (Een integraal vergeleken met zichzelf: 2x ⁣dtlnt\int_2^x \frac{\dd t}{\ln t})

De vergelijkingsgereedschapskist loopt ook op integralen. Zij F(x)=2x ⁣dtlntF(x) = \int_2^x\frac{\dd t}{\ln t} (het integrandum is continu op [2,)\intco2\infty). Partiële integratie:

F(x)=[tlnt]2x+2x ⁣dt(lnt)2=xlnx+O(2x ⁣dt(lnt)2)+O(1),F(x) = \Bigl[\frac{t}{\ln t}\Bigr]_2^x + \int_2^x\frac{\dd t}{(\ln t)^2} = \frac{x}{\ln x} + O\Bigl(\int_2^x\frac{\dd t}{(\ln t)^2}\Bigr) + O(1),

en de restintegraal is o(xlnx)o\bigl(\frac{x}{\ln x}\bigr): splits haar in x\sqrt x, begrensend door

2x ⁣dt(lnt)2xenxx ⁣dt(lnt)2x(lnx)2=4x(lnx)2.\int_2^{\sqrt x}\frac{\dd t}{(\ln t)^2} \leq \sqrt x \qquad\text{en}\qquad \int_{\sqrt x}^{x}\frac{\dd t}{(\ln t)^2} \leq \frac{x}{(\ln\sqrt x)^2} = \frac{4x}{(\ln x)^2} .

Dus F(x)xlnxF(x) \sim \frac{x}{\ln x}. Lezers die de priemgetalstelling in het weekendvraagstuk van dit hoofdstuk ontmoetten herkennen FF: zij is de logaritmische integraal, de betere schatter van π(x)\pi(x), en de berekening toont dat zij tot op eerste orde met xlnx\frac{x}{\ln x} overeenkomt.

Voorbeeld 6.21 (Stirling op een scheve binomiaal)

Dezelfde drie-faculteitsroutine als voor Voorbeeld 6.14 geeft, voor (3nn)=(3n)!n!(2n)!\binom{3n}{n} = \frac{(3n)!}{n!\,(2n)!}:

(3nn)6πn(3n/e)3n2πn(n/e)n4πn(2n/e)2n=34πn(274) ⁣n.\binom{3n}{n} \sim \frac{\sqrt{6\pi n}\,(3n/\eu)^{3n}} {\sqrt{2\pi n}\,(n/\eu)^{n}\cdot\sqrt{4\pi n}\,(2n/\eu)^{2n}} = \sqrt{\frac{3}{4\pi n}}\, \Bigl(\frac{27}{4}\Bigr)^{\!n} .

De exponentiële rate 274=3322\frac{27}4 = \frac{3^3}{2^2} is e3nH(1/3)\eu^{3n\,H(1/3)} in de entropienotatie van informatietheorie: scheve binomialen groeien strikt trager dan de centrale 4n4^n per twee stappen — hier (27/4)1/31.89<2(27/4)^{1/3} \approx 1.89 < 2 per stap. Elke binomiale asymptotiek in combinatoriek en kansrekening (Hoofdstuk 22) is deze ene berekening met andere gewichten.

6.4 Impliciet gedefinieerde rijen

Methode 6.22

Om de asymptotiek van oplossingen xnx_n van een vergelijking F(x,n)=0F(x, n) = 0 te vinden:

  1. Lokaliseer: bewijs bestaan en uniciteit van xnx_n in een bepaald interval (monotoniciteit, tussenwaardestelling), en vind haar ruwe gedrag (limiet, groeiorde).
  2. Bootstrap: substitueer de ruwe vorm xn=(hoofdterm)(1+εn)x_n = (\text{hoofdterm})(1 + \varepsilon_n) in de vergelijking en los de volgende orde van εn\varepsilon_n op; herhaal, elke doorgang verfijnt één orde.

Voorbeeld 6.23

Voor n1n \geq 1 heeft de vergelijking tanx=x\tan x = x precies één oplossing xnx_n in (nππ2,nπ+π2)\intoo{n\pi - \frac\pi2}{n\pi + \frac\pi2} (de functie tanxx\tan x - x stijgt daar van -\infty naar ++\infty, haar afgeleide is tan2x0\tan^2 x \geq 0). Ruw: xn=nπ+π2ynx_n = n\pi + \frac\pi2 - y_n met yn(0,π)y_n \in \intoo{0}{\pi}; aangezien xnx_n \to \infty en tanxn=xn+\tan x_n = x_n \to +\infty, nadert xnx_n de asymptoot van links: yn0y_n \to 0. Bootstrap: tanxn=cotyn=1tanyn1yn\tan x_n = \cot y_n = \frac{1}{\tan y_n} \sim \frac{1}{y_n}, en de vergelijking cotyn=xnnπ\cot y_n = x_n \sim n\pi geeft yn1nπy_n \sim \frac{1}{n\pi}. Dus

xn=nπ+π21nπ+o(1n),x_n = n\pi + \frac\pi2 - \frac{1}{n\pi} + o\Bigl(\frac1n\Bigr),

en het proces gaat door tot elke orde (Oefening 6.6).

Voorbeeld 6.24 (Een tweede run van de methode)

Los x+lnx=nx + \ln x = n asymptotisch op. Lokaliseer: xx+lnxx \mapsto x + \ln x stijgt van -\infty naar ++\infty op (0,+)\intoo{0}{+\infty}: een unieke wortel xnx_n, en xnx_n \to \infty. Ruw: lnxn=o(xn)\ln x_n = o(x_n) geeft xnnx_n \sim n. Bootstrap: uit xn=nlnxnx_n = n - \ln x_n en lnxn=lnn+o(1)\ln x_n = \ln n + o(1) (logaritmen van equivalenten, beide kanten \to \infty):

xn=nlnn+o(1);x_n = n - \ln n + o(1) ;

nog een doorgang, met lnxn=ln(nlnn+o(1))=lnnlnnn+o(lnnn)\ln x_n = \ln\bigl(n - \ln n + o(1)\bigr) = \ln n - \frac{\ln n}{n} + o\bigl(\frac{\ln n}n\bigr):

xn=nlnn+lnnn+o(lnnn).x_n = n - \ln n + \frac{\ln n}{n} + o\Bigl(\frac{\ln n}{n}\Bigr).

(Controle bij n=100n = 100: de wortel is x95.4415x \approx 95.4415; de drietermformule geeft 1004.6052+0.0461=95.4409100 - 4.6052 + 0.0461 = 95.4409, de tweeterm 95.394895.3948 — elke doorgang wint de voorspelde orde.) Zelfde lus, derde landschap: de methode van Methode 6.22 geeft niet om hoe de vergelijking eruitziet, alleen dat elke doorgang de dominante onbekende isoleert.

6.5 Oefeningen

Oefening 6.1

Ontwikkel in ++\infty, twee termen voorbij de leidende:

x2+x+1,ln(x2+x)2lnx,x+sinxxlnx.\sqrt{x^2 + x + 1} , \qquad \ln(x^2 + x) - 2\ln x, \qquad \frac{x + \sin x}{x - \ln x} .
Oplossing

Oplossing van Oefening 6.1.

x2+x+1=x1+1x+1x2=x+12+381x+o(1x)\sqrt{x^2 + x + 1} = x\sqrt{1 + \tfrac1x + \tfrac{1}{x^2}} = x + \frac12 + \frac38\cdot\frac1x + o\bigl(\frac1x\bigr) (binomiale ontwikkeling: 12u18u2\frac12 u - \frac18 u^2 met u=1x+1x2u = \frac1x + \frac{1}{x^2} geeft 12x+12x218x2=12x+38x2\frac{1}{2x} + \frac{1}{2x^2} - \frac{1}{8x^2} = \frac{1}{2x} + \frac{3}{8x^2}, dan vermenigvuldigen met xx).

ln(x2+x)2lnx=ln(1+1x)=1x12x2+o(1x2)\ln(x^2 + x) - 2\ln x = \ln\bigl(1 + \tfrac1x\bigr) = \frac1x - \frac{1}{2x^2} + o\bigl(\frac{1}{x^2}\bigr).

Derde functie: ontwikkel elke factor,

x+sinxxlnx=(1+sinxx)(1+lnxx+(lnx)2x2+O((lnx)3x3)).\frac{x + \sin x}{x - \ln x} = \Bigl(1 + \frac{\sin x}{x}\Bigr) \Bigl(1 + \frac{\ln x}{x} + \frac{(\ln x)^2}{x^2} + O\Bigl(\frac{(\ln x)^3}{x^3}\Bigr)\Bigr).

Orden de bijdragen op de schaal in ++\infty: lnxx1xsinxx(lnx)2x2\frac{\ln x}{x} \gg \frac{1}{x} \geq \bigl|\frac{\sin x}{x}\bigr| \gg \frac{(\ln x)^2}{x^2}. De twee termen na de leidende 11 zijn daarom lnxx\frac{\ln x}{x}, dan de begrensd-oscillerende term sinxx\frac{\sin x}{x}:

x+sinxxlnx=1+lnxx+sinxx+O((lnx)2x2).\frac{x + \sin x}{x - \ln x} = 1 + \frac{\ln x}{x} + \frac{\sin x}{x} + O\Bigl(\frac{(\ln x)^2}{x^2}\Bigr).

Oefening 6.2

Geef de aard (convergentie/divergentie) en, wanneer divergent, de leidende asymptotiek van knkα\sum_{k \leq n} k^\alpha voor α>1\alpha > -1, α=1\alpha = -1, α<1\alpha < -1, via Stelling 6.6.

Oplossing

Oplossing van Oefening 6.2.

f(t)=tαf(t) = t^\alpha (t1t \geq 1).

α>1\alpha > -1: divergentie, en door Stelling 6.6 (2), knkα=nα+1α+1+C+o(1)\sum_{k\leq n} k^\alpha = \frac{n^{\alpha+1}}{\alpha+1} + C + o(1) als α<0\alpha < 0 (waar ff daalt); voor α0\alpha \geq 0 (ff stijgend) geeft dezelfde inkleming met omgekeerde ongelijkheden knkαnα+1α+1\sum_{k \leq n} k^\alpha \sim \frac{n^{\alpha + 1}}{\alpha + 1}.

α=1\alpha = -1: Hn=lnn+γ+o(1)H_n = \ln n + \gamma + o(1) (Voorbeeld 6.7).

α<1\alpha < -1: convergentie, met rest k>nkαnα+1(α+1)\sum_{k > n} k^\alpha \sim \frac{n^{\alpha+1}}{-(\alpha+1)} door de inkleming (1) (beide integraalgrenzen zijn equivalent met die waarde).

Oefening 6.3 ★★

Bewijs Hn=lnn+γ+12n+o(1n)H_n = \ln n + \gamma + \frac{1}{2n} + o\bigl(\frac1n\bigr). (Bestudeer vn=Hnlnnγv_n = H_n - \ln n - \gamma: toon vnvn+1=12n2+O(n3)v_n - v_{n+1} = \frac{1}{2n^2} + O(n^{-3}) en sommeer de staart, vergelijkend met kn12k212n\sum_{k \geq n} \frac{1}{2k^2} \sim \frac{1}{2n}Stelling 6.6 (1).)

Oplossing

Oplossing van Oefening 6.3.

Zij vn=Hnlnnγ0v_n = H_n - \ln n - \gamma \to 0. Dan

vnvn+1=lnn+1n1n+1=(1n12n2)(1n1n2)+O(1n3)=12n2+O(1n3),v_n - v_{n+1} = \ln\frac{n+1}{n} - \frac{1}{n+1} = \Bigl(\frac1n - \frac{1}{2n^2}\Bigr) - \Bigl(\frac1n - \frac{1}{n^2}\Bigr) + O\Bigl(\frac{1}{n^3}\Bigr) = \frac{1}{2n^2} + O\Bigl(\frac{1}{n^3}\Bigr),

met 1n+1=1n1n2+O(n3)\frac{1}{n+1} = \frac1n - \frac{1}{n^2} + O(n^{-3}). Aangezien vn0v_n \to 0, telescoperen van de staart:

vn=kn(vkvk+1)=kn(12k2+O(k3))=12n+O(1n2),v_n = \sum_{k \geq n} (v_k - v_{k+1}) = \sum_{k\geq n} \Bigl(\frac{1}{2k^2} + O(k^{-3})\Bigr) = \frac{1}{2n} + O\Bigl(\frac{1}{n^2}\Bigr),

door Stelling 6.6 (1) toegepast op t2t^{-2} (rest 1n\sim \frac1n, gehalveerd) en op t3t^{-3}. Dus Hn=lnn+γ+12n+o(1n)H_n = \ln n + \gamma + \frac{1}{2n} + o(\frac1n).

Oefening 6.4 ★★

Vind met Stirling equivalenten van: (3n)!(n!)3\dfrac{(3n)!}{(n!)^3};   n!nn\;\dfrac{n!}{n^n};   n!n\;\sqrt[n]{n!} (als ne(1+o(1))\frac n\eu(1 + o(1)), precies gemaakt tot twee termen).

Oplossing

Oplossing van Oefening 6.4.

Stirling driemaal:

(3n)!(n!)36πn(3n/e)3n(2πn)3/2(n/e)3n=6  27n2πn12πn2πn  =327n2πn.\frac{(3n)!}{(n!)^3} \sim \frac{\sqrt{6\pi n}\,(3n/\eu)^{3n}} {(2\pi n)^{3/2}\,(n/\eu)^{3n}} = \frac{\sqrt{6}\; 27^{\,n}}{2\pi n} \cdot \frac{1}{\sqrt{2\pi n}}\cdot\sqrt{2\pi n}\; = \frac{\sqrt3\,27^n}{2\pi n} .

(Zorgvuldig: 6πn(2πn)3/2=6(2πn)2πnπn=32πn\frac{\sqrt{6\pi n}}{(2\pi n)^{3/2}} = \frac{\sqrt6}{(2\pi n)\sqrt{2\pi n}}\sqrt{\pi n} = \frac{\sqrt3}{2\pi n}.)

n!nn2πnen\dfrac{n!}{n^n} \sim \sqrt{2\pi n}\,\eu^{-n}.

n!n=exp(lnn!n)\sqrt[n]{n!} = \exp\bigl(\frac{\ln n!}{n}\bigr) met lnn!=nlnnn+12ln(2πn)+o(1)\ln n! = n\ln n - n + \frac12\ln(2\pi n) + o(1):

n!n=exp(lnn1+ln(2πn)2n+o(lnnn))=ne(1+ln(2πn)2n+o(lnnn)).\sqrt[n]{n!} = \exp\Bigl(\ln n - 1 + \frac{\ln(2\pi n)}{2n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr) = \frac{n}{\eu}\Bigl(1 + \frac{\ln(2\pi n)}{2n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr).

Oefening 6.5 ★★

Bewijs voor n2n \geq 2 dat xn+x=1x^n + x = 1 een unieke oplossing xn(0,1)x_n \in \intoo{0}{1} heeft, dat xn1x_n \to 1, en vestig

xn=1lnnn+o(lnnn).x_n = 1 - \frac{\ln n}{n} + o\Bigl(\frac{\ln n}{n}\Bigr).

(Uit xnn=1xnx_n^n = 1 - x_n: neem logaritmen en bootstrap met xn=1εnx_n = 1 - \varepsilon_n.)

Oplossing

Oplossing van Oefening 6.5.

g(x)=xn+x1g(x) = x^n + x - 1 stijgt strikt op [0,1]\intcc{0}{1} van 1-1 naar 11: unieke wortel xnx_n. Aangezien xnn=1xn(0,1)x_n^n = 1 - x_n \in \intoo{0}{1}: als xnc<1x_n \leq c < 1 langs een deelrij, dan xnncn0x_n^n \leq c^n \to 0, dus 1xn01 - x_n \to 0: tegenspraak met xncx_n \leq c. Dus xn1x_n \to 1.

Schrijf xn=1εnx_n = 1 - \varepsilon_n, εn0+\varepsilon_n \to 0^+. De vergelijking luidt (1εn)n=εn(1 - \varepsilon_n)^n = \varepsilon_n, d.w.z.

nln(1εn)=lnεnnεn(1+o(1))=lnεn.n\ln(1 - \varepsilon_n) = \ln \varepsilon_n \quad\Longrightarrow\quad -n\varepsilon_n\bigl(1 + o(1)\bigr) = \ln\varepsilon_n .

Dus nεn=lnεn(1+o(1))+n\varepsilon_n = -\ln\varepsilon_n\,(1 + o(1)) \to +\infty, en opnieuw logaritmen nemen: lnn+lnεn=ln(lnεn)+o(1)\ln n + \ln\varepsilon_n = \ln(-\ln\varepsilon_n) + o(1). Aangezien ln(lnεn)=o(ln(1/εn))\ln(-\ln \varepsilon_n) = o(\ln(1/\varepsilon_n)), geeft dit lnεnlnn\ln\varepsilon_n \sim -\ln n, vandaar εn=lnεnn(1+o(1))lnnn\varepsilon_n = \frac{-\ln\varepsilon_n}{n}(1 + o(1)) \sim \frac{\ln n}{n}:

xn=1lnnn+o(lnnn).x_n = 1 - \frac{\ln n}{n} + o\Bigl(\frac{\ln n}{n}\Bigr) .

Oefening 6.6 ★★

Duw Voorbeeld 6.23 één orde verder:

xn=nπ+π21nπ+12n2π+o(1n2).x_n = n\pi + \frac\pi2 - \frac{1}{n\pi} + \frac{1}{2n^2\pi} + o\Bigl(\frac{1}{n^2}\Bigr).

(Schrijf cotyn=xn\cot y_n = x_n exact, ontwikkel coty=1yy3+o(y)\cot y = \frac1y - \frac y3 + o(y) en xn=nπ(1+12n)x_n = n\pi(1 + \frac{1}{2n} - \dots), en identificeer.)

Oplossing

Oplossing van Oefening 6.6.

Exacte relatie: cotyn=xn=nπ+π2yn\cot y_n = x_n = n\pi + \frac\pi2 - y_n, met yn1nπy_n \sim \frac{1}{n\pi} (Voorbeeld 6.23). Ontwikkel coty=1yy3+O(y3)\cot y = \frac1y - \frac y3 + O(y^3):

1ynyn3+O(yn3)=nπ+π2yn1yn=nπ+π2+O(1n),\frac{1}{y_n} - \frac{y_n}{3} + O(y_n^3) = n\pi + \frac\pi2 - y_n \quad\Longrightarrow\quad \frac{1}{y_n} = n\pi + \frac\pi2 + O\Bigl(\frac1n\Bigr),

(de termen yn-y_n en yn3-\frac{y_n}{3} zijn O(1n)O(\frac1n)). Inverteer:

yn=1nπ11+12n+O(n2)=1nπ(112n+O(1n2))=1nπ12n2π+O(1n3).y_n = \frac{1}{n\pi}\cdot\frac{1}{1 + \frac{1}{2n} + O(n^{-2})} = \frac{1}{n\pi}\Bigl(1 - \frac{1}{2n} + O\Bigl(\frac{1}{n^2}\Bigr)\Bigr) = \frac{1}{n\pi} - \frac{1}{2n^2\pi} + O\Bigl(\frac{1}{n^3}\Bigr).

Dus

xn=nπ+π2yn=nπ+π21nπ+12n2π+o(1n2).x_n = n\pi + \frac{\pi}{2} - y_n = n\pi + \frac\pi2 - \frac{1}{n\pi} + \frac{1}{2n^2\pi} + o\Bigl(\frac{1}{n^2}\Bigr).

Oefening 6.7 ★★

Bepaal limn1n!k=0nk!\lim_{n\to\infty} \dfrac{1}{n!}\sum_{k=0}^{n} k! (begrens de som van alle termen behalve de laatste twee), en leid de asymptotische ontwikkeling knk!=n!(1+1n+O(n2))\sum_{k \leq n} k! = n!\bigl(1 + \frac1n + O(n^{-2})\bigr) af.

Oplossing

Oplossing van Oefening 6.7.

Splits de twee grootste termen af:

k=0nk!=n!+(n1)!+kn2k!,kn2k!(n1)(n2)!=(n1)!.\sum_{k=0}^{n} k! = n! + (n-1)! + \sum_{k \leq n-2} k! , \qquad \sum_{k\leq n-2} k! \leq (n-1)\,(n-2)! = (n-1)! .

Dus 11n!k!1+2n1 \leq \frac{1}{n!}\sum k! \leq 1 + \frac{2}{n}: de limiet is 11. Verfijnen: (n1)!n!=1n\frac{(n-1)!}{n!} = \frac1n en de ruwe grens kn2k!(n1)!\sum_{k \leq n-2}k! \leq (n-1)! kan op dezelfde manier worden aangescherpt: kn2k!=(n2)!(1+O(1n))=O(n!n2)\sum_{k\leq n-2} k! = (n-2)!\,(1 + O(\frac1n)) = O\bigl(\frac{n!}{n^2}\bigr). Dus

k=0nk!=n!(1+1n+O(1n2)).\sum_{k=0}^{n} k! = n!\Bigl(1 + \frac1n + O\Bigl(\frac{1}{n^2}\Bigr)\Bigr).

Oefening 6.8 ★★★

Zij u0>0u_0 > 0 en un+1=un+1unu_{n+1} = u_n + \dfrac{1}{u_n}. Bewijs dat unu_n \to \infty, dan dat un2nu_n \sim \sqrt{2n} (bestudeer un2u_n^2: haar incrementen zijn 2+un22 + u_n^{-2}; sommeer), en verfijn:

un=2n(1+lnn8n+o(lnnn)).u_n = \sqrt{2n}\Bigl(1 + \frac{\ln n}{8n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr).

(Uit un2=2n+k<nuk2+u02u_n^2 = 2n + \sum_{k<n} u_k^{-2} + u_0^2 en uk22ku_k^2 \sim 2k: de som is 12lnn\sim \frac12\ln n door Stelling 6.6.)

Oplossing

Oplossing van Oefening 6.8.

(un)(u_n) stijgt; als begrensd zou zij naar \ell convergeren met =+1\ell = \ell + \frac1\ell: absurd. Dus unu_n \to \infty.

Kwadraten: un+12=un2+2+un2u_{n+1}^2 = u_n^2 + 2 + u_n^{-2}, dus

un2=u02+2n+k=0n11uk2.u_n^2 = u_0^2 + 2n + \sum_{k=0}^{n-1} \frac{1}{u_k^2} .

De som is o(n)o(n) (termen gaan naar 00, Cesàro), dus un22nu_n^2 \sim 2n en un2nu_n \sim \sqrt{2n}.

Verfijning: 1uk212k\frac{1}{u_k^2} \sim \frac{1}{2k}, dus door vergelijking (Stelling 6.6, of equivalenten van partiële sommen van positieve reeksen) k<nuk212lnn\sum_{k<n} u_k^{-2} \sim \frac12 \ln n. Dus

un2=2n+lnn2(1+o(1))+O(1)un=2n1+lnn4n+o(lnnn)=2n(1+lnn8n+o(lnnn)).u_n^2 = 2n + \frac{\ln n}{2}\,(1 + o(1)) + O(1) \quad\Longrightarrow\quad u_n = \sqrt{2n}\sqrt{1 + \frac{\ln n}{4n} + o\Bigl(\frac{\ln n}{n}\Bigr)} = \sqrt{2n}\Bigl(1 + \frac{\ln n}{8n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr).

Oefening 6.9 ★★★

(Een Riemannsom met een twist) Bepaal het asymptotische gedrag van

Sn=k=1n1n+klnn.S_n = \sum_{k=1}^{n} \frac{1}{n + k\ln n} .

(Factoriseer nn: Sn=1nk(1+klnnn)1S_n = \frac1n\sum_k \bigl(1 + \frac{k\ln n}{n}\bigr)^{-1}; herken een Riemann-achtige som met een langzaam variërende parameter t=lnnt = \ln n, bereken 01 ⁣du1+tu=ln(1+t)t\int_0^1 \frac{\dd u}{1 + tu} = \frac{\ln(1+t)}{t}, en besluit SnlnlnnlnnS_n \sim \frac{\ln\ln n}{\ln n}.)

Oplossing

Oplossing van Oefening 6.9.

Factoriseer nn en stel t=lnnt = \ln n:

Sn=1nk=1n11+tkn.S_n = \frac1n \sum_{k=1}^{n} \frac{1}{1 + t\,\frac kn} .

Voor vaste tt is de som een Riemannsom van u11+tuu \mapsto \frac{1}{1 + tu} op [0,1]\intcc{0}{1}; de functie is monotoon in uu, dus wordt de Riemannsom ingeklemd door de integraal verschoven over één mesh:

01 ⁣du1+tu1nSn01 ⁣du1+tu+1n\int_0^1 \frac{\dd u}{1 + tu} - \frac1n \leq S_n \leq \int_0^1 \frac{\dd u}{1 + tu} + \frac1n

(vergelijking van Riemannsommen van een monotone functie met haar integraal, geldig voor elke nn met haar eigen t=lnnt = \ln n). Nu 01 ⁣du1+tu=ln(1+t)t\int_0^1 \frac{\dd u}{1 + tu} = \frac{\ln(1 + t)}{t}, en 1n=o(lntt)\frac1n = o\bigl(\frac{\ln t}{t}\bigr): dus

Sn=ln(1+lnn)lnn+O(1n)    lnlnnlnn.S_n = \frac{\ln(1 + \ln n)}{\ln n} + O\Bigl(\frac 1n\Bigr) \;\sim\; \frac{\ln\ln n}{\ln n} .

Oefening 6.10

Bewijs de identiteit (lnn)lnn=nlnlnn(\ln n)^{\ln n} = n^{\ln\ln n}, en rangschik dan de volgende in toenemende o()o(\cdot)-orde in oneindig, met bewijzen: n2n^2, (lnn)lnn(\ln n)^{\ln n}, 2n2^n, n!n!, nnn^n.

Oplossing

Oplossing van Oefening 6.10.

Identiteit: (lnn)lnn=elnnlnlnn=(elnn)lnlnn=nlnlnn(\ln n)^{\ln n} = \eu^{\ln n\,\ln\ln n} = \bigl(\eu^{\ln n}\bigr)^{\ln\ln n} = n^{\ln\ln n}. Rangschikking: vergelijk logaritmen. ln(n2)=2lnn\ln(n^2) = 2\ln n; ln((lnn)lnn)=lnnlnlnn\ln\bigl((\ln n)^{\ln n}\bigr) = \ln n\ln\ln n; ln(2n)=nln2\ln(2^n) = n\ln2; ln(n!)=nlnnn+O(lnn)\ln(n!) = n\ln n - n + O(\ln n) (Stirling, of de ruwere inkleming lnn!nlnn\ln n! \sim n\ln n); ln(nn)=nlnn\ln(n^n) = n\ln n. Aangezien 2lnn=o(lnnlnlnn)2\ln n = o(\ln n\ln\ln n), lnnlnlnn=o(n)\ln n\ln\ln n = o(n), nln2=o(nlnnn)n\ln 2 = o(n\ln n - n), en nlnnnnlnnn \ln n - n \sim n\ln n maar n!/nn0n! / n^n \to 0 (het verschil van logs is n+O(lnn)-n + O(\ln n) \to -\infty):

n2=o((lnn)lnn),(lnn)lnn=o(2n),2n=o(n!),n!=o(nn).n^2 = o\bigl((\ln n)^{\ln n}\bigr),\quad (\ln n)^{\ln n} = o(2^n),\quad 2^n = o(n!),\quad n! = o(n^n).

(Voor elke stap: het verschil van logaritmen gaat naar ++\infty, dus de ratio naar 00.)

Oefening 6.11 ★★

(Staart van 1/k2\sum 1/k^2, twee termen) Met de exacte telescopatie k>n1k(k+1)=1n+1\sum_{k > n} \frac{1}{k(k+1)} = \frac{1}{n+1} en de ontbinding 1k2=1k(k+1)+1k2(k+1)\frac1{k^2} = \frac{1}{k(k+1)} + \frac{1}{k^2(k+1)}, bewijs

k>n1k2=1n12n2+O(1n3).\sum_{k > n} \frac{1}{k^2} = \frac1n - \frac{1}{2n^2} + O\Bigl(\frac{1}{n^3}\Bigr).
Oplossing

Oplossing van Oefening 6.11.

Ontbind 1k2=1k(k+1)+1k2(k+1)\frac1{k^2} = \frac1{k(k+1)} + \frac1{k^2(k+1)} en sommeer voor k>nk > n:

k>n1k2=1n+1+k>n1k2(k+1),\sum_{k>n}\frac1{k^2} = \frac1{n+1} + \sum_{k>n}\frac{1}{k^2(k+1)} ,

de eerste som exact telescoperend (1k(k+1)=1k1k+1\frac1{k(k+1)} = \frac1k - \frac1{k+1}). Voor de tweede: 1k2(k+1)=1k3+O(1k4)\frac{1}{k^2(k+1)} = \frac1{k^3} + O\bigl(\frac1{k^4}\bigr) (aangezien 1k2(k+1)1k3=1k3(k+1)\frac{1}{k^2(k+1)} - \frac1{k^3} = \frac{-1}{k^3(k+1)}), en door de integraalvergelijking k>n1k3=12n2+O(1n3)\sum_{k>n}\frac1{k^3} = \frac1{2n^2} + O\bigl(\frac1{n^3}\bigr), k>n1k4=O(1n3)\sum_{k>n}\frac1{k^4} = O\bigl(\frac1{n^3}\bigr). Dus

k>n1k2=1n+1+12n2+O(1n3)=1n1n2+12n2+O(1n3)=1n12n2+O(1n3),\sum_{k>n}\frac1{k^2} = \frac1{n+1} + \frac{1}{2n^2} + O\Bigl(\frac1{n^3}\Bigr) = \frac1n - \frac1{n^2} + \frac{1}{2n^2} + O\Bigl(\frac1{n^3}\Bigr) = \frac1n - \frac{1}{2n^2} + O\Bigl(\frac1{n^3}\Bigr),

met 1n+1=1n1n2+O(1n3)\frac1{n+1} = \frac1n - \frac1{n^2} + O\bigl(\frac1{n^3}\bigr).

Oefening 6.12 ★★★

Zij u0=12u_0 = \frac12 en un+1=un+eunu_{n+1} = u_n + \eu^{-u_n}. Bewijs dat unu_n \to \infty, dan — vn=eunv_n = \eu^{u_n} stellend en tonend vn+1=vn+1+12vn+O(vn2)v_{n+1} = v_n + 1 + \frac{1}{2v_n} + O\bigl(v_n^{-2}\bigr) — vestig

un=lnn+lnn2n+O(1n).u_n = \ln n + \frac{\ln n}{2n} + O\Bigl(\frac1n\Bigr).
Oplossing

Oplossing van Oefening 6.12.

(un)(u_n) stijgt; als zij begrensd was zou zij naar een eindige \ell convergeren met =+e\ell = \ell + \eu^{-\ell}: onmogelijk. Dus unu_n \to \infty. Zij vn=eunv_n = \eu^{u_n} \to \infty: dan

vn+1=eun+eun=vne1/vn=vn(1+1vn+12vn2+O(vn3))=vn+1+12vn+O(vn2).v_{n+1} = \eu^{u_n + \eu^{-u_n}} = v_n\,\eu^{1/v_n} = v_n\Bigl(1 + \frac1{v_n} + \frac1{2v_n^2} + O\bigl(v_n^{-3}\bigr)\Bigr) = v_n + 1 + \frac{1}{2v_n} + O\bigl(v_n^{-2}\bigr).

Sommeren van vk+1vk=1+O(1)v_{k+1} - v_k = 1 + O(1) geeft eerst vn=n+O(n)v_n = n + O(n), dus vncnv_n \geq cn uiteindelijk; opnieuw sommeren met 12vk=O(1k)\frac1{2v_k} = O(\frac1k) geeft vn=n+O(lnn)v_n = n + O(\ln n). Nog een doorgang: 12vk=12k(1+O(lnkk))\frac{1}{2v_k} = \frac{1}{2k}\bigl(1 + O\bigl(\tfrac{\ln k}k\bigr)\bigr), dus

vn=n+k<n12k+O(1)=n+lnn2+O(1).v_n = n + \sum_{k<n}\frac1{2k} + O(1) = n + \frac{\ln n}2 + O(1).

Ten slotte un=lnvn=lnn+ln(1+lnn2n+O(1n))=lnn+lnn2n+O(1n)u_n = \ln v_n = \ln n + \ln\Bigl(1 + \frac{\ln n}{2n} + O\bigl(\tfrac1n\bigr)\Bigr) = \ln n + \frac{\ln n}{2n} + O\bigl(\tfrac1n\bigr).

6.6 Vraagstuk: Bootstrappen, van Euler–Maclaurin tot de priemen

Een impliciete of geaccumuleerde grootheid levert haar asymptotiek zelden in één keer; men extraheert ze in doorgangen, elke doorgang voedt de vorige schatting terug in de definiërende relatie. Dit weekendvraagstuk traint die lus op verse vergelijkingen, bewijst de eerste-orde formule van Euler–Maclaurin (de trapeziumopwaardering van de reeks–integraalvergelijking, met rigoureuze foutenbalken), inverteert xlnx=nx\ln x = n, en cash de beroemdste cheque van de methode: uit de toegegeven priemgetalstelling, de asymptotische wet pnnlnnp_n \sim n\ln n van het nn-de priemgetal.

Probleem 6.1

Weekendvraagstuk — de Euler–Maclaurin-correctie en de asymptotiek van het nn-de priemgetal

Deel I — De bootstrap-lus op een verse vergelijking.

  1. Bewijs de uniciteitsclaim van Definitie 6.2: als f=ikciφi+o(φk)=ikciφi+o(φk)f = \sum_{i\leq k} c_i\varphi_i + o(\varphi_k) = \sum_{i \leq k} c_i'\varphi_i + o(\varphi_k) langs dezelfde schaal, dan ci=cic_i = c_i' voor alle ii. Duw dan het gemengde voorbeeld van de cursus één tree verder:

    1xlnx=1x+lnxx2+(lnx)2x3+o((lnx)2x3)(x+),\frac{1}{x - \ln x} = \frac1x + \frac{\ln x}{x^2} + \frac{(\ln x)^2}{x^3} + o\Bigl(\frac{(\ln x)^2}{x^3}\Bigr) \qquad (x \to +\infty),

    en leg uit waarom geen term cx2\frac{c}{x^2} verschijnt.

  2. Toon dat voor elke n1n \geq 1 de vergelijking ex+x=n\eu^x + x = n precies één reële oplossing xnx_n heeft, en dat xn+x_n \to +\infty met xnlnnx_n \sim \ln n.
  3. Bootstrap tweemaal:

    xn=lnnlnnn(lnn)22n2+o((lnn)2n2).x_n = \ln n - \frac{\ln n}{n} - \frac{(\ln n)^2}{2n^2} + o\Bigl(\frac{(\ln n)^2}{n^2}\Bigr).
  4. Controleer numeriek bij n=1000n = 1000: vergelijk x10006.90083x_{1000} \approx 6.90083 met de één-, twee- en drietermwaarden van vraag 3, tot vijf decimalen.

Deel II — Euler–Maclaurin, orde één.

  1. Bewijs de trapeziumkernidentiteit: voor gg van klasse C2C^2 op [0,1]\intcc{0}{1},

    01g(t) ⁣dt=g(0)+g(1)21201t(1t)g(t) ⁣dt\int_0^1 g(t)\,\dd t = \frac{g(0) + g(1)}{2} - \frac12\int_0^1 t(1 - t)\,g''(t)\,\dd t

    (integreer 12t(1t)g\frac12 t(1-t)g'' tweemaal partieel).

  2. Zij ff C2C^2 op [1,+)\intco{1}{+\infty} met 1f<\int_1^\infty \abs{f''} < \infty. Toon dat

    En=k=1nf(k)1nff(1)+f(n)2E_n = \sum_{k=1}^{n} f(k) - \int_1^n f - \frac{f(1) + f(n)}{2}

    convergeert naar een constante EE, met de staartgrens EEn18nf\abs{E - E_n} \leq \frac18\int_n^\infty\abs{f''}: de formule van Euler–Maclaurin tot eerste orde.

  3. Pas dit toe op f(t)=1tf(t) = \frac1t: bewijs

    Hn=lnn+γ+12n+εn,εn18n2,H_n = \ln n + \gamma + \frac{1}{2n} + \varepsilon_n, \qquad \abs{\varepsilon_n} \leq \frac{1}{8n^2},

    versterkend Oefening 6.3 (identificeer de constante met γ\gamma door te vergelijken met Voorbeeld 6.7).

  4. Extraheer de volgende coëfficiënt: toon εn=112n2+o(1n2)\varepsilon_n = -\frac{1}{12n^2} + o\bigl(\frac1{n^2}\bigr) (de incrementen van EnE_n zijn 1201t(1t)f(n+t) ⁣dt=112f(n)+o(f(n))\frac12\int_0^1t(1-t)f''(n+t)\dd t = \frac1{12}f''(n) + o(f''(n)); sommeer de staart met Stelling 6.6).
  5. Pas vraag 6 toe op f=lnf = \ln: herleid in drie regels de convergentie van dn=lnn!(n+12)lnn+nd_n = \ln n! - (n + \frac12)\ln n + n (Stap 1 van Stelling 6.13), met de bonusfoutsnelheid dn=d+O(1n)d_n = d + O\bigl(\frac1n\bigr).
  6. Pas vraag 6 toe op f(t)=1tf(t) = \frac{1}{\sqrt t}: toon

    k=1n1k=2n+c+12n+O(1n3/2)\sum_{k=1}^{n}\frac1{\sqrt k} = 2\sqrt n + c + \frac{1}{2\sqrt n} + O\Bigl(\frac{1}{n^{3/2}}\Bigr)

    voor zekere constante cc, en evalueer alle termen bij n=104n = 10^4 (de constante is c1.4604c \approx -1.4604).

Deel III — Inversie: de vergelijking xlnx=nx\ln x = n.

  1. Toon dat xlnx=nx\ln x = n precies één oplossing xn[1,+)x_n \in \intco{1}{+\infty} heeft voor n1n \geq 1, dat xnx_n \to \infty, en dat lnxnlnn\ln x_n \sim \ln n.
  2. Leid de eenterminversie xnnlnnx_n \sim \dfrac{n}{\ln n} af, dan bootstrap nog eens:

    lnxn=lnnlnlnn+o(1),xn=nlnn(1+lnlnnlnn+o(lnlnnlnn)).\ln x_n = \ln n - \ln\ln n + o(1), \qquad x_n = \frac{n}{\ln n}\Bigl(1 + \frac{\ln\ln n}{\ln n} + o\Bigl(\frac{\ln\ln n}{\ln n}\Bigr)\Bigr).
  3. Test bij n=106n = 10^6: de ware wortel is x87848x \approx 87\,848; vergelijk met de eenterm (72382\approx 72\,382) en tweeterm (86140\approx 86\,140) waarden, en leg de trage winst uit (de ontwikkelingsparameter is lnlnnlnn\frac{\ln\ln n}{\ln n}, slechts 0.19\approx 0.19 bij n=106n = 10^6).
  4. We geven toe de priemgetalstelling: het aantal π(x)\pi(x) van priemen x\leq x voldoet aan π(x)xlnx\pi(x) \sim \frac{x}{\ln x} als xx \to \infty (eerlijk bewezen in het volume van Jaar 3). Schrijvend pnp_n voor het nn-de priemgetal, rechtvaardig π(pn)=n\pi(p_n) = n, en voer de inversie van vragen 11–12 uit om

    pnnlnn.p_n \sim n \ln n .

    te bewijzen.

  5. Dividenden: (a) toon knpkn2lnn2\sum_{k \leq n} p_k \sim \frac{n^2\ln n}{2} (vergelijk klnk\sum k\ln k met tlnt ⁣dt\int t\ln t\,\dd t); (b) bereken de benaderde kans dat een uniform willekeurig geheel getal met 100100 cijfers priem is (ln10100230.26\ln 10^{100} \approx 230.26: ongeveer één op 230230).

Deel IV — De methode geëxporteerd: xtanx=1x\tan x = 1.

  1. Toon dat voor elke n1n \geq 1 de vergelijking tanx=1x\tan x = \frac1x precies één oplossing xnx_n heeft in (nπ,nπ+π2)\intoo{n\pi}{\,n\pi + \frac\pi2}, en dat zn=xnnπ0+z_n = x_n - n\pi \to 0^+.
  2. Eén term: zn1nπz_n \sim \dfrac{1}{n\pi}.
  3. Toon dat de ontwikkeling van znz_n geen cn2\frac{c}{n^2}-term heeft: zn=1nπ+O(1n3)z_n = \frac1{n\pi} + O\bigl(\frac{1}{n^3}\bigr).
  4. Drie termen: met arctanu=uu33+O(u5)\arctan u = u - \frac{u^3}3 + O(u^5) en 1xn=1nπzn(nπ)2+O(n3zn2)\frac1{x_n} = \frac{1}{n\pi} - \frac{z_n}{(n\pi)^2} + O(n^{-3}\cdot z_n^2), bewijs

    xn=nπ+1nπ43π3n3+o(1n3).x_n = n\pi + \frac{1}{n\pi} - \frac{4}{3\pi^3 n^3} + o\Bigl(\frac{1}{n^3}\Bigr).
  5. Controleer bij n=3n = 3: ware wortel x39.5293344x_3 \approx 9.5293344; vergelijk de één- en drietermwaarden, en contrasteer in één zin met de tanx=x\tan x = x van de cursus (Voorbeeld 6.23): waar elke rij in haar venster zit, en waarom.

Deel V — Een dynamische bootstrap, spelregels, synthese.

  1. Zij u0(0,π)u_0 \in \intoo{0}{\pi} en un+1=sinunu_{n+1} = \sin u_n. Toon un0u_n \to 0 dalend, en bereken de limiet van 1un+121un2\dfrac{1}{u_{n+1}^2} - \dfrac{1}{u_n^2} (ontwikkel sin2\sin^{-2} via sinu=uu36+o(u3)\sin u = u - \frac{u^3}6 + o(u^3)).
  2. Leid via Cesàro-middelen (volume van Jaar 1) de klassieker

    un3n.u_n \sim \sqrt{\frac{3}{n}} .

    af.

  3. (Gecertificeerde numeriek) Toon met de rigoureuze grens van vraag 7 dat evalueren van lnn+γ+12n\ln n + \gamma + \frac1{2n} bij n=106n = 10^6 H106H_{10^6} levert met fout ten hoogste 1.2510131.25\cdot10^{-13} — een miljoen-termen-som berekend tot dertien cijfers door drie termen.
  4. (Spelregels) Bewijs of weerleg, met bewijzen of tegenvoorbeelden: (a) als unvn+u_n \sim v_n \to +\infty dan lnunlnvn\ln u_n \sim \ln v_n; (b) als unvnu_n \sim v_n dan eunevn\eu^{u_n} \sim \eu^{v_n}; (c) als fgf \sim g in ++\infty (f,gf, g differentieerbaar) dan fgf' \sim g'.
  5. (Synthese) In één zin elk: de bootstrap-lus van Methode 6.22 zoals gebruikt in Delen I, III, IV; wat de trapeziumcorrectie toevoegt aan Stelling 6.6; waarom inversie van xlnxx\ln x precies de brug van π(x)\pi(x) naar pnp_n is; en welke van de regels van vraag 24 welke stap beschermde. Noem de twee toppen: de formule van Euler–Maclaurin (eerste orde), en de asymptotische wet van het nn-de priemgetal.
Oplossing

Oplossing van Probleem 6.1.

1. Aftrekken van de twee ontwikkelingen: i(cici)φi=o(φk)\sum_i (c_i - c_i')\varphi_i = o(\varphi_k). Als zekere coëfficiënt verschilt, zij i0i_0 de eerste: delen door φi0\varphi_{i_0} en gebruiken φj=o(φi0)\varphi_j = o(\varphi_{i_0}) voor j>i0j > i_0 geeft ci0ci0=o(1)c_{i_0} - c_{i_0}' = o(1): nul, tegenspraak. Voor de ontwikkeling: met u=lnxx0u = \frac{\ln x}x \to 0,

1xlnx=1x11u=1x(1+u+u2+O(u3))=1x+lnxx2+(lnx)2x3+o((lnx)2x3).\frac{1}{x - \ln x} = \frac1x\cdot\frac{1}{1 - u} = \frac1x\bigl(1 + u + u^2 + O(u^3)\bigr) = \frac1x + \frac{\ln x}{x^2} + \frac{(\ln x)^2}{x^3} + o\Bigl(\frac{(\ln x)^2}{x^3}\Bigr).

Geen cx2\frac c{x^2}-term verschijnt omdat de ontwikkeling een meetkundige reeks in u=lnxxu = \frac{\ln x}{x} is: elke term draagt even veel machten van lnx\ln x als van 1x\frac1x voorbij de eerste; de schaal trede 1x2\frac1{x^2} (coëfficiënt van (lnx)0(\ln x)^0) is simpelweg afwezig, met coëfficiënt 00.

2. f(x)=ex+xf(x) = \eu^x + x is continu, strikt stijgend, met limieten -\infty en ++\infty: een bijectie RR\R \to \R, dus xn=f1(n)x_n = f^{-1}(n) bestaat en is uniek, en xn+x_n \to +\infty (f1f^{-1} stijgt naar ++\infty). Uit exn=nxn\eu^{x_n} = n - x_n: xn=ln(nxn)lnnx_n = \ln(n - x_n) \leq \ln n, dus xn/n0x_n/n \to 0 en xn=lnn+ln(1xn/n)=lnn+o(1)lnnx_n = \ln n + \ln(1 - x_n/n) = \ln n + o(1) \sim \ln n.

3. Schrijf un=xn/nu_n = x_n/n. Tweede doorgang: un=lnn+o(1)nu_n = \frac{\ln n + o(1)}{n}, dus

xn=lnn+ln(1un)=lnnun+O(un2)=lnnlnnn+o(lnnn).x_n = \ln n + \ln(1 - u_n) = \ln n - u_n + O(u_n^2) = \ln n - \frac{\ln n}{n} + o\Bigl(\frac{\ln n}n\Bigr).

Derde doorgang: nu un=lnnnlnnn2+o(lnnn2)u_n = \frac{\ln n}{n} - \frac{\ln n}{n^2} + o\bigl(\frac{\ln n}{n^2}\bigr), en ln(1un)=unun22+O(un3)\ln(1 - u_n) = -u_n - \frac{u_n^2}2 + O(u_n^3):

xn=lnnlnnn+lnnn2(lnn)22n2+o((lnn)2n2)=lnnlnnn(lnn)22n2+o((lnn)2n2),x_n = \ln n - \frac{\ln n}n + \frac{\ln n}{n^2} - \frac{(\ln n)^2}{2n^2} + o\Bigl(\frac{(\ln n)^2}{n^2}\Bigr) = \ln n - \frac{\ln n}{n} - \frac{(\ln n)^2}{2n^2} + o\Bigl(\frac{(\ln n)^2}{n^2}\Bigr),

de term lnnn2\frac{\ln n}{n^2} wordt geabsorbeerd in o((lnn)2n2)o\bigl(\frac{(\ln n)^2}{n^2}\bigr).

4. Bij n=1000n = 1000: ln10006.90776\ln 1000 \approx 6.90776 (fout 71037\cdot10^{-3}); twee termen: 6.900856.90085 (fout 21052\cdot10^{-5}); drie termen: 6.900826.90082 (fout onder 10510^{-5}), tegen x10006.90083x_{1000} \approx 6.90083. Elke doorgang koopt ongeveer de voorspelde factor lnnn\frac{\ln n}{n}.

5. Twee partiële integraties, startend van rechts: met  ⁣d ⁣dt[12t(1t)]=12t\frac{\dd}{\dd t}\bigl[\tfrac12t(1-t)\bigr] = \tfrac12 - t en t(1t)t(1-t) verdwijnend aan beide einden,

1201t(1t)g(t) ⁣dt=01(12t)g(t) ⁣dt=[(12t)g]0101g=g(0)+g(1)201g.\frac12\int_0^1 t(1-t)g''(t)\dd t = -\int_0^1\Bigl(\frac12 - t\Bigr)g'(t)\dd t = -\Bigl[\Bigl(\frac12 - t\Bigr)g\Bigr]_0^1 - \int_0^1 g = \frac{g(0) + g(1)}2 - \int_0^1 g .

Herschikt is dit de gestelde identiteit.

6. Bereken het increment, pas dan vraag 5 toe op g(t)=f(n+t)g(t) = f(n + t):

En+1En=f(n+1)nn+1 ⁣ff(n+1)f(n)2=f(n)+f(n+1)2nn+1 ⁣f=1201t(1t)f(n+t) ⁣dt.\begin{align*} E_{n+1} - E_n &= f(n{+}1) - \int_n^{n+1}\!f - \frac{f(n{+}1) - f(n)}2 \\ &= \frac{f(n) + f(n{+}1)}2 - \int_n^{n+1}\!f = \frac12\int_0^1 t(1-t)f''(n+t)\dd t . \end{align*}

Aangezien 0t(1t)140 \leq t(1-t) \leq \frac14: En+1En18nn+1f\abs{E_{n+1} - E_n} \leq \frac18\int_n^{n+1}\abs{f''}, waarvan de som over nn convergeert door de hypothese: (En)(E_n) convergeert (absoluut sommeerbare incrementen) naar zekere EE, met

EEnknEk+1Ek18nf.\abs{E - E_n} \leq \sum_{k\geq n}\abs{E_{k+1} - E_k} \leq \frac18\int_n^\infty\abs{f''} .

7. f(t)=1tf(t) = \frac1t: f(t)=2t3f''(t) = \frac2{t^3}, 1f=1<\int_1^\infty\abs{f''} = 1 < \infty. Vraag 6:

Hn=lnn+1+1n2+E+(EnE)=lnn+(E+12)+12n+εn,H_n = \ln n + \frac{1 + \frac1n}{2} + E + (E_n - E) = \ln n + \Bigl(E + \frac12\Bigr) + \frac1{2n} + \varepsilon_n,

met εn=EnE18n2 ⁣dtt3=18n2\abs{\varepsilon_n} = \abs{E_n - E} \leq \frac18\int_n^\infty\frac{2\dd t}{t^3} = \frac1{8n^2}. Vergelijken met Hn=lnn+γ+o(1)H_n = \ln n + \gamma + o(1) (Voorbeeld 6.7) identificeert E+12=γE + \frac12 = \gamma.

8. Uit de incrementformule van vraag 6,

εn=EnE=kn1201t(1t)2 ⁣dt(k+t)3=kn(1k301t(1t) ⁣dt+O(1k4)),\varepsilon_n = E_n - E = -\sum_{k\geq n}\frac12\int_0^1 t(1-t)\,\frac{2\,\dd t}{(k+t)^3} = -\sum_{k \geq n}\Bigl(\frac1{k^3}\int_0^1t(1-t)\dd t + O\Bigl(\frac1{k^4}\Bigr)\Bigr),

met 1(k+t)3=1k3+O(1k4)\frac{1}{(k+t)^3} = \frac1{k^3} + O\bigl(\frac1{k^4}\bigr) uniform voor t[0,1]t \in \intcc01. Met 01t(1t)=16\int_0^1 t(1-t) = \frac16 en kn1k312n2\sum_{k\geq n}\frac1{k^3} \sim \frac{1}{2n^2} (Stelling 6.6):

εn=1612n2+o(1n2)=112n2+o(1n2).\varepsilon_n = -\frac16\cdot\frac{1}{2n^2} + o\Bigl(\frac1{n^2}\Bigr) = -\frac{1}{12n^2} + o\Bigl(\frac{1}{n^2}\Bigr).

9. f=lnf = \ln: f(t)=1t2f''(t) = -\frac1{t^2}, absoluut integreerbaar. Vraag 6 geeft

lnn!=1nlnt ⁣dt+lnn2+E+O(18n ⁣dtt2)=(n+12)lnnn+1+E+O(1n),\ln n! = \int_1^n\ln t\,\dd t + \frac{\ln n}2 + E + O\Bigl( \frac1{8}\int_n^\infty\frac{\dd t}{t^2}\Bigr) = \Bigl(n + \frac12\Bigr)\ln n - n + 1 + E + O\Bigl(\frac1n\Bigr),

dus dn=1+E+O(1n)d_n = 1 + E + O\bigl(\frac1n\bigr): convergentie van (dn)(d_n) — Stap 1 van Stelling 6.13 — plus de snelheid O(1/n)O(1/n). (De Stirlingwaarde van de limiet geeft E=ln2π1E = \ln\sqrt{2\pi} - 1.)

10. f(t)=t1/2f(t) = t^{-1/2}: f(t)=34t5/2f''(t) = \frac34 t^{-5/2}, absoluut integreerbaar. Vraag 6:

k=1n1k=2n2+1+1n2+E+O(n3/2)=2n+c+12n+O(n3/2),\sum_{k=1}^n \frac1{\sqrt k} = 2\sqrt n - 2 + \frac{1 + \frac1{\sqrt n}}2 + E + O\bigl(n^{-3/2}\bigr) = 2\sqrt n + c + \frac{1}{2\sqrt n} + O\bigl(n^{-3/2}\bigr),

met c=E32c = E - \frac32. Bij n=104n = 10^4: 2n=2002\sqrt n = 200, c1.46035c \approx -1.46035, 12n=0.005\frac1{2\sqrt n} = 0.005: voorspeld 198.54465198.54465, en inderdaad k104k1/2=198.544645\sum_{k\leq10^4}k^{-1/2} = 198.544645\dots — drie termen, zeven cijfers.

11. ttlntt \mapsto t\ln t is continu en strikt stijgend op [1,)\intco1\infty (afgeleide lnt+11\ln t + 1 \geq 1), van 00 naar ++\infty: een unieke xnx_n bestaat, en xnx_n \to \infty (anders zou xnlnxnx_n\ln x_n begrensd blijven). Logaritmen nemen in xnlnxn=nx_n\ln x_n = n: lnxn+lnlnxn=lnn\ln x_n + \ln\ln x_n = \ln n; aangezien lnlnxn=o(lnxn)\ln\ln x_n = o(\ln x_n), geeft delen door lnxn\ln x_n lnnlnxn1\frac{\ln n}{\ln x_n} \to 1: lnxnlnn\ln x_n \sim \ln n.

12. Uit xn=nlnxnx_n = \frac{n}{\ln x_n} en lnxnlnn\ln x_n \sim \ln n: xnnlnnx_n \sim \frac{n}{\ln n}. Volgende doorgang: lnlnxn=ln(lnn(1+o(1)))=lnlnn+o(1)\ln\ln x_n = \ln\bigl(\ln n\,(1 + o(1))\bigr) = \ln\ln n + o(1), dus lnxn=lnnlnlnn+o(1)\ln x_n = \ln n - \ln\ln n + o(1) en

xn=nlnnlnlnn+o(1)=nlnn11lnlnn+o(1)lnn=nlnn(1+lnlnnlnn+o(lnlnnlnn)).x_n = \frac{n}{\ln n - \ln\ln n + o(1)} = \frac{n}{\ln n}\cdot\frac{1}{1 - \frac{\ln\ln n + o(1)}{\ln n}} = \frac{n}{\ln n}\Bigl(1 + \frac{\ln\ln n}{\ln n} + o\Bigl(\frac{\ln\ln n}{\ln n}\Bigr)\Bigr).

13. Bij n=106n = 10^6: nlnn72382\frac{n}{\ln n} \approx 72\,382 (18%18\% ernaast), twee termen geven 86140\approx 86\,140 (1.9%1.9\% ernaast), tegen de ware x87848x \approx 87\,848. De winst per doorgang is slechts de factor lnlnnlnn2.6313.80.19\frac{\ln\ln n}{\ln n} \approx \frac{2.63}{13.8} \approx 0.19: logaritmische schalen convergeren met waanzinnige traagheid — een feit van het leven overal waar priemen betrokken zijn.

14. Er zijn precies nn priemen pn\leq p_n (namelijk p1,,pnp_1, \dots, p_n): π(pn)=n\pi(p_n) = n. De priemgetalstelling (toegegeven; volume van Jaar 3) geeft n=π(pn)pnlnpnn = \pi(p_n) \sim \frac{p_n}{\ln p_n}, d.w.z. pnnlnpnp_n \sim n\ln p_n: dit is de vergelijking xlnxnx\ln x \approx n achterstevoren gelezen. Logaritmen nemen: lnpn=lnn+lnlnpn+o(1)\ln p_n = \ln n + \ln\ln p_n + o(1), en lnlnpn=o(lnpn)\ln\ln p_n = o(\ln p_n) forceert lnpnlnn\ln p_n \sim \ln n zoals in vraag 11. Terugsubstitueren:

pnnlnpn=nlnnlnpnlnnnlnn.p_n \sim n\ln p_n = n\,\ln n\,\frac{\ln p_n}{\ln n} \sim n\ln n .

15. (a) Fixeer ε>0\varepsilon > 0; voor grote kk, (1ε)klnkpk(1+ε)klnk(1 - \varepsilon)k\ln k \leq p_k \leq (1 + \varepsilon)k\ln k. Door vergelijking met de stijgende tlntt\ln t (Stelling 6.6-achtige inkleming), knklnk=1ntlnt ⁣dt+O(nlnn)=n2lnn2n24+O(nlnn)n2lnn2\sum_{k\leq n}k\ln k = \int_1^n t\ln t\,\dd t + O(n\ln n) = \frac{n^2\ln n}2 - \frac{n^2}4 + O(n\ln n) \sim \frac{n^2\ln n}2. Dus knpk=n2lnn2(1+O(ε)+o(1))\sum_{k\leq n}p_k = \frac{n^2\ln n}{2}(1 + O(\varepsilon) + o(1)) voor elke ε\varepsilon: knpkn2lnn2\sum_{k\leq n}p_k \sim \frac{n^2\ln n}2. (b) Door de priemgetalstelling is onder de gehele getallen tot 1010010^{100} een proportie 1ln10100=1230.26\sim \frac{1}{\ln 10^{100}} = \frac1{230.26\dots} priem: een uniform willekeurig 100100-cijferig geheel getal is priem met kans ongeveer 1230\frac1{230}.

16. Op (nπ,nπ+π2)\intoo{n\pi}{n\pi + \frac\pi2} is g(x)=tanx1xg(x) = \tan x - \frac1x continu en strikt stijgend (g=1+tan2x+1x2>0g' = 1 + \tan^2x + \frac1{x^2} > 0), met g1nπ<0g \to -\frac1{n\pi} < 0 aan het linker einde en g+g \to +\infty aan het rechter: precies één wortel xnx_n. Aangezien tanzn=tanxn=1xn0\tan z_n = \tan x_n = \frac1{x_n} \to 0 met zn(0,π2)z_n \in \intoo{0}{\frac\pi2}: zn0+z_n \to 0^+.

17. tanznzn\tan z_n \sim z_n en 1xn1nπ\frac1{x_n} \sim \frac1{n\pi}: zn1nπz_n \sim \frac1{n\pi}.

18. zn=arctan1xnz_n = \arctan\frac1{x_n} en arctanu=u+O(u3)\arctan u = u + O(u^3). Met zn=O(1n)z_n = O(\frac1n):

1xn=1nπ11+znnπ=1nπznn2π2+O(1n4)=1nπ+O(1n3),\frac1{x_n} = \frac{1}{n\pi}\cdot\frac1{1 + \frac{z_n}{n\pi}} = \frac1{n\pi} - \frac{z_n}{n^2\pi^2} + O\Bigl(\frac1{n^4}\Bigr) = \frac1{n\pi} + O\Bigl(\frac1{n^3}\Bigr),

dus zn=1nπ+O(1n3)z_n = \frac1{n\pi} + O\bigl(\frac1{n^3}\bigr): de cn2\frac{c}{n^2}-trede draagt coëfficiënt 00, omdat de eerste correctie aan 1xn\frac1{x_n} zelf van grootte znn2=O(n3)\frac{z_n}{n^2} = O(n^{-3}) is.

19. Steek zn=1nπ+O(n3)z_n = \frac1{n\pi} + O(n^{-3}) in de vorige weergave:

1xn=1nπ1n3π3+O(1n5),\frac{1}{x_n} = \frac{1}{n\pi} - \frac{1}{n^3\pi^3} + O\Bigl(\frac1{n^5}\Bigr),

dan zn=arctan1xn=1xn13(1xn)3+O(1n5)=1nπ1n3π313n3π3+O(1n5)z_n = \arctan\frac1{x_n} = \frac1{x_n} - \frac{1}{3}\Bigl(\frac1{x_n}\Bigr)^3 + O\Bigl(\frac1{n^5}\Bigr) = \frac1{n\pi} - \frac{1}{n^3\pi^3} - \frac{1}{3n^3\pi^3} + O\Bigl(\frac1{n^5}\Bigr):

xn=nπ+1nπ43π3n3+O(1n5).x_n = n\pi + \frac{1}{n\pi} - \frac{4}{3\pi^3n^3} + O\Bigl(\frac1{n^5}\Bigr).

20. Bij n=3n = 3: één term 9.530889.53088, drie termen 9.529299.52929, ware wortel 9.529339.52933: fouten 1.51031.5\cdot10^{-3} en 51055\cdot10^{-5}. Contrast: voor tanx=x\tan x = x moet de wortel tan\tan enorm maken, dus klemt zij zich aan het rechter einde nπ+π2n\pi + \frac\pi2 van het venster, op afstand 1nπ\sim\frac1{n\pi} vóór de asymptoot; voor xtanx=1x\tan x = 1 moet de wortel tan\tan klein maken, dus zit zij net voorbij het linker einde nπn\pi, op afstand 1nπ\sim\frac1{n\pi} na het nulpunt. Zelfde methode, spiegelgeografie.

21. sinu<u\sin u < u op (0,π)\intoo0\pi en sin\sin beeldt (0,π)\intoo0\pi af in (0,1](0,π)\intoc01 \subseteq \intoo0\pi: na één stap u1(0,1]u_1 \in \intoc{0}{1}, dan daalt (un)(u_n) en is van onderen begrensd door 00: zij convergeert naar een vast punt van sin\sin, d.w.z. naar 00. Ontwikkeling: sinu=u(1u26+o(u2))\sin u = u(1 - \frac{u^2}6 + o(u^2)), dus

1un+121un2=1un2((1un26+o(un2))21)=1un2(un23+o(un2))13.\frac{1}{u_{n+1}^2} - \frac1{u_n^2} = \frac{1}{u_n^2}\Bigl(\bigl(1 - \tfrac{u_n^2}6 + o(u_n^2)\bigr)^{-2} - 1\Bigr) = \frac{1}{u_n^2}\Bigl(\frac{u_n^2}{3} + o(u_n^2)\Bigr) \longrightarrow \frac13 .

22. Door Cesàro (volume van Jaar 1) convergeert het gemiddelde van de incrementen naar dezelfde limiet:

1n1un2=1n(1u02+k=0n1(1uk+121uk2))13,\frac{1}{n}\cdot\frac{1}{u_n^2} = \frac1n\Bigl(\frac1{u_0^2} + \sum_{k=0}^{n-1} \Bigl(\frac1{u_{k+1}^2} - \frac1{u_k^2}\Bigr)\Bigr) \longrightarrow \frac13 ,

dus un23nu_n^2 \sim \frac3n en, alle termen positief, un3/nu_n \sim \sqrt{3/n}.

23. Door vraag 7, Hnlnnγ12n18n2\abs{H_n - \ln n - \gamma - \frac1{2n}} \leq \frac1{8n^2}. Bij n=106n = 10^6 is deze grens 181012=1.251013\frac{1}{8\cdot10^{12}} = 1.25\cdot10^{-13}: drie berekende termen leveren de miljoen-termen-harmonische som tot dertien cijfers, met een volledig rigoureus foutcertificaat — het hele punt van een asymptotische formule met expliciete rest.

24. (a) Waar: lnunlnvn=lnunvn0\ln u_n - \ln v_n = \ln\frac{u_n}{v_n} \to 0 terwijl lnvn+\ln v_n \to +\infty, dus de ratio van logaritmen gaat naar 11. (b) Onwaar: un=n+1vn=nu_n = n + 1 \sim v_n = n, maar eun/evn=e1\eu^{u_n}/\eu^{v_n} = \eu \neq 1. Equivalentie verdraagt additieve fouten o(1)o(1) in de exponent, niet O(1)O(1). (c) Onwaar: f(x)=x+sin(x2)g(x)=xf(x) = x + \sin(x^2) \sim g(x) = x in ++\infty, maar f(x)=1+2xcos(x2)f'(x) = 1 + 2x\cos(x^2) oscilleert onbegrensd terwijl g=1g' = 1: afgeleiden van equivalente functies hoeven helemaal niet vergelijkbaar te zijn.

25. De lus van Methode 6.22 liep identiek driemaal: lokaliseer de wortel, extraheer een ruwe term, voed haar terug voor de volgende orde — op ex+x=n\eu^x + x = n (Deel I), op xlnx=nx\ln x = n (Deel III), op xtanx=1x\tan x = 1 (Deel IV). De trapeziumcorrectie waardeert de reeks–integraalvergelijking op van “het verschil convergeert” tot een expliciete f(1)+f(n)2\frac{f(1) + f(n)}2-term met een gecertificeerde O(nf)O(\int_n^\infty \abs{f''})-rest — constanten en foutenbalken in plaats van louter convergentie. De brug naar priemen is pure inversie: de priemgetalstelling zegt π(x)lnxx\pi(x)\ln x \approx x, dus pnp_n, gedefinieerd door π(pn)=n\pi(p_n) = n, lost een xlnx=nx\ln x = n-vergelijking op — en erft haar asymptotiek. Regel (a) van vraag 24 legitimeerde elke overgang van unvnu_n \sim v_n naar lnunlnvn\ln u_n \sim \ln v_n (vragen 11, 14); de onwaarheid van (b) is waarom we nooit equivalenties exponentiëren. Toppen: de formule van Euler–Maclaurin tot eerste orde (vraag 6), en de asymptotische wet pnnlnnp_n \sim n\ln n van het nn-de priemgetal (vraag 14).