="http://www.w3.org/2000/svg" viewBox="0 0 512 512">

8 Die Ableitung und das Riemann-Integral

In diesem Kapitel wollen wir die zwei wichtigen Themen «Riemann-Integral» aus Kapitel 4 und «Ableitung» aus Kapitel 7 verknüpfen. Wie wir sehen werden, ist dieser Zusammenhang für Anwendungen und für die weitere Theorie von fundamentaler Bedeutung.

8.1 – Der Fundamentalsatz der Integral- und Differentialrechnung

Definition 8.1

Sei [latex][a,b][/latex] ein kompaktes Intervall in [latex]\mathbb {R}[/latex] und sei [latex]f \in R([a,b])[/latex] eine auf [latex][a,b][/latex] Riemann-integrierbare Funktion. Die Funktion

[latex]
\begin{aligned}[]x \in [a,b] \mapsto \int _a^x f(t) \thinspace {\rm {d}} t\end{aligned}
[/latex]

nennt sich das Integral mit veränderlichen oberen Grenze oder das partikuläre Integral von [latex]f[/latex].

Man beachte, dass für [latex]f\in R([a,b])[/latex] und [latex]x \in [a,b][/latex] das Integral [latex]\int _a^x f(t) \thinspace {\rm {d}} t[/latex] wohldefiniert ist (siehe Satz 4.26 über die Intervalladditivität des Riemann-Integrals). Wir erinnern daran, dass wir in Übung 4.28 bereits die Stetigkeit des partikulären Integrals gezeigt haben. Mit etwas stärkeren Annahmen ergibt sich nun folgender Satz von fundamentaler Bedeutung.

Theorem 8.2: Ableitung des Integrals

Sei [latex][a,b][/latex] ein kompaktes Intervall mit Endpunkten [latex]a

Beweis

Sei [latex]\varepsilon > 0[/latex]. Dann existiert ein [latex]\delta >0[/latex] mit

[latex]
\begin{aligned}[]|f(x)-f(x_0)| [/latex]

für alle [latex]x \in (x_0-\delta ,x_0 +\delta )\cap [a,b][/latex]. Wir verwenden dies nun in Kombination mit der Dreiecks-Ungleichung für das Riemann-Integral (Satz 4.24) und der Intervalladditivität des Riemann-Integrals (Satz 4.26), um die Aussage zu zeigen. Für [latex]x \in (x_0,x_0 + \delta )\cap [a,b] \setminus \left \lbrace {x_0} \right \rbrace[/latex] gilt

[latex]
\begin{aligned}[]\left | \frac {F(x)-F(x_0)}{x-x_0} -f(x_0) \right | &= \left |\frac {1}{x-x_0} \left (\int _a^x f(t) \thinspace {\rm {d}} t - \int _a^{x_0}f(t) \thinspace {\rm {d}} t\right ) -f(x_0)\right |\\ &= \left |\frac {1}{x-x_0} \int _{x_0}^x f(t) \thinspace {\rm {d}} t -f(x_0) \right |\\ &= \left |\frac {1}{x-x_0} \int _{x_0}^x f(t) \thinspace {\rm {d}} t -\frac {1}{x-x_0}\int _{x_0}^x f(x_0)\thinspace {\rm {d}} t \right |\\ &= \left |\frac {1}{x-x_0} \int _{x_0}^x (f(t)-f(x_0)) \thinspace {\rm {d}} t \right |\\ &\leq \frac {1}{x-x_0} \int _{x_0}^x |f(t)-f(x_0)| \thinspace {\rm {d}} t \leq \frac {1}{x-x_0} \int _{x_0}^x \varepsilon \thinspace {\rm {d}} t\ =\ \varepsilon .\end{aligned}
[/latex]

Analog gilt für [latex]x \in (x_0 -\delta ,x_0)\cap [a,b] \setminus \left \lbrace {x_0} \right \rbrace[/latex], dass

[latex]
\begin{aligned}[]\left | \frac {F(x)-F(x_0)}{x-x_0} -f(x_0) \right | &= \left | \frac {F(x_0)-F(x)}{x_0-x} -f(x_0) \right | = \left |\frac {1}{x_0-x} \int _{x}^{x_0} (f(t)-f(x_0)) \thinspace {\rm {d}} t \right |\\ &\leq \frac {1}{x_0-x} \int _{x}^{x_0} |f(t)-f(x_0)| \thinspace {\rm {d}} t \leq \varepsilon .\end{aligned}
[/latex]

Da [latex]\varepsilon >0[/latex] beliebig war, beweist dies [latex]\lim _{x\to x_0}\frac {F(x)-F(x_0)}{x-x_0}=f(x_0)[/latex] und damit den Satz. ∎

image

Abbildung 8.1 – Illustration zum Beweis von Theorem 8.2. Der Wert [latex]F(x)-F(x_0)[/latex] lässt sich schreiben als [latex]f(x_0)(x-x_0)[/latex] plus die Fläche in Rot, die kleiner ist als [latex]\varepsilon (x-x_0)[/latex]. Somit ist [latex]\frac {F(x)-F(x_0)}{x-x_0}[/latex], bis auf einen Fehler kleiner als [latex]\varepsilon[/latex], durch [latex]f(x_0)[/latex] gegeben.

Korollar 8.3: Ableitung des Integrals

Sei [latex][a,b][/latex] ein kompaktes Intervall mit Endpunkten [latex]a [latex]
\begin{aligned}[]\label{eq;hauptsatzformel} F(x) = \int _a^x f(t) \thinspace {\rm {d}} t + C\end{aligned}
[/latex]
für alle [latex]x \in [a,b][/latex] und eine Konstante [latex]C \in \mathbb {R}[/latex].

Beweis

Nach Satz 4.42 ist [latex]f[/latex] Riemann-integrierbar. Nach Theorem 8.2 ist

[latex]
\begin{aligned}[]x \in [a,b] \mapsto \int _a^x f(t) \thinspace {\rm {d}} t\end{aligned}
[/latex]

differenzierbar und eine Stammfunktion von [latex]f[/latex]. Nach Lemma 7.59 unterscheidet sich jede weitere Stammfunktion nur um eine Konstante von dieser, was die Formel im Korollar beweist. ∎

Korollar 8.4: Berechnung des Integrals

Sei [latex][a,b][/latex] ein kompaktes Intervall mit Endpunkten [latex]a

[latex]
\begin{aligned}[]\int _a^bf(t)\thinspace {\rm {d}} t=F(b)-F(a).\end{aligned}
[/latex]

Wir werden auch öfter die Abkürzung [latex]\big [F(x)\big ]_a^b=\big [F(x)\big ]_{x=a}^{x=b}=F(b)-F(a)[/latex] verwenden.

Beweis

Nach Korollar 8.3 gilt (8.1). Setzen wir [latex]x=a[/latex] erhalten wir für die Konstante [latex]C=F(a)[/latex] und das Korollar folgt nun indem wir [latex]x=b[/latex] in (8.1) einsetzen. ∎

Korollar 8.5: Integral der Ableitung

Sei [latex][a,b][/latex] ein kompaktes Intervall mit Endpunkten [latex]a

[latex]
\begin{aligned}[]F(x) = F(a) + \int _a^x F'(t) \thinspace {\rm {d}} t\end{aligned}
[/latex]

für alle [latex]x \in [a,b][/latex].

Insbesondere kann der Wert von [latex]F[/latex] an jeder Stelle [latex]x[/latex] vollständig durch den Wert von [latex]F[/latex] bei der Stelle [latex]a[/latex] und die Werte der Funktion [latex]F'[/latex] beschrieben werden.

Beweis

Die Funktion [latex]F_0: x \in [a,b] \mapsto \int _a^x F'(t) \thinspace {\rm {d}} t[/latex] ist nach Korollar 8.3 ebenso wie [latex]F[/latex] eine Stammfunktion. Sei [latex]C \in \mathbb {R}[/latex] mit [latex]F(x)= F_0(x)+C[/latex] für alle [latex]x \in [a,b][/latex] nach Korollar 8.3. Setzen wir [latex]x =a[/latex] ein, so erhalten wir

[latex]
\begin{aligned}[]F(a) = F_0(a) + C = \int _a^a F'(t) \thinspace {\rm {d}} t + C = C.\end{aligned}
[/latex]

Somit ist

[latex]
\begin{aligned}[]F(x) = F(a)+ F_0(x) = F(a) + \int _a^x F'(t) \thinspace {\rm {d}} t\end{aligned}
[/latex]

für alle [latex]x \in [a,b][/latex], wie zu zeigen war. ∎

Die obigen Resultate sind auch als der Fundamentalsatz der Integral- und Differentialrechnung bekannt und gehen auf die Arbeiten von Leibniz, Newton und Barrow zurück, die weitgehend die Ausgangspunkte der Analysis darstellen. Zur Vereinfachung der Diskussion haben wir Stetigkeit von [latex]f[/latex] beziehungsweise [latex]F'[/latex] angenommen. Diese Annahme lässt sich auf verschiedene Arten abschwächen.

Übung 8.6: Theorem 8.2 für «fast überall» stetige Funktionen

Sei [latex][a,b][/latex] ein kompaktes Intervall mit Endpunkten [latex]a

Wichtige Übung 8.7: Fundamentalsatz für komplexe Funktionen

In obigem Satz und dessen Korollaren haben wir für die Aussagen eigentlich nicht wirklich verwendet, dass die betrachteten Funktionen reellwertig sind. Zeigen Sie deswegen, dass alle obigen Resultate für komplexwertige Funktionen zutreffen.

Hinweis.

Sie brauchen nicht nochmals die Beweise durchzuführen. Für eine Riemann-integrierbare Funktion [latex]g:[a,b] \to \mathbb {C}[/latex] lässt sich das Riemann-Integral von [latex]g[/latex] in das Integral über den Realteil und das Integral über den Imaginärteil von [latex]g[/latex] zerlegen.

Zusammenfassend wollen wir noch betonen, dass auf Grund obiger Resultate die Berechnung von Riemann-Integralen und damit insbesondere von Flächeninhalten, Schwerpunkkoordinaten, Arbeitsberechnungen, Bogenlängen und Volumen von Rotationskörpern (siehe Abschnitt 8.3) auf die Berechnung von Stammfunktionen zurückgeführt werden kann.

Übung 8.8: Mittelwertsatz der Integralrechnung

Sei [latex][a,b][/latex] ein kompaktes Intervall mit Endpunkten [latex]a

[latex]
\begin{aligned}[]\int _a^b f(x) \thinspace {\rm {d}} x = f(\xi ) (b-a).\end{aligned}
[/latex]

Können Sie zwei verschiedene Beweise finden?

Hinweis.

Verwenden Sie entweder den Zwischenwertsatz für stetige Funktionen oder den Mittelwertsatz der Differentialrechnung.

Übung 8.9: Riemann-integrierbare Ableitung

Wir möchten hier eine etwas stärkere Version von Korollar 8.5 thematisieren. Sei [latex][a,b][/latex] ein kompaktes Intervall mit Endpunkten [latex]a

Hinweis.

