Mathematics · Book 3 · Bachelor Year 1

Universitaire wiskunde — Bachelor jaar 1

Universitaire wiskunde — Bachelor jaar 1 · Bachelor Year 1

16Taylorformules en asymptotische ontwikkelingen

Nabij een punt is een gladde functie net zo goed als een veelterm — met een beheersbare fout. De Taylorformules maken dit exact in drie gedaanten (integraalrest, Lagrange-rest, Young-rest), en de resulterende asymptotische ontwikkelingen, algebraïsch gemanipuleerd, worden het scherpste gereedschap van de elementaire analyse: limieten, equivalenten, lokaal gedrag, asymptoten.

16.1 Vergelijkingsnotatie

Definitie 16.1 (Landau-notatie)

Laat f,gf, g gedefinieerd zijn nabij x0x_0 (x0Rx_0 \in \R of ±\pm\infty). Men schrijft, als xx0x \to x_0:

  • f=o(g)f = o(g) (“kleine-o”) als f=εgf = \varepsilon g met ε(x)0\varepsilon(x) \to 0;
  • f=O(g)f = O(g) (“grote-O”) als f=ugf = u g met uu begrensd nabij x0x_0;
  • fgf \sim g (“equivalent”) als f=(1+ε)gf = (1 + \varepsilon) g met ε0\varepsilon \to 0 — equivalent 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; fgf \sim g impliceert dat ff en gg dezelfde limieten, tekens (nabij x0x_0) en afwezigheid van nulpunten delen;
  2. equivalenten vermenigvuldigen en delen: f1g1f_1 \sim g_1, f2g2f_2 \sim g_2 impliceren 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 is een korte manipulatie van de definities; 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 onderdelen van (4) verdienen hun eigen regel. uo(g)=o(ug)u\cdot o(g) = o(ug): als f=εgf = \varepsilon g met ε0\varepsilon \to 0, dan uf=ε(ug)uf = \varepsilon\,(ug) met dezelfde ε\varepsilon. o(o(g))=o(g)o(o(g)) = o(g): als f=ε1hf = \varepsilon_1 h en h=ε2gh = \varepsilon_2 g met beide εi0\varepsilon_i \to 0, dan f=(ε1ε2)gf = (\varepsilon_1\varepsilon_2) g en het product van de twee infinitesimalen is er één. De equivalentie fg    f=g+o(g)f \sim g \iff f = g + o(g) is de definitie tweemaal 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 oplopende volgorde van sterkte:

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 gewoontes die het aanleren waard zijn: ten eerste, een O()O(\cdot) die in een kleinere klasse landt wordt stilzwijgend opgewaardeerd (O(lnx)O(\ln x) is ook o(x0.01)o(x^{0.01})); ten tweede, bij x0+x \to 0^+ keert de hele ladder om via de substitutie x1xx \mapsto \frac1xlnx=o(x0.01)\ln x = o(x^{-0.01}) daar, zodat “xαlnx0x^\alpha \ln x \to 0” geldt voor elke α>0\alpha > 0. De schaal recht houden is de helft van elk asymptotisch argument in Hoofdstuk 17.

Voorbeeld 16.4 (Uniciteit van ontwikkelingen, en een pariteitsdividend)

