Mathematics · Boek 3 · Bachelor Year 1

Universitaire wiskunde — Bachelor jaar 1

Universitaire wiskunde — Bachelor jaar 1 · Bachelor Year 1

16Formules van Taylor en asymptotische ontwikkelingen

In de buurt van een punt is een gladde functie even goed als een veelterm — met een beheersbare fout. De formules van Taylor maken dat exact in drie smaken (integraalrestterm, restterm van Lagrange, restterm van Young), en de resulterende asymptotische ontwikkelingen worden, algebraïsch behandeld, het scherpste gereedschap van de elementaire analyse: limieten, equivalenten, lokaal gedrag, asymptoten.

16.1 Vergelijkingsnotatie

Definitie 16.1 (Landau-notatie)

Zij f,gf, g gedefinieerd in de buurt van x0x_0 (x0Rx_0 \in \R of ±\pm\infty). Men schrijft, als xx0x \to x_0:

  • f=o(g)f = o(g) (“kleine o”) wanneer f=εgf = \varepsilon g met ε(x)0\varepsilon(x) \to 0;
  • f=O(g)f = O(g) (“grote O”) wanneer f=ugf = u g met uu begrensd in de buurt van x0x_0;
  • fgf \sim g (“equivalent”) wanneer f=(1+ε)gf = (1 + \varepsilon) g met ε0\varepsilon \to 0 — equivalent hiermee: fg=o(g)f - g = o(g).

Dezelfde notatie geldt voor rijen (nn \to \infty).

Propositie 16.2 (Rekenregels)

Als xx0x \to x_0:

  1. \sim is een equivalentierelatie; uit fgf \sim g volgt dat ff en gg dezelfde limieten, dezelfde tekens (in de buurt van x0x_0) en dezelfde afwezigheid van nulpunten hebben;
  2. equivalenten vermenigvuldigen en delen: uit f1g1f_1 \sim g_1 en f2g2f_2 \sim g_2 volgt f1f2g1g2f_1 f_2 \sim g_1 g_2 en f1f2g1g2\frac{f_1}{f_2} \sim \frac{g_1}{g_2};
  3. equivalenten tellen niet op: x+1xx + 1 \sim x en xx+2-x \sim -x + 2 in ++\infty, en toch zijn de sommen 11 en 22 niet equivalent. Om op te tellen keer je terug naar ontwikkelingen met expliciete o()o(\cdot)-termen;
  4. o(g)+o(g)=o(g)o(g) + o(g) = o(g),   uo(g)=o(ug)\;u \cdot o(g) = o(ug),   o(o(g))=o(g)\;o(o(g)) = o(g), en fg    f=g+o(g)f \sim g \iff f = g + o(g).

Bewijs. Elk punt is een korte manipulatie van de definities; zo is bijvoorbeeld f1f2=(1+ε1)(1+ε2)g1g2f_1 f_2 = (1+\varepsilon_1)(1+\varepsilon_2) g_1 g_2 en (1+ε1)(1+ε2)1(1 + \varepsilon_1)(1+\varepsilon_2) \to 1. Twee punten van (4) verdienen hun eigen regel. uo(g)=o(ug)u\cdot o(g) = o(ug): is f=εgf = \varepsilon g met ε0\varepsilon \to 0, dan is uf=ε(ug)uf = \varepsilon\,(ug) met dezelfde ε\varepsilon. o(o(g))=o(g)o(o(g)) = o(g): is f=ε1hf = \varepsilon_1 h en h=ε2gh = \varepsilon_2 g met beide εi0\varepsilon_i \to 0, dan is f=(ε1ε2)gf = (\varepsilon_1\varepsilon_2) g en is het product van de twee oneindig kleine grootheden er ook een. De equivalentie fg    f=g+o(g)f \sim g \iff f = g + o(g) is de definitie twee keer gelezen: fg=εgf - g = \varepsilon g. Het tegenvoorbeeld in (3) is het bewijs van (3).

Voorbeeld 16.3 (De vergelijkingsschaal)

Als x+x \to +\infty luidt de standaardschaal, in toenemende orde van kracht:

1  =  o(lnx),lnx=o(x0.01),x0.01=o(x),x=o(x10),x10=o(ex),ex=o(e2x),1 \;=\; o(\ln x), \quad \ln x = o(x^{0.01}), \quad x^{0.01} = o(\sqrt x), \quad \sqrt x = o(x^{10}), \quad x^{10} = o(\eu^{x}), \quad \eu^x = o(\eu^{2x}) ,

waarbij elke stap een geval is van de groeivergelijkingen van Propositie 4.6 (machten verslaan logaritmen, exponentiëlen verslaan machten, en binnen één familie beslist de exponent). Twee gewoonten die het vormen waard zijn: ten eerste gaat een O()O(\cdot) die in een kleinere klasse landt stilzwijgend over (O(lnx)O(\ln x) is ook o(x0.01)o(x^{0.01})); ten tweede keert bij x0+x \to 0^+ de hele ladder om via de substitutie x1xx \mapsto \frac1x — daar is lnx=o(x0.01)\ln x = o(x^{-0.01}), zodat “xαlnx0x^\alpha \ln x \to 0” geldt voor elke α>0\alpha > 0. De schaal helder houden is de helft van elk asymptotisch argument in Hoofdstuk 17.

Voorbeeld 16.4 (Eenduidigheid van ontwikkelingen, met een pariteitsdividend)

Laat een functie twee ontwikkelingen in 00 tot dezelfde orde toe,

a0+a1x++anxn+o(xn)=b0+b1x++bnxn+o(xn),a_0 + a_1 x + \dots + a_n x^n + o(x^n) = b_0 + b_1 x + \dots + b_n x^n + o(x^n),

dan is ak=bka_k = b_k voor elke kk: trek af en stel ck=akbkc_k = a_k - b_k; evalueer de identiteit c0+c1x++cnxn=o(xn)c_0 + c_1 x + \dots + c_n x^n = o(x^n) bij x0x \to 0 om c0=0c_0 = 0 te krijgen; deel door xx en herhaal — elke deling is geoorloofd omdat de overblijvende uitdrukking opnieuw o(xnk)o(x^{n-k}) is. De coëfficiënten zijn dus intrinsiek, en men mag ze langs elke weg berekenen (met afgeleiden van Taylor, met algebra op bekende ontwikkelingen, door integratie): alle wegen moeten overeenstemmen. Dividend: een even functie heeft alleen even machten in haar ontwikkeling — vervang xx door x-x en beroep je op de eenduidigheid; evenzo hebben oneven functies oneven machten. Dat is de reden waarom cos\cos in de tabel hieronder o(x2p+1)o(x^{2p+1}) draagt en niet o(x2p)o(x^{2p}): de ontbrekende oneven term is gratis informatie, één orde nauwkeurigheid voor niets.

16.2 De drie formules van Taylor

Stelling 16.5 (Taylor met integraalrestterm)

Zij ff van klasse Cn+1C^{n+1} op een interval dat aa en xx bevat. Dan geldt

f(x)=k=0nf(k)(a)k!(xa)k+ax(xt)nn!f(n+1)(t) ⁣dt.f(x) = \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}\,(x - a)^k + \int_a^x \frac{(x - t)^n}{n!}\, f^{(n+1)}(t)\, \dd t .

Bewijs. Inductie op nn. Voor n=0n = 0 is f(x)=f(a)+axf(t) ⁣dtf(x) = f(a) + \int_a^x f'(t)\dd t de hoofdstelling (Stelling 15.9). Stap: integreer de restterm partieel,

ax(xt)nn!f(n+1)(t) ⁣dt=[(xt)n+1(n+1)!f(n+1)(t)]ax+ax(xt)n+1(n+1)!f(n+2)(t) ⁣dt,\int_a^x \frac{(x-t)^n}{n!} f^{(n+1)}(t)\,\dd t = \Bigl[-\frac{(x-t)^{n+1}}{(n+1)!} f^{(n+1)}(t)\Bigr]_a^x + \int_a^x \frac{(x-t)^{n+1}}{(n+1)!} f^{(n+2)}(t)\,\dd t ,

waarbij de haakterm de bijdrage f(n+1)(a)(n+1)!(xa)n+1\frac{f^{(n+1)}(a)}{(n+1)!}(x - a)^{n+1} levert.

Voorbeeld 16.6 (Een exacte ontwikkeling met haar restterm)

Voor ln(1+x)\ln(1 + x) kan de integraalrestterm volledig expliciet worden gemaakt zonder ook maar iets nn keer af te leiden: integreer de eindige meetkundige identiteit 11+t=k=0n1(t)k+(t)n1+t\frac{1}{1+t} = \sum_{k=0}^{n-1}(-t)^k + \frac{(-t)^n}{1+t} van 00 tot xx:

ln(1+x)=k=1n(1)k1xkk+(1)n0xtn1+t ⁣dt,\ln(1 + x) = \sum_{k=1}^{n} \frac{(-1)^{k-1}x^k}{k} + (-1)^n \int_0^x \frac{t^n}{1 + t}\,\dd t ,

en voor 0x10 \leq x \leq 1 wordt de restterm begrensd door 0xtn ⁣dt=xn+1n+1\int_0^x t^n\,\dd t = \frac{x^{n+1}}{n+1}. Dit is in twee opzichten sterker dan Taylor–Young: het is een identiteit die geldt voor een vaste xx (niet alleen voor x0x \to 0), en de foutgrens is numeriek. De weekendopgave (Probleem 16.1) leeft van zulke exacte vormen; Taylor–Young hieronder is het lichtere gereedschap voor limieten, waar alleen de gedaante van de fout telt.

Stelling 16.7 (Ongelijkheid van Taylor–Lagrange)

Zij ff van klasse Cn+1C^{n+1} met f(n+1)M\abs{f^{(n+1)}} \leq M tussen aa en xx. Dan geldt

f(x)k=0nf(k)(a)k!(xa)kMxan+1(n+1)!.\Bigl| f(x) - \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}(x-a)^k \Bigr| \leq M\, \frac{\abs{x - a}^{n+1}}{(n+1)!} .

Bewijs. Begrens de integraalrestterm: ax(xt)nn!f(n+1)(t) ⁣dtMaxxtnn! ⁣dt=Mxan+1(n+1)!\bigl|\int_a^x \frac{(x-t)^n}{n!} f^{(n+1)}(t)\,\dd t\bigr| \leq M \bigl|\int_a^x \frac{\abs{x-t}^n}{n!}\dd t\bigr| = M\frac{\abs{x-a}^{n+1}}{(n+1)!}.

Voorbeeld 16.8 (Gecertificeerd rekenwerk)

Hoeveel is 1.02\sqrt{1.02}? Pas Taylor–Lagrange toe op f(t)=1+tf(t) = \sqrt{1 + t} in a=0a = 0, orde 22, x=0.02x = 0.02:

1.021+0.0220.0228=1.00995,f(t)=38(1+t)5/238,\sqrt{1.02} \approx 1 + \frac{0.02}{2} - \frac{0.02^2}{8} = 1.00995 , \qquad \abs{f'''(t)} = \frac{3}{8}(1+t)^{-5/2} \leq \frac 38 ,

dus is de fout hoogstens 380.0236=5107\frac38 \cdot \frac{0.02^3}{6} = 5\cdot10^{-7}: 1.02=1.00995\sqrt{1.02} = 1.00995 met zes gecertificeerde decimalen (werkelijke waarde 1.00995049381.0099504938\dots — de grens is bijna scherp). Het afsluitende inzicht: Taylor–Young zegt alleen hoe snel de fout verdwijnt; Taylor–Lagrange maakt van dezelfde veelterm een certificaat, een getal met een bewezen foutmarge. Telkens wanneer in dit boek een uitspraak over decimalen wordt gedaan, staat er een grens van het type Lagrange achter; de weekendopgave (Probleem 16.1) industrialiseert het idee.

Stelling 16.9 (Taylor–Young)

Zij ff nn keer afleidbaar in aa. Dan geldt, als xax \to a:

f(x)=k=0nf(k)(a)k!(xa)k+o((xa)n).f(x) = \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}\,(x-a)^k + o\bigl((x-a)^n\bigr) .

Bewijs. Inductie op nn. Voor n=1n = 1 is dit de definitie van de afgeleide (Definitie 14.1). Neem de uitspraak aan op orde n1n - 1, en zij ff nn keer afleidbaar in aa. Pas de inductiehypothese toe op ff' (die n1n-1 keer afleidbaar is in aa):

f(t)=k=0n1f(k+1)(a)k!(ta)k+r(t),r(t)=o((ta)n1).f'(t) = \sum_{k=0}^{n-1} \frac{f^{(k+1)}(a)}{k!}(t-a)^k + r(t), \qquad r(t) = o\bigl((t-a)^{n-1}\bigr).

Zij g(x)=f(x)k=0nf(k)(a)k!(xa)kg(x) = f(x) - \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}(x-a)^k; dan is g=rg' = r en g(a)=0g(a) = 0. Zij ε>0\varepsilon > 0 gegeven; kies δ\delta met r(t)εtan1\abs{r(t)} \leq \varepsilon\abs{t - a}^{n-1} voor taδ\abs{t-a} \leq \delta; de middelwaarde-ongelijkheid (Stelling 14.9), toegepast op het segment van aa tot xx (waar gεxan1\abs{g'} \leq \varepsilon\abs{x-a}^{n-1}), levert g(x)εxan\abs{g(x)} \leq \varepsilon\abs{x - a}^n: precies g(x)=o((xa)n)g(x) = o((x-a)^n).

Opmerking 16.10 (Drie formules, drie prijzen, drie producten)

De hypothesen lopen precies in de pas met de conclusies. Taylor–Young vraagt het minst (nn afgeleiden in het punt alleen) en levert het minst: een kwalitatieve o((xa)n)o((x-a)^n), perfect voor limieten, nutteloos voor gecertificeerde cijfers. De ongelijkheid van Lagrange vraagt Cn+1C^{n+1} op het interval en daar een grens MM, en geeft een numerieke foutmarge terug. De integraalvorm vraagt dezelfde regelmaat en geeft het meest terug: de fout als een expliciet object dat men kan omvormen (partieel integreren, stuksgewijs begrenzen, van variabele veranderen) — het is de vorm die de irrationaliteitsmachine van Probleem 15.1 aandreef. De zwakste formule kiezen die het doel draagt is geen pedanterie: de vlakke functie van Probleem 16.1 voldoet aan Taylor–Young in elke orde, terwijl elke sterkere conclusie erover buiten 00 vals is.

Propositie 16.11 (Standaardontwikkelingen in 00)

Als x0x \to 0 geldt, voor elke vaste orde nn:

ex=1+x+x22!++xnn!+o(xn),cosx=1x22!+x44!+(1)px2p(2p)!+o(x2p+1),sinx=xx33!++(1)px2p+1(2p+1)!+o(x2p+2),11x=1+x+x2++xn+o(xn),ln(1+x)=xx22+x33+(1)n1xnn+o(xn),(1+x)α=1+αx+α(α1)2!x2++(αn)xn+o(xn),\begin{align*} \eu^x &= 1 + x + \frac{x^2}{2!} + \dots + \frac{x^n}{n!} + o(x^n),\\ \cos x &= 1 - \frac{x^2}{2!} + \frac{x^4}{4!} - \dots + \frac{(-1)^p x^{2p}}{(2p)!} + o(x^{2p+1}),\\ \sin x &= x - \frac{x^3}{3!} + \dots + \frac{(-1)^p x^{2p+1}}{(2p+1)!} + o(x^{2p+2}),\\ \frac{1}{1 - x} &= 1 + x + x^2 + \dots + x^n + o(x^n),\\ \ln(1 + x) &= x - \frac{x^2}{2} + \frac{x^3}{3} - \dots + \frac{(-1)^{n-1} x^n}{n} + o(x^n),\\ (1 + x)^\alpha &= 1 + \alpha x + \frac{\alpha(\alpha-1)}{2!}x^2 + \dots + \binom{\alpha}{n} x^n + o(x^n), \end{align*}

waarbij (αn)=α(α1)(αn+1)n!\binom{\alpha}{n} = \frac{\alpha(\alpha - 1)\cdots(\alpha - n + 1)}{n!} voor reële α\alpha. (cosh\cosh en sinh\sinh: hetzelfde als cos\cos en sin\sin, zonder de wisselende tekens.)

Bewijs. Elke functie is glad in de buurt van 00 met gemakkelijk te evalueren afgeleiden: (ex)(k)=ex(\eu^x)^{(k)} = \eu^x; de afgeleiden van sin\sin en cos\cos doorlopen een cyclus met periode 44; ((1+x)α)(k)=α(α1)(αk+1)(1+x)αk\bigl((1+x)^\alpha\bigr)^{(k)} = \alpha(\alpha - 1)\cdots(\alpha - k + 1)(1+x)^{\alpha - k}; (ln(1+x))(k)=(1)k1(k1)!(1+x)k\bigl(\ln(1+x)\bigr)^{(k)} = \frac{(-1)^{k-1}(k-1)!}{(1+x)^k}. Pas Taylor–Young toe in a=0a = 0. (De meetkundige is exact: 11x0nxk=xn+11x=o(xn)\frac{1}{1-x} - \sum_0^n x^k = \frac{x^{n+1}}{1 - x} = o(x^n).)

Methode 16.12 (Rekenen met ontwikkelingen)

  1. Leg de doelorde nn eerst vast, en kap elk tussenresultaat daar af — hogere termen meeslepen is verspilde moeite, lagere weglaten is een fout.
  2. Sommen en producten: ontwikkel elke factor tot orde nn en vermenigvuldig, waarbij je alles voorbij xnx^n weggooit.
  3. Samenstelling f(u(x))f(u(x)) met u(x)0u(x) \to 0: substitueer de ontwikkeling van uu in die van ff, orde na orde.
  4. Quotiënten: schrijf 11+v\frac{1}{1 + v} met v0v \to 0 en gebruik de meetkundige ontwikkeling.
  5. Integreer een ontwikkeling term voor term (afleiden vraagt meer zorg — rechtvaardiging: de integraal van o(tn)o(t^n) van 00 tot xx is o(xn+1)o(x^{n+1}), door rechtstreeks te begrenzen).

Voorbeeld 16.13 (Samenstelling, met de boekhouding zichtbaar)

Ontwikkel esinx\eu^{\sin x} tot orde 33. Binnenste ontwikkeling: u=sinx=xx36+o(x3)u = \sin x = x - \frac{x^3}{6} + o(x^3), die inderdaad naar 00 nadert. Buitenste: eu=1+u+u22+u36+o(u3)\eu^u = 1 + u + \frac{u^2}{2} + \frac{u^3}{6} + o(u^3), en o(u3)=o(x3)o(u^3) = o(x^3) omdat uxu \sim x. Machten van uu, afgekapt bij x3x^3:

u2=x2+o(x3),u3=x3+o(x3)u^2 = x^2 + o(x^3), \qquad u^3 = x^3 + o(x^3)

(de kruisterm 2x(x36)2x\cdot(-\frac{x^3}{6}) is al x4x^4). Zet in elkaar:

esinx=1+(xx36)+x22+x36+o(x3)=1+x+x22+o(x3):\eu^{\sin x} = 1 + \Bigl(x - \frac{x^3}{6}\Bigr) + \frac{x^2}{2} + \frac{x^3}{6} + o(x^3) = 1 + x + \frac{x^2}{2} + o(x^3) :

de twee bijdragen in x3x^3 heffen elkaar precies op. Het afsluitende inzicht: esinx\eu^{\sin x} en ex\eu^x stemmen overeen tot orde 33 — niet omdat sinxx\sin x \approx x grofweg klopt, maar omdat het eerste meningsverschil van de exponenten (x36-\frac{x^3}{6}) binnenkomt vermenigvuldigd met e0=1\eu^0 = 1 en vervolgens wordt opgevangen door de derdegraadsterm van de buitenste exponentiële; een boekhouding orde na orde spoort zulke samenzweringen op, oogballen nooit. (De volgende term is x48-\frac{x^4}{8}: bij orde 44 eindigt het bestand.)

Voorbeeld 16.14

Ontwikkeling van tan\tan op orde 55. Schrijf tanx=sinx1cosx\tan x = \sin x \cdot \frac{1}{\cos x}:

1cosx=11(x22x424+o(x5))=1+(x22x424)+(x22) ⁣2+o(x5)=1+x22+5x424+o(x5),\frac{1}{\cos x} = \frac{1}{1 - \bigl(\frac{x^2}{2} - \frac{x^4}{24} + o(x^5)\bigr)} = 1 + \Bigl(\frac{x^2}{2} - \frac{x^4}{24}\Bigr) + \Bigl(\frac{x^2}{2}\Bigr)^{\!2} + o(x^5) = 1 + \frac{x^2}{2} + \frac{5x^4}{24} + o(x^5),

en daarna

tanx=(xx36+x5120)(1+x22+5x424)+o(x5)=x+x33+2x515+o(x5).\tan x = \Bigl(x - \frac{x^3}{6} + \frac{x^5}{120}\Bigr) \Bigl(1 + \frac{x^2}{2} + \frac{5x^4}{24}\Bigr) + o(x^5) = x + \frac{x^3}{3} + \frac{2x^5}{15} + o(x^5) .
De sinus (donker) tegenover haar Taylorveeltermen in 0: T_1 = x, T_3 = x - x3/6 (gestreept), T_5 = x - x3/6 + x5/120 (gestippeld). Elk nieuw paar termen omarmt de kromme over een zichtbaar breder venster, maar elke veelterm scheurt uiteindelijk weg: een taylorontwikkeling is een lokaal contract, scherp in 0 en zwijgzaam ver weg. De grens van Taylor–Lagrange |x|n+1(n+1)! kwantificeert het venster; de vlakke functie van  laat zien dat het contract zelfs leeg kan zijn buiten het punt zelf.
De sinus (donker) tegenover haar Taylorveeltermen in 00: T1=xT_1 = x, T3=xx36T_3 = x - \frac{x^3}{6} (gestreept), T5=xx36+x5120T_5 = x - \frac{x^3}{6} + \frac{x^5}{120} (gestippeld). Elk nieuw paar termen omarmt de kromme over een zichtbaar breder venster, maar elke veelterm scheurt uiteindelijk weg: een taylorontwikkeling is een lokaal contract, scherp in 00 en zwijgzaam ver weg. De grens van Taylor–Lagrange xn+1(n+1)!\frac{\abs{x}^{n+1}}{(n+1)!} kwantificeert het venster; de vlakke functie van Probleem 16.1 laat zien dat het contract zelfs leeg kan zijn buiten het punt zelf.

16.3 Toepassingen

Voorbeeld 16.15 (Limieten)

limx0xsinxx3:xsinx=x36+o(x3)x36,dus is de limiet 16\lim_{x \to 0} \frac{x - \sin x}{x^3}: \qquad x - \sin x = \frac{x^3}{6} + o(x^3) \sim \frac{x^3}{6}, \qquad\text{dus is de limiet } \frac16

— waarmee de vraag uit Oefening 4.9 is beslecht. Evenzo limx0(sinxx)1/x2\displaystyle\lim_{x\to0}\Bigl(\frac{\sin x}{x}\Bigr)^{1/x^2}: de logaritme is

1x2ln(1x26+o(x2))=1x2(x26+o(x2))16,limiet e1/6.\frac{1}{x^2}\ln\Bigl(1 - \frac{x^2}{6} + o(x^2)\Bigr) = \frac{1}{x^2}\Bigl(-\frac{x^2}{6} + o(x^2)\Bigr) \longrightarrow -\frac16, \qquad\text{limiet } \eu^{-1/6}.

Opmerking 16.16 (Veelgemaakte fouten met ontwikkelingen)

(i) Tel equivalenten nooit op of af: uit tanxx\tan x \sim x en sinxx\sin x \sim x mag men niet besluiten dat tanxsinx0\tan x - \sin x \sim 0 (betekenisloos) — de eerlijke weg loopt via ontwikkelingen:

tanxsinx=(x+x33)(xx36)+o(x3)=x32+o(x3)x32.\tan x - \sin x = \Bigl(x + \frac{x^3}{3}\Bigr) - \Bigl(x - \frac{x^3}{6}\Bigr) + o(x^3) = \frac{x^3}{2} + o(x^3) \sim \frac{x^3}{2} .

(ii) Ontwikkel voorbij de slachting: in dezelfde berekening ziet orde 11 alleen xx=0x - x = 0; telkens wanneer de leidende termen elkaar opheffen, verhoog je de orde tot een coëfficiënt overleeft die niet nul is, en pas dan zet je terug om naar een equivalent. (iii) Equivalenten gaan niet door exponentiëlen heen: n2+nn2n^2 + n \sim n^2, en toch is en2+n=enen2\eu^{n^2+n} = \eu^{n}\,\eu^{n^2} niet equivalent met en2\eu^{n^2} — neem alleen de exponentiële van ontwikkelingen van de exponent waarvan de fout naar 00 nadert, nooit van equivalenten van de exponent. (Logaritmen zijn veiliger: is unvn1u_n \sim v_n \to \ell \neq 1 met >0\ell > 0, dan is lnunlnvn\ln u_n \sim \ln v_n.) (iv) De rekenkunde van de o()o(\cdot) is eenrichtingsverkeer: o(x2)+o(x2)=o(x2)o(x^2) + o(x^2) = o(x^2), 5o(x2)=o(x2)5\,o(x^2) = o(x^2), xo(x2)=o(x3)x\cdot o(x^2) = o(x^3) — maar een o(x2)o(x^2) is geen bepaalde functie, dus streep er nooit twee tegen elkaar weg: o(x2)o(x2)o(x^2) - o(x^2) is o(x2)o(x^2), niet 00.

Propositie 16.17 (Lokaal gedrag)

Stel dat f(x)=f(a)+c(xa)p+o((xa)p)f(x) = f(a) + c\,(x - a)^p + o\bigl((x-a)^p\bigr) met c0c \neq 0 (de eerste term na de constante die niet nul is; p2p \geq 2 in een kritiek punt).

  • Is pp even, dan heeft ff een lokaal minimum in aa als c>0c > 0, en een lokaal maximum als c<0c < 0.
  • Is pp oneven, dan is er geen extremum (ff(a)f - f(a) wisselt van teken); begint de ontwikkeling bovendien na een lineaire term f(a)(xa)f'(a)(x - a), dan kruist de grafiek haar raaklijn: een buigpunt.

Bewijs. In de buurt van aa heeft f(x)f(a)=(xa)p(c+o(1))f(x) - f(a) = (x-a)^p\bigl(c + o(1)\bigr) het teken van c(xa)pc\,(x-a)^p: een vast teken voor pp even, wisselend voor pp oneven.

Voorbeeld 16.18 (Exponenten moeten tot op o(1)o(1) worden ontwikkeld)

Bepaal een equivalent van un=(1+1n)n2u_n = \bigl(1 + \frac1n\bigr)^{n^2}. Ontwikkel de exponent totdat zijn fout naar 00 nadert:

n2ln(1+1n)=n2(1n12n2+O(1n3))=n12+O(1n),n^2 \ln\Bigl(1 + \frac1n\Bigr) = n^2\Bigl(\frac1n - \frac{1}{2n^2} + O\Bigl(\frac{1}{n^3}\Bigr)\Bigr) = n - \frac12 + O\Bigl(\frac1n\Bigr),

dus un=en1/2eO(1/n)u_n = \eu^{\,n - 1/2}\,\eu^{O(1/n)} met eO(1/n)1\eu^{O(1/n)} \to 1:

un    en12.u_n \;\sim\; \eu^{\,n - \frac12} .

Merk op wat er met minder zorg was misgegaan: de exponent afbreken bij n21n=n+O(1)n^2\cdot\frac1n = n + O(1) laat een factor eO(1)\eu^{O(1)} staan — begrensd, maar niet naar 11 naderend — en dan kan geen enkel equivalent worden beweerd. De regel van de valkuilen hierboven, in positieve vorm: een equivalent van ean\eu^{a_n} vereist de ontwikkeling van ana_n tot en met een term die naar nul nadert, waarbij elke coëfficiënt daarvóór exact wordt bewaard.

Voorbeeld 16.19 (Een vlak kritiek punt classificeren)

Onderzoek f(x)=cosx+x22f(x) = \cos x + \frac{x^2}{2} in de buurt van 00. Zowel f(0)=0f'(0) = 0 als f(0)=cos0+1=0f''(0) = -\cos 0 + 1 = 0: de toets met de tweede afgeleide zwijgt. Ontwikkel in plaats daarvan:

f(x)=(1x22+x424+o(x4))+x22=1+x424+o(x4):f(x) = \Bigl(1 - \frac{x^2}{2} + \frac{x^4}{24} + o(x^4)\Bigr) + \frac{x^2}{2} = 1 + \frac{x^4}{24} + o(x^4) :

de eerste term cxpc\,x^p die niet nul is heeft p=4p = 4 even en c=124>0c = \frac{1}{24} > 0: een lokaal minimum, van ongewone vlakheid (de grafiek verlaat haar minimumwaarde als x4x^4, niet als x2x^2). Het afsluitende inzicht: de ontwikkeling ziet in één regel wat herhaald afleiden versluiert — en Propositie 16.17 is het systematische woordenboek van “eerste overlevende term” naar “lokale gedaante”.

Voorbeeld 16.20 (Ontwikkelingen in het oneindige)

Twee berekeningen waarin de veranderlijke naar ++\infty loopt en de substitutie h=1x0+h = \frac1x \to 0^+ het hele gereedschap importeert. Ten eerste

xx2ln(1+1x)=xx2(1x12x2+O(1x3))=12+O(1x)12.x - x^2\ln\Bigl(1 + \frac1x\Bigr) = x - x^2\Bigl(\frac1x - \frac{1}{2x^2} + O\Bigl(\frac{1}{x^3}\Bigr)\Bigr) = \frac12 + O\Bigl(\frac1x\Bigr) \longrightarrow \frac12 .

Ten tweede de boogtangens in het oneindige: uit arctanx+arctan1x=π2\arctan x + \arctan\frac1x = \frac\pi2 voor x>0x > 0 (Propositie 4.12) en de ontwikkeling van arctan\arctan in 00 (Oefening 16.3) volgt

arctanx=π2arctan1x=π21x+13x3+o(1x3):\arctan x = \frac\pi2 - \arctan\frac1x = \frac\pi2 - \frac1x + \frac{1}{3x^3} + o\Bigl(\frac{1}{x^3}\Bigr) :

de grafiek nadert haar asymptoot y=π2y = \frac\pi2 van onderen met snelheid 1x\frac1x. Het afsluitende inzicht: er is geen afzonderlijke theorie van ontwikkelingen in het oneindige — één omgekeerde substitutie herleidt ze tot ontwikkelingen in 00, mits elke tussenliggende OO en oo eerlijk wordt meegedragen.

Voorbeeld 16.21 (Asymptoot via ontwikkeling)

Als x+x \to +\infty geldt

x2+x=x1+1x=x(1+12x18x2+o(1x2))=x+1218x+o(1x):\sqrt{x^2 + x} = x\sqrt{1 + \tfrac1x} = x\Bigl(1 + \frac{1}{2x} - \frac{1}{8x^2} + o\bigl(\tfrac{1}{x^2}\bigr)\Bigr) = x + \frac12 - \frac{1}{8x} + o\bigl(\tfrac 1x\bigr):

de rechte y=x+12y = x + \frac12 is een asymptoot, benaderd van onderen (de volgende term 18x-\frac{1}{8x} is negatief).

Opmerking 16.22 (Waar ontwikkelingen hierna werken)

Asymptotische ontwikkelingen zijn de vaste taal van de rest van het boek: in Hoofdstuk 17 beslissen zij over convergentie (equivalenten voeden de vergelijkingscriteria, en de studie van 1nα\sum \frac{1}{n^\alpha} is een vermomde ontwikkeling); in het volume van bachelorjaar 2 worden zij machtreeksen, waarin de taylorveelterm oneindig veel termen en een convergentiestraal krijgt; en elke linearisatie in de fysica — de slinger, storing van eerste orde — is een uitspraak van Taylor–Young waarbij de o()o(\cdot) stilzwijgend is weggelaten. De ene waarschuwing die het graveren waard is: een ontwikkeling beschrijft een functie alleen in de buurt van één punt — zie de vlakke functie van de weekendopgave, waarvan de ontwikkeling in 00 identiek nul is zonder dat de functie dat is.

Opmerking 16.23 (Vooruitzichten binnen dit volume)

Ontwikkelingen zijn de werktaal van de resterende analyse en van de meetkunde die komt. Hoofdstuk 17 zet ze om in convergentieoordelen: een equivalent van de algemene term is een ontwikkeling afgekapt bij haar eerste term, en de fijnere criteria (alternerend met foutbeheersing) verbruiken ook de tweede term. Hoofdstuk 24 leest de lokale meetkunde af uit ontwikkelingen van de twee coördinaatfuncties: of een geparametriseerde kromme in een punt kruist, kust of een keerpunt maakt, wordt beslist door welke machten van tt overleven in x(t)x(t) en y(t)y(t) — de vlakke versie van Propositie 16.17. En Hoofdstuk 25 stopt met opzet bij orde één: het raakvlak is een uitspraak van Taylor–Young in twee veranderlijken, terwijl de volledige theorie van de tweede orde (hessianen, zadelpunten) wordt uitgesteld tot het volume van bachelorjaar 2. De rode draad: elke “lokale” vraag in dit boek wordt beantwoord door de eerste overlevende term van een ontwikkeling op te schrijven.

16.4 Oefeningen

Oefening 16.1

Geef de ontwikkelingen in 00: e2x\eu^{2x} tot orde 33;   ln(1x)\;\ln(1 - x) tot orde 44;   1+x\;\sqrt{1 + x} tot orde 33;   11+x2\;\dfrac{1}{1 + x^2} tot orde 66.

Oplossing

Oplossing van Oefening 16.1.

e2x=1+2x+2x2+4x33+o(x3);ln(1x)=xx22x33x44+o(x4);\eu^{2x} = 1 + 2x + 2x^2 + \frac{4x^3}{3} + o(x^3); \qquad \ln(1 - x) = -x - \frac{x^2}{2} - \frac{x^3}{3} - \frac{x^4}{4} + o(x^4);
1+x=1+x2x28+x316+o(x3);11+x2=1x2+x4x6+o(x6),\sqrt{1+x} = 1 + \frac x2 - \frac{x^2}{8} + \frac{x^3}{16} + o(x^3); \qquad \frac{1}{1 + x^2} = 1 - x^2 + x^4 - x^6 + o(x^6),

de laatste door x2-x^2 in de meetkundige ontwikkeling te substitueren.

Oefening 16.2

Bereken de limieten:

limx0ex1xx2,limx0cosx1x2x4,limx0ln(1+x)sinxx2.\lim_{x\to 0} \frac{\eu^x - 1 - x}{x^2}, \qquad \lim_{x\to 0} \frac{\cos x - \sqrt{1 - x^2}}{x^4}, \qquad \lim_{x\to 0} \frac{\ln(1+x) - \sin x}{x^2}.
Oplossing

Oplossing van Oefening 16.2.

ex1x=x22+o(x2)\eu^x - 1 - x = \frac{x^2}{2} + o(x^2): limiet 12\dfrac12.

cosx=1x22+x424+o(x4)\cos x = 1 - \frac{x^2}{2} + \frac{x^4}{24} + o(x^4) en 1x2=1x22x48+o(x4)\sqrt{1 - x^2} = 1 - \frac{x^2}{2} - \frac{x^4}{8} + o(x^4): verschil x424+x48=x46+o(x4)\frac{x^4}{24} + \frac{x^4}{8} = \frac{x^4}{6} + o(x^4): limiet 16\dfrac16.

ln(1+x)sinx=(xx22)x+o(x2)=x22+o(x2)\ln(1+x) - \sin x = \bigl(x - \frac{x^2}{2}\bigr) - x + o(x^2) = -\frac{x^2}{2} + o(x^2): limiet 12-\dfrac12.

Oefening 16.3

Ontwikkel arctanx\arctan x in 00 tot orde 55 door de ontwikkeling van 11+x2\frac{1}{1 + x^2} te integreren, en arcsinx\arcsin x tot orde 55 door die van (1x2)1/2(1 - x^2)^{-1/2} te integreren.

Oplossing

Oplossing van Oefening 16.3.

11+t2=1t2+t4+o(t5)\frac{1}{1+t^2} = 1 - t^2 + t^4 + o(t^5); door van 00 tot xx te integreren (Methode 16.12 (5)):

arctanx=xx33+x55+o(x5)  (zelfs o(x6), wegens de oneven pariteit).\arctan x = x - \frac{x^3}{3} + \frac{x^5}{5} + o(x^5)\ \ (\text{zelfs }o(x^6)\text{, wegens de oneven pariteit}).

(1t2)1/2=1+t22+38t4+o(t5)(1 - t^2)^{-1/2} = 1 + \frac{t^2}{2} + \frac38 t^4 + o(t^5) (binomiaalontwikkeling met α=12\alpha = -\frac12, x=t2x = -t^2: (1/22)=(12)(32)2=38\binom{-1/2}{2} = \frac{(-\frac12)(-\frac32)}{2} = \frac38); na integratie:

arcsinx=x+x36+3x540+o(x5).\arcsin x = x + \frac{x^3}{6} + \frac{3x^5}{40} + o(x^5) .

Oefening 16.4

Bewijs met Taylor–Lagrange voor exp\exp op [0,1]\intcc{0}{1} dat

ek=0n1k!3(n+1)!,\Bigl| \eu - \sum_{k=0}^{n} \frac{1}{k!} \Bigr| \leq \frac{3}{(n+1)!},

en bepaal een nn die 66 exacte decimalen van e\eu garandeert.

Oplossing

Oplossing van Oefening 16.4.

Taylor–Lagrange (Stelling 16.7) voor exp\exp in a=0a = 0, x=1x = 1: de (n+1)(n+1)-de afgeleide is ete<3\eu^t \leq \eu < 3 op [0,1]\intcc{0}{1}, dus

ek=0n1k!3(n+1)!.\Bigl|\eu - \sum_{k=0}^{n} \frac{1}{k!}\Bigr| \leq \frac{3}{(n+1)!} .

Voor 66 exacte decimalen willen we 3(n+1)!<5×107\frac{3}{(n+1)!} < 5\times 10^{-7}, dat wil zeggen (n+1)!>6×106(n+1)! > 6\times 10^{6}: omdat 10!=362880010! = 3\,628\,800 en 11!=3991680011! = 39\,916\,800, volstaat n+1=11n + 1 = 11, dus n=10n = 10.

Oefening 16.5 ★★

Ontwikkel tot orde 22 in 1n\frac1n en leid de limiet en de convergentiesnelheid af:

(1+1n) ⁣n=e(112n+1124n2+o(1n2)).\Bigl(1 + \frac 1n\Bigr)^{\!n} = \eu\Bigl(1 - \frac{1}{2n} + \frac{11}{24n^2} + o\Bigl(\frac{1}{n^2}\Bigr)\Bigr).
Oplossing

Oplossing van Oefening 16.5.

nln(1+1n)=n(1n12n2+13n3+o(1n3))=112n+13n2+o(1n2)n\ln\bigl(1 + \frac1n\bigr) = n\Bigl(\frac1n - \frac{1}{2n^2} + \frac{1}{3n^3} + o\bigl(\frac{1}{n^3}\bigr)\Bigr) = 1 - \frac{1}{2n} + \frac{1}{3n^2} + o\bigl(\frac{1}{n^2}\bigr). Neem de exponentiële, met u=12n+13n2u = -\frac{1}{2n} + \frac{1}{3n^2} en eu=1+u+u22+o(u2)\eu^u = 1 + u + \frac{u^2}2 + o(u^2):

(1+1n)n=eeu=e(112n+13n2+18n2+o(1n2))=e(112n+1124n2+o(1n2)).\Bigl(1 + \frac1n\Bigr)^n = \eu\cdot \eu^{u} = \eu\Bigl(1 - \frac{1}{2n} + \frac{1}{3n^2} + \frac{1}{8n^2} + o\Bigl(\frac{1}{n^2}\Bigr)\Bigr) = \eu\Bigl(1 - \frac{1}{2n} + \frac{11}{24n^2} + o\Bigl(\frac{1}{n^2}\Bigr)\Bigr).

Limiet e\eu; de fout is e2n\sim \dfrac{\eu}{2n}: traag (één cijfer per vertienvoudiging van nn).

Oefening 16.6 ★★

Onderzoek het lokale gedrag in 00 van f(x)=x2x4f(x) = x^2 - x^4 en van g(x)=x3+x5g(x) = x^3 + x^5; en bepaal de ligging van de grafiek van h(x)=exh(x) = \eu^x ten opzichte van haar raaklijn in a=1a = 1, eerst lokaal en daarna globaal.

Oplossing

Oplossing van Oefening 16.6.

f(x)=x2x4=x2(1+o(1))f(x) = x^2 - x^4 = x^2(1 + o(1)): eerste term x2x^2, p=2p = 2 even, coëfficiënt >0> 0: lokaal minimum in 00 (niet globaal: f(2)=12f(2) = -12).

g(x)=x3+x5g(x) = x^3 + x^5: eerste term x3x^3, pp oneven: geen extremum; gg kruist haar (horizontale) raaklijn: buigpunt in 00.

h=exph = \exp in a=1a = 1: h(x)=e+e(x1)+e2(x1)2+o((x1)2)h(x) = \eu + \eu(x-1) + \frac{\eu}{2}(x-1)^2 + o((x-1)^2); het verschil met de raaklijn is e2(x1)2+o()>0\frac{\eu}{2}(x-1)^2 + o(\cdot) > 0 lokaal. Globaal: exex0\eu^x - \eu x \geq 0 voor alle xx wegens de convexiteit (Stelling 14.19 (3)): de grafiek ligt boven elke raaklijn, met gelijkheid alleen in het raakpunt.

Oefening 16.7 ★★

Bepaal de asymptoten in ±\pm\infty van f(x)=x3+x23f(x) = \sqrt[3]{x^3 + x^2} en de ligging van de kromme ten opzichte daarvan.

Oplossing

Oplossing van Oefening 16.7.

Voor x+x \to +\infty:

f(x)=x(1+1x)1/3=x(1+13x19x2+o(1x2))=x+1319x+o(1x):f(x) = x\Bigl(1 + \frac1x\Bigr)^{1/3} = x\Bigl(1 + \frac{1}{3x} - \frac{1}{9x^2} + o\Bigl(\frac{1}{x^2}\Bigr)\Bigr) = x + \frac13 - \frac{1}{9x} + o\Bigl(\frac1x\Bigr):

asymptoot y=x+13y = x + \frac13, met de kromme eronder bij ++\infty. Als xx \to -\infty blijft dezelfde berekening geldig (de derdemachtswortel is voor alle reële getallen gedefinieerd, en 1x0\frac1x \to 0): dezelfde asymptoot y=x+13y = x + \frac13, maar nu is 19x>0-\frac{1}{9x} > 0: de kromme ligt boven de rechte.

Oefening 16.8 ★★

Bepaal het equivalent, als nn \to \infty, van

un=n+1n,vn=ln(n+1)lnn,wn=sin1ntan1n,u_n = \sqrt{n+1} - \sqrt n, \qquad v_n = \ln(n+1) - \ln n, \qquad w_n = \sin\frac{1}{n} - \tan\frac{1}{n},

telkens als een macht van nn maal een constante.

Oplossing

Oplossing van Oefening 16.8.

un=n(1+1n1)=n(12n+o(1n))12nu_n = \sqrt n\bigl(\sqrt{1 + \tfrac1n} - 1\bigr) = \sqrt n\bigl(\frac{1}{2n} + o(\frac1n)\bigr) \sim \dfrac{1}{2\sqrt n}.

vn=ln(1+1n)1nv_n = \ln\bigl(1 + \frac1n\bigr) \sim \dfrac 1n.

wnw_n: met h=1n0h = \frac1n \to 0 is sinhtanh=(hh36)(h+h33)+o(h3)=h32+o(h3)\sin h - \tan h = \bigl(h - \frac{h^3}{6}\bigr) - \bigl(h + \frac{h^3}{3}\bigr) + o(h^3) = -\frac{h^3}{2} + o(h^3), dus wn12n3w_n \sim -\dfrac{1}{2n^3}.

Oefening 16.9 ★★★

Zij ff van klasse C2C^2 op R\R. Bewijs dat voor elke xx en elke h>0h > 0

f(x)f(x+h)f(xh)2h+h2sup[xh,x+h]f,\abs{f'(x)} \leq \frac{\abs{f(x+h) - f(x-h)}}{2h} + \frac{h}{2}\sup_{\intcc{x-h}{x+h}}\abs{f''} ,

en leid de ongelijkheid van het type Landau–Kolmogorov af: is fM0\abs f \leq M_0 en fM2\abs{f''} \leq M_2 op R\R, dan is f2M0M2\abs{f'} \leq \sqrt{2 M_0 M_2} overal. (Optimaliseer over hh.)

Oplossing

Oplossing van Oefening 16.9.

Taylor–Lagrange op orde 11 rond xx, aan beide zijden:

f(x+h)=f(x)+hf(x)+R+,f(xh)=f(x)hf(x)+R,R±h22supf.f(x + h) = f(x) + h f'(x) + R_+,\quad f(x - h) = f(x) - h f'(x) + R_-, \qquad \abs{R_\pm} \leq \frac{h^2}{2} \sup \abs{f''} .

Aftrekken geeft f(x+h)f(xh)=2hf(x)+(R+R)f(x+h) - f(x-h) = 2h f'(x) + (R_+ - R_-), dus

f(x)f(x+h)f(xh)2h+h2sup[xh,x+h]f.\abs{f'(x)} \leq \frac{\abs{f(x+h) - f(x-h)}}{2h} + \frac{h}{2}\sup_{\intcc{x-h}{x+h}}\abs{f''} .

Met globale grenzen: f(x)M0h+M2h2\abs{f'(x)} \leq \frac{M_0}{h} + \frac{M_2 h}{2} voor elke h>0h > 0. Het rechterlid wordt geminimaliseerd bij h=2M0/M2h = \sqrt{2M_0/M_2} (afgeleide nul), met waarde 2M0M2\sqrt{2M_0M_2} — bijgevolg is f2M0M2\abs{f'} \leq \sqrt{2M_0M_2}. (Is M2=0M_2 = 0, laat dan hh \to \infty: f=0f' = 0, consistent.)

Oefening 16.10 ★★★

De rij u0(0,π)u_0 \in \intoo{0}{\pi}, un+1=sinunu_{n+1} = \sin u_n daalt naar 00 (verantwoord dit kort). Om haar snelheid te bepalen, beschouw vn=1un2v_n = \frac{1}{u_n^2}:

  1. bewijs met de ontwikkeling van sin\sin dat vn+1vn13v_{n+1} - v_n \to \frac13;
  2. leid met Cesàro (Oefening 11.10) af dat vnn13\frac{v_n}{n} \to \frac13, en daarna het equivalent un3nu_n \sim \sqrt{\dfrac{3}{n}}.
Oplossing

Oplossing van Oefening 16.10.

Op (0,π)\intoo{0}{\pi} is 0<sinu<u0 < \sin u < u, dus is (un)(u_n) strikt dalend en positief, en dus convergent; de limiet is een vast punt van sin\sin in [0,π]\intcc{0}{\pi}, en sin=\sin \ell = \ell dwingt =0\ell = 0 af (want sinx<x\sin x < x voor x>0x > 0).

  1. Met sinu=uu36+o(u3)\sin u = u - \frac{u^3}{6} + o(u^3) als u0u \to 0:

    vn+1vn=1sin2un1un2=1un2((1un26+o(un2)) ⁣21)=1un2(un23+o(un2))13.v_{n+1} - v_n = \frac{1}{\sin^2 u_n} - \frac{1}{u_n^2} = \frac{1}{u_n^2}\Bigl(\Bigl(1 - \frac{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 .
  2. Volgens Oefening 11.10 (3) (Cesàro voor verschillen) is vnn13\frac{v_n}{n} \to \frac13, dat wil zeggen vnn3v_n \sim \frac n3, dus un23nu_n^2 \sim \frac 3n: omdat un>0u_n > 0,

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

Oefening 16.11 ★★

(Verschillen van oneindigheden) Bereken

limx0(1x21sin2x)enlimx0+(1x1ln(1+x))\lim_{x \to 0} \Bigl(\frac{1}{x^2} - \frac{1}{\sin^2 x}\Bigr) \qquad\text{en}\qquad \lim_{x \to 0^+} \Bigl(\frac 1x - \frac{1}{\ln(1 + x)}\Bigr)

door op één noemer te brengen en teller en noemer afzonderlijk te ontwikkelen.

Oplossing

Oplossing van Oefening 16.11.

Breng op één noemer. Eerste limiet:

1x21sin2x=sin2xx2x2sin2x,sin2x=(xx36+o(x4)) ⁣2=x2x43+o(x5),\frac{1}{x^2} - \frac{1}{\sin^2 x} = \frac{\sin^2 x - x^2}{x^2\sin^2 x}, \qquad \sin^2 x = \Bigl(x - \frac{x^3}{6} + o(x^4)\Bigr)^{\!2} = x^2 - \frac{x^4}{3} + o(x^5) ,

zodat de teller x43+o(x4)-\frac{x^4}{3} + o(x^4) is terwijl de noemer x4\sim x^4 is: de limiet is 13-\dfrac13.

Tweede: 1x1ln(1+x)=ln(1+x)xxln(1+x)\dfrac1x - \dfrac{1}{\ln(1+x)} = \dfrac{\ln(1+x) - x}{x\ln(1+x)}; de teller is x22+o(x2)-\frac{x^2}{2} + o(x^2), de noemer x(x+o(x))x2x\bigl(x + o(x)\bigr) \sim x^2: de limiet is 12-\dfrac12.

Oefening 16.12 ★★★

(Asymptotiek van impliciete wortels) Toon aan dat voor elke kNk \in \N^* de vergelijking tanx=x\tan x = x precies één oplossing xkx_k heeft in (kππ2,kπ+π2)\intoo{k\pi - \frac\pi2}{k\pi + \frac\pi2}, dat xk=kπ+π2εkx_k = k\pi + \frac\pi2 - \varepsilon_k met εk=arctan1xk\varepsilon_k = \arctan\frac{1}{x_k}, en leid de ontwikkeling

xk=kπ+π21kπ+o(1k)(k)x_k = k\pi + \frac\pi2 - \frac{1}{k\pi} + o\Bigl(\frac 1k\Bigr) \qquad (k \to \infty)

af.

Oplossing

Oplossing van Oefening 16.12.

Op Ik=(kππ2,kπ+π2)I_k = \intoo{k\pi - \frac\pi2}{k\pi + \frac\pi2} heeft de functie g(x)=tanxxg(x) = \tan x - x als afgeleide tan2x0\tan^2 x \geq 0, die alleen in het enkele punt kπk\pi nul wordt: gg is strikt stijgend op IkI_k (Gevolg 14.12 (2)), met limieten -\infty en ++\infty in de uiteinden: precies één nulpunt xkx_k. Voor k1k \geq 1 is g(kπ)=kπ<0g(k\pi) = -k\pi < 0, dus xk(kπ,kπ+π2)x_k \in \intoo{k\pi}{k\pi + \frac\pi2}: schrijf xk=kπ+π2εkx_k = k\pi + \frac\pi2 - \varepsilon_k met εk(0,π2)\varepsilon_k \in \intoo{0}{\frac\pi2}. Dan is

xk=tanxk=tan(π2εk)=1tanεkεk=arctan1xk,x_k = \tan x_k = \tan\Bigl(\frac\pi2 - \varepsilon_k\Bigr) = \frac{1}{\tan\varepsilon_k} \quad\Longrightarrow\quad \varepsilon_k = \arctan\frac{1}{x_k} ,

met gebruik van tanεk=1xk\tan\varepsilon_k = \frac{1}{x_k} en εk(0,π2)\varepsilon_k \in \intoo{0}{\frac\pi2}. Omdat xkkπx_k \geq k\pi \to \infty, is εk0\varepsilon_k \to 0, en

εk=arctan1xk=1xk+O(1xk3)=1kπ+O(1)+O(1k3)=1kπ+O(1k2),\varepsilon_k = \arctan\frac{1}{x_k} = \frac{1}{x_k} + O\Bigl(\frac{1}{x_k^3}\Bigr) = \frac{1}{k\pi + O(1)} + O\Bigl(\frac{1}{k^3}\Bigr) = \frac{1}{k\pi} + O\Bigl(\frac{1}{k^2}\Bigr) ,

waaruit xk=kπ+π21kπ+o(1k)x_k = k\pi + \frac\pi2 - \frac{1}{k\pi} + o\bigl(\frac1k\bigr).

16.5 Opgave: alternerende sommen, gecertificeerde cijfers en de irrationaliteit van cos1\cos 1

Probleem 16.1

Weekendopgave — de alternerende schatting SSnan+1\abs{S - S_n} \leq a_{n+1}: ln2\ln 2 en π\pi met bewezen decimalen, de formule van Machin, en cos1Q\cos 1 \notin \Q

Een alternerende som met dalende termen is het vriendelijkste object van de numerieke analyse: haar fout wordt begrensd door de eerste weggelaten term, met bekend teken. Deze opgave bewijst dat beginsel met de stelling over ingesloten rijen, en besteedt het vervolgens op drie manieren: gecertificeerde decimalen voor ln2\ln 2 (drie concurrerende wegen) en voor π\pi (Leibniz, en daarna de formule van Machin uit 1706, eeuwenlang nog steeds het idee achter recordberekeningen), de irrationaliteit van cos1\cos 1, sin1\sin 1 en cosh1\cosh 1, en, als tegenwicht, de gelijkheidsvorm van Taylor–Lagrange en de vlakke functie waarvan de taylorontwikkeling liegt. Overal is de taal van “reeksen” informeel: elke som hier is een rij van partiële sommen, zoals in Voorbeeld 11.12; de eigenlijke theorie opent in Hoofdstuk 17.

Deel I — De alternerende schatting. Zij (ak)k0(a_k)_{k \geq 0} dalend naar 00 en Sn=k=0n(1)kakS_n = \sum_{k=0}^{n} (-1)^k a_k.

  1. Toon aan dat (S2n+1)(S_{2n+1}) stijgend is, (S2n)(S_{2n}) dalend, en dat zij ingesloten zijn (Stelling 11.11): beide convergeren naar een gemeenschappelijke SS met, voor elke nn,

    S2n+1SS2n,SSnan+1,S_{2n+1} \leq S \leq S_{2n}, \qquad \abs{S - S_n} \leq a_{n+1} ,

    waarbij de fout het teken heeft van de eerste weggelaten term. Toon bovendien aan dat al deze ongelijkheden strikt zijn wanneer de daling strikt is.

  2. Eerste dividend: vergelijk voor x=1x = 1 in de exponentiële reeks met Stelling 16.7 in a=0a = 0: toon aan dat Tn=k=0n(1)kk!T_n = \sum_{k=0}^{n} \frac{(-1)^k}{k!} convergeert naar e1\eu^{-1} met e1Tn<1(n+1)!\abs{\eu^{-1} - T_n} < \frac{1}{(n+1)!}.
  3. (Leibniz, 1674) Bewijs uit de exacte eindige identiteit

    11+t2=k=0n(1)kt2k+(1)n+1t2n+21+t2,\frac{1}{1 + t^2} = \sum_{k=0}^{n} (-1)^k t^{2k} + \frac{(-1)^{n+1} t^{2n+2}}{1 + t^2} ,

    geïntegreerd over [0,1]\intcc{0}{1}, dat

    π4=k=0n(1)k2k+1+ρn,ρn12n+3.\frac\pi4 = \sum_{k=0}^{n} \frac{(-1)^k}{2k+1} + \rho_n, \qquad \abs{\rho_n} \leq \frac{1}{2n+3} .
  4. Traagheid: hoeveel termen van Leibniz garanderen zes exacte decimalen van π\pi? (Ongeveer twee miljoen.) Evalueer 4S4=4(113+1517+19)4 S_4 = 4\bigl(1 - \frac13 + \frac15 - \frac17 + \frac19\bigr) en zijn afstand tot π\pi, om de pijn te voelen.

Deel II — ln2\ln 2 langs drie wegen.

  1. (Weg 1: de alternerende harmonische) Uit 11+t=k=0n1(1)ktk+(1)ntn1+t\frac{1}{1+t} = \sum_{k=0}^{n-1}(-1)^k t^k + \frac{(-1)^n t^n}{1+t}, geïntegreerd over [0,1]\intcc{0}{1}:

    ln2=k=1n(1)k1k+(1)nRn,12(n+1)Rn1n+1:\ln 2 = \sum_{k=1}^{n} \frac{(-1)^{k-1}}{k} + (-1)^n R_n, \qquad \frac{1}{2(n+1)} \leq R_n \leq \frac{1}{n+1} :

    de fout is van exacte orde 1n\frac 1n — een miljoen termen voor zes decimalen.

  2. (Weg 2: de snelle reeks) Integreer 11t2=k=0nt2k+t2n+21t2\frac{1}{1 - t^2} = \sum_{k=0}^{n} t^{2k} + \frac{t^{2n+2}}{1-t^2} van 00 tot x(0,1)x \in \intoo{0}{1} en evalueer in x=13x = \frac13 (merk op dat 1+1/311/3=2\frac{1 + 1/3}{1 - 1/3} = 2):

    ln2=2k=0n(1/3)2k+12k+1+ρ~n,0<ρ~n94(1/3)2n+32n+3:\ln 2 = 2\sum_{k=0}^{n} \frac{(1/3)^{2k+1}}{2k+1} + \tilde\rho_n, \qquad 0 < \tilde\rho_n \leq \frac{9}{4}\cdot \frac{(1/3)^{2n+3}}{2n+3} :

    meetkundige convergentie, ruwweg één cijfer per term.

  3. (Weg 3: Riemannsommen en een verborgen identiteit) Bewijs met inductie de identiteit

    k=12n(1)k1k=H2nHn=k=1n1n+k,\sum_{k=1}^{2n} \frac{(-1)^{k-1}}{k} = H_{2n} - H_n = \sum_{k=1}^{n} \frac{1}{n + k} ,

    en vind ln2\ln 2 terug als de limiet van de Riemannsom uit Voorbeeld 15.21: de wegen 1 en 3 zijn heimelijk hetzelfde getal, twee keer gezien.

  4. Krachtmeting bij zes termen: vergelijk k=16(1)k1k=0.6167\sum_{k=1}^{6} \frac{(-1)^{k-1}}{k} = 0.6167 met weg 2 bij n=5n = 5, die al ln2=0.693147\ln 2 = 0.693147 geeft met een fout 1.1107\leq 1.1\cdot10^{-7}. Verklaar de structurele reden (een evaluatiepunt diep binnen het convergentie-interval tegenover een punt op de rand ervan).
  5. Hoeveel termen van weg 2 certificeren tien decimalen van ln2\ln 2? Toon aan dat n=10n = 10 volstaat.

Deel III — De formule van Machin.

  1. Bereken (5+i)4(5 + \iu)^4 en ga de complexe identiteit

    (5+i)4=2(1+i)(239+i)(5 + \iu)^4 = 2\,(1 + \iu)\,(239 + \iu)

    na. Leid, door argumenten te nemen (met de afspraken van Hoofdstuk 3), de formule van Machin

    π4=4arctan15arctan1239\frac\pi4 = 4\arctan\frac15 - \arctan\frac{1}{239}

    af. (Controleer dat geen enkel argument (π2,π2)\intoo{-\frac\pi2}{\frac\pi2} verlaat.)

  2. Stel, zoals in vraag 3, voor 0<x<10 < x < 1 vast dat

    arctanx=k=0n(1)kx2k+12k+1+rn(x),rn(x)x2n+32n+3.\arctan x = \sum_{k=0}^{n} \frac{(-1)^k x^{2k+1}}{2k+1} + r_n(x), \qquad \abs{r_n(x)} \leq \frac{x^{2n+3}}{2n+3} .
  3. Certificeer π\pi tot op zeven decimalen met zes termen: begrens de totale fout van

    π16k=04(1)k(1/5)2k+12k+14(1239132393)\pi \approx 16\sum_{k=0}^{4} \frac{(-1)^k (1/5)^{2k+1}}{2k+1} - 4\Bigl(\frac{1}{239} - \frac{1}{3\cdot239^3}\Bigr)

    door 16(1/5)1111+4(1/239)55<510816\,\frac{(1/5)^{11}}{11} + 4\,\frac{(1/239)^5}{5} < 5\cdot10^{-8}, en geef de resulterende waarde 3.14159263.1415926\dots

  4. Vergelijk de drie wegen naar π\pi die nu beschikbaar zijn — Leibniz (vraag 4), de integralen van Dalzell uit Probleem 15.1 (fout 415m4^{1-5m}), Machin (fout 165(2n+3)\approx 16\cdot 5^{-(2n+3)}) — in cijfers per term, en leg uit waarom het verkleinen van het evaluatiepunt alles verslaat.

Deel IV — De gehele-getallenval, alternerende uitvoering.

  1. Stel dat cos1=pq\cos 1 = \frac pq. Vermenigvuldig de strikte alternerende insluiting van k(1)k(2k)!\sum_{k} \frac{(-1)^k}{(2k)!} (vraag 1) met (2n)!(2n)! waarbij 2nq2n \geq q, en leid een tegenspraak af: cos1\cos 1 is irrationaal.
  2. Pas dit aan voor sin1=k(1)k(2k+1)!\sin 1 = \sum_k \frac{(-1)^k}{(2k+1)!} (vermenigvuldig met (2n+1)!(2n+1)!): sin1Q\sin 1 \notin \Q. Beide irrationaal, en toch is cos21+sin21=1\cos^2 1 + \sin^2 1 = 1: irrationaliteit is niet bestand tegen algebra.
  3. De niet-alternerende neef: cosh1=k1(2k)!\cosh 1 = \sum_k \frac{1}{(2k)!} (in de zin van partiële sommen, met de tweezijdige staartgrens 0<cosh1kn1(2k)!<2(2n+2)!0 < \cosh 1 - \sum_{k \leq n} \frac{1}{(2k)!} < \frac{2}{(2n+2)!}, nog te bewijzen). Besluit met dezelfde val dat cosh1Q\cosh 1 \notin \Q.
  4. Duw door tot cos1m\cos\frac 1m voor elk geheel getal m1m \geq 1: vermenigvuldig met m2n(2n)!m^{2n}(2n)! en besluit dat cos1mQ\cos\frac1m \notin \Q. Waar breekt dezelfde poging voor een algemene cosab\cos\frac ab met b>1b > 1? (Wijs de noemer aan die niet langer wegvalt.)

Deel V — Scherper en donkerder: de gelijkheidsvorm, en een functie die Taylor bedriegt.

  1. (Taylor–Lagrange, gelijkheidsvorm) Zij ff n+1n + 1 keer afleidbaar tussen aa en xx. Definieer g(t)=f(x)k=0nf(k)(t)k!(xt)kA(xt)n+1(n+1)!g(t) = f(x) - \sum_{k=0}^{n} \frac{f^{(k)}(t)}{k!}(x - t)^k - A\,\frac{(x-t)^{n+1}}{(n+1)!} met de constante AA zo gekozen dat g(a)=0g(a) = 0. Bereken gg' (de som telescopeert), pas Rolle toe op [a,x]\intcc{a}{x}, en besluit dat er een cc strikt tussen aa en xx bestaat met

    f(x)=k=0nf(k)(a)k!(xa)k+f(n+1)(c)(n+1)!(xa)n+1.f(x) = \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}(x-a)^k + \frac{f^{(n+1)}(c)}{(n+1)!}\,(x-a)^{n+1} .
  2. Dividend van de gelijkheid: toon voor x>0x > 0 aan dat

    ex>1+x+x22!++xnn!\eu^x > 1 + x + \frac{x^2}{2!} + \dots + \frac{x^n}{n!}

    (strikt, voor elke nn), en lokaliseer waar de ongelijkheid voor x<0x < 0 omklapt naargelang de pariteit van nn.

  3. IJk de resttermen op sin(0.5)\sin(0.5) bij orde 33: Young geeft alleen o(x3)o(x^3) (geen getal); Lagrange geeft sin0.5(0.50.536)0.55120=2.61104\abs{\sin 0.5 - (0.5 - \frac{0.5^3}{6})} \leq \frac{0.5^5}{120} = 2.61\cdot10^{-4}; de alternerende schatting geeft dezelfde grens plus de tekeninformatie sin0.5>0.50.536\sin 0.5 > 0.5 - \frac{0.5^3}{6}. Vergelijk met de werkelijke fout 2.591042.59\cdot10^{-4}: de grens wordt bijna bereikt. Naar welk gereedschap zou je grijpen, en wanneer?
  4. (De vlakke functie) Zij f(x)=e1/x2f(x) = \eu^{-1/x^2} voor x0x \neq 0 en f(0)=0f(0) = 0. Toon aan dat ff continu is in 00, dat f(0)=0f'(0) = 0, en algemener — door te bewijzen dat elke afgeleide de vorm f(k)(x)=Pk(1x)e1/x2f^{(k)}(x) = P_k\bigl(\frac1x\bigr)\eu^{-1/x^2} heeft voor een veelterm PkP_k (inductie) — dat f(k)(0)=0f^{(k)}(0) = 0 voor alle kk (groeivergelijking, Propositie 4.6). Besluit: alle taylorveeltermen van ff in 00 zijn nul, en toch is f(x)>0f(x) > 0 voor x0x \neq 0: Taylor–Young geldt in elke orde en zegt niets over ff buiten 00. Ontwikkelingen beschrijven kiemen, geen functies.

Deel VI — Synthese.

  1. Laat de val nog één keer draaien, op e1=k(1)kk!\eu^{-1} = \sum_k \frac{(-1)^k}{k!}: vermenigvuldig de strikte alternerende insluiting met n!n! en besluit dat e1Q\eu^{-1} \notin \Q, en dus eQ\eu \notin \Q — het derde bewijs van dit feit in het volume. Som de drie op (ingesloten rijen, Oefening 11.9; integralen, Probleem 15.1; alternerende sommen, hier) en wat elk ervan nodig had.
  2. De kleine lettertjes: de monotonie is geen versiering. Zij bk=1kb_k = \frac1k voor oneven kk en bk=1k2b_k = \frac{1}{k^2} voor even kk: de bkb_k zijn positief en naderen tot 00, en toch divergeren de partiële sommen van (1)kbk\sum (-1)^k b_k naar -\infty. Bewijs dit (splits de partiële som in het even deel, begrensd via Voorbeeld 11.22, en het oneven deel, dat de helft van de harmonische reeks domineert, Oefening 11.5), en zeg precies welke stap van vraag 1 de monotonie gebruikte.
  3. Ga de eenvoudiger identiteit van Euler arctan12+arctan13=π4\arctan\frac12 + \arctan\frac13 = \frac\pi4 na via (2+i)(3+i)=5(1+i)(2 + \iu)(3 + \iu) = 5(1 + \iu), schat het aantal termen dat langs deze weg nodig is voor zes decimalen van π\pi (n=10n = 10 volstaat), en plaats haar tussen Leibniz en Machin in de rangschikking van vraag 13.
  4. Synthese, telkens één zin: (i) formuleer de alternerende schatting en haar twee opbrengsten (grens en teken); (ii) waarom exacte eindige identiteiten met expliciete resttermen het winnen van limietuitspraken voor gecertificeerd rekenwerk; (iii) de inventaris van de opgave (π\pi tot op 10710^{-7} met de hand, ln2\ln 2 tot op tien decimalen, vier irrationaliteitsbewijzen, één gelijkheidsstelling, één waarschuwend voorbeeld); (iv) welke van deze draden Hoofdstuk 17 zal oppakken (het criterium voor alternerende reeksen, absolute tegenover voorwaardelijke convergentie, en het herschikkingsdrama van haar weekendopgave).
Oplossing

Oplossing van Probleem 16.1.

1. S2n+1S2n1=a2na2n+10S_{2n+1} - S_{2n-1} = a_{2n} - a_{2n+1} \geq 0 en S2n+2S2n=a2n+2a2n+10S_{2n+2} - S_{2n} = a_{2n+2} - a_{2n+1} \leq 0, terwijl S2nS2n+1=a2n+10S_{2n} - S_{2n+1} = a_{2n+1} \to 0: de rijen (S2n+1)(S_{2n+1}) en (S2n)(S_{2n}) zijn ingesloten en convergeren naar een gemeenschappelijke SS (Stelling 11.11) met S2n+1SS2nS_{2n+1} \leq S \leq S_{2n}. Voor even nn geeft Sn+1SSnS_{n+1} \leq S \leq S_n dat an+1SSn0-a_{n+1} \leq S - S_n \leq 0; voor oneven nn dat 0SSnan+10 \leq S - S_n \leq a_{n+1}. In beide gevallen is SSnan+1\abs{S - S_n} \leq a_{n+1} en heeft SSnS - S_n het teken van (1)n+1(-1)^{n+1}, de eerste weggelaten term. Bij strikte daling wordt elke getoonde ongelijkheid strikt, in het bijzonder 0<SSn<an+10 < \abs{S - S_n} < a_{n+1}.

2. ak=1k!a_k = \frac{1}{k!} daalt strikt naar 00: vraag 1 is van toepassing. Taylor–Lagrange (Stelling 16.7) voor exp\exp tussen 1-1 en 00: e1Tn1(n+1)!\abs{\eu^{-1} - T_n} \leq \frac{1}{(n+1)!} (de afgeleide et\eu^t is daar 1\leq 1), dus Tne1T_n \to \eu^{-1}, en de limiet SS uit vraag 1 is e1\eu^{-1}, met de strikte grenzen 0<e1Tn<1(n+1)!0 < \abs{\eu^{-1} - T_n} < \frac{1}{(n+1)!}.

3. Integratie van de identiteit over [0,1]\intcc{0}{1}: het linkerlid is arctan1=π4\arctan 1 = \frac\pi4 (hoofdstelling), de kk-de term geeft (1)k2k+1\frac{(-1)^k}{2k+1}, en

ρn=(1)n+101t2n+21+t2 ⁣dt,ρn01t2n+2 ⁣dt=12n+3.\rho_n = (-1)^{n+1}\int_0^1 \frac{t^{2n+2}}{1+t^2}\dd t, \qquad \abs{\rho_n} \leq \int_0^1 t^{2n+2}\dd t = \frac{1}{2n+3} .

4. De fout op π\pi is 4ρn42n+34\abs{\rho_n} \leq \frac{4}{2n+3}: onder 10610^{-6} komen vereist 2n+3>41062n + 3 > 4\cdot10^6, ongeveer twee miljoen termen. Ondertussen is 4S4=4(113+1517+19)=4×0.834921=3.3396834S_4 = 4\bigl(1 - \frac13 + \frac15 - \frac17 + \frac19\bigr) = 4 \times 0.834921 = 3.339683, bijna 0.20.2 verwijderd van π\pi: vijf termen, nog geen enkel juist cijfer.

5. Integratie over [0,1]\intcc{0}{1} geeft ln2=k=1n(1)k1k+(1)nRn\ln 2 = \sum_{k=1}^{n} \frac{(-1)^{k-1}}{k} + (-1)^n R_n met Rn=01tn1+t ⁣dtR_n = \int_0^1 \frac{t^n}{1+t}\dd t; uit 1211+t1\frac12 \leq \frac{1}{1+t} \leq 1 volgt 12(n+1)Rn1n+1\frac{1}{2(n+1)} \leq R_n \leq \frac{1}{n+1}. De fout zit gevangen tussen twee veelvouden van 1n\frac1n: zes decimalen kosten ongeveer een miljoen termen.

6. Integratie van 00 tot xx: 12ln1+x1x=k=0nx2k+12k+1+0xt2n+21t2 ⁣dt\frac12\ln\frac{1 + x}{1 - x} = \sum_{k=0}^{n} \frac{x^{2k+1}}{2k+1} + \int_0^x \frac{t^{2n+2}}{1-t^2}\dd t. In x=13x = \frac13 is 1+1/311/3=2\frac{1 + 1/3}{1 - 1/3} = 2, en op [0,13]\intcc{0}{\frac13} is 11t298\frac{1}{1 - t^2} \leq \frac98:

ln2=2k=0n(1/3)2k+12k+1+ρ~n,0<ρ~n94(1/3)2n+32n+3:\ln 2 = 2\sum_{k=0}^{n} \frac{(1/3)^{2k+1}}{2k+1} + \tilde\rho_n, \qquad 0 < \tilde\rho_n \leq \frac{9}{4}\cdot \frac{(1/3)^{2n+3}}{2n+3} :

elke extra term deelt de fout ongeveer door 99.

7. Inductie: voor n=1n = 1 is 112=12=H2H11 - \frac12 = \frac12 = H_2 - H_1. Stap:

H2n+2Hn+1=(H2nHn)+12n+1+12n+21n+1=(H2nHn)+12n+112n+2,H_{2n+2} - H_{n+1} = (H_{2n} - H_n) + \frac{1}{2n+1} + \frac{1}{2n+2} - \frac{1}{n+1} = (H_{2n} - H_n) + \frac{1}{2n+1} - \frac{1}{2n+2},

en dat is precies de aangroei van de alternerende som. En H2nHn=k=1n1n+kH_{2n} - H_n = \sum_{k=1}^{n}\frac{1}{n+k} is de Riemannsom van Voorbeeld 15.21, die naar ln2\ln 2 convergeert: de even partiële sommen van weg 1 zijn de Riemannsommen van weg 3.

8. k=16(1)k1k=0.61667\sum_{k=1}^{6}\frac{(-1)^{k-1}}{k} = 0.61667, met fout 0.07650.0765; weg 2 bij n=5n = 5 geeft 0.69314710.6931471 met fout 94(1/3)1313=1.1107\leq \frac94\cdot\frac{(1/3)^{13}}{13} = 1.1\cdot10^{-7}. De reden: weg 1 evalueert de logaritmereeks in het randpunt x=1x = 1, waar de termen als 1k\frac1k afnemen; weg 2 evalueert in x=13x = \frac13, diep binnenin, waar elke term een verse factor 19\frac19 meedraagt.

9. Tien decimalen: we willen ρ~n51011\tilde\rho_n \leq 5\cdot10^{-11}. Bij n=10n = 10: 94(1/3)2323=941.061011231.01012<51011\frac94 \cdot \frac{(1/3)^{23}}{23} = \frac94\cdot\frac{1.06\cdot10^{-11}}{23} \approx 1.0\cdot10^{-12} < 5\cdot10^{-11}: elf termen volstaan.

10. (5+i)2=24+10i(5+\iu)^2 = 24 + 10\iu, en daarna (5+i)4=(24+10i)2=476+480i(5+\iu)^4 = (24 + 10\iu)^2 = 476 + 480\iu; en 2(1+i)(239+i)=2(238+240i)=476+480i2(1+\iu)(239+\iu) = 2(238 + 240\iu) = 476 + 480\iu: gelijk. Argumenten: arg(5+i)=arctan15\arg(5 + \iu) = \arctan\frac15, dus het linkerlid heeft argument 4arctan150.79(0,π)4\arctan\frac15 \approx 0.79 \in \intoo{0}{\pi}; het rechterlid heeft argument π4+arctan1239(0,π)\frac\pi4 + \arctan\frac{1}{239} \in \intoo{0}{\pi}. Twee gelijke complexe getallen met argumenten in hetzelfde interval van lengte <2π< 2\pi:

4arctan15=π4+arctan1239,4\arctan\frac15 = \frac\pi4 + \arctan\frac{1}{239} ,

en dat is de formule van Machin.

11. Integreer 11+t2=k=0n(1)kt2k+(1)n+1t2n+21+t2\frac{1}{1+t^2} = \sum_{k=0}^n (-1)^k t^{2k} + \frac{(-1)^{n+1}t^{2n+2}}{1+t^2} van 00 tot xx:

arctanx=k=0n(1)kx2k+12k+1+rn(x),rn(x)0xt2n+2 ⁣dt=x2n+32n+3.\arctan x = \sum_{k=0}^{n}\frac{(-1)^k x^{2k+1}}{2k+1} + r_n(x), \qquad \abs{r_n(x)} \leq \int_0^x t^{2n+2}\dd t = \frac{x^{2n+3}}{2n+3} .

12. Fouten: 16(1/5)1111=3.010816\,\frac{(1/5)^{11}}{11} = 3.0\cdot 10^{-8} en 4(1/239)55110124\,\frac{(1/239)^5}{5} \approx 1\cdot10^{-12}: samen <5108< 5\cdot10^{-8}. De getoonde som evalueert tot 3.141592683.14159268\dots, dus π=3.1415926\pi = 3.1415926\dots gecertificeerd tot op 51085\cdot10^{-8}: zeven decimalen uit zes termen (vijf bij 15\frac15, twee bij 1239\frac1{239}, ruim geteld).

13. Leibniz: fout 1n\sim \frac1n, zodat elk nieuw cijfer het werk met tien vermenigvuldigt. Dalzell (Probleem 15.1, vraag 22): fout 415m4^{1-5m}, ongeveer drie cijfers per stap, met per stap een zwaardere veelterm. Machin: foutverhouding 125\frac{1}{25} per term, ongeveer 1.41.4 cijfers per term, met per term één deling. De moraal: de restterm van een ontwikkeling van meetkundig type schaalt als x2nx^{2n}, zodat xx klein maken cijfers koopt tegen vaste kosten per term — de complexe identiteit van Machin is precies een machine om xx te verkleinen.

14. ak=1(2k)!a_k = \frac{1}{(2k)!} daalt strikt naar 00; volgens vraag 1 en Propositie 16.11 (met de grens van Lagrange als in vraag 2) geldt kn(1)k(2k)!cos1\sum_{k \leq n}\frac{(-1)^k}{(2k)!} \to \cos 1 met de strikte insluiting 0<cos1Sn<1(2n+2)!0 < \bigl|\cos 1 - S'_n\bigr| < \frac{1}{(2n+2)!}. Stel cos1=pq\cos 1 = \frac pq en neem 2nq2n \geq q: dan is (2n)!Sn=kn(1)k(2n)!(2k)!Z(2n)!\,S'_n = \sum_{k\leq n} (-1)^k \frac{(2n)!}{(2k)!} \in \Z en (2n)!pqZ(2n)!\,\frac pq \in \Z, terwijl

0<(2n)!cos1(2n)!Sn<(2n)!(2n+2)!=1(2n+1)(2n+2)<1:0 < \Bigl|(2n)!\cos 1 - (2n)!S'_n\Bigr| < \frac{(2n)!}{(2n+2)!} = \frac{1}{(2n+1)(2n+2)} < 1 :

een geheel getal dat niet nul is met absolute waarde <1< 1. Tegenspraak: cos1Q\cos 1 \notin \Q.

15. Identiek met ak=1(2k+1)!a_k = \frac{1}{(2k+1)!}, na vermenigvuldiging met (2n+1)!(2n+1)! waarbij 2n+1q2n + 1 \geq q: sin1Q\sin 1 \notin \Q. En toch is cos21+sin21=1Q\cos^2 1 + \sin^2 1 = 1 \in \Q: producten en sommen van irrationale getallen mogen rationaal zijn — irrationaliteit overleeft geen enkele algebraïsche bewerking gratis.

16. Staartgrens: voor m>nm > n is

k=n+1m1(2k)!1(2n+2)!(1+12+14+)2(2n+2)!,\sum_{k=n+1}^{m} \frac{1}{(2k)!} \leq \frac{1}{(2n+2)!}\Bigl(1 + \frac12 + \frac14 + \dots\Bigr) \leq \frac{2}{(2n+2)!} ,

omdat elke opeenvolgende verhouding 1(2k+1)(2k+2)12\frac{1}{(2k+1)(2k+2)} \leq \frac12 is; de staart is positief (haar eerste term is dat). Dus 0<cosh1kn1(2k)!<2(2n+2)!0 < \cosh 1 - \sum_{k\leq n}\frac{1}{(2k)!} < \frac{2}{(2n+2)!}, en vermenigvuldiging met (2n)!(2n)! waarbij 2nq2n \geq q vangt opnieuw een geheel getal dat niet nul is in (0,1)\intoo{0}{1}: cosh1Q\cosh 1 \notin \Q.

17. cos1m=k(1)km2k(2k)!\cos\frac1m = \sum_k \frac{(-1)^k}{m^{2k}(2k)!}: de termen dalen strikt naar nul, en m2n(2n)!1m2k(2k)!=m2(nk)(2n)!(2k)!Zm^{2n}(2n)!\cdot \frac{1}{m^{2k}(2k)!} = m^{2(n-k)}\frac{(2n)!}{(2k)!} \in \Z voor knk \leq n. Is cos1m=pq\cos\frac1m = \frac pq, vermenigvuldig de strikte insluiting dan met qm2n(2n)!q\,m^{2n}(2n)!: de fout wordt begrensd door qm2(2n+1)(2n+2)<1\frac{q}{m^2(2n+1)(2n+2)} < 1 voor grote nn: tegenspraak. Voor ab\frac ab met a2a \geq 2: het wegwerken van de noemers vermenigvuldigt de staart met b2n(2n)!b^{2n}(2n)!, maar de eerste weggelaten term is a2n+2b2n+2(2n+2)!\frac{a^{2n+2}}{b^{2n+2}(2n+2)!}, en het product a2n+2b2(2n+1)(2n+2)\frac{a^{2n+2}}{b^2(2n+1)(2n+2)} explodeert: de teller a2n+2a^{2n+2} valt niet langer weg, en de val loopt vast. (Het resultaat blijft waar — via machinerie in de stijl van Niven, niet via deze.)

18. In t=xt = x verdwijnt elke term van gg behalve f(x)f(x)=0f(x) - f(x) = 0: g(x)=0g(x) = 0; en AA is zo gekozen dat g(a)=0g(a) = 0. Bij het afleiden telescopeert de som:

g(t)=f(n+1)(t)n!(xt)n+A(xt)nn!=(xt)nn!(Af(n+1)(t)).g'(t) = -\frac{f^{(n+1)}(t)}{n!}(x - t)^n + A\,\frac{(x-t)^n}{n!} = \frac{(x-t)^n}{n!}\bigl(A - f^{(n+1)}(t)\bigr) .

Rolle op het segment van aa tot xx geeft een cc strikt daartussen met g(c)=0g'(c) = 0; omdat (xc)n0(x - c)^n \neq 0, is A=f(n+1)(c)A = f^{(n+1)}(c). Het uitschrijven van g(a)=0g(a) = 0 levert de gelijkheid van Taylor met restterm f(n+1)(c)(n+1)!(xa)n+1\frac{f^{(n+1)}(c)}{(n+1)!}(x-a)^{n+1}.

19. Voor x>0x > 0 is de restterm ec(n+1)!xn+1>0\frac{\eu^{c}} {(n+1)!}x^{n+1} > 0: de exponentiële overtreft elk van haar taylorveeltermen, strikt, in elke orde. Voor x<0x < 0 is het teken van de restterm dat van xn+1x^{n+1}: ex\eu^x ligt boven de veelterm voor nn oneven en eronder voor nn even — afwisselende zijden, zoals de grafieken van 1+x1 + x en 1+x+x221 + x + \frac{x^2}{2} tegenover ex\eu^x al laten zien.

20. Werkelijke fout: sin0.50.4791667=2.59104\sin 0.5 - 0.4791667 = 2.59\cdot10^{-4}, tegenover de grens 0.55120=2.60104\frac{0.5^5}{120} = 2.60\cdot10^{-4}: bijna bereikt (de volgende term domineert de staart). Young: voor limieten en lokale analyse, waar geen constante nodig is. Lagrange: voor gecertificeerde decimalen. Alternerend: waar toepasbaar, dezelfde grens plus de richting van de fout — de beste van de drie, maar de zeldzaamste.

21. Continuïteit in 00: met u=1x2+u = \frac{1}{x^2} \to +\infty is f(x)=eu0=f(0)f(x) = \eu^{-u} \to 0 = f(0). Afgeleide in 00: f(h)h=ueu0\bigl|\frac{f(h)}{h}\bigr| = \sqrt u\,\eu^{-u} \to 0 (Propositie 4.6): f(0)=0f'(0) = 0. Voor x0x \neq 0 is f(x)=2x3e1/x2f'(x) = \frac{2}{x^3}\eu^{-1/x^2}: de vorm P1(1x)e1/x2P_1\bigl(\frac1x\bigr)\eu^{-1/x^2} met P1(X)=2X3P_1(X) = 2X^3; met inductie geeft het afleiden van Pk(1x)e1/x2P_k(\frac1x)\eu^{-1/x^2} dat Pk+1(X)=2X3Pk(X)X2Pk(X)P_{k+1}(X) = 2X^3 P_k(X) - X^2 P_k'(X), een veelterm. Dan is

f(k)(h)0h=vPk(v)ev2v=1/h0\frac{f^{(k)}(h) - 0}{h} = v\,P_k(v)\,\eu^{-v^2} \Big|_{v = 1/h} \longrightarrow 0

(veelterm tegen ev2\eu^{-v^2}, groeivergelijking in ±\pm\infty): met inductie is f(k)(0)=0f^{(k)}(0) = 0 voor alle kk. Alle taylorveeltermen van ff in 00 zijn nul, en toch is f>0f > 0 buiten 00: Taylor–Young is exact in elke orde en blind voorbij de kiem. Een ontwikkeling is uitsluitend lokale informatie.

22. Volgens vraag 2 is 0<e1Tn<1(n+1)!0 < \abs{\eu^{-1} - T_n} < \frac{1}{(n+1)!}, strikt. Is e1=pq\eu^{-1} = \frac pq, neem dan nqn \geq q en vermenigvuldig met n!n!: n!TnZn!\,T_n \in \Z en n!pqZn!\frac pq \in \Z, dus heeft een geheel getal dat niet nul is een absolute waarde <n!(n+1)!=1n+1<1< \frac{n!}{(n+1)!} = \frac{1}{n+1} < 1: tegenspraak. Bijgevolg is e1Q\eu^{-1} \notin \Q, en is e=1e1\eu = \frac{1}{\eu^{-1}} eveneens irrationaal. De drie bewijzen: ingesloten rijen die q!eq!\,\eu insluiten (Oefening 11.9); de integraalrecursie An=enAn1A_n = \eu - nA_{n-1} (Probleem 15.1); de alternerende insluiting (hier). Eén val, drie certificaten van kleinheid.

23. Groepeer de partiële sommen paarsgewijs: k=12n(1)kbk=EnOn\sum_{k=1}^{2n} (-1)^k b_k = E_n - O_n met En=j=1n14j2E_n = \sum_{j=1}^{n} \frac{1}{4j^2}, begrensd (via de telescoperende grens van Voorbeeld 11.22 is En12E_n \leq \frac12), en On=j=1n12j112Hn+O_n = \sum_{j=1}^{n}\frac{1}{2j-1} \geq \frac12 H_n \to +\infty (Oefening 11.5): de partiële sommen naderen tot -\infty. De monotonie werd in vraag 1 precies daar gebruikt waar S2n+1S2n1=a2na2n+1S_{2n+1} - S_{2n-1} = a_{2n} - a_{2n+1} een teken nodig had: zonder daling hoeven de even en de oneven deelrijen niet monotoon te zijn, en stort de insluiting in.

24. (2+i)(3+i)=5+5i=5(1+i)(2+\iu)(3+\iu) = 5 + 5\iu = 5(1+\iu); door argumenten te nemen (alle in (0,π2)\intoo{0}{\frac\pi2}): arctan12+arctan13=π4\arctan\frac12 + \arctan\frac13 = \frac\pi4. Kosten in termen voor zes decimalen: de fout is 4((1/2)2n+32n+3+(1/3)2n+32n+3)\leq 4\bigl(\frac{(1/2)^{2n+3}}{2n+3} + \frac{(1/3)^{2n+3}}{2n+3}\bigr), wat bij n=10n = 10 ongeveer 2108<51072\cdot10^{-8} < 5\cdot10^{-7} is: elf termen. Rangschikking: exponentieel beter dan Leibniz, achter Machin (waarvan het dominante punt 15\frac15 kleiner is dan 12\frac12): ruwweg 0.60.6 cijfers per term tegenover de 1.41.4 van Machin.

25. (i) Voor dalende ak0a_k \to 0 convergeren de alternerende partiële sommen met SSnan+1\abs{S - S_n} \leq a_{n+1}, waarbij de fout het teken draagt van de eerste weggelaten term. (ii) Een eindige identiteit met expliciete restterm kan in een gekozen punt worden geëvalueerd en begrensd, terwijl een limietuitspraak alleen uiteindelijke nabijheid belooft — certificering vraagt het eerste. (iii) Geoogst: π\pi tot op 51085\cdot10^{-8} via Machin, ln2\ln 2 tot op tien decimalen via de 13\frac13-reeks, de irrationaliteit van cos1\cos 1, sin1\sin 1, cosh1\cosh 1, cos1m\cos\frac1m en e1\eu^{-1}, de gelijkheid van Taylor–Lagrange, en de waarschuwing van de vlakke functie. (iv) Hoofdstuk 17 tilt vraag 1 op tot het criterium voor alternerende reeksen, scheidt absolute van voorwaardelijke convergentie, en haar weekendopgave voert het herschikkingsdrama op waarvoor de alternerende harmonische reeks van weg 1 de kroongetuige is.