Betrachten Sie eine Zerlegung [latex]\mathfrak {Z} = \left \lbrace {a=x_0

Applet 8.10: Fundamentalsatz

Wir sehen hier nochmals die Idee des Beweises des Fundamentalsatzes der Integral- und Differentialrechnung (Theorem 8.2), wobei unten das partikuläre Integral der Funktion im oberen Fenster dargestellt wird.

8.1.1 – Differentiation von Potenzreihen

Wir wenden nun den Fundamentalsatz an, um zu zeigen, dass sich Potenzreihen nicht nur integrieren (siehe Satz 6.85), sondern auch differenzieren lassen.

Korollar 8.11: Differentiation von Potenzreihen

Sei [latex]f(x) = \sum _{n=0}^\infty c_n x^n[/latex] eine Potenzreihe mit Konvergenzradius [latex]R > 0[/latex]. Dann ist [latex]f:(-R,R) \to \mathbb {R}[/latex] differenzierbar und es gilt

[latex]
\begin{aligned}[]f'(x) = \sum _{n=1}^\infty nc_n x^{n-1}\end{aligned}
[/latex]

für alle [latex]x \in (-R,R)[/latex], wobei die Potenzreihe rechts ebenfalls Konvergenzradius [latex]R[/latex] hat.

Beweis

Sei [latex]g(x)[/latex] die Potenzreihe [latex]\sum _{n=1}^\infty nc_n x^{n-1} = \sum _{k=0}^\infty (k+1)c_{k+1} x^{k}[/latex] und sei [latex]S[/latex] der zugehörige Konvergenzradius. Nach Satz 6.56 ist [latex]g[/latex] auf [latex](-S,S)[/latex] stetig und nach Satz 6.85 darf [latex]g(x)[/latex] gliedweise integriert werden. Genauer gesagt besagt Satz 6.85, dass [latex]G: (-S,S) \to \mathbb {R}[/latex] gegeben durch

[latex]
\begin{aligned}[]G(x) = \int _0^x g(t) \thinspace {\rm {d}} t = \sum _{k=0}^\infty \frac {(k+1)c_{k+1}}{k+1} x^{k+1} = \sum _{k=0}^\infty c_{k+1} x^{k+1}\end{aligned}
[/latex]

wieder eine Potenzreihe mit Konvergenzradius [latex]S[/latex] darstellt. Nach Korollar 8.3 gilt weiters, dass [latex]G[/latex] eine Stammfunktion von [latex]g=G'[/latex] darstellt. Da dies aber abgesehen vom ersten Glied der Reihe genau die Potenzreihe [latex]f[/latex] ist, stimmen die Konvergenzradien [latex]R=S[/latex] überein, [latex]f[/latex] ist differenzierbar und [latex]f'(x) = G'(x) = g(x) = \sum _{n=1}^\infty nc_n x^{n-1}[/latex] für alle [latex]x \in (-R,R)[/latex]. ∎

Korollar 8.11 lässt sich leicht anpassen, um zu zeigen, dass durch Potenzreihen gegebene Funktionen glatt sind.

Übung 8.12: Verallgemeinerung von Korollar 8.11

Sei [latex]f(x) = \sum _{n=0}^\infty c_n x^n[/latex] eine Potenzreihe mit Konvergenzradius [latex]R > 0[/latex]. Zeigen Sie, dass [latex]f:(-R,R) \to \mathbb {R}[/latex] glatt ist und finden Sie eine Darstellung von [latex]f^{(n)}[/latex] durch eine Potenzreihe für jedes [latex]n \in \mathbb {N}[/latex].

Übung 8.13: Koeffizientenvergleich für Potenzreihen

Seien [latex]f(x) = \sum _{n=0}^\infty c_n x^n[/latex] und [latex]g(x) = \sum _{n=0}^\infty d_n x^n[/latex] Potenzreihen mit reellen Koeffizienten und positiven Konvergenzradien [latex]R_f,R_g[/latex]. Sei [latex]R = \min \left \lbrace {R_f,R_g} \right \rbrace[/latex] und angenommen [latex]f(x) = g(x)[/latex] für alle [latex]x \in (-R,R)[/latex]. Zeigen Sie, dass dann [latex]c_n = d_n[/latex] für alle [latex]n \in \mathbb {N}_0[/latex] und damit [latex]R_f=R_g[/latex] gilt.

Hinweis.

Verwenden Sie Korollar 8.11 und Übung 8.12.

Übung 8.14: Potenzreihenentwicklung für Wurzeln

Sei [latex]\alpha \in \mathbb {C}[/latex]. Wir wollen hier zeigen, dass

[latex]
\begin{aligned}[](1+x)^\alpha = \sum _{n=0}^\infty \binom {\alpha }{n} x^n\end{aligned}
[/latex]

für alle [latex]x \in (-1,1)[/latex], wobei die verallgemeinerten Binomialkoeffizienten für [latex]n \in \mathbb {N}_0[/latex] durch

[latex]
\begin{aligned}[]\binom {\alpha }{n} = \frac {\prod _{j=0}^{n-1} (\alpha -j)}{n!} = \frac {\alpha (\alpha -1)\cdots (\alpha -n+1)}{n!}\end{aligned}
[/latex]

definiert sind.

  1. Zeigen Sie, dass für [latex]\alpha \not \in \mathbb {N}_0[/latex] die Potenzreihe
    [latex]
    \begin{aligned}[]g(x) = \sum _{n=0}^\infty \binom {\alpha }{n} x^n\end{aligned}
    [/latex]

    Konvergenzradius [latex]1[/latex] hat.

  2. Berechnen Sie die Ableitung von [latex]g[/latex] und zeigen Sie, dass [latex]f(x) = (1+x)^\alpha[/latex] und [latex]g(x)[/latex] die Differentialgleichung
    [latex]
    \begin{aligned}[]y' = \alpha \frac {y}{1+x}\end{aligned}
    [/latex]

    erfüllen.

  3. Berechnen Sie die Ableitung von [latex]\frac {g(x)}{f(x)}[/latex] und schliessen Sie die Behauptung.

Hinweis.

Für (b) müssen Sie die Gleichung

[latex]
\begin{aligned}[](n+1) \binom {\alpha }{n+1} + n \binom {\alpha }{n} = \alpha \binom {\alpha }{n}\end{aligned}
[/latex]

zeigen (was mit der richtigen Interpretation des Produkts [latex]\prod _{j=0}^{n-1} (\alpha -j)[/latex] auch für [latex]n=0[/latex] gilt).

8.1.2 – Die alternierende harmonische Reihe

Wir haben bereits in Beispiel 6.20 gesehen, dass die alternierende harmonische Reihe [latex]\sum _{n=1}^\infty \frac {(-1)^{n+1}}{n}[/latex] konvergiert, was auch aus dem Leibniz-Kriterium (Proposition 6.25) folgt. Mit den Resultaten von Kapitel 6 konnten wir den Wert der Reihe aber nicht bestimmen. Nun können wir mit Hilfe des Fundamentalsatzes der Integral- und Differentialrechnung die Identität

[latex]
\begin{aligned}[]\sum _{n=1}^\infty \frac {(-1)^{n+1}}{n} = \log (2)\end{aligned}
[/latex]

beweisen.

Es ist in diesem Fall (vielleicht überraschenderweise) einfacher, eine allgemeinere Aussage zu zeigen. Wir beginnen hierfür mit

[latex]
\begin{aligned}[](\log (1+x))' = \frac {1}{1+x} = \frac {1}{1-(-x)} = \sum _{n=0}^\infty (-1)^n x^n\end{aligned}
[/latex]

für alle [latex]x\in (-1,1)[/latex], wobei die Reihe rechts Konvergenzradius [latex]1[/latex] hat (und an den Endpunkten divergiert). Nach Korollar 8.5 und Satz 6.85 folgt daraus, dass

[latex]
\begin{aligned}[]\log (1+x) = \log (1) + \int _0^x \frac {1}{1+t} \thinspace {\rm {d}} t = \sum _{n=0}^\infty \frac {(-1)^{n}}{n+1}x^{n+1} = \sum _{k=1}^\infty \frac {(-1)^{k+1}}{k}x^k\end{aligned}
[/latex]

für alle [latex]x \in (-1,1)[/latex]. Da die Potenzreihe rechts auch für [latex]x=1[/latex] konvergiert, ist die Funktion

[latex]
\begin{aligned}[]f: x \in (-1,1] \mapsto \sum _{k=1}^\infty \frac {(-1)^{k+1}}{k}x^k\end{aligned}
[/latex]

nach dem Abelschen Grenzwertsatz (Satz 6.65) auch bei [latex]x=1[/latex] stetig. Da die Funktion mit Definitionsbereich [latex](-1,1][/latex] definiert durch [latex]\log (1+x)[/latex] für [latex]x\in (-1,1][/latex] ebenfalls stetig ist und für [latex]x \in (-1,1)[/latex] mit [latex]f(x)[/latex] übereinstimmt, ist

[latex]
\begin{aligned}[]\log (2) = \lim _{x \nearrow 1}\log (1+x) = \lim _{x \nearrow 1} f(x) = \sum _{n=1}^\infty \frac {(-1)^{n+1}}{n}.\end{aligned}
[/latex]

8.1.3 – Die Leibniz-Reihe

Wir verwenden obige Methode nochmals, um
[latex]
\begin{aligned}[]\label{eq:fund-leibnizreihe} \sum _{n=0}^\infty \frac {(-1)^n}{2n+1} = \frac {\pi }{4}\end{aligned}
[/latex]
zu beweisen.

Nach dem Leibniz-Kriterium (Satz 6.25) ist die Reihe [latex]\sum _{n=0}^\infty \frac {(-1)^n}{2n+1}[/latex] konvergent. Wir beginnen die Berechnung ihres Wertes mit

[latex]
\begin{aligned}[]\arctan '(x) = \frac {1}{1+x^2} = \sum _{n=0}^\infty (-1)^n x^{2n}.\end{aligned}
[/latex]

für [latex]x \in (-1,1)[/latex]. Nach Korollar 8.5 gilt

[latex]
\begin{aligned}[]\arctan (x) = \arctan (0) + \int _0^x \frac {1}{1+t^2} \thinspace {\rm {d}} t = \sum _{n=0}^\infty \frac {(-1)^n}{2n+1}x^{2n+1}\end{aligned}
[/latex]

Stetigkeit des Arkustangens bei [latex]x=1[/latex], der Abelsche Grenzwertsatz (Satz 6.65) und die Identität [latex]\arctan (1) = \frac {\pi }{4}[/latex] beweisen nun (8.2).

8.2 – Integrationsmethoden

Wir erinnern daran, dass das unbestimmte Integral einer Funktion [latex]f[/latex] in der Variablen [latex]x[/latex] der Ausdruck

[latex]
\begin{aligned}[]\int f(x) \thinspace {\rm {d}} x = F(x)+C\end{aligned}
[/latex]

ist, wobei [latex]F[/latex] eine Stammfunktion von [latex]f[/latex] ist. In den Abschnitten 7.5.2 und 7.4 haben wir bereits einige Regeln zur Berechnung konkreter unbestimmter Integrale kennengelernt: Für [latex]s \in \mathbb {R}[/latex] (oder sogar [latex]s\in \mathbb {C}[/latex]) ist

[latex]
\begin{aligned}[]\int x^s \thinspace {\rm {d}} x = \left \lbrace \begin{array}{cl} \frac {1}{s+1}x^{s+1} + C & \text {falls } s \neq {-1} \\ \log |x|+ C & \text {falls } s=-1\end{array} \right .\end{aligned}
[/latex]

und

[latex]
\begin{aligned}[]\int \exp (x) \thinspace {\rm {d}} x &= \exp (x) + C\\ \int \cos (x) \thinspace {\rm {d}} x &= \sin (x) + C\\ \int \sin (x) \thinspace {\rm {d}} x &= - \cos (x) + C\\ \int \sinh (x) \thinspace {\rm {d}} x &= \cosh (x) + C\\ \int \cosh (x) \thinspace {\rm {d}} x &= \sinh (x) + C\\ \int \frac {1}{\sqrt {1-x^2}} \thinspace {\rm {d}} x &= \arcsin (x) + C\\ \int \frac {1}{1+x^2} \thinspace {\rm {d}} x &= \arctan (x) + C \\ \int \frac {1}{\sqrt {1+x^2}}\thinspace {\rm {d}} x&=\operatorname {arsinh}(x)+C\\ \int \frac {1}{\sqrt {x^2-1}}\thinspace {\rm {d}} x&=\operatorname {arcosh}(x)+C.\end{aligned}
[/latex]

Des Weiteren gilt für Funktionen [latex]f_1,f_2[/latex] in der Variable [latex]x[/latex] und Zahlen [latex]\alpha _1,\alpha _2 \in \mathbb {R}[/latex]

[latex]
\begin{aligned}[]\int \alpha _1 f_1(x) + \alpha _2 f_2(x) \thinspace {\rm {d}} x = \alpha _1 \int f_1(x) \thinspace {\rm {d}} x + \alpha _2 \int f_2(x) \thinspace {\rm {d}} x.\end{aligned}
[/latex]

Denn falls [latex]F_1[/latex] eine Stammfunktion von [latex]f_1[/latex] ist und [latex]F_2[/latex] eine Stammfunktion von [latex]f_2[/latex] ist, so muss die Funktion [latex]\alpha _1 F_1 + \alpha _2 F_2[/latex] auf Grund der Linearität der Ableitung (Proposition 7.5) eine Stammfunktion von [latex]\alpha _1 f_1 + \alpha _2 f_2[/latex] sein.

Auf ähnliche Weise lassen sich die anderen Regeln der Differentiation als Identitäten für unbestimmte Integrale auffassen, wie wir nun ausführen wollen.

8.2.1 – Partielle Integration

Die Produktregel in Proposition 7.5

[latex]
\begin{aligned}[](uv)' = u'v+uv'\end{aligned}
[/latex]

für zwei differenzierbare Funktionen [latex]u,v[/latex] führt ebenso zu einer Integrationsregel, nämlich der partiellen Integration
[latex]
\begin{aligned}[]\label{eq:abl-partial integration} uv + C &= \int (u'v+uv') \thinspace {\rm {d}} x \nonumber \\ \int uv' \thinspace {\rm {d}} x &= uv - \int u'v \thinspace {\rm {d}} x + C.\end{aligned}
[/latex]
In der Leibniz-Notation ist [latex]v' = \frac {\thinspace {\rm {d}} v}{\thinspace {\rm {d}} x}[/latex] und [latex]u' = \frac {\thinspace {\rm {d}} u}{\thinspace {\rm {d}} x}[/latex]. Deswegen schreibt man die partielle Integration oft auch als

[latex]
\begin{aligned}[]\int u \thinspace {\rm {d}} v = uv - \int v \thinspace {\rm {d}} u + C,\end{aligned}
[/latex]

was formal bloss als Kurzform der Formel in (8.3) verstanden werden sollte. Die Regel der partiellen Integration ist bereits ein Beispiel, wo eine einfache Regel des Differenzierens eine komplexere Regel des Integrierens als Entsprechung hat. Die Produktregel erlaubt uns, die Ableitung jedes Produkts mittels der Ableitung dessen Faktoren auszudrücken. Die partielle Integration hingegen erlaubt uns, das unbestimmte Integral eines Produkts mittels dem Integral eines Faktors und eines weiteren Integrals auszudrücken. Mit etwas Glück (und Geschick) ist das zweite Integral einfacher und kann anschliessend berechnet werden. Wir demonstrieren dies anhand zweier Beispiele.

Beispiel 8.15: Beispiele partieller Integration

  1. Wir berechnen das unbestimmte Integral [latex]\int x \exp (x) \thinspace {\rm {d}} x[/latex]. Dafür setzen wir [latex]u(x) = x[/latex] und [latex]v'(x) = \exp (x)[/latex]. Eine Stammfunktion von [latex]v'[/latex] ist [latex]v(x) = \exp (x)[/latex]. Damit erhalten wir
    [latex]
    \begin{aligned}[]\int x \exp (x) \thinspace {\rm {d}} x = x \exp (x) - \int 1 \cdot \exp (x) \thinspace {\rm {d}} x + C = x\exp (x) - \exp (x) + C.\end{aligned}
    [/latex]

    Wir bemerken, dass es genügt, in solchen Berechnungen immer bloss eine unbekannte Integrationskonstante [latex]C[/latex] zu verwenden, da mehrere solche einfach zusammengefasst werden können. (Kontrollieren Sie diese Rechnung durch Ableiten.) Dieselbe Berechnungsmethode führt auch für unbestimmte Integrale der Form [latex]\int x^n \exp (x) \thinspace {\rm {d}} x[/latex], [latex]\int x^n \sin (x) \thinspace {\rm {d}} x[/latex] und [latex]\int x^n \cos (x) \thinspace {\rm {d}} x[/latex] für [latex]n \in \mathbb {N}[/latex] zum Erfolg.

  2. Wir wollen das unbestimmte Integral [latex]\int \log (x) \thinspace {\rm {d}} x[/latex] berechnen. Es mag zuerst etwas überra­schend sein, dass wir dazu partielle Integration verwenden wollen. Sei [latex]u(x) = \log (x)[/latex] und [latex]v' = 1[/latex]. Dann ist [latex]v(x) = x[/latex] eine Stammfunktion von [latex]v'[/latex], womit
    [latex]
    \begin{aligned}[]\int \log (x) \thinspace {\rm {d}} x &= \int \log (x) \cdot 1 \thinspace {\rm {d}} x = \log (x)\cdot x - \int \frac {1}{x} x \thinspace {\rm {d}} x + C = x\log (x) - \int 1 \thinspace {\rm {d}} x + C\\ &= x\log (x)-x + C.\end{aligned}
    [/latex]

    Dies kann man wiederum durch Ableiten verifizieren (was nicht notwendig ist, aber einen sehr einfachen Test darstellt).

Übung 8.16

  1. Berechnen Sie [latex]\int x^2\sin (x)\thinspace {\rm {d}} x[/latex].
  2. Geben Sie eine rekursive Formel zur Berechnung von [latex]\int x^n \exp (x) \thinspace {\rm {d}} x[/latex], [latex]\int x^n \sin (x) \thinspace {\rm {d}} x[/latex] und [latex]\int x^n \cos (x) \thinspace {\rm {d}} x[/latex] für [latex]n \in \mathbb {N}[/latex] an.
  3. Berechnen Sie [latex]\int x^s\log (x)\thinspace {\rm {d}} x[/latex] für jedes [latex]s\in \mathbb {R}[/latex]. Beachten Sie hierbei, dass der Fall [latex]s=-1[/latex] getrennt zu behandeln ist.
  4. Berechnen Sie das unbestimmte Integral [latex]\int e^{ax}\sin (bx) \thinspace {\rm {d}} x[/latex] für [latex]a,b \in \mathbb {R} \setminus \left \lbrace {0} \right \rbrace[/latex].

    Hinweis.

    Bei (iv) ergibt sich nach zweifacher partieller Integration ein Gleichungssystem.

8.2.2 – Substitution

Falls eine Funktion [latex]g[/latex] auf einem Intervall [latex]I_u[/latex] das unbestimmte Integral [latex]\int g(u) \thinspace {\rm {d}} u = G(u) + C[/latex] besitzt und [latex]f: I_x \to I_u[/latex] eine stetig differenzierbare Abbildung auf dem Intervall [latex]I_x[/latex] ist, dann gilt

[latex]
\begin{aligned}[]\int (g \circ f)(x) f'(x) \thinspace {\rm {d}} x = (G \circ f)(x) + C\end{aligned}
[/latex]

auf [latex]I_x[/latex]. Dies folgt unmittelbar aus der Kettenregel in Satz 7.8 und wird oft auch geschrieben als
[latex]
\begin{aligned}[]\label{eq:fund-substitution} \int (g \circ f)(x) f'(x) \thinspace {\rm {d}} x = \int g(u) \thinspace {\rm {d}} u\end{aligned}
[/latex]
für die «neue Variable» [latex]u = f(x)[/latex]. Alternativ werden wir die obige Substitutionsregel gemeinsam mit der Leibniz-Notation auch in folgender informellen Schreibweise

[latex]
\begin{aligned}[]\int (g \circ f)(x) f'(x) \thinspace {\rm {d}} x = \int (g \circ f)(x) \frac {\thinspace {\rm {d}} u}{\thinspace {\rm {d}} x} \thinspace {\rm {d}} x = \int g(u) \thinspace {\rm {d}} u = G(u) + C = G(f(x)) + C,\end{aligned}
[/latex]

verwenden, wobei [latex]u = f(x)[/latex] und [latex]\thinspace {\rm {d}} u = f'(x) \thinspace {\rm {d}} x[/latex].

Beispiel 8.17

  1. Es gilt
    [latex]
    \begin{aligned}[]\int \frac {x}{1+x^2}\thinspace {\rm {d}} x =\int \underbrace {\frac {1}{1+x^2}}_{=\frac 1u}\underbrace {x\thinspace {\rm {d}} x}_{=\frac 12\thinspace {\rm {d}} u} = \frac {1}{2} \int \frac {1}{u}\thinspace {\rm {d}} u = \frac {1}{2} \log |u| = \frac 12 \log (1+x^2) + C,\end{aligned}
    [/latex]

    wobei [latex]u = 1+x^2[/latex] gesetzt wurde, womit [latex]\thinspace {\rm {d}} u = 2x \thinspace {\rm {d}} x[/latex].

  2. Es gilt
    [latex]
    \begin{aligned}[]\int \frac {1}{\sin (x) }\thinspace {\rm {d}} x &= \int \frac {1}{2\sin \left (\frac {x}{2}\right )\cos \left (\frac {x}{2}\right )}\thinspace {\rm {d}} x = \int \frac {1}{\tan (u) \cos ^2(u)} \thinspace {\rm {d}} u= \int \frac {1}{v}\thinspace {\rm {d}} v = \log |v| + C\\ &= \log \left | \tan \left (\frac {x}{2}\right ) \right |+ C,\end{aligned}
    [/latex]

    wobei [latex]u = \frac {x}{2}[/latex], [latex]\thinspace {\rm {d}} u = \frac {1}{2}\thinspace {\rm {d}} x[/latex], und [latex]v = \tan (u)[/latex], [latex]\thinspace {\rm {d}} v = \frac {1}{\cos ^2(u)}\thinspace {\rm {d}} u[/latex].

Wie bereits erwähnt, benötigt die Integration mehr Übung und Vorraussicht als die Differentiation. Obige Substitutionen benötigen zum Beispiel den Blick ob gewisse Faktoren vielleicht die gewünschte Ableitung [latex]f'[/latex] einer inneren Funktionen [latex]f[/latex] darstellen könnte. Manchmal ist dies naheliegend wie in Beispiel 8.17(i), doch manchmal erfordert dies Erfahrung und eine längere Suche wie in Beispiel 8.17(ii).

8.2.3 – Integration rationaler Funktionen

Wir erinnern daran, dass eine rationale Funktion eine Funktion der Form [latex]x \mapsto \frac {p(x)}{q(x)}[/latex] für Polynome [latex]p(t), q(t) \in \mathbb {R}[t][/latex] und [latex]q(t) \neq 0[/latex] ist, wobei der Definitionsbereich [latex]\mathbb {R}[/latex] ohne die Nullstellen von [latex]q(t)[/latex] ist. Wir wollen hier ein Verfahren zur Berechnung des unbestimmten Integrals einer rationalen Funktion besprechen. Nach Division mit Rest für Polynome (siehe Übung 3.17) können wir als erstes ein Polynom abspalten, so dass die verbleibende rationale Funktion von der Form [latex]\frac {p_1(x)}{q(x)}[/latex] für [latex]\deg (p_1)

Da Polynome mittels der Formel [latex]\int x^n \thinspace {\rm {d}} x = \frac {1}{n+1} x^{n+1} + C[/latex] integriert werden können, nehmen wir nun an, dass der Grad von [latex]p[/latex] kleiner als der Grad von [latex]q[/latex] ist. Wir betrachten zuerst einige Spezialfälle.

Beispiel 8.18: Integration von elementaren rationalen Funktionen

Sei [latex]a \in \mathbb {R}[/latex] beliebig und [latex]n \geq 2[/latex] eine natürliche Zahl.

  1. Es gilt
    [latex]
    \begin{aligned}[]\int \frac {1}{x-a} \thinspace {\rm {d}} x = \int \frac {1}{u}\thinspace {\rm {d}} u = \log |u|+C = \log |x-a| + C,\end{aligned}
    [/latex]

    wobei [latex]u=x-a[/latex] gesetzt wurde und [latex]\thinspace {\rm {d}} u = \thinspace {\rm {d}} x[/latex] ist.

  2. Für [latex]n\geq 2[/latex] gilt
    [latex]
    \begin{aligned}[]\int \frac {1}{(x-a)^n} \thinspace {\rm {d}} x = \int u^{-n}\thinspace {\rm {d}} u = \frac {1}{-n+1} u^{-n+1}+C = \frac {1}{-n+1} (x-a)^{-n+1}+C,\end{aligned}
    [/latex]

    wobei wieder [latex]u=x-a[/latex] gesetzt wurde und [latex]\thinspace {\rm {d}} u = \thinspace {\rm {d}} x[/latex] ist.

  3. Falls [latex]a \neq 0[/latex] ist, so gilt
    [latex]
    \begin{aligned}[]\int \frac {1}{a^2+x^2} \thinspace {\rm {d}} x &= \frac {1}{a^2} \int \frac {1}{1+\left (\frac {x}{a}\right )^2}\thinspace {\rm {d}} x = \frac {1}{a} \int \frac {1}{1+u^2}\thinspace {\rm {d}} u = \frac {1}{a} \arctan (u) +C\\ &= \frac {1}{a} \arctan \left (\frac {x}{a}\right ) +C,\end{aligned}
    [/latex]

    wobei [latex]u = \frac {x}{a}[/latex] gesetzt wurde und [latex]\thinspace {\rm {d}} u = \frac {1}{a}\thinspace {\rm {d}} x[/latex] ist.

  4. Es gilt
    [latex]
    \begin{aligned}[]\int \frac {x}{a^2+x^2} \thinspace {\rm {d}} x = \frac {1}{2} \int \frac {1}{u} \thinspace {\rm {d}} u = \frac 12\log |u| +C = \frac 12 \log (a^2+x^2) + C,\end{aligned}
    [/latex]

    wobei [latex]u = a^2+x^2[/latex] und [latex]\thinspace {\rm {d}} u = 2x \thinspace {\rm {d}} x[/latex].

  5. Für [latex]n\geq 2[/latex] gilt
    [latex]
    \begin{aligned}[]\int \frac {x}{(a^2+x^2)^n} \thinspace {\rm {d}} x = \frac {1}{2}\int \frac {1}{u^n} \thinspace {\rm {d}} u = \frac {1}{2(1-n)} u^{1-n} + C = \frac {1}{2(1-n)} (a^2+x^2)^{1-n} + C,\end{aligned}
    [/latex]

    wobei [latex]u = a^2+x^2[/latex] und [latex]\thinspace {\rm {d}} u = 2x \thinspace {\rm {d}} x[/latex].

Im Allgemeinen verwenden wir die sogenannte Partialbruchzerlegung von [latex]\frac {p(x)}{q(x)}[/latex], um die Integration rationaler Funktionen auf obige Beispiele zurückzuführen. In der Tat lässt sich [latex]\frac {p(x)}{q(x)}[/latex] als Linearkombination von einfacheren rationalen Funktionen darstellen. Diese sind von der Form

[latex]
\begin{aligned}[]\frac {1}{(x-a)},\quad \frac {1}{(x-a)^2},\quad \ldots \quad ,\quad \frac {1}{(x-a)^k},\end{aligned}
[/latex]

oder von der Form

[latex]
\begin{aligned}[]\frac {A_1 x + B_1}{(x-\lambda )(x-\overline {\lambda })},\quad \ldots \quad ,\quad \frac {A_\ell x + B_\ell }{\big ((x-\lambda )(x-\overline {\lambda })\big )^\ell },\end{aligned}
[/latex]

wobei [latex]a\in \mathbb {R}[/latex] eine Nullstelle von [latex]q[/latex] mit Vielfachheit [latex]k[/latex] und [latex]\lambda \in \mathbb {C} \setminus \mathbb {R}[/latex] eine Nullstelle von [latex]q[/latex] mit Vielfachheit [latex]\ell \in \mathbb {N}[/latex] ist und [latex]A_1,B_1,\ldots ,A_\ell ,B_\ell \in \mathbb {R}[/latex].

Beispiel 8.19: Integration rationaler Funktionen und die Partialbruchzerlegung

  1. Wir wollen das unbestimmte Integral [latex]\int \frac {x^4+1}{x^2(x+1)}\thinspace {\rm {d}} x[/latex] bestimmen. Als erstes führen wir Division mit Rest
    [latex]
    \begin{aligned}[](&x^4+1):(x^3+x^2) = x-1\\ -&\underline{x^4-x^3}\\ & \quad -x^3+1\\ & \quad \, \quad \underline{x^3+x^2}\\ & \qquad \ \qquad x^2+1 \quad \mathrm {Rest}\end{aligned}
    [/latex]

    durch, womit

    [latex]
    \begin{aligned}[]\int \frac {x^4+1}{x^2(x+1)}\thinspace {\rm {d}} x = \int \left (x-1 + \frac {x^2+1}{x^2(x+1)}\right )\thinspace {\rm {d}} x = \frac {x^2}{2}-x + \int \frac {x^2+1}{x^2(x+1)}\thinspace {\rm {d}} x.\end{aligned}
    [/latex]

    Um die Partialbruchzerlegung von [latex]\frac {x^2+1}{x^2(x+1)}[/latex] zu erhalten, setzen wir

    [latex]
    \begin{aligned}[]\frac {x^2+1}{x^2(x+1)} = \frac {A}{x^2} + \frac {B}{x} + \frac {C}{x+1}\end{aligned}
    [/latex]

    für noch unbekannte Zahlen [latex]A,B,C \in \mathbb {R}[/latex], multiplizieren mit [latex]x^2(x+1)[/latex] und erhalten

    [latex]
    \begin{aligned}[]x^2+1 = A(x+1) + Bx(x+1)+Cx^2.\end{aligned}
    [/latex]

    Nun setzen wir in diesem [latex]x=0[/latex] um [latex]A=1[/latex] zu erhalten und [latex]x=-1[/latex], um [latex]C=2[/latex] zu erhalten. Für [latex]x=1[/latex] ergibt sich nun [latex]2 = 1 \cdot 2 + B \cdot 2 + 2 \cdot 1[/latex] und somit [latex]B=-1[/latex]. (Alternativ kann man auch beide Seiten ausmultiplizieren, die Koeffizienten links und rechts vergleichen, und auf diese Weise drei Gleichungen in den unbekannten Variablen [latex]A,B,C[/latex] erhalten.) Daher ist

    [latex]
    \begin{aligned}[]\int \frac {x^2+1}{x^2(x+1)}\thinspace {\rm {d}} x &= \int \frac {1}{x^2}\thinspace {\rm {d}} x - \int \frac {1}{x} \thinspace {\rm {d}} x + 2 \int \frac {1}{x+1}\thinspace {\rm {d}} x \\ &= - \frac {1}{x} - \log |x| + 2\log |x+1| + D\end{aligned}
    [/latex]
  2. Wir berechnen das unbestimmte Integral [latex]\int \frac {1}{x(x^2+2x+2)} \thinspace {\rm {d}} x[/latex]. Man beachte dabei, dass das Polynom [latex]x^2+2x+2[/latex] keine reellen Nullstellen hat. Für die Partialbruchzerlegung machen wir den Ansatz
    [latex]
    \begin{aligned}[]\frac {1}{x(x^2+2x+2)} = \frac {A}{x} + \frac {Bx+C}{x^2+2x+2}.\end{aligned}
    [/latex]

    Nun multiplizieren wir mit [latex]x(x^2+2x+2)[/latex] und erhalten

    [latex]
    \begin{aligned}[]1 = A(x^2+2x+2) + (Bx+C)x.\end{aligned}
    [/latex]

    Für [latex]x = 0[/latex] ergibt sich [latex]A = \frac 12[/latex]. Daher ist

    [latex]
    \begin{aligned}[]1 = (\tfrac 12 + B) x^2 + (1+C)x + 1\end{aligned}
    [/latex]

    und [latex]B=-\frac 12[/latex] und [latex]C = -1[/latex]. Es folgt

    [latex]
    \begin{aligned}[]\int \frac {1}{x(x^2+2x+2)} \thinspace {\rm {d}} x &= \tfrac 12 \int \frac {1}{x} \thinspace {\rm {d}} x - \tfrac 12 \int \frac {x+2}{x^2+2x+2}\thinspace {\rm {d}} x\\ &= \tfrac 12 \log |x|- \tfrac {1}{2} \int \frac {x+2}{(x+1)^2 +1} \thinspace {\rm {d}} x \\ &= \tfrac 12 \log |x|- \tfrac {1}{2} \int \frac {u+1}{u^2+1}\thinspace {\rm {d}} x\\ &= \tfrac 12 \log |x|- \tfrac {1}{4} \log |u^2+1| - \tfrac 12 \arctan (u) + D\\ &= \tfrac 12 \log |x|- \tfrac {1}{4} \log ((x+1)^2+1) - \tfrac 12 \arctan (x+1) + D,\end{aligned}
    [/latex]

    wobei wir [latex]u=x+1[/latex] gesetzt haben und Beispiele 8.18 (c) und (d) verwendet haben.

In manchen Fällen kann obiges Verfahren auch auf das Integral [latex]\int \frac {1}{(a^2+x^2)^n}\thinspace {\rm {d}} x[/latex] für ein [latex]a \in \mathbb {R}[/latex] und [latex]n \geq 2[/latex] führen, was wir mit der trigonometrischen Substitution [latex]\tan (u) = \frac {x}{a}[/latex] (siehe unten) behandeln können. Eine andere, allgemeinere Herangehensweise möchten wir in folgender Bemerkung für Interessierte behandeln.

Bemerkung: Integration rationaler Funktionen mit mehrfachen komplexen Nullstellen

Wie oben schon bemerkt, kann man nach der Partialbruchzerlegung ein Integral einer rationalen Funktion auf die Integration von Ausdrücken der Form [latex]\frac {1}{(x-a)^k}[/latex] oder von [latex]\frac {Ax+B}{(x^2+bx+c)^k}[/latex] für [latex]k \in \mathbb {N}[/latex] und für Konstanten [latex]a,A,B,b,c\in \mathbb {R}[/latex] zurückführen, wobei die Polynome der Form [latex]x^2+bx+c[/latex] keine rellen Nullstellen haben. Für die Berechnung eines Integrals des zweiten Typs mit [latex]k>1[/latex] möchten wir hier einen Algorithmus erläutern, wobei wir uns auf den Fall [latex]x^2+bx+c=x^2+1[/latex] beschränken (auf welchen man den allgemeinen Fall mit quadratischem Ergänzen zurückführen kann).

Seien also [latex]k>1[/latex] und ein Polynom [latex]q[/latex] von Grad kleiner als [latex]2k[/latex] gegeben. Dann ist das unbestimmte Integral [latex]\int \frac {q(x)}{(x^2+1)^k} \thinspace {\rm {d}} x[/latex] immer von der Form
[latex]
\begin{aligned}[]\label{eq:fund-hoehereOrdkomplNST} \frac {p(x)}{(x^2+1)^{k-1}} + \alpha \arctan (x) + \beta \log (x^2+1)+C.\end{aligned}
[/latex]
für ein Polynom [latex]p[/latex] von Grad kleiner [latex]2k-2[/latex] und Konstanten [latex]\alpha ,\beta[/latex]. Durch Ableiten, auf den gemeinsamen Nenner bringen und Vergleich der Koeffizienten lässt sich somit die Stammfunktion ermitteln.

Übung 8.20

Wir möchten in dieser Übung den oben erklärten Algorithmus genauer erklären und beginnen mit einem konkreten Beispiel.

  1. Berechnen Sie das Integral [latex]\int \frac {1}{(x^2+1)^2} \thinspace {\rm {d}} x[/latex].

Sei nun [latex]k>1[/latex] und [latex]q[/latex] ein Polynom von Grad kleiner als [latex]2k[/latex].

  1. Zeigen Sie, dass die Ableitung von [latex]\frac {p(x)}{(x^2+1)^{k-1}}[/latex] für ein beliebiges Polynom [latex]p[/latex] von Grad kleiner als [latex]2k-2[/latex] durch
    [latex]
    \begin{aligned}[]\frac {(x^2+1)p'(x) - 2(k-1)xp(x)}{(x^2+1)^{k}}\end{aligned}
    [/latex]

    gegeben ist.

  2. Berechnen Sie die Matrixdarstellung [latex]M[/latex] der Abbildung
    [latex]
    \begin{aligned}[]\phi : p(x) \mapsto (x^2+1)p'(x) - 2(k-1)xp(x)\end{aligned}
    [/latex]

    bezüglich der Basis der Monome.

  3. Schliessen Sie auf die Darstellung in (8.5), indem Sie zeigen, dass das Bild von [latex]\phi[/latex] zusammen mit [latex](x^2+1)^{k-1}[/latex] und [latex]2x(x^2+1)^{k-1}[/latex] den Vektorraum der Polynome von Grad kleiner gleich [latex]2k-1[/latex] aufspannt.

8.2.4 – Trigonometrische Substitution

In allen bisherigen Beispielen der Substitutionsregel in Abschnitt 8.2.2 hatten wir das Glück, dass das vorhandene Integral (vielleicht nach etwas Arbeit) bereits die richtige Struktur besass. Man verwendet die Substitutionsregel aber oft auch bevor wir noch wissen welches Integral sich eigentlich nach der Substitution ergibt, wobei es gewisse Funktionentypen gibt bei denen eine gewisse Substitution erfahrungsgemäss erfolgreich sein könnte. Wir wenden uns nun einem konkreten Beispiel dessen zu.

Beispiel 8.21: Kreisfläche

Wir möchten für [latex]r>0[/latex] das unbestimmte Integral [latex]\int \sqrt {r^2-x^2} \thinspace {\rm {d}} x[/latex] berechnen. Auf Grund der trigonometrischen Identitäten [latex]\sqrt {r^2 - r^2 \sin ^2(\theta )} = r \cos (\theta )[/latex] bietet es sich nun an, die Funktion

[latex]
\begin{aligned}[]f:I_\theta = \left (- \tfrac {\pi }{2},\tfrac {\pi }{2} \right )\to I_x=(-r,r),\ \theta \mapsto x=r \sin (\theta )\end{aligned}
[/latex]

für die Substitution zu verwenden. Denn mit dieser Substitution haben wir die Hoffnung, die Wurzel in einen anderen Ausdruck zu verwandeln.

Allerdings ist dies umgekehrt zu der Substitution in Abschnitt 8.2.2, da wir hier die «neue Variable»  [latex]\theta[/latex] verwenden um die «alte Variable»  [latex]x=r\sin (\theta )[/latex] auszudrücken. (Anstatt wie in Abschnitt 8.2.2 wo wir die neue Variable [latex]u=f(x)[/latex] als Funktion der alten Variable [latex]x[/latex] gesehen haben). Da [latex]f[/latex] bijektiv ist, ist dies kein Problem: denn [latex]x=r\sin (\theta )\in (-r,r)[/latex] ist zu [latex]\theta =\arcsin (\frac {x}r)\in (-\frac \pi 2,\frac \pi 2)[/latex] äquivalent. Weiters ist die Ableitung von [latex]f=\frac {\thinspace {\rm {d}} x}{\thinspace {\rm {d}}\theta }[/latex] gleich [latex]r\cos \theta[/latex] und damit auf ganz [latex](- \frac {\pi }{2},\frac {\pi }{2})[/latex] ungleich [latex]0[/latex]. Gemeinsam mit dem Satz über die Ableitung der inversen Funktion (Satz 7.14) erhalten wir daher

[latex]
\begin{aligned}[]\int \sqrt {r^2-x^2} \thinspace {\rm {d}} x & = \int r\cos (\theta ) \left (\tfrac {\thinspace {\rm {d}} x}{\thinspace {\rm {d}} \theta }\right )\left (\tfrac {\thinspace {\rm {d}} \theta }{\thinspace {\rm {d}} x}\right )\thinspace {\rm {d}} x\\ &=\int r^2 \cos ^2(\theta ) \thinspace {\rm {d}}\theta \\ &= r^2 \int \frac {1+\cos (2\theta )}{2} \thinspace {\rm {d}} \theta \\ &= \tfrac {r^2}{2} \left ( \theta + \tfrac {1}{2} \sin (2\theta ) \right ) + C \\ &= \tfrac {r^2}2 \arcsin \left (\tfrac {x}{r}\right ) + \tfrac {1}{2}x \sqrt {r^2-x^2} + C,\end{aligned}
[/latex]

wobei wir eben [latex]x = r \sin (\theta )[/latex], [latex]\sqrt {r^2-x^2} = r \cos (\theta )[/latex], [latex]\thinspace {\rm {d}} x = r \cos (\theta )\thinspace {\rm {d}} \theta[/latex] und die trigonometrischen Identitäten

[latex]
\begin{aligned}[]\cos (2 \theta ) &= \cos ^2(\theta ) - \sin ^2(\theta ) = 2\cos ^2(\theta ) -1\\ \cos ^2(\theta ) &= \frac {\cos (2\theta )+1}{2}\\ \sin (2\theta ) &= 2\sin (\theta ) \cos (\theta )\end{aligned}
[/latex]

für [latex]\theta \in (-\frac {\pi }2,\frac {\pi }2)[/latex] verwendet haben.

Veranschaulichen Sie sich die Substitution und die wichtigsten der obigen Identitäten in einem rechtwinkeligen Dreieck. Geben Sie weiters eine geometrische Interpretation der beiden Terme des unbestimmten Integrals bei der Berechnung des bestimmten Integrals [latex]\int _0^b\sqrt {r^2-x^2} \thinspace {\rm {d}} x[/latex] für [latex]0

Dies zeigt, dass

[latex]
\begin{aligned}[]G_0:x\in (-r,r)\mapsto \tfrac 12r^2 \arcsin \left (\tfrac {x}{r}\right ) + \tfrac {1}{2}x \sqrt {r^2-x^2}\end{aligned}
[/latex]

eine Stammfunktion von [latex]g_0:x\in (-r,r)\mapsto \sqrt {r^2-x^2}[/latex] ist (was wie immer viel einfacher zu überprüfen ist). Da aber sogar die Funktion [latex]g:x\in [-r,r]\mapsto \sqrt {r^2-x^2}[/latex] stetig ist, besitzt [latex]g[/latex] nach Korollar 8.3 auch auf ganz [latex][-r,r][/latex] eine Stammfunktion [latex]G_1[/latex], welche auf [latex](-r,r)[/latex] mit [latex]G_0+C[/latex] übereinstimmt. Da aber

[latex]
\begin{aligned}[]G:x\in [-r,r]\mapsto \tfrac 12r^2 \arcsin \left (\tfrac {x}{r}\right ) + \tfrac {1}{2}x \sqrt {r^2-x^2}\end{aligned}
[/latex]

auch eine auf ganz [latex][-r,r][/latex] stetige Funktion definiert, folgt aus Stetigkeit von [latex]G[/latex] und [latex]G_1[/latex], dass [latex]G_1=G+C[/latex] und damit ist [latex]G[/latex] auf ganz [latex][-r,r][/latex] eine Stammfunktion von [latex]g[/latex]. Wir bemerken allerdings, dass [latex]\arcsin \left (\frac {x}r\right )[/latex] keine Ableitung in den Punkten [latex]-r[/latex] und [latex]r[/latex] besitzt. (Wieso ist dies kein Widerspruch zu obiger Diskussion?)

Substitutionen wie obige nennen sich vielfach trigonometrische Substitutionen. Wir werden bei diesen Berechnungen nicht immer so sorgfältig argumentieren und vielmehr der Leibniz Notation vertrauen, doch muss immer Invertierbarkeit der Funktion gegeben sein wenn wir die alte Variable durch die neue Variable ausdrücken. Für die folgende Auflistung der trigonometrischen Substitutionen sei [latex]n\in \mathbb {Z}[/latex].

  • In Ausdrücken der Form [latex](a^2-x^2)^{\frac {n}{2}}[/latex] für [latex]a>0[/latex] führt wie bereits im obigen Beispiel oft die Substitution [latex]x=a\sin (\theta )[/latex] mit [latex]\theta \in (-\frac \pi 2,\frac \pi 2)[/latex] zum Ziel, wobei sich damit [latex]\thinspace {\rm {d}} x = a \cos (\theta ) \thinspace {\rm {d}} \theta[/latex] und [latex](a^2-x^2)^{\frac {1}{2}} = a \cos (\theta )[/latex] ergibt.
  • In Ausdrücken der Form [latex](a^2+x^2)^{\frac {n}{2}}[/latex] für [latex]a > 0[/latex] führt oft die Substitution [latex]x = a\tan (\theta )[/latex] mit [latex]\theta \in (-\frac \pi 2,\frac \pi 2)[/latex] zum Ziel, wobei sich damit [latex]\thinspace {\rm {d}} x = \frac {a}{\cos ^2(\theta )} \thinspace {\rm {d}} \theta[/latex] und [latex](a^2+x^2)^{\frac 12} = \frac {a}{\cos (\theta )}[/latex] ergibt.
  • Obwohl dies keine trigonometrische Substitution darstellt, bemerken wir noch Folgendes. Falls ein «einzelnes» [latex]x[/latex] vor dem Ausdruck [latex](a^2-x^2)^{\frac {n}{2}}[/latex] oder dem Ausdruck [latex](a^2+x^2)^{\frac {n}{2}}[/latex] steht, ist die Substitution [latex]u=a^2-x^2[/latex] respektive [latex]u=a^2+x^2[/latex] teilweise viel einfacher.

Als Merkhilfe kann es helfen für die beiden trigonometrischen Substitution ein rechtwinkeliges Dreieck zu skizzieren und abhängig von der Substitution die Seiten mit Hilfe von Pythagoras entsprechend zu beschriften.

Beispiel 8.22: Trigonometrische Substitution

  1. Es gilt für [latex]a >0[/latex]
    [latex]
    \begin{aligned}[]\int \frac {1}{(a^2+x^2)^{\frac 32}} \thinspace {\rm {d}} x &= \int \frac {\cos ^3(\theta )}{a^3} a\frac {1}{\cos ^2(\theta )} \thinspace {\rm {d}} \theta = \tfrac {1}{a^2} \int \cos (\theta ) d\theta = \tfrac {1}{a^2} \sin (\theta ) + C \\ &= \frac {x}{a^2\sqrt {a^2+x^2}} + C,\end{aligned}
    [/latex]

    wobei wir [latex]x = a \tan (\theta )[/latex], [latex]\sqrt {a^2+x^2} = a \frac {1}{\cos (\theta )}[/latex], [latex]\thinspace {\rm {d}} x = a \frac {1}{\cos ^2(\theta )} \thinspace {\rm {d}} \theta[/latex] verwendet haben. (Veranschaulichen Sie sich die Substitution und obige Identitäten in einem Bild.)

  2. Es ist
    [latex]
    \begin{aligned}[]\int x \sqrt {1-x^2} \thinspace {\rm {d}} x = - \tfrac {1}{2} \int u^{\frac 12} \thinspace {\rm {d}} u = - \tfrac {1}{2} \tfrac {2}{3} u^{\frac {3}{2}} + C = - \tfrac {1}{3} (1-x^2)^{\frac {3}{2}} + C,\end{aligned}
    [/latex]

    wobei [latex]u = 1-x^2[/latex], [latex]\thinspace {\rm {d}} u = -2x \thinspace {\rm {d}} x[/latex].

8.2.5 – Weitere Integrationsmethoden

Es gibt viele weitere Methoden zur Integration; viele davon beruhen auf spezielle Substitutionen.

Beispielsweise lassen sich gewisse unbestimmte Integrale mit hyperbolischen Substitutionen berechnen. Sei [latex]n\geq -1[/latex]. In Ausdrücken der Form [latex](x^2-a^2)^{\frac {n}{2}}[/latex] für [latex]a\in \mathbb {R}[/latex] führt oft die Substitution [latex]x=a\cosh (u)[/latex] zum Ziel, wobei sich damit [latex]\thinspace {\rm {d}} x = a \sinh (u) \thinspace {\rm {d}} u[/latex] und [latex](x^2-a^2)^{\frac {1}{2}} = a \sinh (u)[/latex] ergibt.

Beispiel 8.23

Wir berechnen

[latex]
\begin{aligned}[]\int \sqrt {x^2-1} \thinspace {\rm {d}} x &= \int \sinh ^2(u) \thinspace {\rm {d}} u = \cosh (u)\sinh (u) - \int \cosh ^2(u) \thinspace {\rm {d}} u\\ &= \cosh (u)\sinh (u) - \int \sinh ^2(u) + 1 \thinspace {\rm {d}} u+C\\ &= \cosh (u)\sinh (u) - u - \int \sinh ^2(u) \thinspace {\rm {d}} u+C,\end{aligned}
[/latex]

wobei [latex]x = \cosh (u)[/latex] und [latex]\thinspace {\rm {d}} x = \sinh (u)\thinspace {\rm {d}} u[/latex]. Nach Auflösen ergibt sich somit

[latex]
\begin{aligned}[]\int \sqrt {x^2-1} \thinspace {\rm {d}} x = \int \sinh ^2(u) \thinspace {\rm {d}} u = \frac {\cosh (u)\sinh (u)-u}{2} + C = \frac {x \sqrt {x^2-1}-\operatorname {arcosh}(x)}{2} + C.\end{aligned}
[/latex]

Eine andere Methode, die wir hier kurz erwähnen möchten, ist die sogenannte Halbwinkelmethode (oder auch Weierstrass-Substitution). Diese ist dann nützlich, wenn man das Integral einer rationalen Funktion in [latex]\cos (x)[/latex] und [latex]\sin (x)[/latex] wie zum Beispiel [latex]\frac {\cos ^2(x)}{\sin (x) + 2017}[/latex] in die Integration einer rationalen Funktion in [latex]u=\tan \left (\frac {x}{2}\right )[/latex] umwandeln möchte (siehe auch Beispiel 8.17 (b)).

Übung 8.24: Halbwinkelmethode

Wir möchten das unbestimmte Integral [latex]\int \frac {\cos (x)}{2+\sin (x)}\thinspace {\rm {d}} x[/latex] mit der Substitution [latex]u = \tan \left (\frac {x}{2}\right )[/latex] berechnen. Zeigen Sie dafür zuerst die Identitäten

[latex]
\begin{aligned}[]\sin (x)=\frac {2u}{1+u^2},\quad \cos (x)=\frac {1-u^2}{1+u^2}.\end{aligned}
[/latex]

Zeigen Sie anschliessend, dass das obige Integral nach Substitution zu einem Integral einer rationalen Funktion in [latex]u[/latex] wird und berechnen Sie es.

Manchmal führt man auch die eine oder die andere Substitution durch, weil in der zu integrierenden Funktion eine verschachtelte Funktion vorliegt und man einfach keine andere Methode zur Verfügung hat. Zum Beispiel bei dem Integral [latex]\int \sin (\sqrt {x})\thinspace {\rm {d}} x[/latex] steht keine der erwähnten Methoden zur Verfügung, doch ist man versucht [latex]u=\sqrt {x}[/latex] zu setzen um zu sehen was sich daraus ergibt. Dies führt in der Tat zum Erfolg (wieso?). Ebenso in dem Integral der Form [latex]\int \frac 1{1+\exp (x)}\thinspace {\rm {d}} x[/latex] führt der Ansatz [latex]u=\exp (x)[/latex] zu einem unbestimmten Integral einer rationalen Funktion (wieso?).

8.2.6 – Neue Funktionen

Manchmal führen obige Methoden zur Bestimmung eines unbestimmten Integrals einer Funktion zu keinem Ergebnis. Dies kann daran liegen, dass die gesuchte Stammfunktion sich nicht mit den bisher bekannten Funktionen ausdrücken lässt.

Beispiel 8.25: Integralsinus

Der Integralsinus ist die Stammfunktion [latex]\operatorname {Si}: \mathbb {R} \to \mathbb {R}[/latex] der stetigen Funktion

[latex]
\begin{aligned}[]x \in \mathbb {R} \mapsto \left \lbrace \begin{array}{cc} \frac {\sin (x)}{x} & \text {falls } x \neq 0 \\ 1 & \text {falls } x= 0\end{array} \right .\end{aligned}
[/latex]

mit der Normalisierung [latex]\operatorname {Si}(0) = 0[/latex]. Er lässt sich als Potenzreihe schreiben, denn nach Satz 6.85 gilt

[latex]
\begin{aligned}[]\operatorname {Si}(x) = \int _0^x \frac {\sin (t)}{t} \thinspace {\rm {d}} t = \int _0^x \sum _{n=0}^\infty \frac {(-1)^n}{(2n+1)!} t^{2n} \thinspace {\rm {d}} t = \sum _{n=0}^\infty \frac {(-1)^n}{(2n+1)! (2n+1)} x^{2n+1}\end{aligned}
[/latex]

für alle [latex]x \in \mathbb {R}[/latex].

Beispiel 8.26: Integralkosinus

Der Integralkosinus [latex]\operatorname {Ci}: (0,\infty ) \to \mathbb {R}[/latex] ist definiert als die Stammfunktion von [latex]x \in (0,\infty ) \mapsto \frac {\cos (x)}{x} \in \mathbb {R}[/latex] mit der Normalisierung [latex]\lim _{x \to \infty } \operatorname {Ci}(x) = 0[/latex].

Dabei möchten wir auf folgende Übung verweisen, die zeigt, dass der Integralkosinus so wohldefiniert ist.

Übung 8.27

Sei [latex]F[/latex] eine Stammfunktion von [latex]x \in (0,\infty ) \mapsto \frac {\cos (x)}{x} \in \mathbb {R}[/latex]. Zeigen Sie, dass der Grenzwert [latex]\lim _{x \to \infty } F(x)[/latex] existiert. Drücken Sie [latex]\operatorname {Ci}[/latex] als Summe einer Konstanten (der sogenannten Euler-Mascheroni Konstanten), der Logarithmusfunktion und einer Potenzreihe aus.

Unter Verwendung uneigentlicher Integrale werden wir später weitere wichtige Funktionen kennenlernen, die sich nicht in Termen bekannter Funktionen ausdrücken lassen — siehe zum Beispiel 8.45.

8.2.7 – Das bestimmte Integral

Alle obigen Regeln zur Berechnung des unbestimmten Integrals lassen sich nach dem Fundamentalsatz der Integral- und Differentialrechnung eins zu eins auch für das Riemann-Integral, welches im Gegensatz zum unbestimmten Integral auch das bestimmte Integral genannt wird, anwenden. Dabei haben wir zwei Möglichkeiten.

  • Eine erste Möglichkeit ist mit obigen Methoden zuerst das unbestimmte Integral zu berechnen und dann Korollar 8.4 zur Berechnung des Riemann-Integrals zu verwenden.
  • Falls wir aber nur an einem einzigen Riemann-Integral interessiert sind, ist es oft einfacher, die Ausdrücke ausserhalb des Integrals so früh wie möglich zu berechnen. Wir erklären dies im Folgenden für die partielle Integration und die Substitution.

Sind [latex]u,v[/latex] zwei stetig differenzierbare Funktionen auf einem kompakten Intervall [latex][a,b][/latex] mit Endpunkten [latex]a

[latex]
\begin{aligned}[]\int _a^b u v' \thinspace {\rm {d}} x = \left [ uv \right ]_a^b - \int _a^b u'v \thinspace {\rm {d}} x.\end{aligned}
[/latex]

Denn falls [latex]F[/latex] eine Stammfunktion von [latex]uv'[/latex] und [latex]G[/latex] eine Stammfunktion von [latex]u'v[/latex] ist, dann gilt für alle [latex]x \in [a,b][/latex]

[latex]
\begin{aligned}[]F(x) + C_1 = \int uv' \thinspace {\rm {d}} x = u(x)v(x) - \int u'v \thinspace {\rm {d}} x + C_2 = u(x)v(x) - G(x)+ C_3\end{aligned}
[/latex]

für gewisse Integrationskonstanten [latex]C_1,C_2,C_3[/latex]. Somit ist nach Korollar 8.4

[latex]
\begin{aligned}[]\int _a^b uv' \thinspace {\rm {d}} x = \left [F(x)\right ]_a^b &= F(b) - F(a) = u(b)v(b)-G(b)- (u(a)v(a)-G(a))\\ &= \left [u(x)v(x)\right ]_a^b - \left [G(x)\right ]_a^b= \left [u(x)v(x)\right ]_a^b - \int _a^b u'v \thinspace {\rm {d}} x.\end{aligned}
[/latex]

Ebenso können wir bei einer Substitution in Abschnitt 8.2.2 die Grenzen für ein Riemann-Integral enstsprechend der Substitution neu berechnen. Sei [latex]I_x[/latex] ein Intervall mit Endpunkten [latex]a_x

[latex]
\begin{aligned}[]\int _a^b g \circ f(x) f'(x)\thinspace {\rm {d}} x = \int _{f(a)}^{f(b)} g(u) \thinspace {\rm {d}} u.\end{aligned}
[/latex]

In der Tat, wenn [latex]G[/latex] eine Stammfunktion von [latex]g[/latex] auf [latex]I_u[/latex] ist, dann ist nach der Kettenregel [latex]G \circ f[/latex] eine Stammfunktion von [latex]x\in I_x \mapsto g \circ f(x) f'(x)[/latex]. Nach Korollar 8.4 gilt also

[latex]
\begin{aligned}[]\int _a^b g \circ f(x) f'(x) \thinspace {\rm {d}} x = \left [G \circ f\right ]_a^b = G(f(b))-G(f(a)) = \left [G\right ]_{f(a)}^{f(b)} = \int _{f(a)}^{f(b)} g(u) \thinspace {\rm {d}} u.\end{aligned}
[/latex]

Die Annahme der Stetigkeit an [latex]g[/latex] kann abgeschwächt werden — siehe die entsprechende Übung im Abschnitt 8.8.2.

8.2.8 – Leibniz-Notation

Wir werden die Leibniz-Notation in der Berechnung von unbestimmten und bestimmten Integralen wie bereits oben im Folgenden immer wieder verwenden. Diese Notation verpackt in einem natürlichen Formalismus die partielle Integration

[latex]
\begin{aligned}[]\int u \thinspace {\rm {d}} v = uv - \int v \thinspace {\rm {d}} u + C\end{aligned}
[/latex]

und die Substitutionsregeln

[latex]
\begin{aligned}[]\int g(u(x)) \frac {\thinspace {\rm {d}} u}{\thinspace {\rm {d}} x} \thinspace {\rm {d}} x &= \int g(u) \thinspace {\rm {d}} u\\ \int g(u(x))\thinspace {\rm {d}} x &= \int g(u) \frac {\thinspace {\rm {d}} x}{\thinspace {\rm {d}} u}(u) \thinspace {\rm {d}} u,\end{aligned}
[/latex]

wobei wir in der zweiten Formulierung der Substitution vorraussetzen, dass [latex]u:I_x\to I_u[/latex] bijektiv mit nicht verschwindender Ableitung ist und dadurch im linken Integral mit [latex]1=\frac {\thinspace {\rm {d}} x}{\thinspace {\rm {d}} u} \frac {\thinspace {\rm {d}} u}{\thinspace {\rm {d}} x}[/latex] multiplizieren konnten und die erste Formulierung der Substitutionsregel anwenden konnten. Wie wir gesehen haben, sind diese Regeln Umformulierungen der Produktregel für die Ableitung und der Kettenregel für die Ableitung (gemeinsam mit der Ableitungsregel für die inverse Abbildung).

Bei konkreten Integralberechnungen verwenden wir mitunter auch Gleichungen, die [latex]\thinspace {\rm {d}} x[/latex] und [latex]\thinspace {\rm {d}} u[/latex] miteinander verbinden. Zum Beispiel bei der trigonometrischen Substitution [latex]x=a\sin \theta[/latex] (für [latex]x\in (-a,a)[/latex] und [latex]\theta \in (-\frac \pi 2,\frac \pi 2)[/latex]) verwenden wir auch die Formel [latex]\thinspace {\rm {d}} x=a\cos \theta \thinspace {\rm {d}}\theta[/latex], die formal gesehen keine Bedeutung hat (und deswegen auf keinen Fall in dieser Form in Beweisen auftreten sollte), doch eben im Zuge der Substitution in der Formulierung der Leibniz-Notation einen bequemen Zwischenschritt darstellt.

Informell taucht in Anwendungen das Symbol [latex]\thinspace {\rm {d}} x[/latex] auch oft in Diskussionen auf, die zu einem Riemann-Integral führen, wobei [latex]\thinspace {\rm {d}} x[/latex] dann für ein (sehr) kleines [latex]\Delta x[/latex] stehen sollte. In Anwendungen werden häufig die Begriffe der Riemann-Summe oder der additiven Intervallfunktion vermieden, wobei es genau diese Begriffe sind, die diese Verwendung von [latex]\thinspace {\rm {d}} x[/latex] genau und formal korrekt machen würden (siehe Abschnitte 4.4 und 5.5). Auf jeden Fall hat in diesem Zusammenhang eine Formel der Gestalt [latex]\thinspace {\rm {d}} x=a\cos \theta \thinspace {\rm {d}}\theta[/latex] auch eine Interpretation: Da [latex]\Delta x[/latex] die Länge eines kleinen Teilintervalls von [latex](-a,a)[/latex] angibt und [latex]\Delta \theta[/latex] die Länge des entsprechenden Teilintervalls in [latex](-\frac \pi {2},\frac \pi 2)[/latex] so gibt die Ableitung [latex]a\cos \theta[/latex] (bis auf einen kleinen und wie sich herausstellt vernachlässigbaren Fehler) den Grössenunterschied [latex]\frac {\Delta x}{\Delta \theta }[/latex] an, der bei Betrachtung von etwaigen Riemann-Summen in der Variable [latex]x[/latex] und der Variable [latex]\theta[/latex] als zusätzlicher Faktor auftreten würde. Wir müssen dies nicht genauer ausführen oder die Substitution auf diese Art und Weise beweisen, da wir ja mittels der Kettenregel und dem Fundamentalsatz der Integral- und Differentialrechnung bereits die Substitutionsregel für Riemann-Integrale bewiesen haben und obiger Formalismus diese nur auf eine andere Art präsentiert. Dieser Beweis über den Fundamentalsatz verwendet allerdings etwas stärkere Annahmen als notwendig (siehe folgende Übung für den direkten Beweis mit schwächeren Annahmen).

Übung 8.28: Substitution für Riemann-integrierbare Funktionen

Sei [latex][a,b][/latex] ein kompaktes Intervall in [latex]\mathbb {R}[/latex] mit Endpunkten [latex]a

[latex]
\begin{aligned}[]\int _a^b g \circ f(t) f'(t) \thinspace {\rm {d}} t = \int _{f(a)}^{f(b)} g(x) \thinspace {\rm {d}} x.\end{aligned}
[/latex]

8.3 – Anwendungen

8.3.1 – Flächeninhalte

Wir wollen hier nochmals Beispiele für Flächenberechnungen besprechen, welche unter anderem den Namen der Umkehrfunktionen der hyperbolischen Funktionen erklären.

Beispiel 8.29

Wir berechnen den Flächeninhalt des Kreises mit Radius [latex]r> 0[/latex]. Dieser ist durch [latex]2 \int _{-r}^r \sqrt {r^2-x^2} \thinspace {\rm {d}} x[/latex] definiert (wieso?), und gemeinsam mit Bespiel 8.21 ergibt sich daraus

[latex]
\begin{aligned}[]2 \int _{-r}^r \sqrt {r^2-x^2} \thinspace {\rm {d}} x = 2\left [\tfrac 12r^2 \arcsin \left (\tfrac {x}{r}\right ) + \tfrac {1}{2}x \sqrt {r^2-x^2}\right ]_{-r}^r = r^2\pi\end{aligned}
[/latex]

Beispiel 8.30: Hyperbolische Umkehrfunktionen

Wir verwenden die Funktion
[latex]
\begin{aligned}[]\label{eq:fund-param hyperbel} t\in \mathbb {R} \mapsto (x,y) = (\cosh (t), \sinh (t)) \in \mathbb {R}^2\end{aligned}
[/latex]
um die «positive Hälfte» der Hyperbel [latex]\left \lbrace {(x,y)} \mid {x^2-y^2 = 1}\right \rbrace[/latex] zu parametrisieren. Wir stellen uns den Parameter [latex]t \in \mathbb {R}[/latex] vorerst als Zeit vor. In diesem Sinne beschreibt (8.6) eine Bewegung im [latex]\mathbb {R}^2[/latex]. Wir wollen den Flächeninhalt des folgenden Gebietes in Rosa zwischen dem Ursprung und einem Teil der Hyperbel berechnen.

image

Dieser ist der Flächeninhalt [latex]\frac {1}{2}x_0y_0 = \frac {1}{2}\cosh (t_0)\sinh (t_0)[/latex] des eingezeichneten Dreiecks minus dem Flächeninhalt unterhalb der Hyperbel zwischen der [latex]1[/latex] und [latex]x_0 = \cosh (t_0)[/latex] in Blau. Letztere Fläche ist durch [latex]\int _1^{x_0} \sqrt {x^2-1}\thinspace {\rm {d}} x[/latex] gegeben. Um dieses Integral zu berechnen, verwenden wir die hyperbolische Substitution [latex]x = \cosh (t)[/latex], [latex]\thinspace {\rm {d}} x = \sinh (t) \thinspace {\rm {d}} t[/latex] und erhalten

[latex]
\begin{aligned}[]\int _1^{x_0} \sqrt {x^2-1}\thinspace {\rm {d}} x &= \int _0^{t_0} \sinh ^2(t) \thinspace {\rm {d}} t = \int _0^{t_0} \frac {\mathrm {e}^{2t}-2 + \mathrm {e}^{-2t}}{4} \thinspace {\rm {d}} t = \left [ \frac {\mathrm {e}^{2t}-\mathrm {e}^{-2t}}{8} - \frac {1}{2}t \right ]_0^{t_0}\\ &= \left [ \frac {1}{4} \sinh (2t) - \frac {1}{2}t \right ]_0^{t_0} = \frac {1}{2}x_0y_0 - \frac {1}{2}t_0.\end{aligned}
[/latex]

Somit ist der Flächeninhalt des gesuchten Gebiets

[latex]
\begin{aligned}[]\frac {1}{2}t_0 = \frac 12\operatorname {arcosh}(x_0)=\frac {1}{2} \operatorname {arsinh}(y_0).\end{aligned}
[/latex]

Dies erklärt die Namen «Areasinus Hyperbolicus» und «Areakosinus Hyperbolicus» der Umkehrfunktionen der hyperbolischen Funktionen (wieso?).

8.3.2 – Bogenlänge

Im Folgenden möchten wir einen stetige Funktion [latex]\gamma : [a,b] \to \mathbb {R}^d[/latex] für [latex]d \geq 2[/latex], auch Weg oder Kurve von [latex]\gamma (a)[/latex] nach [latex]\gamma (b)[/latex] genannt, betrachten. Dabei fassen wir [latex]t \in [a,b][/latex] als Zeitparameter und [latex]\gamma (t)[/latex] als die Position zum Zeitpunkt [latex]t[/latex] auf.

Falls alle Komponenten [latex]\gamma _1,\ldots , \gamma _d[/latex] von [latex]\gamma = (\gamma _1,\ldots ,\gamma _d)^t[/latex] stetig differenzierbar sind, interpretieren wir für einen Zeitpunkt [latex]t \in [a,b][/latex] den Ausdruck [latex]\| \dot {\gamma }(t)\| _2=\sqrt {\dot {\gamma }_1(t)^2 +\ldots + \dot {\gamma }_d(t)^2}[/latex] als die Geschwindigkeit zum Zeitpunkt [latex]t \in [a,b][/latex].

Wir möchten nun die Bogenlänge des Weges [latex]\gamma[/latex] als die gesamte Strecke, die zwischen den Zeiten [latex]a[/latex] und [latex]b[/latex] zurückgelegt wurde, definieren. Dabei soll gelten, dass die zurückgelegte Strecke zwischen gleichen Zeiten [latex]\alpha[/latex] und [latex]\alpha[/latex] Null ist und dass sich Strecken additiv verhalten, also dass die zwischen den Zeiten [latex]\alpha 4.29 ist die zurückgelegte Strecke also eine additive Intervallfunktion auf [latex][a,b][/latex].

Des Weiteren möchten wir natürlich verlangen, dass die in einem Teilintervall [latex][\alpha ,\beta ] \subseteq [a,b][/latex] mit [latex]\alpha 4.30 ist daher die einzig vernünftige Definition der Bogenlänge des Weges [latex]\gamma[/latex] der Ausdruck

[latex]
\begin{aligned}[]L(\gamma ) = \int _a^b \| \dot {\gamma }(t)\| _2\thinspace {\rm {d}} t=\int _a^b\sqrt {\dot {\gamma }_1(t)^2 +\ldots + \dot {\gamma }_d(t)^2} \thinspace {\rm {d}} t.\end{aligned}
[/latex]

Anders formuliert ist also die Länge des zurückgelegten Weges das Integral über die Geschwindigkeitsfunktion.

Beispiel 8.31: Umfang des Kreises

Wir betrachten den Weg

[latex]
\begin{aligned}[]\gamma : t \in [0,2\pi ] \mapsto (\cos (t),\sin (t))^t \in \mathbb {R}^2.\end{aligned}
[/latex]

Wegen [latex]\gamma (0) = \gamma (2\pi ) = (1,0)^t[/latex] sind der Start- und der Endpunkt von [latex]\gamma[/latex] gleich (wir sagen auch, dass der Weg [latex]\gamma[/latex] geschlossen ist). Auch gilt für die Geschwindigkeit zu jedem Zeitpunkt [latex]t\in [0,2\pi ][/latex]

[latex]
\begin{aligned}[]\| {\dot {\gamma }(t)}\| _2 = \sqrt {\dot {\gamma }_1(t)^2 +\dot {\gamma }_2(t)^2} = \sqrt {\sin ^2(t)+ \cos ^2(t)} = 1.\end{aligned}
[/latex]

Der Weg (oder die Kurve) [latex]\gamma[/latex] durchläuft (wegen [latex]\cos ^2(t) + \sin ^2(t) = 1[/latex] für alle [latex]t\in \mathbb {R}[/latex]) den Einheitskreis also mit konstanter Geschwindigkeit Eins. Deswegen gilt

[latex]
\begin{aligned}[]L(\gamma ) = \int _0^{2\pi } \sqrt {\dot {\gamma }_1(t)^2 +\dot {\gamma }_2(t)^2} \thinspace {\rm {d}} t = \int _0^{2\pi } 1 \thinspace {\rm {d}} t = 2\pi .\end{aligned}
[/latex]

Des Weiteren besucht [latex]\gamma[/latex] jeden Punkt (bis auf den Endpunkt) genau einmal (siehe auch Abschnitt 6.6.4). Einen solchen Weg nennen wir auch einfach. Deswegen lässt sich die Bogenlänge von [latex]\gamma[/latex] auch als den Umfang des Einheitskreises auffassen, der somit [latex]2\pi[/latex] ist. Dies gilt analog für Teilstrecken und definiert den Begriff Winkel als Bogenlänge am Einheitskreis.

Sei [latex]\gamma : [a,b] \to \mathbb {R}^d[/latex] ein stetig differenzierbarer Weg ausgehend von einem Intervall [latex][a,b][/latex] mit Endpunkten [latex]aReparametrisierung von [latex]\gamma[/latex] ist ein Weg der Form [latex]\gamma \circ \psi :[\tilde {a},\tilde {b}] \to \mathbb {R}^d[/latex], wobei [latex][\tilde {a},\tilde {b}][/latex] ein kompaktes Intervall mit Endpunkten [latex]\tilde {a}

Intuitiv ausgedrückt ist eine Reparametrisierung eines Weges also ein Weg mit denselben Endpunkten (da [latex]\psi (\tilde {a}) = a[/latex] und [latex]\psi (\tilde {b}) = b[/latex]) und der immer in dieselbe Richtung geht (wegen Monotonie). Anschaulich kann man deswegen erwarten, dass jede Reparametrisierung eines Weges dieselbe Bogenlänge hat. Auch wollen wir zeigen, dass ein nie anhaltender Weg so reparametrisiert werden kann, dass der neue Weg Einheitsgeschwindigkeit hat. Falls [latex]\gamma : [a,b] \to \mathbb {R}^d[/latex] nie anhält oder genauer falls [latex]\| {\dot {\gamma }(t)}\| _2 > 0[/latex] für alle [latex]t \in [a,b][/latex], so nennen wir [latex]\gamma[/latex] regulär.

Lemma 8.32: Reparametrisierungen eines Weges

Sei [latex]\gamma : [a,b] \to \mathbb {R}^d[/latex] ein stetig differenzierbarer Weg für [latex]a Parametrisierung nach Bogenlänge genannt wird.

In Beispiel 8.31 ist der betrachtete Weg bereits nach Bogenlänge parametrisiert.

Beweis

Sei [latex][\tilde {a},\tilde {b}][/latex] ein kompaktes Intervall mit Endpunkten [latex]\tilde {a}

[latex]
\begin{aligned}[]L(\gamma \circ \psi ) &= \int _{\tilde {a}}^{\tilde {b}} \sqrt {(\gamma _1 \circ \psi )'(s)^2 +\ldots + (\gamma _d \circ \psi )'(s)^2} \thinspace {\rm {d}} s \\ &= \int _{\tilde {a}}^{\tilde {b}} \sqrt { \dot {\gamma }_1(\psi (s))^2 \psi '(s)^2 +\ldots + \dot {\gamma }_d(\psi (s))^2 \psi '(s)^2} \thinspace {\rm {d}} s \\ &= \int _{\tilde {a}}^{\tilde {b}} \sqrt { \dot {\gamma }_1(\psi (s))^2 +\ldots + \dot {\gamma }_d(\psi (s))^2 } \psi '(s)\thinspace {\rm {d}} s \\ &= \int _a^b \sqrt {\dot {\gamma }_1(t)^2 +\ldots + \dot {\gamma }_d(t)^2} \thinspace {\rm {d}} t = L(\gamma ).\end{aligned}
[/latex]

Für die zweite Aussage konstruieren wir nun eine geeignete Funktion [latex]\psi[/latex] wie oben. Sei

[latex]
\begin{aligned}[]\phi :[a,b] \to [0,L(\gamma )],\ t \mapsto \int _{a}^t \sqrt {\dot {\gamma }_1(s)^2 +\ldots + \dot {\gamma }_d(s)^2} \thinspace {\rm {d}} s.\end{aligned}
[/latex]

Wegen [latex]\dot \phi (t) = \sqrt {\dot {\gamma }_1(t)^2 +\ldots + \dot {\gamma }_d(t)^2} > 0[/latex] für alle [latex]t \in [a,b][/latex] sowie [latex]\phi (a) = 0[/latex] und [latex]\phi (b) = L(\gamma )[/latex] ist [latex]\phi :[a,b] \to [0,L(\gamma )][/latex] eine streng monoton wachsende, stetig differenzierbare Bijektion. Insbesondere ist [latex]\psi = \phi ^{-1}:[0,L(\gamma )] \to [a,b][/latex] ebenfalls streng monoton wachsend und stetig differenzierbar. Zur Zeit [latex]s \in [0,L(\gamma )][/latex] berechnen wir nun die Geschwindigkeit von [latex]\gamma \circ \psi[/latex]. Ist [latex]t = \psi (s)[/latex], so gilt wegen [latex]\psi '(s)>0[/latex] auch

[latex]
\begin{aligned}[]\sqrt {(\gamma _1 \circ \psi )'(s)^2 +\ldots + (\gamma _d \circ \psi )'(s)^2} &= \sqrt { \dot {\gamma }_1(t)^2 \psi '(s)^2 +\ldots + \dot {\gamma }_d(t)^2 \psi '(s)^2}\\ &= \psi '(s) \sqrt { \dot {\gamma }_1(t)^2 +\ldots + \dot {\gamma }_d(t)^2 }\\ &= \frac {1}{\dot \phi (t)}\sqrt { \dot {\gamma }_1(t)^2 +\ldots + \dot {\gamma }_d(t)^2 } = 1.\end{aligned}
[/latex]

Somit hat die Reparametrisierung [latex]\gamma \circ \psi[/latex] die gewünschte Eigenschaft. ∎

Übung 8.33: Eindeutigkeit der Parametrisierung

In Lemma 8.32 wird bereits von der Parametrisierung nach Bogenlänge gesprochen. Wir wollen dies hier begründen. Sei [latex]\gamma :[a,b] \to \mathbb {R}^d[/latex] ein stetig differenzierbarer, regulärer Weg. Nach Lemma 8.32 dürfen wir annehmen, dass [latex]\gamma[/latex] Einheitsgeschwindigkeit hat. Zeigen Sie, dass es keine weitere Reparametrisierung von [latex]\gamma[/latex] mit Einheitsgeschwindigkeit gibt.

Übung 8.34: Totale Variation des Weges

In dieser Übung wollen wir noch eine weitere Begründung für die Definition der Bogenlänge eines Weges [latex]\gamma :[a,b]\to \mathbb {R}^d[/latex] geben. Hierfür interpretieren wir [latex]d(v,w)=\| v-w\| _2[/latex] als den Abstand zweier Punkte [latex]v,w\in \mathbb {R}^d[/latex]. Die totale Variation von [latex]\gamma :[a,b]\to \mathbb {R}^d[/latex] ist definiert als

[latex]
\begin{aligned}[]V(\gamma )=\sup _{\mathfrak {Z}}\sum _{k=1}^n\| \gamma (x_k)-\gamma (x_{k-1})\| ,\end{aligned}
[/latex]

wobei das Supremum über alle Zerlegungen [latex]\mathfrak {Z}=\left \lbrace {x_0=a

Hinweis.

Verwenden Sie den Mittelwertsatz für jede Komponente von [latex]\gamma _j[/latex] in jedem Intervall [latex][x_{k-1},x_k][/latex] für [latex]j=1,\ldots ,d[/latex] und [latex]k=1,\ldots ,n[/latex] gemeinsam mit gleichmässiger Stetigkeit der Funktion [latex]t\in [a,b]\mapsto (\gamma _1'(t),\ldots ,\gamma _d'(t))^t[/latex].

Für einen Weg [latex]\gamma :[a,b]\to \mathbb {R}^d[/latex] und eine stetige Funktion [latex]f:\mathbb {R}^d\to \mathbb {R}[/latex] kann ein Integral der Form

[latex]
\begin{aligned}[]\int _a^b f(\gamma (t)) \| \dot {\gamma }(t)\| _2\thinspace {\rm {d}} t\end{aligned}
[/latex]

auch physikalische Bedeutung haben. Zum Beispiel kann der Weg einen verbogenen Draht (mit konstanter Dichte [latex]1kg/m[/latex]) beschreiben. In diesem Fall gibt

[latex]
\begin{aligned}[]\frac 1{L(\gamma )} \int _a^b \gamma _j(t) \| \dot {\gamma }(t)\| _2\thinspace {\rm {d}} t\end{aligned}
[/latex]

die [latex]j[/latex]-te Koordinate des Schwerpunktes des Drahtes an, wobei [latex]j\in \{ 1,\ldots ,d\}[/latex].

8.3.3 – Wegintegrale von Vektorfeldern

Wir kommen nun zu einem weiteren Typ von Wegintegralen, der sowohl für die Physik als auch für die weitere Analysis wichtig sein wird. Hierfür betrachten wir nochmals reelle Zahlen [latex]a

Sei [latex]\mathbf {{f}}:\mathbb {R}^d\to \mathbb {R}^d[/latex] eine stetige Funktion (siehe Definition 5.96), welche wir als ein Kraftfeld interpretieren und bei jedem Punkt [latex]\mathbf {{v}}\in \mathbb {R}^d[/latex] die Richtung und Stärke einer Krafteinwirkung zum Beispiel auf Grund von Wind angibt (mit jeder Koordinate in [latex]N[/latex]). Wir nennen in diesem Zusammenhang [latex]\mathbf {{f}}[/latex] auch ein Vektorfeld und visualisieren für [latex]d=2[/latex] (und etwas schwieriger auch für [latex]d=3[/latex]) dieses durch eine Ansammlung von Vektoren bei mehreren Punkten im Definitionsbereich, siehe folgendes Bild.

image

Das innere Produkt [latex]\left \langle {\mathbf {{f}}(\gamma (t))}, {\dot {\gamma }(t)} \right \rangle[/latex] gibt damit die Leistung (in [latex]W=Nm/s[/latex]) an, die bei Bewegung mit vorgeschriebener Geschwindigkeit von der Krafteinwirkung zum Zeitpunkt [latex]t[/latex] geleistet wird. Hierbei kann es vorkommen, dass Krafteinwirkung und Geschwindigkeit ähnliche Richtungen haben und das innere Produkt positiv ist. Ebenso kann es aber vorkommen, dass Krafteinwirkung und Geschwindigkeit entgegengesetzt sind und das innere Produkt negativ ist. In diesem Sinne (siehe auch Abschnitt 4.4.4) berechnet das sogenannte Wegintegral

[latex]
\begin{aligned}[]\int _\gamma \mathbf {{f}} \cdot \thinspace {\rm {d}} \mathbf {{s}} = \int _a^b \left \langle {\mathbf {{f}}(\gamma (t))}, {\dot {\gamma }(t)} \right \rangle \thinspace {\rm {d}} t\end{aligned}
[/latex]

die Arbeit, die von der Krafteinwirkung insgesamt geleistet wurde.

Wir werden im zweiten Semester derartige Integrale nochmals genauer untersuchen und dann zum Beispiel folgende Frage beantworten können: Wie kann man einem Kraftfeld [latex]\mathbf {{f}}[/latex] ansehen, ob das Wegintegral nur von Anfangspunkt [latex]\gamma (a)[/latex] und Endpunkt [latex]\gamma (b)[/latex] abhängt und nicht von der Wahl des konkreten Weges von [latex]\gamma (a)[/latex] nach [latex]\gamma (b)[/latex]?

Beispiel 8.35: Abhängigkeit von der Wahl des Weges

Sei [latex]\mathbf {{f}}:\mathbb {R}^2\to \mathbb {R}^2[/latex] definiert durch [latex]\mathbf {{f}}(x,y)= \Big ({\scriptsize \arraycolsep =0.3\arraycolsep \ensuremath {\begin{matrix}y\\ x^2\end{matrix}}}\Big )[/latex]. Wir betrachten den Weg [latex]\gamma :[0,1]\to \mathbb {R}^2[/latex] definiert durch [latex]\gamma (t)=\Big ({\scriptsize \arraycolsep =0.3\arraycolsep \ensuremath {\begin{matrix}t\\ t^2\end{matrix}}}\Big )[/latex] für [latex]t\in [0,1][/latex]. Dann ist das Wegintegral von [latex]\mathbf {{f}}[/latex] über den Weg [latex]\gamma[/latex] von [latex]\gamma (0)=\Big ({\scriptsize \arraycolsep =0.3\arraycolsep \ensuremath {\begin{matrix}0\\ 0\end{matrix}}}\Big )[/latex] nach [latex]\gamma (1)=\Big ({\scriptsize \arraycolsep =0.3\arraycolsep \ensuremath {\begin{matrix}1\\ 1\end{matrix}}}\Big )[/latex] durch

[latex]
\begin{aligned}[]\int _0^1 \left \langle {\begin{pmatrix}t^2\\ t^2\end{pmatrix}}, {\begin{pmatrix}1\\ 2t\end{pmatrix}} \right \rangle \thinspace {\rm {d}} t =\int _0^1 (t^2+2t^3)\thinspace {\rm {d}} t=\frac 13+\frac 24=\frac 56\end{aligned}
[/latex]

gegeben. Verwenden wir allerdings den Weg [latex]\eta :[0,1]\to \mathbb {R}^2[/latex] definert durch [latex]\eta (t)=\Big ({\scriptsize \arraycolsep =0.3\arraycolsep \ensuremath {\begin{matrix}t^2\\ t\end{matrix}}}\Big )[/latex] für [latex]t\in [0,1][/latex], so sind zwar Anfangs- und Endpunkte unverändert, doch ist das Wegintegral durch

[latex]
\begin{aligned}[]\int _0^1 \left \langle {\begin{pmatrix}t\\ t^4\end{pmatrix}}, {\begin{pmatrix}2t\\ 1\end{pmatrix}} \right \rangle \thinspace {\rm {d}} t =\int _0^1 (2t^2+t^4)\thinspace {\rm {d}} t=\frac 23+\frac 15=\frac {13}{15}\end{aligned}
[/latex]

gegeben.

Applet 8.36: Wegintegral

Wir stellen sowohl das Vektorfeld [latex]\mathbf {{f}}[/latex], einen verschiebbaren Weg [latex]\gamma[/latex] mit animiertem Punkt [latex]\gamma (t)[/latex], die Ableitung [latex]\gamma '(t)[/latex] und darunter den Graph der Funktion [latex]t\in [0,1]\mapsto \left \langle {\mathbf {{f}}(\gamma (t))}, {\gamma '(t)} \right \rangle[/latex] dar.

8.3.4 – Volumen von Rotationskörpern

Sei [latex][a,b] \subseteq \mathbb {R}[/latex] ein kompaktes Intervall mit Endpunkten [latex]a

[latex]
\begin{aligned}[]G = \left \lbrace {(x,y) \in \mathbb {R}^2} \mid {a \leq x \leq b,\ 0 \leq y \leq f(x)}\right \rbrace\end{aligned}
[/latex]

und den zugehörigen Körper

[latex]
\begin{aligned}[]K = \left \lbrace {(x,y,z) \in \mathbb {R}^3} \mid {a \leq x \leq b,\ 0 \leq \sqrt {y^2+z^2} \leq f(x)}\right \rbrace ,\end{aligned}
[/latex]

der sich aus Rotation von [latex]G[/latex] um die [latex]x[/latex]-Achse ergibt. Sind die beiden Zylinder [latex]Z_1,Z_2[/latex] mit Radius [latex]\min _{x \in [a,b]} f(x)[/latex] respektive [latex]\max _{x \in [a,b]} f(x)[/latex] um die [latex]x[/latex]-Achse gegeben, so will man wegen den Enthaltungen [latex]Z_1 \subseteq K \subseteq Z_2[/latex], dass das Volumen von [latex]K[/latex] zwischen [latex]\pi \left (\min _{x \in [a,b]} f(x)\right )^2 (b-a)[/latex] und [latex]\pi \left (\max _{x \in [a,b]} f(x)\right )^2 (b-a)[/latex] liegt. Wir halten dies in folgendem Bild fest, wo gemeinsam mit dem Rotationskörper eine von vielen «Scheiben» , die zusammen den Körper approximieren, dargestellt werden.

image

Deswegen (siehe auch Übung 8.38) definieren wir das Volumen des Rotationskörpers [latex]K[/latex] durch
[latex]
\begin{aligned}[]\pi \int _a^b f(x)^2 \thinspace {\rm {d}} x.\end{aligned}
[/latex]

Beispiel 8.37: Volumen der Kugel

Sei [latex]K = \left \lbrace {(x,y,z) \in \mathbb {R}^3} \mid {x^2+y^2+z^2 \leq r^2}\right \rbrace[/latex] für [latex]r>0[/latex] die Kugel mit Radius [latex]r[/latex]. Die Kugel [latex]K[/latex] lässt sich auch als Rotationskörper mittels der Funktion [latex]f:x \in [-r,r] \mapsto \sqrt {r^2-x^2}[/latex] auffassen. Ihr Volumen ist deswegen durch

[latex]
\begin{aligned}[]\pi \int _{-r}^r \left (\sqrt {r^2-x^2}\right )^2 \thinspace {\rm {d}} x = \pi \int _{-r}^r (r^2-x^2) \thinspace {\rm {d}} x = \pi \left [r^2x-\tfrac {x^3}{3} \right ]^r_{-r} &= \pi \left (r^3- \tfrac {r^3}{3}+r^3-\tfrac {r^3}{3}\right )\\ &= \tfrac {4\pi }{3}r^3\end{aligned}
[/latex]

gegeben.

Übung 8.38

Motivieren Sie die Definition des Volumen eines Rotationskörpers mit mehr Details in Analogie zu Abschnitt 8.3.2 unter Verwendung von Proposition 4.30.

8.3.5 – Oberflächen von Rotationskörpern

Obwohl Proposition 4.30 oft ein guter Wegweiser für das Auffinden einer geeigneten Definition darstellt, müssen oder können wir diese nicht immer als Grundlage wählen. Manchmal begnügen wir uns mit geometrischer Intuition als Motivation der Definition.[1]

Wir betrachten [latex]a

[latex]
\begin{aligned}[]K = \left \lbrace {(x,y,z) \in \mathbb {R}^3} \mid {a \leq x \leq b,\ 0 \leq \sqrt {y^2+z^2} \leq f(x)}\right \rbrace\end{aligned}
[/latex]

wie im letzten Abschnitt.

image

In einem kleinen Teilintervall [latex][x_{k-1},x_k] \subseteq [a,b][/latex] der Länge [latex]\triangle x_k[/latex] ist die Funktion [latex]y = f(x)[/latex] der Tangente [latex]y-f(\xi ) = f'(\xi ) (x-\xi )[/latex] für ein [latex]\xi \in [x_{k-1},x_k][/latex] sehr nahe. Ausserdem wird die Oberfläche des Anteils von [latex]K[/latex], der dem Intervall [latex][x_{k-1},x_k][/latex] entspricht, sehr gut durch die Aussenoberfläche des Kegelstumpfs beschrieben, der entsteht, wenn man obiges Tangentenstück zwischen [latex]x_{k-1}[/latex] und [latex]x_k[/latex] um die [latex]x[/latex]-Achse rotiert. Die Aussenoberfläche eines Kegelstumpfs ist näherungsweise [latex]\ell \cdot U[/latex], wobei [latex]\ell[/latex] die Länge der Aussenkante des Kegelstumpfs und [latex]U[/latex] der Umfang einer der beiden Kreise darstellt (wieso?).

Die Oberfläche sollte also näherungsweise durch

[latex]
\begin{aligned}[]\sum _{k=1}^n \sqrt {(\triangle x_k)^2 + (\triangle y_k)^2}\cdot 2 \pi f(x_k) &= 2\pi \sum _{k=1}^n \sqrt {1 + \left (\frac {\triangle y_k}{\triangle x_k}\right )^2} f(x_k) \triangle x_k\\ &= 2\pi \sum _{k=1}^n \sqrt {1 + f'(\xi _k)^2} f(x_k) (x_k-x_{k-1})\end{aligned}
[/latex]

gegeben sein, wobei [latex]\xi _k \in [x_{k-1},x_k][/latex] für jedes [latex]k\in \left \lbrace {1,\ldots ,n} \right \rbrace[/latex] einen Zwischenpunkt darstellt. Deswegen definieren wir nun die Oberfläche des Rotationskörpers [latex]K[/latex] als

[latex]
\begin{aligned}[]2\pi \int _a^b \sqrt {1+f'(x)^2}\ f(x) \thinspace {\rm {d}} x.\end{aligned}
[/latex]

Beispiel 8.39: Kugeloberfläche

Wie in Beispiel 8.37 betrachten wir zu [latex]r > 0[/latex] die Funktion [latex]x \in [-r,r] \mapsto \sqrt {r^2-x^2}[/latex], deren Rotationskörper gerade die Kugel von Radius [latex]r[/latex] ist. Für alle [latex]x \in [-r,r][/latex] ist

[latex]
\begin{aligned}[]f'(x) = -\frac {x}{\sqrt {r^2-x^2}}.\end{aligned}
[/latex]

Damit ist die Kugeloberfläche gleich

[latex]
\begin{aligned}[]2\pi \int _{-r}^r \sqrt {1+ \frac {x^2}{r^2-x^2}} \sqrt {r^2-x^2} \thinspace {\rm {d}} x = 2\pi \int _{-r}^r \sqrt {r^2} \thinspace {\rm {d}} x = 2\pi r [x]_{-r}^r = 4 \pi r^2.\end{aligned}
[/latex]

8.4 – Das uneigentliche Integral

Wir wollen nun den Begriff des Riemann-Integrals auf mehrere Arten erweitern.

8.4.1 – Uneigentliche Integrationsgrenzen

Für [latex]a \in \mathbb {R}[/latex] und eine komplexwertige Funktion [latex]f:[a,\infty ) \to \mathbb {C}[/latex] mit [latex]f|_{[a,b]} \in R([a,b])[/latex] für alle [latex]b > a[/latex] definieren wir das uneigentliche Integral

[latex]
\begin{aligned}[]\int _a^\infty f(x) \thinspace {\rm {d}} x = \lim _{b \to \infty } \int _a^b f(x) \thinspace {\rm {d}} x,\end{aligned}
[/latex]

falls der Grenzwert existiert. Weiter sagen wir, dass das uneigentliche Integral konvergiert, falls der obige Grenzwert in [latex]\mathbb {C}[/latex] existiert. Ansonsten nennen wir das uneigentliche Integral [latex]\int _a^\infty f(x) \thinspace {\rm {d}} x[/latex] divergent.

Beispiel 8.40

Es gilt

[latex]
\begin{aligned}[]\int _0^\infty \frac {1}{1+x^2} \thinspace {\rm {d}} x = \lim _{b \to \infty } \int _0^b \frac {1}{1+x^2} \thinspace {\rm {d}} x = \lim _{b \to \infty } \arctan (b) = \frac {\pi }{2}.\end{aligned}
[/latex]

Beispiel 8.41

Es gilt für [latex]\alpha \in \mathbb {R}[/latex]

[latex]
\begin{aligned}[]\int _1^\infty x^{-\alpha } \thinspace {\rm {d}} x = \left \lbrace \begin{array}{cc} \frac {1}{\alpha -1} & \text {falls } \alpha {\; > \; } 1 \\ +\infty & \text {falls } \alpha \leq 1.\end{array} \right .\end{aligned}
[/latex]

Insbesondere ist das obige uneigentliche Integral genau dann konvergent, wenn [latex]\alpha >1[/latex].

In der Tat ist

[latex]
\begin{aligned}[]\int _1^b x^{-\alpha } \thinspace {\rm {d}} x = \left \lbrace \begin{array}{cl} \left [\frac {1}{- \alpha +1} x^{-\alpha +1}\right ]_1^b = \frac {1}{- \alpha +1} b^{-\alpha +1} - \frac {1}{- \alpha +1} & \text {falls } \alpha \neq 1 \\ \left [\log (x)\right ]_1^b = \log (b) & \text {falls } \alpha = 1\end{array} \right .\end{aligned}
[/latex]

und

[latex]
\begin{aligned}[]\lim _{b \to \infty } \frac {1}{- \alpha +1} b^{-\alpha +1} &= \left \lbrace \begin{array}{cc} +\infty & \text {falls }\alpha {\; \; } 1\end{array} \right .\\ \lim _{b \to \infty } \log (b) &= +\infty .\end{aligned}
[/latex]

Übung 8.42

Berechnen Sie [latex]\int _1^\infty x^{-\alpha } \thinspace {\rm {d}} x[/latex] auch für [latex]\alpha \in \mathbb {C}[/latex].

Uneigentliche Integrale der Form [latex]\int _{-\infty }^b f(x) \thinspace {\rm {d}} x[/latex] sind ähnlich definiert. Ebenso definieren wir für eine komplexwertige Funktion [latex]f: \mathbb {R} \to \mathbb {C}[/latex] mit [latex]f|_{[a,b]} \in R([a,b])[/latex] für alle [latex]a

[latex]
\begin{aligned}[]\int _{-\infty }^\infty f(x) \thinspace {\rm {d}} x = \int _{-\infty }^0 f(x) \thinspace {\rm {d}} x + \int _{0 }^\infty f(x) \thinspace {\rm {d}} x = \lim _{a \to -\infty } \int _{a}^0 f(x) \thinspace {\rm {d}} x + \lim _{b \to \infty } \int _{0}^b f(x) \thinspace {\rm {d}} x,\end{aligned}
[/latex]

falls beide Grenzwerte existieren. Wir möchten dazu anmerken, dass man sich bewusst dazu entscheidet, die Bewegungen gegen [latex]-\infty[/latex] respektive [latex]+\infty[/latex] komplett getrennt zu behandeln. Alles andere würde zu komischen Phänomenen führen, wie folgendes Beispiel zeigt.

Beispiel 8.43

Das uneigentliche Integral [latex]\int _{-\infty }^\infty x \thinspace {\rm {d}} x[/latex] existiert nicht, da [latex]\int _{-\infty }^0 x \thinspace {\rm {d}} x[/latex] sowie [latex]\int _{0}^\infty x \thinspace {\rm {d}} x[/latex] nicht existieren. Wir bemerken aber, dass der Grenzwert [latex]\lim _{c \to \infty } \int _{-c}^c x \thinspace {\rm {d}} x = \lim _{c \to \infty }0 = 0[/latex] existieren würde aber [latex]\lim _{c \to \infty } \int _{-c}^{c+1} x \thinspace {\rm {d}} x = \lim _{c \to \infty }\frac 12((c+1)^2-c^2) = \infty[/latex] wäre.

Wie wir nun besprechen wollen, haben uneigentlichen Integrale oft sehr enge Beziehungen zu Reihen. Genau wie bei Folgen und Reihen (siehe Satz 5.34 und Proposition 6.11) ist es bei uneigentlichen Integralen nicht-negativer Funktionen einfacher über Konvergenz zu entscheiden.

Lemma 8.44

Sei [latex]a \in \mathbb {R}[/latex] und [latex]f:[a,\infty )\to \mathbb {R}_{\geq 0}[/latex] eine nicht-negative Funktion mit [latex]f \in R([a,b])[/latex] für alle [latex]b > a[/latex]. Entweder konvergiert das uneigentliche Integral über [latex]f[/latex] oder es divergiert gegen Unendlich. In beiden Fällen gilt

[latex]
\begin{aligned}[]\int _a^\infty f(x)\thinspace {\rm {d}} x=\sup \left \lbrace {\int _a^bf(x)\thinspace {\rm {d}} x} \mid {b> a}\right \rbrace .\end{aligned}
[/latex]
Beweis

Die Funktion [latex]b\in [a,\infty )\mapsto \int _a^bf(x)\thinspace {\rm {d}} x[/latex] ist monoton wachsend. Wenn das Supremum [latex]S=\sup \left \lbrace {\int _a^bf(x)\thinspace {\rm {d}} x} \mid {b> a}\right \rbrace[/latex] Unendlich ist, dann divergiert das uneigentliche Integral auf Grund der Monotonie gegen Unendlich. Wenn [latex]S 0[/latex] ein [latex]B > a[/latex] mit

[latex]
\begin{aligned}[]S - \varepsilon \leq \int _a^B f(x)\thinspace {\rm {d}} x \leq S\end{aligned}
[/latex]

Insbesondere gilt für [latex]b > B[/latex] auf Grund der Monotonie und der Definition von [latex]S[/latex] dieselbe Ungleichung auch für [latex]\int _a^b f(x)\thinspace {\rm {d}} x[/latex]. Dies beweist die Konvergenz des uneigentlichen Integrals. ∎

Beispiel 8.45: Gaussche Glockenkurve

Wir wollen das uneigentliche Integral

[latex]
\begin{aligned}[]\int _{-\infty }^\infty \mathrm {e}^{-x^2}\thinspace {\rm {d}} x = \int _{-\infty }^{-1} \mathrm {e}^{-x^2}\thinspace {\rm {d}} x+ \int _{-1}^1 \mathrm {e}^{-x^2}\thinspace {\rm {d}} x+ \int _1^\infty \mathrm {e}^{-x^2}\thinspace {\rm {d}} x\end{aligned}
[/latex]

besprechen, wobei die Funktion [latex]x\in \mathbb {R}\mapsto \mathrm {e}^{-x^2}[/latex] die Gaussche Glockenkurve genannt wird. Auf Grund von Lemma 8.44 reicht es aus eine «Majorantenfunktion» zu finden, die ein konvergentes uneigentliches Integral definiert. Für [latex]x\in [1,\infty )[/latex] gilt zum Beispiel [latex]x^2\geq x[/latex] und daher [latex]\mathrm {e}^{-x^2}\leq \mathrm {e}^{-x}[/latex], woraus

[latex]
\begin{aligned}[]\int _1^\infty \mathrm {e}^{-x^2}\thinspace {\rm {d}} x\leq \int _1^\infty \mathrm {e}^{-x}\thinspace {\rm {d}} x[/latex]

folgt. Dies zeigt die Konvergenz des zweiten uneigentlichen Integrals, auf Grund der Symmetrie der Funktion ist daher auch [latex]\int _{-\infty }^\infty \mathrm {e}^{-x^2}\thinspace {\rm {d}} x[/latex] konvergent. Wir werden den Wert [latex]I[/latex] dieses Integral erst im zweiten Semester berechnen können. Doch wollen wir noch erwähnen, dass die streng monoton wachsende Funktion

[latex]
\begin{aligned}[]\Phi :x\in \mathbb {R}\mapsto I^{-1}\int _{-\infty }^x\mathrm {e}^{-t^2}\thinspace {\rm {d}} t\end{aligned}
[/latex]

die Verteilungsfunktion der Normalverteilung mit Erwartungswert [latex]0[/latex] und Standardabweichung [latex]\frac 1{\sqrt 2}[/latex] genannt wird. Diese Funktion lässt sich ebenso wie die Funktionen aus Abschnitt 8.2.6 nicht durch die sonst üblichen Funktionen ausdrücken und ist in der Wahrscheinlichkeitsrechnung, der Statstik und in vielen Anwendungen von fundamentaler Bedeutung.

Der folgende Satz charakterisiert nun Konvergenz uneigentlicher Integrale wie in obigem Lemma durch Konvergenz von Reihen (auf hinreichende und notwendige Weise).

Satz 8.46: Integraltest für Reihen

Sei [latex]f:[1,\infty ) \to \mathbb {R}_{\geq 0}[/latex] eine monoton fallende Funktion. Dann gilt

[latex]
\begin{aligned}[]\sum _{n=2}^\infty f(n) \leq \int _1^\infty f(x) \thinspace {\rm {d}} x \leq \sum _{n=1}^\infty f(n).\end{aligned}
[/latex]

Insbesondere konvergiert die Reihe [latex]\sum _{n=1}^\infty f(n)[/latex] genau dann, wenn das uneigentliche Integral [latex]\int _1^\infty f(x) \thinspace {\rm {d}} x[/latex] konvergiert. Dies gilt analog für Integrale der Form [latex]\int _N^\infty f(x)\thinspace {\rm {d}} x[/latex] für [latex]N \in \mathbb {N}[/latex].

Wir bemerken, dass auf Grund der Monotonieannahme an [latex]f[/latex] in obigem Satz die Eigenschaft [latex]f|_{[1,b]} \in R([1,b])[/latex] für alle [latex]b > 1[/latex] erfüllt ist nach Satz 4.31.

Beweis

Für [latex]n \in \mathbb {N}[/latex] und einen beliebigen Zwischenpunkt [latex]x \in [n,n+1][/latex] gilt nach Monotonie von [latex]f[/latex] die Ungleichung [latex]f(n+1) \leq f(x) \leq f(n)[/latex] und somit

[latex]
\begin{aligned}[]f(n+1) \leq \int _n^{n+1} f(x) \thinspace {\rm {d}} x \leq f(n),\end{aligned}
[/latex]

was auch in folgendem Bild ersichtlich ist.

image

Nach Summation von [latex]1[/latex] bis [latex]n[/latex] erhält man mit Intervalladditivität des Riemann-Integrals

[latex]
\begin{aligned}[]\sum _{\ell =2}^{n+1} f(\ell ) = \sum _{k=1}^n f(k+1) \leq \int _1^{n+1} f(x) \thinspace {\rm {d}} x \leq \sum _{k=1}^n f(k).\end{aligned}
[/latex]

Falls das uneigentliche Integral [latex]\int _1^{\infty } f(x) \thinspace {\rm {d}} x[/latex] existiert, dann folgt

[latex]
\begin{aligned}[]\sum _{\ell =2}^{n+1} f(\ell ) \leq \int _1^{n+1}f(x) \thinspace {\rm {d}} x \leq \int _1^{\infty }f(x) \thinspace {\rm {d}} x.\end{aligned}
[/latex]

Daher ist die monoton wachsende Folge [latex]\left (\sum _{\ell =2}^{n+1} f(\ell )\right )_n[/latex] nach oben beschränkt und konvergiert somit nach Satz 5.34. Insbesondere gilt auch [latex]\sum _{\ell =2}^{\infty } f(\ell )\leq \int _1^{\infty }f(x) \thinspace {\rm {d}} x[/latex].

Falls [latex]\sum _{k=1}^\infty f(k)[/latex] konvergiert, dann ist für [latex]b > 1[/latex] und [latex]n = \lfloor {b} \rfloor[/latex]

[latex]
\begin{aligned}[]\int _1^b f(x) \thinspace {\rm {d}} x \leq \int _1^{n+1} f(x) \thinspace {\rm {d}} x \leq \sum _{k=1}^\infty f(k).\end{aligned}
[/latex]

Nach Lemma 8.44 ist somit das uneigentliche Integral [latex]\int _1^\infty f(x) \thinspace {\rm {d}} x[/latex] konvergent und durch die Zahl [latex]\sum _{k=1}^\infty f(k)[/latex] beschränkt. ∎

Übung 8.47: Divergenzrate der harmonischen Reihe

Verwenden Sie obigen Satz, um den [latex]p[/latex]-Test in Beispiel 6.17 zu erhalten. Imitieren Sie des Weiteren die Methodik im obigen Beweis von Satz 8.46, um die Divergenzrate

[latex]
\begin{aligned}[]\sum _{n=1}^N \frac {1}{n} = \log (N) + O(1)\end{aligned}
[/latex]

für [latex]N \to \infty[/latex] für die harmonische Reihe zu beweisen.

Übung 8.48: Eine lange Nadel

Berechnen Sie das Volumen und die Oberfläche des «un­eig­ent­lichen Rotations­körpers» , der entsteht, wenn man das Gebiet unter dem Graphen der Funktion [latex]x \in [1,\infty ) \mapsto \frac {1}{x}[/latex] um die [latex]x[/latex]-Achse rotiert.

Übung 8.49: Ein oszillierendes Integral

Entscheiden Sie für welche [latex]p \in \mathbb {R}_{\geq 0}[/latex] das uneigentliche Integral [latex]\int _0^\infty x \sin (x^p) \thinspace {\rm {d}} x[/latex] konvergiert.

Hinweis.

Für [latex]p > 0[/latex] verwenden Sie am besten die Substitution [latex]u = x^p[/latex] für das Integral [latex]\int _1^b x\sin (x^p)\thinspace {\rm {d}} x[/latex] und [latex]b>1[/latex] und ziehen je nach Fall auf geeignete Weise das Leibniz-Kriterium hinzu.

8.4.2 – Das Integral über unbeschränkte Funktionen

Für [latex]a uneigentliche Integral

[latex]
\begin{aligned}[]\int _a^B f(x)\thinspace {\rm {d}} x = \lim _{b \nearrow B} \int _a^b f(x) \thinspace {\rm {d}} x,\end{aligned}
[/latex]

falls der Grenzwert existiert.

Wie folgende Übung zeigt, steht diese Notation nicht im Widerspruch zum Riemann-Integral.

Übung 8.50: Kompatibilität

Sei [latex]f:[a,B) \to \mathbb {C}[/latex] wie oben. Angenommen [latex]f[/latex] ist beschränkt. Zeigen Sie, dass das uneigentliche Integral [latex]\int _a^B f(x)\thinspace {\rm {d}} x[/latex] existiert und gleich dem Riemann-Integral [latex]\int _a^B f(x) \thinspace {\rm {d}} x[/latex] ist, wobei man [latex]f[/latex] auf beliebige Weise auf den Punkt [latex]B[/latex] erweitert.

Für [latex]A

[latex]
\begin{aligned}[]\int _A^b f(x) \thinspace {\rm {d}} x = \lim _{a \searrow A} \int _a^b f(x) \thinspace {\rm {d}} x.\end{aligned}
[/latex]

Beispiel 8.51

Wir berechnen [latex]\int _0^1 \log (x) \thinspace {\rm {d}} x[/latex] mittels

[latex]
\begin{aligned}[]\int _0^1\log (x) \thinspace {\rm {d}} x = \lim _{a \searrow 0} \int _a^1 \log (x) \thinspace {\rm {d}} x &= \lim _{a\searrow 0} \left [x \log (x) -x \right ]_a^1 \\ &= \lim _{a\searrow 0} (\log (1)-1-a \log (a)+a) = -1\end{aligned}
[/latex]

nach Beispiel 8.15 und Beispiel 5.76.

Weitere uneigentliche Integrale führen wir mittels Intervalladditivität auf obige uneigentliche Integrale zurück. Wir überlassen es Interessierten, sich hier einige Möglichkeiten auszudenken, und führen stattdessen ein Beispiel vor.

Beispiel 8.52

Wir betrachten das Integral [latex]\int _{-1}^1 \frac {1}{x} \thinspace {\rm {d}} x[/latex]. Dieses ist uneigentlich, da [latex]x \mapsto \frac {1}{x}[/latex] auf jeder Umgebung von [latex]0[/latex] unbeschränkt ist. Es gilt daher auf Grund der Definition in diesem Fall, dass

[latex]
\begin{aligned}[]\int _{-1}^1 \frac {1}{x} \thinspace {\rm {d}} x = \int _{-1}^0 \frac {1}{x} \thinspace {\rm {d}} x + \int _{0}^1 \frac {1}{x} \thinspace {\rm {d}} x,\end{aligned}
[/latex]

wobei beide Integrale rechts uneigentlich sind und das uneigentliche Integral [latex]\int _{-1}^1 \frac {1}{x} \thinspace {\rm {d}} x[/latex] per Definition genau dann existiert, wenn die beiden Integrale rechts existieren. Des Weiteren gilt

[latex]
\begin{aligned}[]\int _{-1}^0 \frac {1}{x} \thinspace {\rm {d}} x &= \lim _{b \nearrow 0} \int _{-1}^b \frac {1}{x} \thinspace {\rm {d}} x = \lim _{b \nearrow 0} \log |b| - \log |-1| = \lim _{b \nearrow 0} \log |b| = - \infty \\ \int _0^1 \frac {1}{x}\thinspace {\rm {d}} x &= \lim _{a \searrow 0} \int _a^1 \frac {1}{x} \thinspace {\rm {d}} x = \lim _{a \searrow 0} \log (1) - \log (a) = +\infty ,\end{aligned}
[/latex]

wodurch [latex]\int _{-1}^1 \frac {1}{x} \thinspace {\rm {d}} x[/latex] nicht existiert (und wir diesem auch nicht das Symbol [latex]\infty[/latex] oder [latex]-\infty[/latex] zuweisen).

Beispiel 8.53: Bogenlänge des Kreises

Wir wollen nochmals die Bogenlänge des Kreises berechnen. Doch verwenden wir diesmal die Gleichung [latex]y=f(x)=\sqrt {1-x^2}[/latex] als Definition des oberen Halbkreises. Die Bogenlänge des Kreises ist demnach gegeben durch das Integral

[latex]
\begin{aligned}[]2\int _{-1}^1\sqrt {1+f'(x)^2}\thinspace {\rm {d}} x&=4\int _{0}^1\sqrt {1+\frac {x^2}{1-x^2}}\thinspace {\rm {d}} x\\ &=4\lim _{b\nearrow 1}\int _{0}^b\frac 1{\sqrt {1-x^2}}\thinspace {\rm {d}} x\\ &=4\lim _{b\nearrow 1}\arcsin (b)-0=4\tfrac \pi 2=2\pi .\end{aligned}
[/latex]

Übung 8.54

Berechnen Sie [latex]\int _0^1 \frac {1}{\sqrt {x}} \thinspace {\rm {d}} x[/latex] und [latex]\int _0^{\frac {\pi }{2}} \tan (x) \thinspace {\rm {d}} x[/latex].

Übung 8.55: Absolute Konvergenz

Sei [latex]a \in \mathbb {R}[/latex] und [latex]f:[a,\infty )\to \mathbb {C}[/latex] eine komplexwertige Funktion mit [latex]f \in R([a,b])[/latex] für alle [latex]b > a[/latex]. Wir nennen das uneigentliche Integral [latex]\int _a^\infty f(x) \thinspace {\rm {d}} x[/latex] absolut konvergent, falls [latex]\int _a^\infty |f(x)| \thinspace {\rm {d}} x[/latex] konvergent ist. Zeigen Sie, dass absolute Konvergenz des uneigentlichen Integrals [latex]\int _a^\infty f(x) \thinspace {\rm {d}} x[/latex] auch die Konvergenz dieses Integrals impliziert.

Bemerkung

Zusammenfassend haben wir bei den Definitionen in diesem Abschnitt bei jedem Problempunkt eines möglichen Riemann-Integrals einen Grenzwert verwendet, um den Integralbegriff zu erweiteren. Dies wirft nochmals die Frage auf, ob es nicht vielleicht einen Integralbegriff gibt, der diese und auch andere bereits erwähnte Probleme des Riemann-Integrals auf natürliche Art und Weise löst. Diese Frage wird im zweiten Studienjahr des Mathematikstudiums mit der Theorie des Lebesgue-Integrals in der Vorlesung «Mass und Integral» positiv beantwortet.

8.4.3 – Die Gamma-Funktion*

Die Gamma-Funktion [latex]\Gamma[/latex] ist bei [latex]s\in (0,\infty )[/latex] durch das konvergente uneigentliche Integral
[latex]
\begin{aligned}[]\label{eq:fund-def gammafct} \Gamma (s) = \int _0^\infty x^{s-1} \mathrm {e}^{-x} \thinspace {\rm {d}} x\end{aligned}
[/latex]
definiert. Für [latex]s\in (0,1)[/latex] ist dies aus zwei Gründen ein uneigentliches Integral und wir müssen die Integrationsgrenzen [latex]A=0[/latex] und [latex]B=\infty[/latex] getrennt untersuchen. Für [latex]a > 0[/latex] und [latex]b > a[/latex] gilt jedoch

[latex]
\begin{aligned}[]\int _a^b x^{s-1} \mathrm {e}^{-x} \thinspace {\rm {d}} x = \tfrac {1}{s} \left [ x^{s} \mathrm {e}^{-x} \right ]_{a}^b + \tfrac {1}{s} \int _{a}^b x^{s} \mathrm {e}^{-x} \thinspace {\rm {d}} x.\end{aligned}
[/latex]

Wir setzen [latex]b=1[/latex] und erhalten

[latex]
\begin{aligned}[]\int _0^1 x^{s-1} \mathrm {e}^{-x} \thinspace {\rm {d}} x =\lim _{a\searrow 0}\left (\tfrac {1}{s} \left [ x^{s} \mathrm {e}^{-x} \right ]_{a}^1 + \tfrac {1}{s} \int _{a}^1 x^{s} \mathrm {e}^{-x} \thinspace {\rm {d}} x\right ) =\tfrac {1}{s\mathrm {e}}+\tfrac {1}{s} \int _{0}^1 x^{s} \mathrm {e}^{-x} \thinspace {\rm {d}} x,\end{aligned}
[/latex]

wobei das Integral rechts (für alle [latex]s\in (0,\infty )[/latex]) ein eigentliches Riemann-Integral darstellt. Für [latex]a=1[/latex] erhalten wir

[latex]
\begin{aligned}[]\int _1^\infty x^{s-1} \mathrm {e}^{-x} \thinspace {\rm {d}} x = \lim _{b\to \infty }\left (\tfrac {1}{s} \left [ x^{s} \mathrm {e}^{-x} \right ]_{1}^b + \tfrac {1}{s} \int _{1}^b x^{s} \mathrm {e}^{-x} \thinspace {\rm {d}} x\right ) =-\tfrac {1}{s\mathrm {e}}+\tfrac {1}{s} \int _{1}^\infty x^{s} \mathrm {e}^{-x} \thinspace {\rm {d}} x.\end{aligned}
[/latex]

Um die Konvergenz von [latex]\int _1^\infty x^{s} \mathrm {e}^{-x} \thinspace {\rm {d}} x[/latex] zu zeigen, wollen wir den Integraltest für Reihen in Satz 8.46 verwenden. Die erste Vorraussetzung des Integraltests ist erfüllt, da die Funktion [latex]f(x)=x^s\mathrm {e}^{-x}[/latex] nichtnegativ ist. Die zweite Vorraussetzung ist, dass [latex]f(x)[/latex] monoton abnehmend sein soll. Wir berechnen daher die Ableitung und sehen, dass

[latex]
\begin{aligned}[]f'(x)=sx^{s-1}\mathrm {e}^{-x}-x^s\mathrm {e}^{-x}.\end{aligned}
[/latex]

Da [latex]f'(x)s[/latex], sehen wir, dass diese Vorraussetzung zumindest für [latex]x\geq N=\lfloor s\rfloor +1[/latex] erfüllt ist. Es folgt daher, dass das Integral [latex]\int _N^\infty x^{s} \mathrm {e}^{-x} \thinspace {\rm {d}} x[/latex] genau dann konvergiert wenn

[latex]
\begin{aligned}[]\sum _{n=N}^\infty n^{s}\mathrm {e}^{-n}[/latex]

konvergiert, was aber nach dem Quotientenkriterium in Korollar 6.32 (oder dem Wurzelkriterium in Korollar 6.30) für Reihen in der Tat zutrifft.

Addieren wir die beiden Integrale wieder und verwenden wir die Definition in (8.8) so erhalten wir, dass [latex]\Gamma (s)[/latex] für alle [latex]s\in (0,\infty )[/latex] wohldefiniert ist und
[latex]
\begin{aligned}[]\label{eq:fund-prop1 gammafct} \Gamma (s+1) = s\Gamma (s)\end{aligned}
[/latex]
erfüllt.

Oft wird die Gamma-Funktion als eine Erweiterung der Fakultätsfunktion auf [latex]\mathbb {N}[/latex] bezeichnet. In der Tat gilt für alle [latex]n \in \mathbb {N}[/latex]

[latex]
\begin{aligned}[]\Gamma (n+1) = n!\ ,\end{aligned}
[/latex]

was wir hier noch beweisen möchten. Für [latex]n = 0[/latex] gilt

[latex]
\begin{aligned}[]\Gamma (1) = \int _0^\infty \mathrm {e}^{-x} \thinspace {\rm {d}} x = \lim _{b \to \infty } -\left [\mathrm {e}^{-x}\right ]_{0}^b = 1\end{aligned}
[/latex]

und somit folgt nach (8.9) und Induktion

[latex]
\begin{aligned}[]\Gamma (n+1) = n \Gamma (n) = n (n-1) \Gamma (n-1) = \ldots = n!\ \Gamma (1) = n!\ .\end{aligned}
[/latex]

Wir werden in der Fortsetzung dieser Vorlesung weitere Eigenschaften der Gamma-Funktion nachweisen können. Beispielsweise stellt sich heraus, dass die Gamma-Funktion glatt ist. Wir können dies hier aber nicht zeigen, da [latex]\Gamma[/latex] in (8.8) durch ein sogenanntes Parameterintegral definiert ist. Auch können wir den Wert

[latex]
\begin{aligned}[]\Gamma (\tfrac {1}{2}) &= \int _0^\infty \frac {1}{\sqrt {x}} \mathrm {e}^{-x} \thinspace {\rm {d}} x = \lim _{\varepsilon \to 0} \lim _{b \to \infty } \int _\varepsilon ^b \frac {1}{\sqrt {x}} \mathrm {e}^{-x} \thinspace {\rm {d}} x = \lim _{\varepsilon \to 0} \lim _{b \to \infty } 2\int _{\sqrt {\varepsilon }}^{\sqrt {b}} \mathrm {e}^{-u^2} \thinspace {\rm {d}} u\\ &= 2\int _0^\infty \mathrm {e}^{-u^2} \thinspace {\rm {d}} u = \int _{-\infty }^\infty \mathrm {e}^{-u^2} \thinspace {\rm {d}} u\end{aligned}
[/latex]

mit den uns bis jetzt bekannten Integrationsmethoden nicht berechnen, werden jedoch später mittels einem zweidimensionalen Integral sehen, dass dieser [latex]\sqrt {\pi }[/latex] ist.

Die Gamma-Funktion enthüllt ihre wahre Schönheit erst wenn man komplexe Parameter [latex]z\in \mathbb {C}\setminus \{ 0,-1,-2,\ldots \}[/latex] erlaubt. Wir laden Interessierte ein, diese Funktion in folgender Übung zu konstruieren.

Übung 8.56: Challenge

Für [latex]z\in \mathbb {C}[/latex] mit [latex]\Re (z)>0[/latex] definiert man
[latex]
\begin{aligned}[]\label{eq:fund-def gammafct-komplex} \Gamma (z) = \int _0^\infty x^{z-1} \mathrm {e}^{-x} \thinspace {\rm {d}} x.\end{aligned}
[/latex]

(a) Zeigen Sie, dass [latex]\int _0^1 x^{z-1} \mathrm {e}^{-x} \thinspace {\rm {d}} x[/latex] für alle [latex]z\in \mathbb {C}[/latex] mit [latex]\Re (z)>0[/latex] konvergiert.

(b) Zeigen Sie, dass [latex]\int _1^\infty x^{z-1} \mathrm {e}^{-x} \thinspace {\rm {d}} x[/latex] für alle [latex]z\in \mathbb {C}[/latex] mit [latex]\Re (z)>0[/latex] konvergiert.

(c) Zeigen Sie (8.9) für alle [latex]z\in \mathbb {C}[/latex] mit [latex]\Re (z)>0[/latex].

(d) Verwenden Sie (8.9) um rekursiv [latex]\Gamma (z)[/latex] für [latex]z\in \mathbb {C}\setminus \{ 0,-1,-2,\ldots \}[/latex] mit [latex]\Re (z)>-1[/latex], oder [latex]\Re (z)>-2[/latex], …, zu definieren, so dass anschliessend [latex]\Gamma (z)[/latex] für alle [latex]z\in \mathbb {C}\setminus \{ 0,-1,-2,\ldots \}[/latex] definiert ist und (8.9) auf dem ganzen Definitionsbereich erfüllt.

Hinweis: Sie können in (a) und (b) Übung 8.55 verwenden.

Hilbert (1862–1943) verwendete in seinem Artikel [2]von 1893 uneigentliche Integrale im Stile der Gamma-Funktion, um zu beweisen, dass [latex]\mathrm {e}[/latex] (wie erstmals von Hermite in 1873 bewiesen) und [latex]\pi[/latex] (wie erstmals von Lindemann 1882 bewiesen) transzendent sind. Wir bemerken dabei, dass sich die blosse Irrationalität dieser Zahlen deutlich einfacher beweisen lässt — für [latex]\mathrm {e}[/latex] gibt es hierzu eine Übung in Abschnitt 6.9.2 und für [latex]\pi[/latex] eine Übung in Abschnitt 8.8.2. Transzendenzbeweise sind jedoch im Allgemeinen deutlich schwieriger. Wie schwierige derartige Aussagen tatsächlich sind, illustriert vielleicht die Tatsache, dass immer noch nicht bekannt ist, ob [latex]\mathrm {e}+\pi[/latex] eine transzendente Zahl ist oder nicht. Hilbert’s Beweis der Transzendenz von [latex]\mathrm {e}[/latex] und [latex]\pi[/latex] ist mit den uns bisher bekannten Hilfsmitteln allerdings gut lesbar, weswegen wir Ihnen einen Blick auf diese Lektüre und die damit verbundene Zeitreise empfehlen möchten.

8.5 – Taylor Approximation

Wir erinnern daran, dass die Ableitung [latex]f'(x_0)[/latex] einer reellwertigen differenzierbaren Funktion [latex]f[/latex] auf einem Intervall die Steigung der Tangente

[latex]
\begin{aligned}[]y(x)=f(x_0)+f'(x_0)(x-x_0)\end{aligned}
[/latex]

des Graphen von [latex]f[/latex] bei [latex]x_0[/latex] angibt, wobei die Tangente den Graphen gut approximiert (im Sinne von [latex]f(x)=y(x)+o(x-x_0)[/latex] für [latex]x\to x_0[/latex]). Wir wollen hier die Güte der Approximation erhöhen indem wir statt linearen Approximationen auch noch höhere polynomiale Approximationen erlauben.

Sei also [latex]f: (a,b) \to \mathbb {C}[/latex] eine [latex]n[/latex]-mal differenzierbare Funktion auf einem offenen, nicht-leeren (möglicherweise unbeschränkten) Intervall [latex](a,b) \subseteq \mathbb {R}[/latex]. Die [latex]\mathbf {n}[/latex]-te Taylor-Approximation von [latex]f[/latex] um einen Punkt [latex]x_0 \in (a,b)[/latex] ist das Polynom

[latex]
\begin{aligned}[]P_{x_0,n}^{f}(x) = \sum _{k=0}^n \frac {f^{(k)}(x_0)}{k!} (x-x_0)^k.\end{aligned}
[/latex]

Die Koeffizienten wurden dabei so gewählt, dass [latex]P_{x_0,n}^{f}(x_0) = f(x_0)[/latex] und allgemeiner

[latex]
\begin{aligned}[](P_{x_0,n}^{f})^{(k)}(x_0) = f^{(k)}(x_0)\end{aligned}
[/latex]

für [latex]k \in \left \lbrace {0,\ldots ,n} \right \rbrace[/latex] gilt (wieso?). Die vielleicht naive Hoffnung ist hierbei, dass die Taylor-Approximation, wie der Name sagt, die Funktion [latex]f[/latex] approximiert. Falls [latex]f[/latex] glatt ist, dann ist die Taylorreihe von [latex]f[/latex] um [latex]x_0 \in (a,b)[/latex] definiert als die Potenzreihe

[latex]
\begin{aligned}[]T_{x_0}^f(x) = \sum _{k=0}^\infty \frac {f^{(k)}(x_0)}{k!} (x-x_0)^k.\end{aligned}
[/latex]

Damit wir hier von einer Potenzreihe sprechen dürfen, erweitern wir die Definition 6.54 wie folgt. Eine Potenzreihe um einen Punkt [latex]z_0 \in \mathbb {C}[/latex] in der Variable [latex]z \in \mathbb {C}[/latex] ist ein formaler Ausdruck der Form [latex]\sum _{n=0}^\infty a_n (z-z_0)^n[/latex], wobei [latex]a_n \in \mathbb {C}[/latex] für alle [latex]n \in \mathbb {N}[/latex] die Koeffizienten der Potenzreihe sind. Der Konvergenzradius

[latex]
\begin{aligned}[]R = \left (\limsup _{n \to \infty } \sqrt [n]{|a_n|}\right )^{-1}\end{aligned}
[/latex]

ist wie in Abschnitt 6.4.1 definiert und die Potenzreihe [latex]\sum _{n=0}^\infty a_n (z-z_0)^n[/latex] konvergiert für alle [latex]z \in \mathbb {C}[/latex] mit Abstand kleiner [latex]R[/latex] von [latex]z_0[/latex] und divergiert für alle [latex]z \in \mathbb {C}[/latex] mit [latex]|z-z_0| > R[/latex] (was aus Satz 6.56 folgt, wie?).

Soweit ist nicht klar, was der Konvergenzradius der Taylor-Reihe einer glatten Funktion [latex]f[/latex] um einen Punkt [latex]x_0[/latex] im Definitionsbereich ist und ob er positiv ist. Die naive Hoffnung ist, dass die Taylor-Reihe, wo definiert, gleich der Funktion [latex]f[/latex] ist, da die Taylorreihe bei [latex]x_0[/latex] die gleichen Ableitungen wie [latex]f[/latex] hat. Dies ist in der Tat für viele Funktionen der Fall, doch, wie folgendes Beispiel zeigt, nicht immer.

Beispiel 8.57: Verschwindende Taylorreihe

Wir betrachten die Funktion

[latex]
\begin{aligned}[]\psi : x \in \mathbb {R} \mapsto \left \lbrace \begin{array}{cc} \exp \left (-\frac 1x\right ) & \text {falls } x > 0 \\ 0 & \text {falls } x \leq 0\end{array} \right . ,\end{aligned}
[/latex]

die nach Beispiel 7.23 glatt ist und [latex]\psi ^{(n)}(0) = 0[/latex] für alle [latex]n \in \mathbb {N}_0[/latex] erfüllt. Wir zeigen, dass sich [latex]\psi[/latex] nicht durch eine Potenzreihe um [latex]0[/latex] darstellen lässt.

Wir nehmen also indirekt an, dass es eine Potenzreihe [latex]f(x) = \sum _{n=0}^\infty c_n x^n[/latex] mit reellen Koeffizienten und Konvergenzradius [latex]R > 0[/latex] gibt, so dass [latex]f(x) = \psi (x)[/latex] für alle [latex]x \in (-R,R)[/latex]. Nach Korollar 8.11 (oder genauer Übung 8.12) ist [latex]f^{(n)}(0) = c_n n![/latex]. Aber da [latex]f(x) = \psi (x)[/latex] für alle [latex]x \in (-R,R)[/latex] angenommen wurde und da Ableiten eine lokale Operation ist, schliessen wir, dass [latex]f^{(n)}(0) = 0[/latex] und somit [latex]c_n = 0[/latex] für alle [latex]n \in \mathbb {N}_0[/latex]. Dies widerspricht jedoch [latex]\psi (x) > 0[/latex] für alle [latex]x > 0[/latex]. Also kann [latex]\psi[/latex] in keiner Umgebung von [latex]0[/latex] durch eine Potenzreihe dargestellt werden.

Folgender Satz liefert nun einen direkten Vergleich zwischen einer Funktion und ihrer Taylor-Approximationen. Er impliziert für viele Funktionen, dass sie mit ihrer Taylorreihe übereinstimmen.

Theorem 8.58: Taylor-Approximation

Sei [latex](a,b) \subseteq \mathbb {R}[/latex] ein offenes, nicht-leeres Intervall und sei [latex]f:(a,b) \to \mathbb {C}[/latex] eine [latex](n+1)[/latex]-mal stetig differenzierbare Funktion. Sei [latex]x_0 \in (a,b)[/latex]. Dann gilt für alle [latex]x \in (a,b)[/latex]

[latex]
\begin{aligned}[]f(x) = P_{x_0,n}^{f}(x) + R_{x_0,n}^{f}(x),\end{aligned}
[/latex]

wobei [latex]P_{x_0,n}^{f}[/latex] die [latex]n[/latex]-te Taylor-Approximation ist und wir den Fehlerterm [latex]R_{x_0,n}^{f}[/latex] durch das sogenannte Integral-Restglied

[latex]
\begin{aligned}[]x \in (a,b) \mapsto R_{x_0,n}^{f}(x) = \int _{x_0}^x f^{(n+1)}(t) \frac {(x-t)^n}{n!} \thinspace {\rm {d}} t\end{aligned}
[/latex]

darstellen können. Dies gilt auch für Funktionen auf [latex][x_0,b)[/latex] und Punkte [latex]x \in [x_0,b)[/latex] (beziehungsweise [latex](a,x_0][/latex] und [latex]x \in (a,x_0][/latex]).

Die Annahme im Theorem, dass [latex]f[/latex] eine [latex](n+1)[/latex]-mal stetig differenzierbare Funktion ist, ist essentiell für obige Formulierung. In der Tat ist damit [latex]f^{(n+1)}[/latex] eine stetige Funktion und das Integral der stetigen Funktion [latex]t \mapsto f^{(n+1)}(t) \frac {(x-t)^n}{n!}[/latex] im Integral-Restglied existiert.

Beweis von Satz 8.58

Das Theorem ergibt sich mit Induktion über [latex]n[/latex] und partieller Integration (die ebenso für komplexwertige Funktionen auf [latex](a,b)[/latex] gilt, wieso?). Ist [latex]n = 0[/latex] und [latex]f:(a,b) \to \mathbb {R}[/latex] stetig differenzierbar, so gilt nach dem Fundamentalsatz der Differential- und Integralrechnung (genauer Korollar 8.5)

[latex]
\begin{aligned}[]f(x) = f(x_0) + \int _{x_0}^x f'(t) \thinspace {\rm {d}} t = P_{x_0,0}^{f}(x) + R_{x_0,0}^{f}(x).\end{aligned}
[/latex]

Ist [latex]f[/latex] nun zweimal stetig differenzierbar (also [latex]n=1[/latex]), so können wir auf obiges Integral partielle Integration mit [latex]u(t) = f'(t)[/latex] und [latex]v(t) = t-x[/latex] anwenden und erhalten

[latex]
\begin{aligned}[]f(x)&= f(x_0) + \left [f'(t)(t-x)\right ]_{t=x_0}^{t=x}- \int _{x_0}^x f''(t) (t-x) \thinspace {\rm {d}} t\\ &= f(x_0) + f'(x_0) (x-x_0) + \int _{x_0}^x f''(t) \frac {(x-t)^1}{1!}\thinspace {\rm {d}} t\\ &= P_{x_0,1}^{f}(x) + R_{x_0,1}^{f}(x).\end{aligned}
[/latex]

Wir sehen also wie sich ausgehend vom Fundamentalsatz mittels partieller Integration der nächste Fall des Satzes ergibt.

Angenommen die Aussage des Satzes stimmt für [latex]n-1\geq 0[/latex] und sei [latex]f:(a,b) \to \mathbb {R}[/latex] eine [latex](n+1)[/latex]-mal stetig differenzierbare Funktion. Dann gilt auf Grund der Induktionsvorraussetzung

[latex]
\begin{aligned}[]f(x) = \sum _{k=0}^{n-1} \frac {f^{(k)}(x_0)}{k!} (x-x_0)^k + \int _{x_0}^x f^{(n)}(t)\frac {(x-t)^{n-1}}{(n-1)!} \thinspace {\rm {d}} t\end{aligned}
[/latex]

für alle [latex]x \in (a,b)[/latex].

Wir setzen [latex]u(t) = f^{(n)}(t)[/latex] und [latex]v(t) = - \frac {(x-t)^n}{n!}[/latex], bemerken [latex]v'(t) = \frac {(x-t)^{n-1}}{(n-1)!}[/latex] und wenden partielle Integration an, um

[latex]
\begin{aligned}[]f(x) &= \sum _{k=0}^{n-1} \frac {f^{(k)}(x_0)}{k!} (x-x_0)^k - \left [f^{(n)}(t)\frac {(x-t)^n}{n!}\right ]_{t=x_0}^{t=x} + \int _{x_0}^x f^{(n+1)}(t) \frac {(x-t)^n}{n!} \thinspace {\rm {d}} t\\ &= \sum _{k=0}^{n} \frac {f^{(k)}(x_0)}{k!} (x-x_0)^k + \int _{x_0}^x f^{(n+1)}(t) \frac {(x-t)^n}{n!} \thinspace {\rm {d}} t\end{aligned}
[/latex]

zu erhalten. Dies beweist den Induktionsschritt und damit den Satz. ∎

Oft werden wir das Theorem von Taylor (Theorem 8.58) in folgender Form verwenden.

Korollar 8.59: Taylor-Abschätzung

Sei [latex](a,b) \subseteq \mathbb {R}[/latex] ein offenes, nicht-leeres Intervall, sei [latex]f:(a,b) \to \mathbb {C}[/latex] eine [latex](n+1)[/latex]-mal stetig differenzierbare Funktion und seien [latex]x_0,x \in (a,b)[/latex] zwei Punkte. Wir setzen [latex]M_{n+1} = \max \left \lbrace {|f^{(n+1)}(t)|} \mid {t \text { zwischen } x_0 \text { und } x}\right \rbrace[/latex]. Dann gilt

[latex]
\begin{aligned}[]|f(x) - P_{x_0,n}^{f}(x)| = |R_{x_0,n}^{f}(x)| \leq \frac {M_{n+1} |x-x_0|^{n+1}}{(n+1)!}.\end{aligned}
[/latex]

Insbesondere ist [latex]f(x) = P_{x_0,n}^{f}(x) + O((x-x_0)^{n+1})[/latex] für [latex]x \to x_0[/latex].

Beweis

Angenommen [latex]x \geq x_0[/latex]. Dann gilt

[latex]
\begin{aligned}[]|R_{x_0,n}^{f}(x)| \leq \int _{x_0}^{x} |f^{(n+1)}(t)| \frac {(x-t)^{n}}{n!} \thinspace {\rm {d}} t \leq \frac {M_{n+1}}{n!} \int _{x_0}^{x} (x-t)^{n} \thinspace {\rm {d}} t\\ = \frac {M_{n+1}}{n!} \left [ - \frac {1}{n+1} (x-t)^{n+1} \right ]_{x_0}^x = \frac {M_{n+1} (x-x_0)^{n+1}}{(n+1)!}.\end{aligned}
[/latex]

Der Beweis für [latex]x

[latex]
\begin{aligned}[]|R_{x_0,n}^{f}(x)| = \left | \int _{x_0}^{x} f^{(n+1)}(t) \frac {(x-t)^{n}}{n!} \thinspace {\rm {d}} t \right | \leq \int _{x}^{x_0} |f^{(n+1)}(t)| \frac {(t-x)^{n}}{n!} \thinspace {\rm {d}} t.\end{aligned}
[/latex]

Bemerkung

Wir hätten andere Versionen des obigen Satz zu Taylor-Approximationen bereits in Abschnitt 7.2 für reellwertige Funktionen beweisen können. Dies sogar unter der etwas schwächeren Annahme an die [latex](n+1)[/latex]-te Ableitung, dass diese bloss zwischen [latex]x_0[/latex] und [latex]x[/latex] existieren soll und nicht unbedingt stetig sein muss. Unter dieser Voraussetzung gibt es ein [latex]\xi _C[/latex] zwischen [latex]x_0[/latex] und [latex]x[/latex], so dass das sogenannte Restglied nach Cauchy durch

[latex]
\begin{aligned}[]R_{x_0,n}^{f}(x) = \frac {1}{n!}f^{(n+1)}(\xi _C) (x-\xi _C)^n (x-x_0)\end{aligned}
[/latex]

gegeben ist. Es gibt unter denselben Voraussetzungen auch ein [latex]\xi _L[/latex] zwischen [latex]x_0[/latex] und [latex]x[/latex], so dass das Restglied nach Lagrange durch

[latex]
\begin{aligned}[]R_{x_0,n}^{f}(x) = \frac {1}{(n+1)!}f^{(n+1)}(\xi _L)(x-x_0)^{n+1}\end{aligned}
[/latex]

gegeben ist. Wir verweisen dafür auf folgende Übung.

Übung 8.60: Andere Formeln für das Restglied

Gegeben sei ein nicht-leeres Intervall [latex](a,b) \subseteq \mathbb {R}[/latex], Zahlen [latex]x,x_0 \in (a,b)[/latex] und eine [latex](n+1)[/latex]-mal differenzierbare Funktion [latex]f:(a,b) \to \mathbb {R}[/latex].

  1. Zeigen Sie, dass die Ableitung der Funktion
    [latex]
    \begin{aligned}[]F:t \in (a,b) \mapsto \sum _{k=0}^n \frac {f^{(k)}(t)}{k!} (x-t)^k,\end{aligned}
    [/latex]

    durch [latex]t \in (a,b) \mapsto \frac {f^{(n+1)}(t)}{n!}(x-t)^n[/latex] gegeben ist.

  2. Wenden Sie den Mittelwertsatz (Theorem 7.29) auf die obige Funktion [latex]F[/latex] an, um die Formel [latex]R_{x_0,n}^{f}(x) = \frac {1}{n!}f^{(n+1)}(\xi _C) (x-\xi _C)^n (x-x_0)[/latex] für das Restglied nach Cauchy zu beweisen.
  3. Verwenden Sie den verallgemeinerten Mittelwertsatz (Satz 7.48) für obiges [latex]F[/latex] und die Funktion [latex]g: t \mapsto (x-t)^{n+1}[/latex], um die Formel [latex]R_{x_0,n}^{f}(x) = \frac {1}{(n+1)!}f^{(n+1)}(\xi _L)(x-x_0)^{n+1}[/latex] für das Restglied nach Lagrange zu beweisen.

Beispiel 8.61: Extremwerte

Wir wollen die Taylor-Approximation verwenden um die Diskussion in Abschnitt 7.1.3 und Korollar 7.37 zu verfeinern. Sei [latex](a,b) \subseteq \mathbb {R}[/latex] ein offenes, nicht-leeres Intervall und sei [latex]f:(a,b) \to \mathbb {C}[/latex] eine [latex]n[/latex]-mal stetig differenzierbare Funktion. Angenommen [latex]x_0\in (a,b)[/latex] erfüllt

[latex]
\begin{aligned}[]f'(x_0)=\ldots =f^{(n-1)}(x_0)=0.\end{aligned}
[/latex]

Dann gelten folgende Implikationen.

  • Falls [latex]f^{(n)}(x_0)
  • Falls [latex]f^{(n)}(x_0)>0[/latex] ist und [latex]n[/latex] gerade ist, so nimmt [latex]f[/latex] in [latex]x_0[/latex] ein lokales Minimum an.
  • Falls [latex]f^{(n)}(x_0)\neq 0[/latex] und [latex]n[/latex] ungerade ist, so ist [latex]x_0[/latex] kein lokales Extremum von [latex]f[/latex].

Alle drei Aussagen folgen aus der Taylor-Approximation in Theorem 8.58, die in diesem Fall für [latex]x\in (a,b)[/latex] die Form

[latex]
\begin{aligned}[]f(x)=f(x_0)+\int _{x_0}^x f^{(n)}(t) \frac {(x-t)^{n-1}}{(n-1)!} \thinspace {\rm {d}} t\end{aligned}
[/latex]

annimmt. Falls [latex]f^{(n)}(x_0)>0[/latex] ist, existiert ein [latex]\delta >0[/latex] mit [latex]f^{(n)}(t)>0[/latex] für alle [latex]t\in (x_0-\delta ,x_0+\delta )[/latex]. Wir betrachten nun mehrere Fälle. Ist [latex]x\in (x_0,x_0+\delta )[/latex], so ist obiges Integral positiv und damit [latex]f(x)>f(x_0)[/latex]. Ist [latex]x\in (x_0-\delta ,x_0)[/latex] und [latex]n-1[/latex] gerade, so ist obiges Integral (auf Grund der umgekehrten Reihenfolge der Integrationsgrenzen) negativ, womit [latex]f(x)f(x_0)[/latex], womit [latex]f[/latex] ein lokales Minimum in [latex]x_0[/latex] annimmt. Für [latex]f^{(n)}(x_0)

8.5.1 – Analytische Funktionen

Applet 8.62: Taylor-Approximationen

Wir stellen einige Taylor-Approximationen bei verschiebbaren Fusspunkten zu bekannten Funktionen dar.

Betrachtet man eine glatte Funktion [latex]f[/latex] in Theorem 8.58, für die der Wert der Ableitung von [latex]f[/latex] nicht zu wild wird für hohe Ableitungen, dann geht das Restglied zu einem fest gewählten [latex]x[/latex] für [latex]n \to \infty[/latex] gegen Null und die Funktion [latex]f[/latex] wird durch ihre Taylorreihe beschrieben. Wie bereits gesehen, muss dies jedoch nicht unbedingt der Fall sein. Eine mögliche Abschätzung, die ausreichen würde, ist

[latex]
\begin{aligned}[]\max _{t \in [x_0,x]} |f^{(n+1)}(t)|\leq cA^n\end{aligned}
[/latex]

für alle [latex]n\in \mathbb {N}[/latex] und zwei Konstanten [latex]c,A\geq 1[/latex]. Eine Abschätzung dieser Art findet man beispielsweise für [latex]\exp ,\sin ,\cos[/latex] und Kombinationen dieser Funktionen.

Übung 8.63

Sei [latex][a,b] \subseteq \mathbb {R}[/latex] ein kompaktes Intervall mit Endpunkten [latex]a

  1. Sei [latex]f:[a,b] \to \mathbb {R}[/latex] ein Polynom oder eine der Funktionen [latex]\exp ,\sin ,\cos[/latex]. Zeigen Sie, dass Konstanten [latex]c,A \geq 1[/latex] mit
    [latex]
    \begin{aligned}[]\label{eq:fund-polygrowthofderivatives} \max _{t \in [a,b]} |f^{(n+1)}(t)|\leq cA^n\end{aligned}
    [/latex]
    existieren.
  2. Zeigen Sie, dass falls [latex]f,g:[a,b] \to \mathbb {R}[/latex] Funktionen mit der Eigenschaft (8.11) auch [latex]f+g[/latex] und [latex]f \cdot g[/latex] diese Eigenschaften besitzen.

Eine Funktion, die sich um jeden Punkt im Definitionsbereich als Potenzreihe darstellen lässt, nennen wir analytisch. Genauer sei [latex](a,b) \subseteq \mathbb {R}[/latex] ein offenes, nicht-leeres Intervall und [latex]f:(a,b) \to \mathbb {C}[/latex] eine glatte Funktion. Die Funktion [latex]f[/latex] heisst analytisch, falls die Taylorreihe von [latex]f[/latex] um jeden Punkt in [latex]x_0 \in (a,b)[/latex] positiven Konvergenzradius [latex]R[/latex] hat und auf [latex](x_0-R,x_0+R) \cap (a,b)[/latex] mit [latex]f[/latex] übereinstimmt.

Wir wollen ein weiteres Beispiel, wo sich eine glatte Funktion mit ihrer Taylorreihe vergleichen lässt, Interessierten als Übung überlassen (vergleichen Sie dies mit Übung 8.14).

Übung 8.64

Sei [latex]\alpha \in \mathbb {C}[/latex]. Zeigen Sie, dass

[latex]
\begin{aligned}[](1+x)^\alpha = \sum _{n=0}^\infty \binom {\alpha }{n} x^n\end{aligned}
[/latex]

für alle [latex]x \in (-1,1)[/latex], wobei für [latex]n \in \mathbb {N}_0[/latex]

[latex]
\begin{aligned}[]\binom {\alpha }{n} = \frac {1}{n!}\prod _{k=0}^{n-1} (\alpha -k) = \frac {\alpha (\alpha -1)\cdots (\alpha -n+1)}{n!}.\end{aligned}
[/latex]

Hinweis.

Sei [latex]f: x\in (-1,1) \mapsto (1+x)^\alpha[/latex]. Berechnen Sie zuerst [latex]f^{(n)}[/latex] für alle [latex]n \in \mathbb {N}[/latex] und geben Sie die Taylorreihe von [latex]f[/latex] um Null an. Zeigen Sie dann, dass das Restglied «klein» wird.

Bemerkung

Analytische Funktionen haben im Zusammenhang mit «holomorphen Funktionen» auf offenen Teilmengen von [latex]\mathbb {C}[/latex] nach [latex]\mathbb {C}[/latex] eine schöne, geschlossene Behandlung, die in der Vorlesung «Funktionentheorie» im zweiten Studienjahr thematisiert wird. Der komplexe Blickwinkel erklärt zum Beispiel, warum die analytische Funktion [latex]x\in \mathbb {R}\to \frac 1{1+x^2}[/latex] bei [latex]x_0=0[/latex] eine Taylorreihe mit Konvergenzradius [latex]1[/latex] besitzt, obwohl die Funktion auf ganz [latex]\mathbb {R}[/latex] definiert ist.

8.5.2 – Konvergenzgeschwindigkeit des Newton-Verfahrens*

Als Anwendung des Satzes von Taylor (Theorem 8.58) möchten wir in diesem Teilabschnitt das Newton-Verfahren zur approximativen Berechnung einer Nullstelle einer gegebenen Funktion diskutieren. Weitere Anwendung werden in den nächsten Abschnitten dieses Kapitels folgen.

Beispiel 8.65

Sei [latex][a,b] \subseteq \mathbb {R}[/latex] ein kompaktes Intervall mit Endpunkten [latex]a0[/latex] für alle [latex]x \in [a,b][/latex]. (Wir erinnern daran, dass dies im Allgemeinen nicht dasselbe wie strenge Monotonie von [latex]f[/latex] ist, aber strenge Monotonie impliziert.) Angenommen es gilt [latex]f(a) 0[/latex]. Dann gibt es auf Grund des Zwischenwertsatzes (Satz 3.59) ein [latex]z \in [a,b][/latex] mit [latex]f(z) = 0[/latex], das wegen der strengen Monotonie von [latex]f[/latex] eindeutig bestimmt ist. Beginnend mit einem gewählten Punkt [latex]x_0 \in [a,b][/latex] definieren wir rekursiv
[latex]
\begin{aligned}[]x_{n+1} = x_{n} - \frac {f(x_n)}{f'(x_n)}\end{aligned}
[/latex]
für alle [latex]n \in \mathbb {N}_0[/latex] (falls [latex]x_n \in [a,b][/latex]). Die Idee hinter diesem Verfahren möchten wir in folgendem Bild erklären.

image

Abbildung 8.2 – In einem ersten Schritt approximiert man [latex]f[/latex] mit der Tangenten [latex]g_0: t \mapsto f(x_0) + f'(x_0) (t-x_0)[/latex] bei [latex]x_0[/latex]. Von dieser berechnet man nun die eindeutige Nullstelle [latex]x_1 = x_{0} - \frac {f(x_0)}{f'(x_0)}[/latex]. In einem zweiten Schritt verwendet man nun diesen neuen Punkt [latex]x_1[/latex] und berechnet die eindeutige Nullstelle [latex]x_2[/latex] der Tangenten [latex]g_1[/latex] bei [latex]x_1[/latex]. Dies führt man iterativ so fort.

Falls die Folge [latex](x_n)_n[/latex] definiert ist (das heisst, [latex]x_n \in [a,b][/latex] für alle [latex]n \in \mathbb {N}[/latex]) und konvergiert, dann ist ihr Grenzwert eine Nullstelle von [latex]f[/latex] und somit unter unseren Annahmen gleich der Nullstelle [latex]z[/latex]. In der Tat gilt für [latex]z' = \lim _{n \to \infty }x_n[/latex]

[latex]
\begin{aligned}[]z' = \lim _{n \to \infty }x_n = \lim _{n \to \infty }x_{n+1} = \lim _{n \to \infty }\left (x_{n} - \frac {f(x_n)}{f'(x_n)}\right ) = z' - \frac {f(z')}{f'(z')}\end{aligned}
[/latex]

und somit [latex]f(z') = 0[/latex] und [latex]z'=z[/latex].

Das Konvergenzverhalten der Folge [latex](x_n)_n[/latex] ist im Allgemeinen «sehr chaotisch» . Unter etwas stärkeren Annahmen möchten wir nun aber zeigen, dass die Folge [latex](x_n)_n[/latex] konvergiert und die Konvergenzgeschwindigkeit untersuchen.

Wir nehmen zusätzlich an, dass [latex]f[/latex] zweimal stetig differenzierbar ist und definieren

[latex]
\begin{aligned}[]M &= \max \left \lbrace {|f''(t)|} \mid {t \in [a,b]}\right \rbrace \\ m &= \min \left \lbrace {|f'(t)|} \mid {t \in [a,b]}\right \rbrace .\end{aligned}
[/latex]

Nun wenden wir Korollar 8.59 an, um den Fehlerterm [latex]R_{x_0,1}^f[/latex] bei der Approximation von [latex]f[/latex] durch die Tangente um [latex]x_0[/latex] abzuschätzen. Bei der Nullstelle [latex]z[/latex] ergibt dies
[latex]
\begin{aligned}[]\label{eq:fund-newton1} |R_{x_0,1}^f(z)| \leq \frac {M |z-x_0|^2}{2}\end{aligned}
[/latex]
Des Weiteren gilt per Definition
[latex]
\begin{aligned}[]\label{eq:fund-newton2} 0 = f(z) = f(x_0) + f'(x_0) (z-x_0) + R_{x_0,1}^f(z).\end{aligned}
[/latex]
Wir dividieren Gleichung (8.14) durch [latex]f'(x_0)[/latex] und erhalten

[latex]
\begin{aligned}[]0 = z-\underset {=x_1}{\underbrace {\left (x_0-\frac {f(x_0)}{f'(x_0)}\right )}} + \frac {R_{x_0,1}^f(z)}{f'(x_0)}.\end{aligned}
[/latex]

Mit der Abschätzung (8.13) ergibt dies
[latex]
\begin{aligned}[]\label{eq:fund-newton3} |x_1-z| = \left | \frac {R_{x_0,1}^f(z)}{f'(x_0)}\right | \leq \tfrac {M}{2m}\, |x_0-z|^2\end{aligned}
[/latex]
Falls nun [latex]\varepsilon \in (0,\frac {m}{M})[/latex] klein genug ist, so dass [latex]B_\varepsilon (z) \subseteq [a,b][/latex] gilt, und falls [latex]x_0 \in B_\varepsilon (z)[/latex] ist, dann folgt aus (8.15), dass

[latex]
\begin{aligned}[]|x_1-z| \leq \tfrac {M}{2m}\, |x_0-z|^2 \leq \tfrac {M}{2m}\, |x_0-z|\varepsilon \leq \tfrac {1}{2} |x_0-z| [/latex]

Daher liegt [latex]x_1[/latex] auch in [latex][a,b][/latex] und ist bereits näher an [latex]z[/latex] als [latex]x_0[/latex]. Komplett analog beweist man, dass für alle [latex]n \in \mathbb {N}_0[/latex]

[latex]
\begin{aligned}[]|x_{n+1}-z| \leq \tfrac {1}{2}\, |x_n-z|\end{aligned}
[/latex]

gilt, was mit vollständiger Induktion auch

[latex]
\begin{aligned}[]|x_n-z| \leq 2^{-n} |x_0-z|\end{aligned}
[/latex]

für alle [latex]n \in \mathbb {N}_0[/latex] zeigt. Insbesondere konvergiert die Folge [latex](x_n)_n[/latex] und wir können das Newton-Verfahren (unter geeigneter Wahl eines Startpunktes [latex]x_0[/latex]) verwenden, um die Nullstelle mit hoher Genauigkeit zu approximieren. Die Konvergenz ist noch schneller als diese Abschätzung beweist; man spricht von quadratischer Konvergenz.

Übung 8.66: Quadratische Konvergenz

Wir möchten die Konvergenzgeschwindigkeit des Newton-Verfahrens hier genauer analysieren. In obigem Beispiel wurde gezeigt, dass sich der Fehler in jedem Schritt mindestens halbiert. Wir betrachten nun das Argument etwas genauer (und behalten dementsprechend die Notation). Sei [latex]\beta = \frac {M}{2m}[/latex]. Zeigen Sie für den mit [latex]\beta[/latex] gewichteten Abstand [latex]\beta |x_n-z|[/latex] die Abschätzung

[latex]
\begin{aligned}[]\beta |x_n-z| \leq (\beta |x_0-z|)^{2^{n}}\end{aligned}
[/latex]

für jedes [latex]n \in \mathbb {N}[/latex] gilt.

Übung 8.67

Verwenden Sie das Newton-Verfahren, um [latex]\sqrt {2}[/latex] (oder alternativ eine andere irrationale, reelle algebraische Zahl) auf drei Dezimalstellen genau zu berechnen. Für die genaue, numerische Berechnung des Resultats dürfen Sie einen Rechner oder ein Computer-Programm benutzen.

8.6 – Numerische Integration

Wie bereits erwähnt wurde, gibt es Integrale, die sich nicht in Ausdrücken der üblichen (den uns bisher bekannten) Funktionen darstellen lassen. Besonderes in diesen Fällen sind folgende Abschätzungen zur approximativen Berechnung von Integralen sehr nützlich. Im nächsten Abschnitt werden wir weitere Beispiele sehen, die die zugrundeliegende Idee des folgenden Satzes in anderen Überlegungen gewinnbringend einsetzen.

Satz 8.68

Seien [latex]a

  1. (Rechtecksregel*) Falls [latex]f[/latex] stetig differenzierbar ist, dann gilt
    [latex]
    \begin{aligned}[]\int _a^{b} f(x) \thinspace {\rm {d}} x = h \big ( f(x_0) + \ldots + f(x_{n-1}) \big ) + F_1,\end{aligned}
    [/latex]

    wobei der Fehler [latex]F_1[/latex] durch [latex]|F_1| \leq \frac {(b-a)^{2}}{2n} \max _{x \in [a,b]} |f'(x)|[/latex] beschränkt ist.

  2. (Sehnentrapezregel) Falls [latex]f[/latex] zweimal stetig differenzierbar ist, dann gilt
    [latex]
    \begin{aligned}[]\int _a^{b} f(x) \thinspace {\rm {d}} x &= h \left ( \frac {f(x_0)+f(x_1)}{2} + \ldots + \frac {f(x_{n-1})+f(x_n)}{2} \right ) + F_2 \\ &= \frac {h}{2} \big ( f(x_0) + 2f(x_1) + \ldots + 2f(x_{n-1}) + f(x_n) \big ) + F_2,\end{aligned}
    [/latex]

    wobei der Fehler [latex]F_2[/latex] durch [latex]|F_2| \leq \frac {(b-a)^{3}}{6n^{2}} \max _{x \in [a,b]} |f''(x)|[/latex] beschränkt ist.

  3. (Simpson-Regel*) Falls [latex]f[/latex] viermal stetig differenzierbar ist und [latex]n[/latex] gerade ist, dann gilt
    [latex]
    \begin{aligned}[]\int _a^{b} f(x) \thinspace {\rm {d}} x = \frac {h}{3} \big ( f(x_0) + 4 f(x_1) + 2 f(x_2) &+ 4 f(x_3) + 2 f(x_4) + \ldots \\ &+ 2 f(x_{n-2}) + 4f(x_{n-1}) + f(x_n)\big ) + F_3,\end{aligned}
    [/latex]

    wobei der Fehler [latex]F_3[/latex] durch [latex]|F_3| \leq \frac {(b-a)^{5}}{45 n^{4}} \max _{x \in [a,b]}|f^{(4)}(x)|[/latex] beschränkt ist.

*Wir haben zwar die Simpson-Regel ausführlicher besprochen. Doch in Hinblick auf eine mögliche Testfrage eignet sich die Sehnentrapezregel, wo bereits alle Ideen im Beweis auftauchen, besser.

Insbesondere verhält sich der Fehler für das Rechtecksverfahren wie [latex]O_{f}(n^{-1})[/latex] für [latex]n \to \infty[/latex], für das Sehnentrapezverfahren wie [latex]O_{f}(n^{-2})[/latex] für [latex]n \to \infty[/latex] und für das Simpson-Verfahren wie [latex]O_{f}(n^{-4})[/latex] für [latex]n \to \infty[/latex].

Wir möchten anmerken, dass die Konstanten in obigen Abschätzungen nicht optimal sind. Alle drei obigen Approximationsverfahren sind sogenannte Newton-Cotes-Verfahren. Die wesentliche Idee eines solchen Verfahrens ist die folgende: zuerst schreibt man nach Intervalladditivität

[latex]
\begin{aligned}[]\int _a^b f(x) = \sum _{\ell =0}^{n-1} \int _{x_\ell }^{x_{\ell +1}} f(x) \thinspace {\rm {d}} x.\end{aligned}
[/latex]

Nun approximiert man jedes obige Integralstück [latex]\int _{x_\ell }^{x_{\ell +1}} f(x) \thinspace {\rm {d}} x[/latex] durch einen Ausdruck der Form [latex]\sum _{k=1}^{K} w_k f(z_k)[/latex] für Gewichte [latex]w_1,\ldots ,w_K \in (0,1)[/latex] mit [latex]\sum _{k=1}^K w_k = 1[/latex] und Stützpunkte [latex]z_1,\ldots ,z_K \in [x_\ell ,x_{\ell +1}][/latex]. Beispielsweise nimmt man für die Sehnentrapezregel zwei Stützpunkte ([latex]K=2[/latex]), nämlich die beiden Endpunkte des Intervalls [latex][x_{\ell },x_{\ell +1}][/latex], mit den Gewichten [latex]w_1 = w_2 = \frac {1}{2}[/latex].

Die Summe der Fehler, die auf den Stücken [latex]\int _{x_\ell }^{x_{\ell +1}} f(x) \thinspace {\rm {d}} x[/latex] zustandekommen, ergeben dann den Gesamtfehler der Approximation. Obiger Satz gibt demnach an, wie dieser Fehler kontrolliert werden kann.

Vor dem Beweis des obigen Satzes möchten wir kurz erklären, wie sich die Simpson-Regel als Newton-Cotes-Verfahren auffassen lässt. Für die äquidistante Zerlegung des Intervalles [latex][a,b][/latex] mit den Punkten [latex]y_\ell = a + \ell \frac {b-a}{m}[/latex] für [latex]\ell = 0, \ldots ,m[/latex] betrachtet man auf [latex][y_\ell ,y_{\ell +1}][/latex] die Gewichte [latex]w_1 = \frac {1}{6}[/latex], [latex]w_2 = \frac {4}{6}[/latex] und [latex]w_3 = \frac {1}{6}[/latex] und die Stützpunkte [latex]y_\ell[/latex], [latex]\frac {y_\ell +y_{\ell +1}}{2}[/latex] und [latex]y_{\ell +1}[/latex]. Das dazugehörige Newton-Cotes-Verfahren ist genau das Simpson-Verfahren (wieso?). Wir wenden uns nun dem Beweis des obigen Satzes zu.

Beweis

Für (a) verwenden wir den Mittelwertsatz, wonach es zu [latex]\ell \in \left \lbrace {0,\ldots ,n-1} \right \rbrace[/latex] und [latex]x \in [x_\ell ,x_{\ell +1}][/latex] ein [latex]\xi _x \in (x_\ell ,x_{\ell +1})[/latex] gibt mit [latex]f(x) - f(x_\ell ) = f'(\xi _x) (x-x_\ell )[/latex]. Insbesondere gilt

[latex]
\begin{aligned}[]|f(x) - f(x_\ell )| = |f'(\xi _x)| (x-x_\ell ) \leq (x-x_\ell ) \max _{t \in [a,b]}|f'(t)|.\end{aligned}
[/latex]

Damit erhalten wir

[latex]
\begin{aligned}[]\left | \int _{x_\ell }^{x_{\ell +1}} f(x) \thinspace {\rm {d}} x - f(x_\ell ) h \right | &\leq \int _{x_\ell }^{x_{\ell +1}} |f(x)-f(x_\ell )| \thinspace {\rm {d}} x \leq \max _{t \in [a,b]}|f'(t)| \int _{x_\ell }^{x_{\ell +1}} (x-x_\ell ) \thinspace {\rm {d}} x \\ &\leq \frac {h^{2}}{2} \max _{t \in [a,b]}|f'(t)| .\end{aligned}
[/latex]

Durch Summation, Intervalladditivität des Integrals und die Dreiecksungleichung erhalten wir

[latex]
\begin{aligned}[]\left | \int _{a}^{b} f(x) \thinspace {\rm {d}} x - \sum _{\ell =0}^{n-1}f(x_\ell ) h \right | \leq \frac {nh^{2}}{2} \max _{t \in [a,b]}|f'(t)| = \frac {(b-a)^{2}}{2n} \max _{t \in [a,b]}|f'(t)| .\end{aligned}
[/latex]

Für (b) betrachten wir zuerst zu [latex]\ell \in \left \lbrace {0,\ldots ,n-1} \right \rbrace[/latex] die Endpunkte [latex]x_- = x_{\ell }[/latex] und [latex]x_+ = x_{\ell +1}[/latex] und den Mittelpunkt [latex]\tilde {x} = \frac {x_{-}+x_+}{2}[/latex] des Intervalls [latex][x_-,x_+][/latex]. Des Weiteren definieren wir den Wert [latex]M_2 = \max _{x \in [a,b]} |f''(x)|[/latex]. Nach Korollar 8.59 gilt für die Approximation durch das erste Taylor-Polynom um [latex]\tilde {x}[/latex]

[latex]
\begin{aligned}[]\left |f(t) - \big (f(\tilde {x}) + f'(\tilde {x}) (t-\tilde {x}) \big )\right | \leq \frac {M_2}{2!}\ |t-\tilde {x}|^{2} \leq \frac {M_2}{2} \left (\frac {h}{2}\right )^{2} \leq \frac {M_2}{8} h^{2}\end{aligned}
[/latex]

für alle [latex]t \in [x_-,x_+][/latex]. Wir verwenden dies für die Endpunkte [latex]t = x_-[/latex] und [latex]t = x_+[/latex] des Intervalls [latex][x_-,x_+][/latex] und erhalten aus der Dreiecksungleichung

[latex]
\begin{aligned}[]\bigg | \frac {f(x_-) + f(x_+)}{2} - f(\tilde {x}) \bigg | \\ =\tfrac {1}{2} \big | \big (f(x_{-})-f(\tilde {x}) - f'(\tilde {x}) (x_{-}-\tilde {x})\big ) + \big (f(x_+)-f(\tilde {x}) - f'(\tilde {x}) (x_{+}-\tilde {x})\big ) \big | \leq \frac {M_2 h^{2}}{8},\end{aligned}
[/latex]

da sich der lineare Term wegen [latex]x_- - \tilde {x} = -(x_+-\tilde {x})[/latex] aufhebt.

Aus demselben Grund erhalten wir

[latex]
\begin{aligned}[]\left | \int _{x_{-}}^{x_+} f(t) \thinspace {\rm {d}} t - f(\tilde {x})h \right | &=\left | \int _{x_{-}}^{x_+} \big (f(t) - f(\tilde {x})\big ) \thinspace {\rm {d}} t \right | =\left | \int _{x_{-}}^{x_+} \big (f(t)- (f(\tilde {x}) + f'(\tilde {x})(t-\tilde {x}))\big ) \thinspace {\rm {d}} t \right | \\ &\leq \tfrac {M_2}{2} \int _{x_{-}}^{x_+} (t-\tilde {x})^{2} \thinspace {\rm {d}} t = \tfrac {M_2}{2} \int _{-\frac {h}{2}}^{\frac {h}{2}} s^{2} \thinspace {\rm {d}} s = \frac {M_2h^{3}}{3 \cdot 8}\end{aligned}
[/latex]

Zusammenfassend gilt also für [latex]\ell \in \left \lbrace {0,\ldots ,n-1} \right \rbrace[/latex] und [latex]\tilde {x}_\ell = \frac {x_\ell + x_{\ell +1}}{2}[/latex]
[latex]
\begin{aligned}[]\bigg |f(\tilde {x}_\ell ) -\frac {f(x_{\ell -1}) + f(x_\ell )}{2} \bigg | &\leq \frac {M_2 h^{2}}{8}\label{eq:fund-bewtrapez1} \\ \left | \int _{x_{\ell -1}}^{x_\ell } f(t) \thinspace {\rm {d}} t - f(\tilde {x}_\ell )h \right | &\leq \frac {M_2h^{3}}{3 \cdot 8} \label{eq:fund-bewtrapez2}\end{aligned}
[/latex]
Wir multiplizieren (8.16) mit [latex]h[/latex] und summieren sowohl (8.16) als auch () über [latex]\ell[/latex] in [latex]\left \lbrace {1,\ldots , n} \right \rbrace[/latex]. Daraus folgt mit Intervalladditivität des Riemann-Integrals

[latex]
\begin{aligned}[]\left | \int _a^{b}f(x) \thinspace {\rm {d}} x -h \sum _{\ell = 1}^{n} \frac {f(x_{\ell -1})+f(x_\ell ) }{2} \right | \leq n \left ( \frac {M_2 h^{3}}{8} + \frac {M_2 h^{3}}{3 \cdot 8}\right ) = \frac {(b-a)^{3} M_2}{6 n^{2}}.\end{aligned}
[/latex]

Für (c) betrachten wir wieder zuerst zu [latex]k \in \left \lbrace {0,\ldots , \frac {n}{2}-1} \right \rbrace[/latex] die Endpunkte [latex]x_- = x_{2k}[/latex] und [latex]x_+ = x_{2k+2}[/latex] und den Mittelpunkt [latex]\tilde {x} = x_{2k+1}[/latex] des Intervalls [latex][x_-,x_+][/latex] und verwenden Korollar 8.59 bei [latex]\tilde {x}[/latex]. Dies ergibt für [latex]t \in [x_-,x_+][/latex]
[latex]
\begin{aligned}[]\label{eq:fund-bewsimpson1} \bigg | f(t) - \bigg (f(\tilde {x}) &+ f'(\tilde {x}) (t-\tilde {x}) + \frac {f''(\tilde {x})}{2} (t-\tilde {x})^{2} + \frac {f'''(\tilde {x})}{6} (t-\tilde {x})^3\bigg ) \bigg |\nonumber \\ &\leq \frac {M_4 (t-\tilde {x})^4}{4!} \leq \frac {M_4h^4}{4!},\end{aligned}
[/latex]
wobei [latex]M_4 = \max _{x \in [a,b]} |f^{(4)}(x)|[/latex]. Durch Integration über [latex][x_-,x_+][/latex] erhalten wir

[latex]
\begin{aligned}[]\left |\int _{x_{-}}^{x_{+}} f(t) \thinspace {\rm {d}} t - \left ( 2h f(\tilde {x}) + \frac {f''(\tilde {x})}{2} \int _{-h}^h s^2 \thinspace {\rm {d}} s \right ) \right | \leq \frac {M_4}{4!} \int _{-h}^h s^4 \thinspace {\rm {d}} s\end{aligned}
[/latex]

oder auch

[latex]
\begin{aligned}[]\left |\int _{x_{-}}^{x_{+}} f(t) \thinspace {\rm {d}} t - \left ( 2h f(\tilde {x}) + \frac {f''(\tilde {x})}{3 } h^3 \right ) \right | \leq \frac {M_4h^5}{60}\end{aligned}
[/latex]

Setzen wir [latex]t = x_{-}[/latex] und [latex]t = x_{+}[/latex] in Gleichung (8.17), so erhalten wir

[latex]
\begin{aligned}[]\frac {h}{3} \left |\big (f(x_{-}) + f(x_{+})\big ) - \left (2 f(\tilde {x}) + f''(\tilde {x})h^2 \right )\right | \leq \frac {h \cdot 2 \cdot M_4 h^4}{3\cdot 24} = \frac {M_4 h^5}{36},\end{aligned}
[/latex]

da sich die linearen und kubischen Terme gegenseitig aufheben. Daher gilt

[latex]
\begin{aligned}[]\bigg |\int _{x_{-}}^{x_{+}} f(t) \thinspace {\rm {d}} t - \tfrac {h}{3} \big (f(x_{-}) + 4 f(\tilde {x}) + f(x_{+}) \big ) \bigg |\\ = \left | \int _{x_{-}}^{x_{+}} f(t) \thinspace {\rm {d}} t - \tfrac {h}{3}\bigl (6 f(\tilde {x}) + f''(\tilde {x}) h^2 \bigr ) + \tfrac {h}{3} \big (2 f(\tilde {x}) + f''(\tilde {x})h^2\big ) -\tfrac {h}{3}\big ( f(x_{-}) + f(x_{+})\big ) \right |\\ \leq \tfrac 1{60}M_4h^5+ \tfrac 1{36}M_4 h^5= \tfrac 2{45}M_4h^5.\end{aligned}
[/latex]

Nach Summation über [latex]k \in \left \lbrace {0,\ldots , \frac {n}{2}-1} \right \rbrace[/latex] ergibt sich die Folge [latex]1,4,2,4,2,\ldots ,2,4,1[/latex] der Gewichte für die Funktionswerte in der Simpson-Regel wie im Satz und die Abschätzung genau wie im Beweis von (b) oben. ∎

Übung 8.69

Erklären Sie unter Verwendung der Simpson-Regel, wie man [latex]\pi = 4 \int _0^1 \frac {1}{1+x^2} \thinspace {\rm {d}} x[/latex] auf beliebig viele Dezimalstellen genau bestimmen kann.

8.6.1 – Landau-Notation II

Wie in obigem Beweis der Simpson-Regel ersichtlich wurde, ist das genaue Buchführen der Konstanten relativ anstrengend und der eigentliche Wert, den man zum Schluss erhält, steckt nicht so sehr in der konkreten Konstante sondern in den anderen Ausdrücken. Wir möchten nun deswegen ein weiteres Stück Notation einführen, welche uns erlaubt bei Abschätzungen unsere Denkleistung auf das Wesentliche zu fokussieren.

Sei [latex]X[/latex] eine Menge und seien [latex]f,g:X\to \mathbb {C}[/latex] zwei Funktion. Wir schreiben

[latex]
\begin{aligned}[]f(x) = O(g(x)) \quad \text {für } x \in X,\end{aligned}
[/latex]

falls eine Konstante [latex]C > 0[/latex] existiert mit [latex]|f(x)| \leq C |g(x)|[/latex] für alle [latex]x \in X[/latex].

Die obige Notation kann leicht mit der Landau-Notation aus Abschnitt 5.6 verwirrt werden, weswegen wir uns Mühe geben werden, den Zusatz «für [latex]x \in X[/latex]» auch immer anzugeben. In einem gewissen Sinne ist die Aussage [latex]f(x) = O(g(x)) \ \text {für } x \in [a,\infty )[/latex] eine «globale Aussage» , da alle Zahlen [latex]x[/latex] in [latex][a,\infty )[/latex] betroffen sind. Hingegen ist [latex]f(x) = O(g(x)) \ \text {für } x \to \infty[/latex] eine «lokale Aussage» , da nur Zahlen [latex]x\in [a,\infty )[/latex] betroffen sind, die gross genug sind (das heisst nahe genug bei unendlich sind, siehe Abschnitt 5.6).

In gewissen Situation bedeutet die Notation aber dasselbe, wie folgende Übung zeigt.

Übung 8.70

Sei [latex]a \in \mathbb {R}[/latex] und seien [latex]f,g:[a,\infty ) \to \mathbb {R}[/latex] zwei stetige Funktionen mit [latex]g(x)>0[/latex] für alle [latex]x\in [a,\infty )[/latex]. Zeigen Sie, dass die Aussagen

[latex]
\begin{aligned}[]f(x) &= O(g(x)) \quad \text {für } x \in [a,\infty )\\ f(x) &= O(g(x)) \quad \text {für } x \to \infty\end{aligned}
[/latex]

äquivalent sind.

Die Gross-O-Notation ist per Definition also insbesondere dann nützlich, wenn wir Konstanten «verstecken» wollen. Beispielsweise ist [latex]100! = O(1)[/latex]. Damit die Notation auch in Rechnungen nützlich ist, erlauben wir auch arithmetische Operationen mit ihr (vergleiche auch Abschnitt 5.6): Ist [latex]X[/latex] eine Menge und sind [latex]f_1,f_2,g:X\to \mathbb {C}[/latex] drei Funktionen, dann bedeutet

[latex]
\begin{aligned}[]f_1(x) = f_2(x) + O(g(x)) \quad \text {für } x \in X,\end{aligned}
[/latex]

dass [latex]f_1(x) - f_2(x) = O(g(x))[/latex] für [latex]x \in X[/latex] oder intuitiv ausgedrückt, dass die Differenz von [latex]f_1[/latex] und [latex]f_2[/latex] durch [latex]g[/latex] kontrolliert ist. Des Weiteren folgt aus [latex]f_1=O(g)[/latex] für [latex]x\in X[/latex], auch [latex]f_1f_2=O(gf_2)[/latex] für [latex]x\in X[/latex].

Die Gross-O-Notation wird auch dazu verwendet, Unabhängigkeiten von gewissen Parametern zum Ausdruck zu bringen. Wir möchten dies an einem Beispiel illustrieren.

Beispiel 8.71: Parameterabhängigkeit bei Landau-Notation

Sei [latex]k \in \mathbb {N}[/latex] und [latex]\alpha \in (-\infty ,2][/latex].

  1. Nach Taylorapproximation im Sinne von Korollar 8.59 gilt
    [latex]
    \begin{aligned}[]\left |x^{\frac {3}{2}} - \left ( k^{\frac {3}{2}} + \tfrac {3}{2} k^{\frac {1}{2}}(x-k) \right )\right | \leq \tfrac {1}{2}\max _{t\in [1,\infty )}\bigl |\tfrac 34t^{-\frac 12}\bigr | (x-k)^2 = \tfrac {3}{8} (x-k)^2\end{aligned}
    [/latex]

    für [latex]x\in [1,\infty )[/latex] . Insbesondere gilt

    [latex]
    \begin{aligned}[]x^{\frac {3}{2}} = k^{\frac {3}{2}} + \tfrac {3}{2} k^{\frac {1}{2}} (x-k) + O((x-k)^2)\quad \text {für } x \in [1,\infty ),\end{aligned}
    [/latex]

    wobei die in [latex]O(\cdot )[/latex] versteckte implizite Konstante also nicht vom Parameter [latex]k[/latex] abhängt.

  2. Hängt die implizite Konstante, nicht wie in (a) oben, vom Parameter ab, so indizieren wir den Parameter bei [latex]O(\cdot )[/latex]. Ein konkretes Beispiel: Wenn [latex]\alpha \in (-\infty ,2][/latex], dann gilt
    [latex]
    \begin{aligned}[]x^\alpha = 1+ \alpha (x-1) + O_\alpha ((x-1)^2)\end{aligned}
    [/latex]

    für [latex]x \in [\frac {1}{2},\infty )[/latex] auf Grund der Taylorapproximation in Korollar 8.59.

  3. Auch eine Abhängigkeit des Definitionsbereichs vom Parameter ist denkbar und oft nützlich. Für alle Zahlen [latex]t \in \left [k-\frac 12,k+\frac 12\right ][/latex] gilt [latex](\log )''(t)=-\frac 1{t^{2}}=O\left (\frac 1{k^2}\right )[/latex] und daher
    [latex]
    \begin{aligned}[]\log (t) = \log (k) + \tfrac {1}{k} (t-k) + O\left (\tfrac {(t-k)^2}{k^2} \right ) = \log (k) + \tfrac {1}{k} (t-k) + O\left (\tfrac {1}{k^2} \right )\end{aligned}
    [/latex]

    auf Grund der Taylorapproximation in Korollar 8.59.

Im Allgemeinen sollte man die Landau-Notation sorgfältig verwenden, da sich hier oft Fehler einschleichen insbesondere bei der Frage der Abhängigkeit der impliziten Konstanten.

8.7 – Drei asymptotische Formeln

8.7.1 – Das Wallissche Produkt

Die Formel
[latex]
\begin{aligned}[]\label{eq:fund-wallisprod} \pi = \lim _{n \to \infty } \frac {1}{n} \frac {(2^n n!)^4}{((2n)!)^2}\end{aligned}
[/latex]
nach John Wallis (1616-1703) wird das Wallissche Produkt genannt.

Um Gleichung (8.18) zu beweisen, wollen wir zuerst für jedes [latex]k \in \mathbb {N}_0[/latex] das Integral

[latex]
\begin{aligned}[]I_k = \int _0^{{\pi }} \sin ^k(x) \thinspace {\rm {d}} x\end{aligned}
[/latex]

mit vollständiger Induktion berechnen.

Für [latex]k = 0[/latex] gilt [latex]I_0 = {\pi }[/latex] und für [latex]k=1[/latex] gilt [latex]I_1 = \int _0^{{\pi }} \sin (x) \thinspace {\rm {d}} x = 2[/latex]. Mit partieller Integration erhalten wir für alle [latex]k\geq 2[/latex]

[latex]
\begin{aligned}[]I_{k} &= \int _0^{\pi } \sin ^{k}(x) \thinspace {\rm {d}} x = \int _0^{\pi } \sin ^{k-1}(x) \sin (x) \thinspace {\rm {d}} x\\ &= \left [\sin ^{k-1}(x) \big (-\cos (x)\big )\right ]_{0}^{\pi } + \int _0^{\pi } (k-1) \sin ^{k-2}(x) \cos ^2(x) \thinspace {\rm {d}} x\\ &= (k-1) \int _0^{\pi } \sin ^{k-2}(x) (1-\sin ^2(x)) \thinspace {\rm {d}} x\\ &= (k-1)(I_{k-2}-I_{k}),\end{aligned}
[/latex]

woraus sich die Rekursionsformel
[latex]
\begin{aligned}[]\label{eq:wallisrecursion} I_{k} = \frac {k-1}{k} I_{k-2}\end{aligned}
[/latex]
für jedes ganze [latex]k\geq 2[/latex] ergibt.

Somit ist für [latex]n \in \mathbb {N}[/latex]

[latex]
\begin{aligned}[]I_{2n} = \frac {2n-1}{2n}I_{2n-2} &= \frac {2n-1}{2n}\cdot \frac {2n-3}{2n-2}\cdots \frac {1}{2} I_0\\ &= \frac {(2n)(2n-1)(2n-2)\cdots 1}{\big ((2n)(2n-2)\cdots 2\big )^2} I_0 = \frac {(2n)!}{(2^n n!)^2}\cdot \pi ,\end{aligned}
[/latex]

wobei wir den Bruch mit [latex]2n(2n-2)\cdots 2[/latex] erweitert haben um im Zähler [latex](2n)![/latex] zu erhalten. Auf dieselbe Weise ergibt sich

[latex]
\begin{aligned}[]I_{2n+1} = \frac {2n}{2n+1}\cdot \frac {2n-2}{2n-1} \cdots \frac {2}{3} I_1 = \frac {\big ((2n)(2n-2)\cdots 2\big )^2}{(2n+1)(2n)\cdots 1} \cdot 2= \frac {(2^n n!)^2}{(2n+1)!}\cdot 2.\end{aligned}
[/latex]

Die Folge [latex](I_n)_n[/latex] erfüllt noch eine weitere Eigenschaft: Denn für [latex]x \in \left [0,\pi \right ][/latex] ist [latex]\sin (x) \in [0,1][/latex] und somit

[latex]
\begin{aligned}[]\sin ^{2n+2}(x) \leq \sin ^{2n+1}(x) \leq \sin ^{2n}(x)\end{aligned}
[/latex]

für alle [latex]x\in [0,\pi ][/latex] and [latex]n \in \mathbb {N}[/latex]. Durch Integration folgt die Monotonieungleichung

[latex]
\begin{aligned}[]I_{2n+2} = \int _0^{\pi } \sin ^{2n+2}(x) \thinspace {\rm {d}} x \leq I_{2n+1} = \int _0^{\pi } \sin ^{2n+1}(x) \thinspace {\rm {d}} x \leq I_{2n} = \int _0^{\pi } \sin ^{2n}(x) \thinspace {\rm {d}} x\end{aligned}
[/latex]

und mittels (8.19) auch

[latex]
\begin{aligned}[]I_{2n+2}=\frac {2n+1}{2n+2}\underbrace {\frac {(2n)!}{(2^n n!)^2} \pi }_{I_{2n}} \leq \frac {(2^n n!)^2}{(2n+1)!} 2 \leq I_{2n} =\frac {(2n)!}{(2^n n!)^2} \pi .\end{aligned}
[/latex]

Wir multiplizieren diese Ungleichungskette mit [latex]\frac {(2^n n!)^2}{(2n)!}\frac {2n+1}{2n}[/latex] und erhalten

[latex]
\begin{aligned}[]\frac {(2n+1)^2}{2n(2n+2)}\pi \leq \frac {1}{n} \frac {(2^n n!)^4}{((2n)!)^2} \leq \frac {2n+1}{2n} \pi .\end{aligned}
[/latex]

Dies beweist die Formel (8.18) nach Wallis unter Anwendung des Sandwich-Lemmas (Lemma 5.31).

8.7.2 – Stirling-Formel

Wir verwenden nun das Wallissche Produkt und Taylor-Approximationen, um das asymptotische Verhalten von [latex]n![/latex] zu bestimmen. Genauer formuliert wollen wir die Stirling-Formel
[latex]
\begin{aligned}[]\label{eq:fund-stirling} \lim _{n \to \infty } \frac {n!}{\sqrt {2\pi n} \left (\frac {n}{\mathrm {e}}\right )^n} = 1\end{aligned}
[/latex]
beweisen. Diese ist nach James Stirling (1692-1770) benannt und wird auch als [latex]n! \approx \sqrt {2\pi n} \left (\frac {n}{\mathrm {e}}\right )^n[/latex] geschrieben.

Der Grundgedanke für den Beweis von (8.20) ist die Gleichung

[latex]
\begin{aligned}[]\log (n!) = \log (1) + \log (2) + \ldots + \log (n)\end{aligned}
[/latex]

für [latex]n \in \mathbb {N}[/latex] genauer zu betrachten, wobei wir die Summe aber durch das Integral [latex]\int _{\frac {1}{2}}^{n + \frac {1}{2}} \log (x) \thinspace {\rm {d}} x[/latex] ersetzen wollen, da dieses einfacher zu berechnen ist.

Für [latex]k \in \mathbb {N}[/latex] und [latex]t \in \left [k-\frac 12,k+\frac 12\right ][/latex] gilt nach Beispiel 8.71

[latex]
\begin{aligned}[]\log (t) = \log (k) + \tfrac {1}{k} (t-k) + O\left (\tfrac {1}{k^2} \right ).\end{aligned}
[/latex]

Wir erhalten

[latex]
\begin{aligned}[]\int _{k-\frac 12}^{k+\frac 12} \log (t) \thinspace {\rm {d}} t = \log (k) + O\left (\tfrac {1}{k^2}\right )\end{aligned}
[/latex]

für alle [latex]k \in \mathbb {N}[/latex], da sich der lineare Term bei der Integration weglöscht.

Für ein festes [latex]n\in \mathbb {N}[/latex] und die Summe über [latex]k \in \left \lbrace {1,\ldots ,n} \right \rbrace[/latex] folgt daraus

[latex]
\begin{aligned}[]\log (n!) &= \sum _{k=1}^n \log (k) = \sum _{k=1}^n \int _{k-\frac 12}^{k+\frac 12} \log (t) \thinspace {\rm {d}} t + \sum _{k=1}^n O\left (\tfrac {1}{k^2}\right ) = \int _{\frac 12}^{n+\frac 12} \log (t) \thinspace {\rm {d}} t + \sum _{k=1}^n O\left (\tfrac {1}{k^2}\right )\\ &= (n+\tfrac {1}{2}) \log (n+\tfrac {1}{2}) - (n+\tfrac {1}{2}) - c + \sum _{k=1}^n O\left (\tfrac {1}{k^2}\right ),\end{aligned}
[/latex]

wobei sich die Konstante [latex]c\in \mathbb {R}[/latex] durch Auswerten der Stammfunktion bei der unteren Grenze [latex]\frac 12[/latex] ergibt. Wiederum nach Taylorapproximation gilt [latex]\log (n+\tfrac {1}{2}) = \log (n) + \tfrac {1}{n}\cdot \tfrac {1}{2} + O(\tfrac {1}{n^2})[/latex]. Wir multiplizieren diese Approximation mit [latex](n+\frac 12)[/latex] und erhalten (mit [latex](n+\frac 12)O(\frac 1{n^2})=O(\frac 1n)[/latex]) somit

[latex]
\begin{aligned}[]\log (n!) &= (n+\tfrac {1}{2}) \left (\log (n) + \tfrac {1}{2n}+ O(\tfrac {1}{n^2}) \right )- (n+\tfrac {1}{2}) -c+ \sum _{k=1}^nO\left (\tfrac {1}{k^2}\right )\\ &= (n+\tfrac {1}{2})\log (n) + \tfrac {1}{2}+ O\left (\tfrac {1}{n}\right ) -n -\tfrac 12 -c+ \sum _{k=1}^n O\left (\tfrac {1}{k^2}\right )\\ &= \frac 12\log (n)+n\log (n)-n+ a_n\end{aligned}
[/latex]

für eine konvergente Folge [latex](a_n)_n[/latex] in [latex]\mathbb {R}[/latex] (wobei wir Proposition 6.28 verwendet haben um die absolut konvergenten Fehlerterme [latex]\varepsilon _n=O(\frac 1{n^2})[/latex] aufzusummieren). Wenden wir die Exponential­abbildung an, so erhalten wir für alle [latex]n \in \mathbb {N}[/latex]
[latex]
\begin{aligned}[]\label{eq:fund-proofsterling1} n! = \sqrt {n}n^n\mathrm {e}^{-n} b_n,\end{aligned}
[/latex]
wobei [latex]b_n = \exp (a_n)[/latex] ist. Da [latex](a_n)_n[/latex] konvergent ist, ist auch [latex](b_n)_n[/latex] konvergent und der Grenzwert [latex]B = \lim _{n \to \infty } \exp (a_n) = \exp ( \lim _{n \to \infty } a_n)[/latex] ist positiv.

Um [latex]B[/latex] zu berechnen, verwenden wir (8.21) gemeinsam mit dem Wallisschen Produkt (8.18). Es gilt

[latex]
\begin{aligned}[]\pi = \lim _{n \to \infty } \frac {1}{n} \frac {(2^n n!)^4}{((2n)!)^2} = \lim _{n \to \infty } \frac {1}{n} \frac { \bigl (2^n \sqrt {n}n^n\mathrm {e}^{-n} b_n \bigr )^4}{\bigl (\sqrt {2n}(2n)^{2n}e^{-2n}b_{2n}\bigr )^2} = \lim _{n \to \infty } \frac {b_n^4 }{2b_{2n}^2},\end{aligned}
[/latex]

da die Potenzen von [latex]2[/latex], von [latex]n[/latex] und von [latex]e[/latex] sich jeweils direkt kürzen, und die Wurzelterme gemeinsam mit dem [latex]\frac 1n[/latex] vor dem Bruch sich zu einem [latex]\frac 12[/latex] vereinfachen. Daraus folgt

[latex]
\begin{aligned}[]\pi = \frac {B^4}{2B^2} = \tfrac {1}{2}B^2\end{aligned}
[/latex]

und somit muss [latex]B = \sqrt {2\pi }[/latex] gelten. Aus Gleichung (8.21) erhalten wir

[latex]
\begin{aligned}[]\lim _{n \to \infty } \frac {n!}{\sqrt { n} \left (\frac {n}{\mathrm {e}}\right )^n} = \lim _{n \to \infty }b_n = \sqrt {2\pi },\end{aligned}
[/latex]

was zur Stirling-Formel (8.20) äquivalent ist.

Auf Grund der Bedeutung von der Fakultätsfunktion, zum Beispiel in kombinatorischen Überlegungen, ist es nicht verwunderlich, dass die Sterling-Formel vielfache Anwendungen besitzt. In den folgenden Übungen wollen wir einige wenige solcher Anwendungen präsentieren.

Übung 8.72: Besuche bei der Irrfahrt auf [latex]\mathbb {Z}[/latex]

Wir betrachten eine (diskrete) zufällige Bewegung auf den ganzen Zahlen ausgehend von Null. Also beginnen wir zur Zeit [latex]t = 0[/latex] im Ursprung, hüpfen dann mit Wahrscheinlichkeit [latex]\tfrac {1}{2}[/latex] entweder nach links auf die Position [latex]-1[/latex] oder nach rechts auf die Position [latex]1[/latex]. Bezeichnet im Allgemeinen [latex]x_t \in \mathbb {Z}[/latex] die Position zur Zeit [latex]t\in \mathbb {N}[/latex], so bewegen wir uns auf die Zeit [latex]t+1[/latex] hin wieder entweder links oder rechts mit gleicher Wahrscheinlichkeit von [latex]\frac {1}{2}[/latex].

Berechnen Sie die Wahrscheinlichkeit, dass man sich zur Zeit [latex]t \in \mathbb {N}[/latex] bei Null befindet und approximieren Sie diese mit der Stirling-Formel.

Übung 8.73

Eine Münze wird [latex]n[/latex]-mal geworfen. Sei [latex]\varepsilon >0[/latex]. Zeigen Sie, dass die Wahrscheinlichkeit, dass die Münze höchstens [latex](\frac 12-\varepsilon )n[/latex]-mal auf Kopf fällt, für [latex]n \to \infty[/latex] gegen Null geht.

Hinweis: Sie dürfen dabei annehmen, dass [latex]2^{-n}\binom {n}{k}[/latex] die Wahrscheinlichkeit ist, dass bei [latex]n[/latex]-maligem Wurf die Münze [latex]k[/latex]-mal auf Kopf fällt. Zeigen Sie auch, dass die Funktion [latex]x\in (0,1)\mapsto x\log (x)+(1-x)\log (1-x)[/latex] ein striktes lokales Maximum bei [latex]x=\frac 12[/latex] hat.

Übung 8.74: Ein Ausblick auf Algorithmik

In dieser Übung wollen wir die worst case Geschwindigkeit eines Sortieralgorithmus abschätzen. Sei also gegeben ein Algorithmus, der zu einem Datensatz [latex]a_1,\ldots ,a_n \in \mathbb {Z}[/latex] diesen umordnet, das heisst, eine Abbildung [latex]j: \left \lbrace {1,\ldots ,n} \right \rbrace \to \left \lbrace {1,\ldots ,n} \right \rbrace[/latex] findet, so dass [latex]a_{j(1)}\leq \ldots \leq a_{j(n)}[/latex]. Wir interessieren uns nun für die Zahl [latex]W(n)[/latex] der Anzahl Vergleiche zweier Zahlen des Datensatzes, die der Algorithmus schlimmstenfalls durchführen muss, um den Datensatz zu sortieren. Zeigen Sie, dass eine von [latex]n[/latex] unabhängige Konstante [latex]C> 0[/latex] existiert, so dass [latex]W(n) \geq C n \log (n)[/latex].

Hinweis.

Zählen Sie die Anzahl Datensätze, die mit [latex]A[/latex] Umsortierungen aus einem fixen Datensatz konstruiert werden können.

8.7.3 – Asymptotik der harmonischen Reihe*

Für den Beweis der Stirling-Formel wurde eine Kombination von Taylor-Approximation und Integration verwendet, um das asymptotische Verhalten der Reihe [latex]\log (n!) = \sum _{k=1}^n \log (k)[/latex] zu bestimmen. Diese Methode kann auch allgemeiner zur Untersuchung des Divergenzverhaltens von Reihen eingesetzt werden. Zum Beispiel gilt die «asymptotische Entwicklung»
[latex]
\begin{aligned}[]\label{eq:fund-asym sum 1/n} \sum _{k=1}^n \frac {1}{k} = \log (n) + \gamma + \frac {1}{2n} + O(\tfrac {1}{n^2})\end{aligned}
[/latex]
für [latex]n \to \infty[/latex], wobei

[latex]
\begin{aligned}[]\gamma = \int _1^\infty \left (\frac {1}{\lfloor {x} \rfloor }-\frac {1}{x}\right ) \thinspace {\rm {d}} x = 0.577215664901532\ldots\end{aligned}
[/latex]

die Euler-Mascheroni-Konstante genannt wird.

Übung 8.75

Wir wollen die Asymptotik (8.22) in dieser Übung beweisen.

  1. Beweisen Sie, dass das Integral in der Definition von [latex]\gamma[/latex] konvergiert und dass die Gleichung
    [latex]
    \begin{aligned}[]\sum _{k=1}^n \frac {1}{k} = \log (n) + \gamma + \frac {1}{n} - \int _{n}^\infty \left (\frac {1}{\lfloor {x} \rfloor }-\frac {1}{x}\right ) \thinspace {\rm {d}} x\end{aligned}
    [/latex]

    für alle [latex]n \in \mathbb {N}[/latex] gilt.

  2. Verwenden Sie Taylor-Approximation, um
    • [latex]\log (m+1) = \log (m) + \tfrac {1}{m}- \tfrac {1}{2m^2} + O(\tfrac {1}{m^3})[/latex]
    • [latex]\frac {1}{m} - \int _m^{m+1} \frac {1}{t} \thinspace {\rm {d}} t = \frac {1}{2m^2} + O(\tfrac {1}{m^3}) = \tfrac {1}{2} \int _m^{m+1} t^{-2} \thinspace {\rm {d}} t + O(\tfrac {1}{m^3})[/latex]

    für [latex]m \in \mathbb {N}[/latex] zu beweisen.

  3. Verwenden Sie das Integralkriterium (Satz 8.46) um zu zeigen, dass
    [latex]
    \begin{aligned}[]\sum _{m=n}^\infty \frac {1}{m^3} = O(\tfrac {1}{n^2})\end{aligned}
    [/latex]

    für [latex]n \to \infty[/latex].

  4. Verbinden Sie die obigen Aussagen zu einem Beweis von (8.22).
  5. Wiederholen Sie Beispiel 6.5 und approximieren Sie mit der Formel (8.22) die notwendige Anzahl Bausteine [latex]n[/latex] für einen Sprungturm, der [latex]10[/latex]m in den See hinreicht. Vergleichen Sie dies dann zum exakten Resultat [latex]n=12367[/latex].

8.8 – Weitere Lernmaterialien

8.8.1 – Verwendung des Kapitels

Dieses Kapitel verbindet mit Hilfe des Fundamentalsatz der Analysis die zentralen Begriffe der Ableitung und des Riemann-Integrals. Damit haben wir das vollständige Arsenal an Ableitungsregeln auch für die Berechnung von Integralen verwenden können, womit die Berechnung von Integralen mitunter deutlich einfacher wurde — auch wenn dies zusätzliche Übung erfordert. Umgekehrt haben wir aber auch mit dem Integralrestglied im Satz zur Taylor-Approximation gesehen, dass das Riemann-Integral nützlich sein kann um den Zusammenhang zwischen den Ableitungen und der ursprünglichen Funktion besser zu verstehen. Falls die symbolische Integration sich als nicht machbar erweist, so ist wiederum die Taylor-Approximation nützlich um das Riemann-Integral numerisch — zum Beispiel mit der Simpson-Methode — mit überraschend hoher Genauigkeit zu berechnen. Ebenso ist aber der Satz zur Taylor-Approximation auch von theoretischer Wichtigkeit, da wir mit diesen asymptotische Formeln wie zum Beispiel die Sterling-Formel beweisen können.

Zusammenfassend können wir also sagen, dass dieses Kapitel den Aufbau der eindimensionalen Analysis vollendet. Die Inhalte dieses Kapitels bilden einen zentralen Bestandteil der Analysis I/II Vorlesung und ihrer Anwendungen.

8.8.2 – Übungen

Da die Schwierigkeit beim Integrieren vor allem in der Auswahl der richtigen Methode liegt, wollen wir in folgender Übung noch einige weitere Aufgaben ohne Angabe der richtigen Technik auflisten. Dazu wollen wir noch erwähnen, dass es oft auch mehr als eine Methode gibt, die zum Erfolg führen kann.

Übung

Berechnen Sie die unbestimmten Integrale

[latex]
\begin{aligned}[]\begin{array}{ll} \int \frac {1}{\sqrt {x^2+2017^2}} \thinspace {\rm {d}} x, &\int \frac {1}{x \sqrt {1+x^2}} \thinspace {\rm {d}} x,\\ \int \exp (\sqrt [3]x)\thinspace {\rm {d}} x, &\int \frac {1}{1+\exp (-x)}\thinspace {\rm {d}} x,\\ \int \frac {x}{1+\sqrt {x}}\thinspace {\rm {d}} x, &\int \cos ^5(x) \sin ^3(x) \thinspace {\rm {d}} x, \\ \int \cos ^2(x) \sin ^2(x) \thinspace {\rm {d}} x, &\int \frac {1}{\cos ^2(x) \sin ^2(x)} \thinspace {\rm {d}} x, \\ \int \frac {\log (x)+\sin (\log (x))x^2}{x^3} \thinspace {\rm {d}} x, & \int \frac {1}{\cos (x)}\thinspace {\rm {d}} x\end{array}\end{aligned}
[/latex]

Übung

Zeigen Sie, dass jede stetig differenzierbare Funktion [latex]f: [a,b] \to \mathbb {R}[/latex] auf einem kompakten Intervall [latex][a,b][/latex] mit Endpunkten [latex]a 4.8.2) hat und dass [latex]V(f)(x) = \int _a^x |f'(t)| \thinspace {\rm {d}} t[/latex] gilt für alle [latex]x \in [a,b][/latex].

Hinweis.

Gegeben eine Zerlegung [latex]\mathfrak {Z} = \left \lbrace {a = x_0

Übung: Abel-Summation und Partielle Integration

Wir wollen in dieser Übung den Zusammenhang zwischen der Abel-Summation und der partiellen Integration erklären. Sei also [latex]a3.3, um die partielle Integration [latex]\int _a^b uv'\thinspace {\rm {d}} x=[uv]_a^b-\int _a^bu'v\thinspace {\rm {d}} x[/latex] zu beweisen.

Hinweis.

Erklären Sie zuerst, wieso Sie ohne Beschränkung der Allgemeinheit [latex]u(b) = v(a) = 0[/latex] annehmen können. Sei [latex]\mathfrak {Z} = \left \lbrace {a= x_0

[latex]
\begin{aligned}[]\sum _{k=1}^nu(x_k)(v(x_k)-v(x_{k-1}))\end{aligned}
[/latex]

mit [latex]a_k=(v(x_k)-v(x_{k-1}))[/latex] und [latex]b_k=u(x_k)[/latex]. Verwenden Sie anschliessend den Mittelwertsatz auf [latex]v(x_k)-v(x_{k-1})[/latex] an und vergleichen Sie die resultierende Summen mit der Riemann-Summen für das Integral [latex]\int _a^buv'\thinspace {\rm {d}} x[/latex].

Übung: Ein alternativer Beweis von Proposition 4.30

Wir möchten hier einen Beweis von Proposition 4.30 unter der Annahme, dass [latex]f[/latex] stetig ist, durchführen. Gehen Sie wie folgt vor:

  1. Zeigen Sie, dass die Abbildung [latex]F:x\in [a,b] \mapsto \mathcal {I}(a,x)[/latex] differenzierbar ist und dass [latex]F' = f[/latex].
  2. Verwenden Sie nun den Fundamentalsatz der Integral- und Differentialrechnung, um auf Proposition 4.30 zu schliessen.

Übung: Volumen einer Vase

Berechnen Sie das Volumen der Vase

[latex]
\begin{aligned}[]\left \lbrace {(x,y,z)\in \mathbb {R}^3} \mid {x \in [-\pi ,2\pi ],\ 0 \leq \sqrt {y^2+z^2} \leq \sin (x)+2}\right \rbrace .\end{aligned}
[/latex]

Übung: Gleichmässige Konvergenz der Ableitungen

Sei [latex](f_n)_n[/latex] eine Folge stetig differenzierbarer reellwertiger Funktionen auf einem kompakten Intervall [latex][a,b][/latex] mit Endpunkten [latex]a

Übung

Wir verwenden obige Übung, um eine glatte Funktion zu konstruieren, deren Taylorreihe um einen Punkt Konvergenzradius Null hat.

  1. Zeigen Sie, dass die Reihe [latex]\sum _{n=0}^\infty \mathrm {e}^{-n} \cos (n^2 x)[/latex] für alle [latex]x\in \mathbb {R}[/latex] absolut konvergiert.
  2. Nach (i) definieren wir die Funktion
    [latex]
    \begin{aligned}[]f: x \in \mathbb {R} \mapsto \sum _{n=0}^\infty \mathrm {e}^{-n} \cos (n^2 x).\end{aligned}
    [/latex]

    Zeigen Sie, dass [latex]f[/latex] glatt ist.

  3. Berechnen Sie die Taylorreihe von [latex]f[/latex] um Null und deren Konvergenzradius.

Übung: Krümmung ebener Kurven

In dieser Übung möchten wir die Krümmung ebener Kurven betrachten. Alle Kurven, die wir dabei betrachten wollen, sollen zweimal differenzierbar, regulär und einfach sein, hier der Einfachheit vorerst inklusive den Endpunkten. Für eine solche Kurve [latex]\gamma :[a,b]\to \mathbb {R}^2[/latex] und ein [latex]p\in \mathbb {R}^2[/latex] auf der Kurve sei [latex]t\in [a,b][/latex] der eindeutige Zeitpunkt mit [latex]\gamma (t) = p[/latex]. Dann ist die Krümmung von [latex]\gamma[/latex] bei [latex]p[/latex] definiert als

[latex]
\begin{aligned}[]\kappa _\gamma (p) = \frac {\langle \dot {\gamma }(t),R \ddot {\gamma }(t) \rangle }{\| {\dot {\gamma }(t)}\| ^3},\end{aligned}
[/latex]

wobei [latex]R = \Big ({\scriptsize \arraycolsep =0.3\arraycolsep \ensuremath {\begin{matrix}0 & 1 \\ -1 & 0\end{matrix}}} \Big )[/latex]. Ist [latex]\gamma[/latex] einfach, aber erfüllt [latex]\gamma (a) = \gamma (b)[/latex], so reicht es anzunehmen, dass [latex]\dot {\gamma }(a) = \dot {\gamma }(b)[/latex] und [latex]\ddot {\gamma }(a) = \ddot {\gamma }(b)[/latex] gilt, damit obiger Ausdruck für die Krümmung Sinn ergibt.

  1. Berechnen Sie die Krümmung der Kurven [latex]t \mapsto (\cos (t),\sin (t))[/latex] und [latex]t \mapsto (\cosh (t),\sinh (t))[/latex] (definiert auf geeigneten Intervallen).
  2. Zeigen Sie, dass die Krümmung unabhängig ist von der Parametrisierung, das heisst, dass für jede Reparametrisierung [latex]\gamma \circ \psi[/latex] einer Kurve [latex]\gamma[/latex] wie oben gilt [latex]\kappa _{\gamma \circ \psi }(p) = \kappa _\gamma (p)[/latex] für alle Punkte [latex]p[/latex] auf der Kurve.

Übung: Existenz von Kurven vorgegebener Krümmung

Wie in vorheriger Übung möchten wir hier die Krümmung ebener Kurven betrachten, aber dabei zulassen, dass die betrachteten Kurven nicht einfach sind, womit für eine reguläre Kurve [latex]\gamma :[a,b]\to \mathbb {R}^2[/latex] die Krümmung [latex]\kappa _\gamma[/latex] definiert ist als Funktion auf [latex][a,b][/latex] via [latex]\kappa _\gamma (t) = \frac {\langle \dot {\gamma }(t),R \ddot {\gamma }(t) \rangle }{\| {\dot {\gamma }(t)}\| ^3}[/latex].

Sei nun [latex]\kappa :[a,b] \to \mathbb {R}[/latex] eine beliebige stetige Funktion auf einem Intervall [latex][a,b][/latex] mit Endpunkten [latex]a

[latex]
\begin{aligned}[]\theta :[a,b]\to \mathbb {R}^2,\ s \mapsto \int _a^s \kappa (r) \thinspace {\rm {d}} r\end{aligned}
[/latex]

und setzen

[latex]
\begin{aligned}[]\gamma : [a,b] \to \mathbb {R}^2,\ t \mapsto \int _a^t (\cos (\theta (s)),\sin (\theta (s)) \thinspace {\rm {d}} s.\end{aligned}
[/latex]

Zeigen Sie, dass die Komponenten von [latex]\gamma[/latex] jeweils zweimal stetig differenzierbar ist und dass [latex]\gamma[/latex] eine reguläre, nach Bogenlänge parametrisierte Kurve ist. Verifizieren Sie anschliessend, dass [latex]\kappa _\gamma = \kappa[/latex] gilt.

Übung: Irrationalität der Kreiszahl

In dieser Übung möchten wir zeigen, dass [latex]\pi[/latex] irrational ist, wobei wir dem Beweis von Niven [3]folgen werden.

Per Widerspruch wollen wir annehmen, dass [latex]\pi = \frac {a}{b}[/latex] ist für [latex]a,b \in \mathbb {N}[/latex]. Nun betrachten wir die Polynome

[latex]
\begin{aligned}[]f(x) &= \frac {x^n(a-bx)^n}{n!},\\ F(x) &= f(x) - f^{(2)}(x) + f^{(4)}(x) - \ldots + (-1)^n f^{(2n)}(x)\end{aligned}
[/latex]

für ein [latex]n \in \mathbb {N}[/latex], welches wir später wählen werden.

  1. Begründen Sie, wieso [latex]f[/latex], alle Ableitungen von [latex]f[/latex] und [latex]F[/latex] bei [latex]0[/latex] und bei [latex]\pi[/latex] ganzzahlige Werte annehmen. Verifizieren Sie des Weiteren, dass [latex]f|_{[0,\pi ]}[/latex] eine nicht-negative Funktion ist, welche genau bei [latex]0[/latex] und [latex]\pi[/latex] verschwindet.
  2. Zeigen Sie, dass
    [latex]
    \begin{aligned}[](F'(x) \sin (x) - F(x) \cos (x))' = f(x) \sin (x)\end{aligned}
    [/latex]

    und [latex]\int _0^\pi f(x) \sin (x) \thinspace {\rm {d}} x = F(\pi ) + F(0)[/latex].

  3. Schliessen Sie auf einen Widerspruch.

Übung: Summe der Reziproken der Primzahlen

In dieser Übung möchten wir für natürliche Zahlen [latex]n \in \mathbb {N}[/latex] die Summe [latex]\sum _{p \in \mathbb {P}: p \leq n} \frac {1}{p}[/latex] betrachten, wobei [latex]\mathbb {P} \subseteq \mathbb {N}[/latex] die Menge der Primzahlen bezeichnet. Dabei möchten wir zeigen, dass

[latex]
\begin{aligned}[]\sum _{p \in \mathbb {P}: p \leq n} \frac {1}{p} \geq \log (\log (n)) - \log (C)\end{aligned}
[/latex]

für alle [latex]n\geq 3[/latex] und eine Konstante [latex]C >1[/latex]. Insbesondere gibt es unendlich viele Primzahlen (wieso?). Die obige Ungleichung stellt eine (stark) abgeschwächte Version des zweiten Theorems von Mertens (und damit einen Vorreiter des Primzahlsatzes) dar.

  1. Verifizieren Sie für alle [latex]n\in \mathbb {N}[/latex] die Ungleichung
    [latex]
    \begin{aligned}[]\exp \bigg (\sum _{p \in \mathbb {P}: p \leq n} \frac {1}{p}\bigg ) \geq \prod _{p \in \mathbb {P}: p \leq n} \left ( 1 + \frac {1}{p}\right ).\end{aligned}
    [/latex]
  2. Zeigen Sie für alle [latex]n \in \mathbb {N}[/latex]
    [latex]
    \begin{aligned}[]\sum _{k=1}^n \frac {1}{k} \leq \bigg ( \sum _{\ell =1}^n \frac {1}{\ell ^2} \bigg ) \prod _{p \in \mathbb {P}: p \leq n} \left ( 1 + \frac {1}{p}\right )\end{aligned}
    [/latex]
  3. Schliessen Sie auf die Aussage.

Hinweis.

Für den zweiten Teil können Sie die Existenz und Eindeutigkeit einer Primfaktorzerlegung verwenden, wonach sich insbesondere jede natürliche Zahl zwischen [latex]1[/latex] und [latex]n\in \mathbb {N}[/latex] als Produkt einer Quadratzahl mit einer Zahl schreiben lässt, in deren Faktorisierung jede Primzahl nur einmal vorkommt. Für den letzten Teil können Sie Übung 8.47 benutzen.

8.8.3 – Lernkarten

Sie können wiederum die Lernkarten oder den Graphen für Ihre Wiederholung der Themen des Kapitels verwenden.

<!– post meta –>


  1. Man kann die Sinnhaftigkeit einer Definition zwar hinterfragen, doch kann man eine Definition ohnehin nicht beweisen.
  2. D. Hilbert: Über die Transzendenz der Zahlen [latex]\mathrm {e}[/latex] und [latex]\pi[/latex] (Mathematische Annalen, 1893)
  3. I. Niven: A simple proof that [latex]\pi[/latex] is irrational (Bull. Amer. Math. Soc., 1947)

License

Analysis-Skript CHAB MATH PHYS: 18/19 Copyright © by Manfred Einsiedler and Andreas Wieser. All Rights Reserved.

}