Als een functie twee ontwikkelingen in 00 toelaat tot dezelfde orde,

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 ak=bka_k = b_k voor elke kk: door af te trekken en ck=akbkc_k = a_k - b_k te stellen, evalueer je de identiteit c0+c1x++cnxn=o(xn)c_0 + c_1 x + \dots + c_n x^n = o(x^n) voor x0x \to 0 om c0=0c_0 = 0 te krijgen; deel door xx en herhaal — elke deling is legitiem omdat de resterende uitdrukking opnieuw o(xnk)o(x^{n-k}) is. De coëfficiënten zijn dus intrinsiek, en men kan ze langs elke weg berekenen (Taylor-afgeleiden, algebra op bekende ontwikkelingen, 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 uniciteit; evenzo hebben oneven functies oneven machten. Daarom draagt cos\cos o(x2p+1)o(x^{2p+1}) in plaats van o(x2p)o(x^{2p}) in de tabel hieronder: de afwezige oneven term is gratis informatie, één orde precisie voor niets.

16.2 De drie Taylorformules

Stelling 16.5 (Taylor met integraalrest)

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

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: f(x)=f(a)+axf(t) ⁣dtf(x) = f(a) + \int_a^x f'(t)\dd t is de hoofdstelling (Stelling 15.9). Stap: integreer de rest 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 haakjes de term f(n+1)(a)(n+1)!(xa)n+1\frac{f^{(n+1)}(a)}{(n+1)!}(x - a)^{n+1} bijdragen.

Voorbeeld 16.6 (Een exacte ontwikkeling met haar rest)

Voor ln(1+x)\ln(1 + x) kan de integraalrest volledig expliciet gemaakt worden 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 is de rest begrensd door 0xtn ⁣dt=xn+1n+1\int_0^x t^n\,\dd t = \frac{x^{n+1}}{n+1}. Dit is sterker dan Taylor–Young op twee manieren: het is een identiteit die geldig is voor een vaste xx (niet alleen x0x \to 0), en de foutgrens is numeriek. Het weekendprobleem (Probleem 16.1) leeft op zulke exacte vormen; Taylor–Young hieronder is het lichtere gereedschap voor limieten, waar alleen de vorm van de fout ertoe doet.

Stelling 16.7 (Taylor–Lagrange-ongelijkheid)

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

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 integraalrest: 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 (Gecertificeerde numeriek)

Wat 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 de fout is 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 (echte waarde 1.00995049381.0099504938\dots — de grens is bijna scherp). Het afsluitende inzicht: Taylor–Young zegt alleen hoe snel de fout verdwijnt; Taylor–Lagrange zet dezelfde veelterm om in een certificaat, een getal plus een bewezen foutmarge. Telkens wanneer in dit boek een decimale bewering wordt gedaan, staat er een Lagrange-achtige grens achter; het weekendprobleem (Probleem 16.1) industrialiseert het idee.

Stelling 16.9 (Taylor–Young)

Laat ff nn keer afleidbaar zijn in aa. Dan, 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 bij orde n1n - 1, en laat ff nn keer afleidbaar zijn 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).

Laat 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 g=rg' = r en g(a)=0g(a) = 0. Gegeven ε>0\varepsilon > 0, kies δ\delta met r(t)εtan1\abs{r(t)} \leq \varepsilon\abs{t - a}^{n-1} voor taδ\abs{t-a} \leq \delta; de middelwaardeongelijkheid (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 gelijk op met de conclusies. Taylor–Young vraagt het minst (nn afgeleiden alleen in het punt) en levert het minst: een kwalitatieve o((xa)n)o((x-a)^n), perfect voor limieten, nutteloos voor gecertificeerde cijfers. De Lagrange-ongelijkheid vraagt Cn+1C^{n+1} op het interval en een grens MM daar, en geeft een numerieke foutmarge terug. De integraalvorm vraagt dezelfde regulariteit en geeft het meest terug: de fout als een expliciet object dat men kan transformeren (partieel integreren, stuksgewijs begrenzen, van variabele wisselen) — het is de vorm die de irrationaliteitsmachine van Probleem 15.1 aandreef. Het kiezen van de zwakste formule die het doel ondersteunt is geen muggenzifterij: de vlakke functie van Probleem 16.1 voldoet aan Taylor–Young in elke orde terwijl elke sterkere conclusie over haar onwaar is weg van 00.

Propositie 16.11 (Standaardontwikkelingen in 00)

Als x0x \to 0, 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, sin\sin zonder de alternerende tekens.)

Bewijs. Elke functie is glad nabij 00 met afgeleiden die eenvoudig te evalueren zijn: (ex)(k)=ex(\eu^x)^{(k)} = \eu^x; de afgeleiden van sin\sin en cos\cos lopen rond 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 laten vallen is een fout.
  2. Sommen, 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 vergt meer voorzichtigheid — rechtvaardiging: de integraal van o(tn)o(t^n) van 00 tot xx is o(xn+1)o(x^{n+1}), door directe begrenzing).

Voorbeeld 16.13 (Samenstelling, met de boekhouding getoond)

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 gaat. 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) aangezien 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 samen:

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 x3x^3-bijdragen 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, maar omdat het eerste verschil van de exponenten (x36-\frac{x^3}{6}) binnenkomt vermenigvuldigd met e0=1\eu^0 = 1 en dan wordt beantwoord door de kubische term van de buitenste exponentiële; orde-na-orde-boekhouding detecteert zulke samenzweringen, het oog nooit. (De volgende term is x48-\frac{x^4}{8}: het bestand eindigt bij orde 44.)

Voorbeeld 16.14

Ontwikkeling van tan\tan tot 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),

dan

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) tegen zijn 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, aangescherpt in 0 en stil ver weg. De Taylor–Lagrange-grens |x|n+1(n+1)! kwantificeert het venster; de vlakke functie van  toont dat het contract zelfs leeg kan zijn voorbij het punt zelf.
De sinus (donker) tegen zijn 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, aangescherpt in 00 en stil ver weg. De Taylor–Lagrange-grens xn+1(n+1)!\frac{\abs{x}^{n+1}}{(n+1)!} kwantificeert het venster; de vlakke functie van Probleem 16.1 toont dat het contract zelfs leeg kan zijn voorbij het punt zelf.

16.3 Toepassingen

Voorbeeld 16.15 (Limieten)

limx0xsinxx3:xsinx=x36+o(x3)x36,dus de limiet is 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 de limiet is } \frac16

— wat de vraag beslecht die in Oefening 4.9 werd opgeworpen. 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 (Veelvoorkomende valkuilen met ontwikkelingen)

(i) Tel of trek equivalenten nooit op of af: uit tanxx\tan x \sim x en sinxx\sin x \sim x mag men niet concluderen tanxsinx0\tan x - \sin x \sim 0 (betekenisloos) — de eerlijke weg zijn 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 het bloedbad: in dezelfde berekening ziet orde 11 alleen xx=0x - x = 0; telkens wanneer de leidende termen elkaar opheffen, verhoog je de orde totdat een coëfficiënt ongelijk aan nul overleeft, 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 aan en2\eu^{n^2} — exponentieer alleen ontwikkelingen van de exponent waarvan de fout naar 00 gaat, nooit equivalenten van de exponent. (Logaritmen zijn veiliger: als unvn1u_n \sim v_n \to \ell \neq 1, >0\ell > 0, dan lnunlnvn\ln u_n \sim \ln v_n.) (iv) De o()o(\cdot)-rekenkunde 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 specifieke functie, dus hef nooit twee ervan 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 (eerste term ongelijk aan nul na de constante; p2p \geq 2 in een kritiek punt).

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

Bewijs. Nabij 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: constant teken voor pp even, wisselend voor pp oneven.

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

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

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 misgegaan zou zijn met minder zorg: de exponent stoppen bij n21n=n+O(1)n^2\cdot\frac1n = n + O(1) laat een eO(1)\eu^{O(1)}-factor over — begrensd maar niet naar 11 gaand — en geen enkel equivalent kan 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 een term die naar nul gaat, waarbij elke coëfficiënt daarvoor exact behouden blijft.

Voorbeeld 16.19 (Een vlak kritiek punt classificeren)

Bestudeer f(x)=cosx+x22f(x) = \cos x + \frac{x^2}{2} nabij 00. Zowel f(0)=0f'(0) = 0 als f(0)=cos0+1=0f''(0) = -\cos 0 + 1 = 0: de tweede-afgeleidetest is stom. 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) :

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

Voorbeeld 16.20 (Ontwikkelingen op oneindig)

Twee berekeningen waar de variabele naar ++\infty loopt en de substitutie h=1x0+h = \frac1x \to 0^+ de hele gereedschapskist 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 op oneindig: 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),

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 zijn asymptoot y=π2y = \frac\pi2 van onderaf met snelheid 1x\frac1x. Het afsluitende inzicht: er is geen aparte theorie van ontwikkelingen op oneindig — één omgekeerde substitutie herleidt ze tot ontwikkelingen in 00, mits elke tussenliggende OO en oo eerlijk wordt meegedragen.

Voorbeeld 16.21 (Asymptoot door ontwikkeling)

Als x+x \to +\infty,

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 onderaf (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 ze convergentie (equivalenten voeden de vergelijkingstesten, en de studie van 1nα\sum \frac{1}{n^\alpha} is een vermomde ontwikkeling); in het Jaar 2-volume worden ze machtreeksen, waar de Taylorveelterm oneindig veel termen en een convergentiestraal verwerft; en elke linearisatie in de fysica — de slinger, eerste-orde-verstoring — is een Taylor–Young-uitspraak met de o()o(\cdot) stilzwijgend weggelaten. De ene waarschuwing die het graveren waard is: een ontwikkeling beschrijft een functie alleen nabij één punt — zie de vlakke functie van het weekendprobleem, waarvan de ontwikkeling in 00 identiek nul is zonder dat de functie dat is.

Opmerking 16.23 (Perspectieven binnen dit volume)

Ontwikkelingen zijn de werktaal van de resterende analyse en van de meetkunde die komt. Hoofdstuk 17 zet ze om in convergentieverdicten: een equivalent van de algemene term is een ontwikkeling afgekapt bij zijn eerste term, en de fijnere testen (alternerend met foutcontrole) verbruiken ook de tweede term. Hoofdstuk 24 leest lokale meetkunde af uit ontwikkelingen van de twee coördinaatfuncties: of een geparametriseerde kromme kruist, kust of een keerpunt heeft in een punt 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 bij orde één met opzet: het raakvlak is een tweevariabelen- Taylor–Young-uitspraak, waarbij de volledige tweede-orde-theorie (Hessianen, zadelpunten) wordt uitgesteld tot het Jaar 2-volume. De gemeenschappelijke 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); integreren van 00 tot xx (Methode 16.12 (5)):

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

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

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}, d.w.z. (n+1)!>6×106(n+1)! > 6\times 10^{6}: aangezien 10!=362880010! = 3\,628\,800 en 11!=3991680011! = 39\,916\,800, volstaat n+1=11n + 1 = 11, d.w.z. 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). Exponentiëren, 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 ★★

Bestudeer 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 vind de positie van de grafiek van h(x)=exh(x) = \eu^x ten opzichte van zijn raaklijn in a=1a = 1, lokaal en dan 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, oneven pp: geen extremum; gg kruist zijn (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 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 positie van de kromme ten opzichte van hen.

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

Oefening 16.8 ★★

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

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

Laat ff van klasse C2C^2 zijn op R\R. Bewijs dat voor elke xx en 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 Landau–Kolmogorov-achtige ongelijkheid af: als fM0\abs f \leq M_0 en fM2\abs{f''} \leq M_2 op R\R, dan f2M0M2\abs{f'} \leq \sqrt{2 M_0 M_2} overal. (Optimaliseer over hh.)

Oplossing

Oplossing van Oefening 16.9.

Taylor–Lagrange bij 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: 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 is geminimaliseerd bij h=2M0/M2h = \sqrt{2M_0/M_2} (afgeleide nul), met waarde 2M0M2\sqrt{2M_0M_2} — vandaar f2M0M2\abs{f'} \leq \sqrt{2M_0M_2}. (Als M2=0M_2 = 0, laat 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 (rechtvaardig kort). Om zijn snelheid te vinden, 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 dan het equivalent un3nu_n \sim \sqrt{\dfrac{3}{n}}.
Oplossing

Oplossing van Oefening 16.10.

Op (0,π)\intoo{0}{\pi}: 0<sinu<u0 < \sin u < u, dus (un)(u_n) is strikt dalend, positief, 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 (aangezien 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), vnn13\frac{v_n}{n} \to \frac13, d.w.z. vnn3v_n \sim \frac n3, d.w.z. un23nu_n^2 \sim \frac 3n: aangezien 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 tot een gemeenschappelijke noemer te herleiden en teller en noemer afzonderlijk te ontwikkelen.

Oplossing

Oplossing van Oefening 16.11.

Gemeenschappelijke noemers. 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) ,

dus de teller is x43+o(x4)-\frac{x^4}{3} + o(x^4) 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 af

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).
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 afgeleide tan2x0\tan^2 x \geq 0, die alleen in het enkele punt kπk\pi verdwijnt: 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, 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

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 tanεk=1xk\tan\varepsilon_k = \frac{1}{x_k} en εk(0,π2)\varepsilon_k \in \intoo{0}{\frac\pi2}. Aangezien xkkπx_k \geq k\pi \to \infty: ε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 Probleem: Alternerende sommen, gecertificeerde cijfers, en de irrationaliteit van cos1\cos 1

Probleem 16.1

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

Een alternerende som met dalende termen is het vriendelijkste object in de numerieke analyse: zijn fout is begrensd door de eerste weggelaten term, met bekend teken. Dit probleem bewijst dat principe met de stelling van de ingesloten rijen, en besteedt het dan op drie manieren: gecertificeerde decimalen voor ln2\ln 2 (drie concurrerende wegen) en voor π\pi (Leibniz, dan de formule van Machin uit 1706, nog steeds het idee achter recordberekeningen gedurende eeuwen), de irrationaliteit van cos1\cos 1, sin1\sin 1 en cosh1\cosh 1, en, als tegengewicht, de Taylor–Lagrange-gelijkheid en de vlakke functie waarvan de Taylorontwikkeling liegt. Overal is de “reeks”-taal 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. Laat (ak)k0(a_k)_{k \geq 0} dalen 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}) niet-dalend is, (S2n)(S_{2n}) niet-stijgend, en dat ze 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 als de daling strikt is, al deze ongelijkheden strikt zijn.

  2. Eerste dividend: voor x=1x = 1 in de exponentiële reeks, vergelijk 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!} naar e1\eu^{-1} convergeert met e1Tn<1(n+1)!\abs{\eu^{-1} - T_n} < \frac{1}{(n+1)!}.
  3. (Leibniz, 1674) 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}, bewijs

    π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 op drie manieren.

  1. (Manier 1: alternerend harmonisch) 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. (Manier 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 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. (Manier 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 herwin ln2\ln 2 als de Riemannsom-limiet van Voorbeeld 15.21: Manieren 1 en 3 zijn stiekem hetzelfde getal tweemaal gezien.

  4. Schietwedstrijd bij zes termen: vergelijk k=16(1)k1k=0.6167\sum_{k=1}^{6} \frac{(-1)^{k-1}}{k} = 0.6167 met Manier 2 bij n=5n = 5, die al ln2=0.693147\ln 2 = 0.693147 geeft met fout 1.1107\leq 1.1\cdot10^{-7}. Verklaar de structurele reden (evaluatiepunt diep binnen het convergentie-interval versus op de rand ervan).
  5. Hoeveel termen van Manier 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 verifieer de complexe identiteit

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

    Door argumenten te nemen (Hoofdstuk 3-conventies), leid de formule van Machin af

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

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

  2. Zoals in vraag 3, stel voor 0<x<10 < x < 1 vast:

    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 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 Dalzell-integralen van Probleem 15.1 (fout 415m4^{1-5m}), Machin (fout 165(2n+3)\approx 16\cdot 5^{-(2n+3)}) — in cijfers per term, en verklaar waarom het verkleinen van het evaluatiepunt alles verslaat.

Deel IV — De gehele-getallenval, alternerende editie.

  1. Stel 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)! met 2nq2n \geq q, en leid een tegenspraak af: cos1\cos 1 is irrationaal.
  2. Pas aan op 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 cos21+sin21=1\cos^2 1 + \sin^2 1 = 1: irrationaliteit is niet stabiel onder algebra.
  3. De niet-alternerende neef: cosh1=k1(2k)!\cosh 1 = \sum_k \frac{1}{(2k)!} (in de partiële-som-zin, met de tweezijdige staartgrens 0<cosh1kn1(2k)!<2(2n+2)!0 < \cosh 1 - \sum_{k \leq n} \frac{1}{(2k)!} < \frac{2}{(2n+2)!}, te bewijzen). Concludeer cosh1Q\cosh 1 \notin \Q met dezelfde val.
  4. Duw door naar cos1m\cos\frac 1m voor elk geheel getal m1m \geq 1: vermenigvuldig met m2n(2n)!m^{2n}(2n)! en concludeer cos1mQ\cos\frac1m \notin \Q. Waar breekt dezelfde poging voor cosab\cos\frac ab met b>1b > 1 algemeen? (Identificeer de noemer die niet langer wordt weggewerkt.)

Deel V — Scherper en donkerder: de gelijkheidsvorm, en een functie die Taylor voor de gek houdt.

  1. (Taylor–Lagrange, gelijkheidsvorm) Laat ff n+1n + 1 keer afleidbaar zijn 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 concludeer 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: voor x>0x > 0 toon aan

    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 omkeert voor x<0x < 0 afhankelijk van de pariteit van nn.

  3. Vergelijk de resten 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 teken- informatie sin0.5>0.50.536\sin 0.5 > 0.5 - \frac{0.5^3}{6}. Vergelijk met de echte fout 2.591042.59\cdot10^{-4}: de grens wordt bijna bereikt. Naar welk gereedschap zou je grijpen, en wanneer?
  4. (De vlakke functie) Laat f(x)=e1/x2f(x) = \eu^{-1/x^2} voor x0x \neq 0, 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). Concludeer: alle Taylorveeltermen van ff in 00 verdwijnen, en toch f(x)>0f(x) > 0 voor x0x \neq 0: Taylor–Young geldt in elke orde, en zegt niets over ff weg van 00. Ontwikkelingen beschrijven kiemen, geen functies.

Deel VI — Synthese.

  1. Loop de val nog één keer, op e1=k(1)kk!\eu^{-1} = \sum_k \frac{(-1)^k}{k!}: vermenigvuldig de strikte alternerende insluiting met n!n! en concludeer e1Q\eu^{-1} \notin \Q, 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 nodig had.
  2. De kleine lettertjes: monotonie is niet decoratief. Laat bk=1kb_k = \frac1k voor oneven kk en bk=1k2b_k = \frac{1}{k^2} voor even kk: de bkb_k zijn positief en gaan naar 00, en toch divergeren de partiële sommen van (1)kbk\sum (-1)^k b_k naar -\infty. Bewijs het (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 monotonie gebruikte.
  3. Verifieer de eenvoudigere identiteit van Euler arctan12+arctan13=π4\arctan\frac12 + \arctan\frac13 = \frac\pi4 via (2+i)(3+i)=5(1+i)(2 + \iu)(3 + \iu) = 5(1 + \iu), schat de termen nodig voor zes decimalen van π\pi langs deze weg (n=10n = 10 volstaat), en plaats het tussen Leibniz en Machin in de rangschikking van vraag 13.
  4. Synthese, één zin elk: (i) formuleer de alternerende schatting en haar twee uitkomsten (grens en teken); (ii) waarom exacte eindige identiteiten met expliciete resten limietuitspraken verslaan voor gecertificeerde numeriek; (iii) inventaris van het probleem (π\pi tot 10710^{-7} met de hand, ln2\ln 2 tot tien decimalen, vier irrationaliteitsbewijzen, één gelijkheidsstelling, één waarschuwend voorbeeld); (iv) welke van deze draden Hoofdstuk 17 zal oppikken (de alternerende-reekstest, absolute versus voorwaardelijke convergentie, en het herschikkingsdrama van zijn weekendprobleem).
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}), (S2n)(S_{2n}) zijn ingesloten, convergerend naar een gemeenschappelijke SS (Stelling 11.11) met S2n+1SS2nS_{2n+1} \leq S \leq S_{2n}. Voor even nn: Sn+1SSnS_{n+1} \leq S \leq S_n geeft an+1SSn0-a_{n+1} \leq S - S_n \leq 0; voor oneven nn: 0SSnan+10 \leq S - S_n \leq a_{n+1}. In beide gevallen SSnan+1\abs{S - S_n} \leq a_{n+1} en SSnS - S_n heeft het teken van (1)n+1(-1)^{n+1}, de eerste weggelaten term. Strikte daling maakt elke weergegeven 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 1\leq 1 daar), dus Tne1T_n \to \eu^{-1}, en de limiet SS van 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. De identiteit integreren 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} vereist 2n+3>41062n + 3 > 4\cdot10^6, ongeveer twee miljoen termen. Ondertussen 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 cijfer.

5. Integreren over [0,1]\intcc{0}{1}: 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: 12(n+1)Rn1n+1\frac{1}{2(n+1)} \leq R_n \leq \frac{1}{n+1}. De fout is ingesloten tussen twee veelvouden van 1n\frac1n: zes decimalen kosten ongeveer een miljoen termen.

6. Integreren 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: 1+1/311/3=2\frac{1 + 1/3}{1 - 1/3} = 2, en op [0,13]\intcc{0}{\frac13}, 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 door ongeveer 99.

7. Inductie: voor n=1n = 1: 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},

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

8. k=16(1)k1k=0.61667\sum_{k=1}^{6}\frac{(-1)^{k-1}}{k} = 0.61667, fout 0.07650.0765; Manier 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: Manier 1 evalueert de logaritmereeks in het randpunt x=1x = 1, waar de termen afnemen als 1k\frac1k; Manier 2 evalueert in x=13x = \frac13, diep binnenin, waar elke term een verse factor 19\frac19 draagt.

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

wat de formule van Machin is.

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}: totaal <5108< 5\cdot10^{-8}. De weergegeven som evalueert tot 3.141592683.14159268\dots, vandaar π=3.1415926\pi = 3.1415926\dots gecertificeerd tot 51085\cdot10^{-8}: zeven decimalen uit zes termen (vijf in 15\frac15, twee in 1239\frac1{239}, royaal geteld).

13. Leibniz: fout 1n\sim \frac1n, dus elk nieuw cijfer vermenigvuldigt de werklast met tien. Dalzell (Probleem 15.1, vraag 22): fout 415m4^{1-5m}, ongeveer drie cijfers per stap, elke stap een zwaardere veelterm. Machin: foutverhouding 125\frac{1}{25} per term, ongeveer 1.41.4 cijfers per term, elke term één deling. De moraal: de rest van een expansie van meetkundig type schaalt als x2nx^{2n}, dus xx klein maken koopt cijfers tegen een vaste kost per term — Machins complexe identiteit 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 (Lagrange- grens zoals in vraag 2), 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 (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 ongelijk aan nul met absolute waarde <1< 1. Tegenspraak: cos1Q\cos 1 \notin \Q.

15. Identiek met ak=1(2k+1)!a_k = \frac{1}{(2k+1)!}, vermenigvuldigend met (2n+1)!(2n+1)! met 2n+1q2n + 1 \geq q: sin1Q\sin 1 \notin \Q. En toch cos21+sin21=1Q\cos^2 1 + \sin^2 1 = 1 \in \Q: producten en sommen van irrationalen kunnen rationaal zijn — irrationaliteit gaat door geen enkele algebraïsche bewerking gratis.

16. Staartgrens: voor m>nm > n,

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

aangezien elke opeenvolgende verhouding 1(2k+1)(2k+2)12\frac{1}{(2k+1)(2k+2)} \leq \frac12 is; de staart is positief (zijn 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 vermenigvuldigen met (2n)!(2n)! met 2nq2n \geq q vangt opnieuw een geheel getal ongelijk aan nul 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. Als cos1m=pq\cos\frac1m = \frac pq, vermenigvuldig de strikte insluiting met qm2n(2n)!q\,m^{2n}(2n)!: de fout is 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 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} wordt niet langer weggewerkt, en de val klemt. (Het resultaat is nog steeds waar — via Niven-achtige machinerie, niet 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; AA is zo gekozen dat g(a)=0g(a) = 0. Afleiden, de som telescopeert:

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 cc strikt ertussen met g(c)=0g'(c) = 0; aangezien (xc)n0(x - c)^n \neq 0: A=f(n+1)(c)A = f^{(n+1)}(c). Het ontvouwen van g(a)=0g(a) = 0 levert de Taylor-gelijkheid met rest 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 rest ec(n+1)!xn+1>0\frac{\eu^{c}} {(n+1)!}x^{n+1} > 0: de exponentiële overtreft elk van zijn Taylorveeltermen, strikt, in elke orde. Voor x<0x < 0 is het teken van de rest dat van xn+1x^{n+1}: ex\eu^x ligt boven de veelterm voor nn oneven, onder haar voor nn even — afwisselende zijden, zoals de grafieken van 1+x1 + x en 1+x+x221 + x + \frac{x^2}{2} tegen ex\eu^x al tonen.

20. Echte fout: sin0.50.4791667=2.59104\sin 0.5 - 0.4791667 = 2.59\cdot10^{-4}, tegen 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: wanneer van toepassing, 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, 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, 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; inductief, het afleiden van Pk(1x)e1/x2P_k(\frac1x)\eu^{-1/x^2} geeft Pk+1(X)=2X3Pk(X)X2Pk(X)P_{k+1}(X) = 2X^3 P_k(X) - X^2 P_k'(X), een veelterm. Dan

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): door inductie f(k)(0)=0f^{(k)}(0) = 0 voor alle kk. Alle Taylorveeltermen van ff in 00 verdwijnen, en toch f>0f > 0 buiten 00: Taylor–Young is exact in elke orde en blind voorbij de kiem. Een ontwikkeling is enkel lokale informatie.

22. Volgens vraag 2, 0<e1Tn<1(n+1)!0 < \abs{\eu^{-1} - T_n} < \frac{1}{(n+1)!}, strikt. Als e1=pq\eu^{-1} = \frac pq, neem nqn \geq q en vermenigvuldig met n!n!: n!TnZn!\,T_n \in \Z en n!pqZn!\frac pq \in \Z, dus een geheel getal ongelijk aan nul heeft absolute waarde <n!(n+1)!=1n+1<1< \frac{n!}{(n+1)!} = \frac{1}{n+1} < 1: tegenspraak. Vandaar e1Q\eu^{-1} \notin \Q, en e=1e1\eu = \frac{1}{\eu^{-1}} is eveneens irrationaal. De drie bewijzen: ingesloten rijen die q!eq!\,\eu insluiten (Oefening 11.9); de integraalrecurrentie 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 in paren: 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 (door de telescoperende grens van Voorbeeld 11.22, 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 gaan naar -\infty. Monotonie werd in vraag 1 precies 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 oneven deelrijen niet monotoon te zijn, en de insluiting stort in.

24. (2+i)(3+i)=5+5i=5(1+i)(2+\iu)(3+\iu) = 5 + 5\iu = 5(1+\iu); argumenten nemend (alle in (0,π2)\intoo{0}{\frac\pi2}): arctan12+arctan13=π4\arctan\frac12 + \arctan\frac13 = \frac\pi4. Reekskost voor zes decimalen: fout 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 2108<5107\approx 2\cdot10^{-8} < 5\cdot10^{-7} is: elf termen. Rangschikking: beter dan Leibniz met een exponentiële marge, achter Machin (wiens dominante punt 15\frac15 kleiner is dan 12\frac12): ruwweg 0.60.6 cijfers per term tegen Machins 1.41.4.

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} en draagt de fout het teken van de eerste weggelaten term. (ii) Een eindige identiteit met expliciete rest kan geëvalueerd en begrensd worden in een gekozen punt, terwijl een limietuitspraak alleen eventuele nabijheid belooft — certificering vereist de eerste. (iii) Uitgehaald: π\pi tot 51085\cdot10^{-8} met Machin, ln2\ln 2 tot tien decimalen met de 13\frac13-reeks, irrationaliteit van cos1\cos 1, sin1\sin 1, cosh1\cosh 1, cos1m\cos\frac1m en e1\eu^{-1}, de Taylor–Lagrange-gelijkheid, en de vlakke-functiewaarschuwing. (iv) Hoofdstuk 17 waardeert vraag 1 op tot de alternerende- reekstest, scheidt absolute van voorwaardelijke convergentie, en zijn weekendprobleem ensceneert het herschikkingsdrama waarvoor de alternerende harmonische reeks van Manier 1 de sterkroongetuige is.