fix all the Bessel stuff

13 files changed, 555 insertions, 110 deletions

diff --git a/buch/chapters/040-rekursion/gamma.tex b/buch/chapters/040-rekursion/gamma.tex
index 45acf9f..2b0700e 100644
--- a/buch/chapters/040-rekursion/gamma.tex
+++ b/buch/chapters/040-rekursion/gamma.tex
@@ -73,6 +73,9 @@ gilt.
 Der Plan ist, dies so umzuformen, dass man für $x$ eine beliebige
 komplexe Zahl einsetzen kann.
 
+%
+% Pochhammer-Symbol
+%
 \subsubsection{Pochhammer-Symbol}
 Die spezielle Form des Nenners und des zweiten Faktors im Zähler
 von \eqref{buch:rekursion:gamma:eqn:fakultaet}
@@ -115,6 +118,9 @@ x!
 Der erste Faktor in diesem Ausdruck enthält jetzt nur noch Dinge,
 die für beliebige $x\in\mathbb{C}$ definiert sind.
 
+%
+% Grenwertdefinition
+%
 \subsubsection{Grenzwertdefinition}
 Der zweite Bruch in \eqref{buch:rekursion:gamma:eqn:produkt3}
 besteht aus Termen, die zwar nur für natürliches $x$ definiert sind,
@@ -147,6 +153,9 @@ $x\in\mathbb{C}\setminus\{0,-1,-2,-3,\dots\}$ ist der Grenzwert
 \index{Gamma-Funktion!Grenzwertdefinition}%
 \end{definition}
 
+%
+% Rekursionsgleichung für Gamma(x)
+%
 \subsubsection{Rekursionsgleichung für $\Gamma(x)$}
 Es ist aus der Herleitung klar, dass $\Gamma(n)=(n-1)!$ sein muss.
 Wir sollten dies aber auch direkt aus der
@@ -199,15 +208,50 @@ x\lim_{n\to\infty}
 \frac{n^{x-1}}{(n+1)^{x-1}}
 \\
 &=
+x
 \Gamma(x)
 \lim_{n\to\infty} \biggl(\frac{n}{n+1}\biggr)^{x-1}
 =
-\Gamma(x),
+x\Gamma(x),
 \end{align*}
 Weil $n/(n+1)\to 1$ ist und die Funktion $z\mapsto z^{x-1}$ für alle
 nach der Definition zulässigen Werte von $x$ eine stetige Funktion ist.
 
 %
+% Gamma-Funktion und Pochhammer-Symbol
+%
+\subsubsection{Gamma-Funktion und Pochhammer-Symbol}
+Durch Iteration der Rekursionsformel für $\Gamma(x)$ folgt jetzt
+\begin{align*}
+\Gamma(x+n)
+&=
+(x+n-1) \Gamma(x+n-1)
+\\
+&=
+(x+n-1)(x+n-2)\Gamma(x+n-2)
+\\
+&=
+\underbrace{
+(x+n-1)(x+n-2)\cdots(x-1)(x)
+}_{\text{$n$ Faktoren}} \Gamma(x)
+\\
+&=(x)_n \Gamma(x).
+\end{align*}
+Damit folgt
+
+\begin{satz}
+\label{buch:rekursion:gamma:satz:gamma-pochhammer}
+Die Rekursionsformel für die Gamma-Funktion kann geschrieben werden als
+\[
+\Gamma(x+n) = (x)_n \Gamma(x).
+\]
+Das Pochhammer-Symbol $(x)_n$ ist für alle natürlichen $n$ gegeben durch
+\[
+(x)_n = \frac{\Gamma(x+n)}{\Gamma(x)}.
+\]
+\end{satz}
+
+%
 % Numerische Unzulänglichkeit der Grenzwertdefinition
 %
 \subsubsection{Numerische Unzulänglichkeiten der Grenzwertdefinition}
diff --git a/buch/chapters/050-differential/bessel.tex b/buch/chapters/050-differential/bessel.tex
index 383c360..a3237fe 100644
--- a/buch/chapters/050-differential/bessel.tex
+++ b/buch/chapters/050-differential/bessel.tex
@@ -18,6 +18,9 @@ die sich durch bekannte Funktionen ausdrücken lassen, es ist also
 nötig, eine neue Familie von speziellen Funktionen zu definieren,
 die Bessel-Funktionen.
 
+%
+% Besselsche Differentialgleichung
+%
 \subsection{Die Besselsche Differentialgleichung}
 % XXX Wo taucht diese Gleichung auf
 Die Besselsche Differentialgleichung ist die Differentialgleichung
@@ -30,6 +33,9 @@ für eine auf dem Interval $[0,\infty)$ definierte Funktion $y(x)$.
 Der Parameter $\alpha$ ist eine beliebige komplexe Zahl $\alpha\in \mathbb{C}$,
 die Lösungsfunktionen hängen daher von $\alpha$ ab.
 
+%
+% Eigenwertproblem
+%
 \subsubsection{Eigenwertproblem}
 Die Besselsche Differentialgleichung
 \eqref{buch:differentialgleichungen:eqn:bessel}
@@ -46,12 +52,15 @@ erfüllt
 \[
 By
 =
-x^2y''+xy+x^2y
+x^2y''+xy'+x^2y
 =\alpha^2 y,
 \]
 ist also eine Eigenfunktion des Bessel-Operators zum Eigenwert
 $\alpha^2$.
 
+%
+% Indexgleichung
+%
 \subsubsection{Indexgleichung}
 Die Besselsche Differentialgleichung ist eine Differentialgleichung
 der Art~\eqref{buch:differentialgleichungen:eqn:dglverallg} mit
@@ -117,8 +126,9 @@ Nur eine Lösung kann man finden, wenn
 \]
 ist.
 
-
-
+%
+% Bessel-Funktionen erster Art
+%
 \subsection{Bessel-Funktionen erster Art}
 Für $\alpha \ge 0$ gibt es immer mindestens eine Lösung der Besselgleichung
 als verallgemeinerte Potenzreihe mit $\varrho=\alpha$.
@@ -138,6 +148,9 @@ Da $F(\varrho+1)\ne 0$ ist, folgt dass $a_1=0$ sein muss.
 
 % Fall n=1 gesondert behandeln
 
+%
+% Der allgemeine Fall
+%
 \subsubsection{Der allgemeine Fall}
 Für die höheren Potenzen $n>1$ wird die Rekursionsformel für die
 Koeffizienten $a_n$ der verallgemeinerten Potenzreihe
@@ -201,10 +214,11 @@ x^\varrho\biggl(
 x^\varrho \sum_{k=0}^\infty
 \frac{1}{(\varrho+1)_k} \frac{(-x^2/4)}{k!}
 =
+x^\varrho
+\cdot
 \mathstrut_0F_1\biggl(;\varrho+1;-\frac{x^2}{4}\biggr)
 \end{align*}
-Falls also $\alpha$ kein ganzzahliges Vielfaches von $\frac12$ ist, finden
-wir zwei Lösungsfunktionen
+Wir finden also zwei Lösungsfunktionen
 \begin{align}
 y_1(x)
 %J_\alpha(x)
@@ -214,8 +228,10 @@ x^{\alpha\phantom{-}}
 \frac{1}{(\alpha+1)_k}
 \frac{(-x^2/4)^k}{k!}
 =
+x^\alpha
+\cdot
 \mathstrut_0F_1\biggl(;\alpha+1;-\frac{x^2}{4}\biggr),
-\label{buch:differentialgleichunge:bessel:erste}
+\label{buch:differentialgleichunge:bessel:eqn:erste}
 \\
 y_2(x)
 %J_{-\alpha}(x)
@@ -223,32 +239,50 @@ y_2(x)
 x^{-\alpha} \sum_{k=0}^\infty
 \frac{1}{(-\alpha+1)_k} \frac{(-x^2/4)^k}{k!}
 =
+x^{-\alpha}
+\cdot
 \mathstrut_0F_1\biggl(;-\alpha+1;-\frac{x^2}{4}\biggr).
-\label{buch:differentialgleichunge:bessel:zweite}
+\label{buch:differentialgleichunge:bessel:eqn:zweite}
 \end{align}
+Man beachte, dass die zweite Lösung für ganzzahliges $\alpha>0$ nicht
+definiert ist.
+Man kann auch direkt nachrechnen, dass diese Funktionen Lösungen
+der Besselschen Differentialgleichung sind.
 
+%
+% Bessel-Funktionen
+%
 \subsubsection{Bessel-Funktionen}
 Da die Besselsche Differentialgleichung linear ist, ist auch
 jede Linearkombination der Funktionen
-\eqref{buch:differentialgleichunge:bessel:erste}
+\eqref{buch:differentialgleichunge:bessel:eqn:erste}
 und
-\eqref{buch:differentialgleichunge:bessel:zweite}
+\eqref{buch:differentialgleichunge:bessel:eqn:zweite}
 eine Lösung.
-Man kann zum Beispiel das Pochhammer-Symbol im Nenner loswerden,
-wenn man im Nenner mit $\Gamma(\alpha+1)$
-multipliziert:
+Satz~\ref{buch:rekursion:gamma:satz:gamma-pochhammer}
+ermöglicht, das Pochhammer-Symbol durch Werte der Gamma-Funktion
+wie in
 \[
-\frac{(1/2)^\alpha}{\Gamma(\alpha+1)}
+(\alpha+1)_n = \frac{\Gamma(\alpha+k+1)}{\Gamma(\alpha+1)}
+\]
+auszudrücken.
+Damit wird
+\begin{align}
 y_1(x)
+&=
+x^\alpha
+\sum_{k=0}^\infty
+\frac{\Gamma(\alpha+1)}{\Gamma(\alpha+k+1)}
+\frac{(-x^2/4)^k}{k!}
 =
+\Gamma(\alpha+1) 2^{\alpha}
 \biggl(\frac{x}{2}\biggr)^\alpha
 \sum_{k=0}^\infty
-\frac{(-1)^k}{k!\,\Gamma(\alpha+k+1)}
-\biggl(\frac{x}{2}\biggr)^{2k}.
-\]
-Dabei haben wir es durch
-Multiplikation mit $(\frac12)^\alpha$ auch geschafft, die Funktion
-einheitlich als Funktion von $x/2$ auszudrücken.
+\frac{(-1)^k}{k!\,\Gamma(\alpha+k+1)} \biggl(\frac{x}{2}\biggr)^{2k}
+\label{buch:differentialgleichungen:bessel:normierungsgleichung}
+\end{align}
+Nur gerade der Faktor $2^\alpha\Gamma(\alpha+1)$ ist von $k$ und $x$ 
+unabhängig, daher ist die folgende Definition sinnvoll:
 
 \begin{definition}
 \label{buch:differentialgleichungen:bessel:definition}
@@ -262,8 +296,26 @@ J_{\alpha}(x)
 \biggl(\frac{x}{2}\biggr)^{2k}
 \]
 heisst {\em Bessel-Funktion erster Art der Ordnung $\alpha$}.
+\index{Bessel-Funktion!erster Art}%
 \end{definition}
 
+Die Bessel-Funktion $J_\alpha(x)$ der Ordnung $\alpha$ unterscheidet sich
+nur durch einen Normierungsfaktor von der Lösung $y_1(x)$.
+Dasselbe gilt für $J_{-\alpha}(x)$ und $y_2(x)$:
+\begin{align*}
+J_{\alpha}(x)
+&=
+\frac{1}{2^\alpha\Gamma(\alpha+1)}
+\cdot
+y_1(x)
+\\
+J_{-\alpha}(x)
+&=
+\frac{1}{2^{-\alpha}\Gamma(-\alpha+1)}
+\cdot
+y_2(x).
+\end{align*}
+
 Man beachte, dass diese Definition für beliebige ganzzahlige 
 $\alpha$ funktioniert.
 Ist $\alpha=-n<0$, $n\in\mathbb{N}$, dann hat der Nenner Pole 
@@ -285,6 +337,8 @@ J_{-n}(x)
 (-1)^n
 J_{n}(x).
 \end{align*}
+Insbesondere unterscheiden sich $J_n(x)$ und $J_{-n}(x)$ nur durch
+ein Vorzeichen.
 
 \subsubsection{Erzeugende Funktion}
 \begin{figure}
@@ -388,6 +442,9 @@ Die beiden Exponentialreihen sind
 \notag
 \end{align}
 
+%
+% Additionstheorem
+%
 \subsubsection{Additionstheorem}
 Die erzeugende Funktion kann dazu verwendet werden, das Additionstheorem
 für die Besselfunktionen zu beweisen.
@@ -438,7 +495,9 @@ J_l(x+y) &= \sum_{m=-\infty}^\infty J_m(x)J_{l-m}(y)
 für alle $l$.
 \end{proof}
 
-
+%
+% Der Fall \alpha=0
+% 
 \subsubsection{Der Fall $\alpha=0$}
 Im Fall $\alpha=0$ hat das Indexpolynom eine doppelte Nullstelle, wir
 können daher nur eine Lösung erwarten.
@@ -453,8 +512,19 @@ J_0(x)
 \]
 geschrieben werden kann.
 
-% XXX Zweite Lösung explizit durchrechnen
+Als lineare Differentialgleichung zweiter Ordnung erwarten wir noch
+eine zweite, linear unabhängige Lösung.
+Diese kann jedoch nicht allein mit der Potenzreihenmethode,
+dazu sind die Methoden der Funktionentheorie nötig.
+Im Abschnitt~\ref{buch:funktionentheorie:subsection:dglsing}
+wird gezeigt, wie dies möglich ist und auf
+Seite~\pageref{buch:funktionentheorie:subsubsection:bessel2art}
+werden die damit zu findenden Bessel-Funktionen 0-ter Ordnung und
+zweiter Art vorgestellt.
 
+%
+% Der Fall \alpha=p, p\in \mathbb{N}
+%
 \subsubsection{Der Fall $\alpha=p$, $p\in\mathbb{N}, p > 0$}
 In diesem Fall kann nur die erste
 Lösung~\eqref{buch:differentialgleichunge:bessel:erste}
@@ -467,8 +537,9 @@ J_p(x)
 \frac{(-1)^k}{k!(p+k)!}\biggl(\frac{x}{2}\biggr)^{p+2k}.
 \]
 
-TODO: Lösung für $\alpha=-n$
-
+%
+% Der Fall $\alpha=n+\frac12$
+%
 \subsubsection{Der Fall $\alpha=n+\frac12$, $n\in\mathbb{N}$}
 Obwohl $2\alpha$ eine Ganzzahl ist, sind die beiden Lösungen
 \label{buch:differentialgleichunge:bessel:erste}
@@ -491,7 +562,7 @@ Es ist
 =
 \frac{1}{2^k}\bigl(3\cdot 5\cdot\ldots\cdot (2k+1)\bigr)
 =
-\frac{(2k+1)!}{2^{2k+1}\cdot k!}
+\frac{(2k+1)!}{2^{2k}\cdot k!}
 \\
 \biggl(-\frac12 + 1\biggr)_k
 &=
@@ -508,63 +579,181 @@ Es ist
 =
 \frac{1}{2^k}\bigl(1\cdot 3 \cdot\ldots\cdot (2(k-1)+1)\bigr)
 =
-\frac{(2k-1)!}{2^{2k}\cdot (k-1)!}
+\frac{(2k-1)!}{2^{2k-1}\cdot (k-1)!}
 \end{align*}
 Damit können jetzt die Reihenentwicklungen der Lösung wie folgt
 umgeformt werden
 \begin{align*}
 y_1(x)
 &=
-\sqrt{x}
+x^\alpha
 \sum_{k=0}^\infty
 \frac{1}{(\alpha+1)_k}
 \frac{(-x^2/4)^k}{k!}
 =
 \sqrt{x}
 \sum_{k=0}^\infty
-\frac{2^{2k+1}k!}{(2k+1)!}
+\frac{2^{2k}k!}{(2k+1)!}
 \frac{(-x^2/4)^k}{k!}
 =
 \sqrt{x}
 \sum_{k=0}^\infty
 (-1)^k
-\frac{2\cdot x^{2k}}{(2k+1)!}
+\frac{x^{2k}}{(2k+1)!}
 \\
 &=
-\frac{1}{2\sqrt{x}}
+\frac{1}{\sqrt{x}}
 \sum_{k=0}^\infty
 (-1)^k
 \frac{x^{2k+1}}{(2k+1)!}
 =
-\frac{1}{2\sqrt{x}} \sin x
+\frac{1}{\sqrt{x}} \sin x
 \\
 y_2(x)
 &=
-\frac{1}{\sqrt{x}}
+x^{-\alpha}
 \sum_{k=0}^\infty
-\frac{2^{2k}\cdot (k-1)!}{(2k-1)!}
+\frac{1}{(-\alpha+1)_k}
 \frac{(-x^2/4)^k}{k!}
 =
-\frac{1}{\sqrt{x}}
+x^{-\frac12}
 \sum_{k=0}^\infty
-(-1)^k
-\frac{x^{2k}}{(2k-1)!\cdot k}
+\frac{2^{2k-1}\cdot (k-1)!}{(2k-1)!}
+\frac{(-x^2/4)^k}{k!}
 \\
 &=
-\frac{2}{\sqrt{x}}
+\frac{1}{\sqrt{x}}
 \sum_{k=0}^\infty
 (-1)^k
 \frac{x^{2k}}{(2k-1)!\cdot 2k}
 =
-\frac{2}{\sqrt{x}} \cos x.
+\frac{1}{\sqrt{x}} \cos x.
 \end{align*}
 
-% XXX Nachrechnen, dass diese Funktionen
-% XXX Lösungen der Differentialgleichung sind
-
-\subsection{Analytische Fortsetzung und Bessel-Funktionen zweiter Art}
-
-
-
+Die Bessel-Funktionen verwenden aber eine andere Normierung. 
+Die Gleichung~\eqref{buch:differentialgleichungen:bessel:normierungsgleichung}
+zeigt, dass die Bessel-Funktionen durch Division
+der Funktion $y_1(x)$ und $y_2(x)$ durch $2^\alpha \Gamma(\alpha+1)$ 
+erhalten werden können.
+Dies ergibt
+\begin{equation*}
+\renewcommand{\arraycolsep}{1pt}
+\begin{array}{rclclclcl}
+J_{\frac12}(x)
+&=&
+\displaystyle\frac{1}{2^{\frac12}\Gamma(\frac12+1)}
+y_1(x)
+&=&
+\displaystyle\frac{1}{2^{\frac12}\frac12\Gamma(\frac12)}
+y_1(x)
+&=&
+\displaystyle\frac{\sqrt{2}}{\Gamma(\frac12)}
+y_1(x)
+&=&
+\displaystyle\frac{1}{\Gamma(\frac12)}
+\sqrt{ \frac{2}{x}}
+\sin x,
+\\
+J_{-\frac12}(x)
+&=&
+\displaystyle\frac{1}{2^{-\frac12}\Gamma(-\frac12+1)}
+y_2(x)
+&=&
+\displaystyle\frac{2^{\frac12}}{\Gamma(\frac12)}
+y_2(x)
+&=&
+\displaystyle\frac{\sqrt{2}}{\Gamma(\frac12)}
+y_2(x)
+&=&
+\displaystyle\frac{1}{\Gamma(\frac12)}
+\sqrt{\frac{2}{x}}
+\cos x.
+\end{array}
+\end{equation*}
+Wegen $\Gamma(\frac12)=\sqrt{\pi}$ sind die
+halbzahligen Bessel-Funktionen daher
+\begin{align*}
+J_{\frac12}(x)
+&=
+\sqrt{\frac{2}{\pi x}} \sin x
+=
+\sqrt{\frac{2}{\pi}} x^{-\frac12}\sin x
+&
+&\text{und}&
+J_{-\frac12}(x)
+&=
+\sqrt{\frac{2}{\pi x}} \cos x
+=
+\sqrt{\frac{2}{\pi}} x^{-\frac12}\cos x.
+\end{align*}
 
+%
+% Direkte Verifikation der Lösungen
+%
+\subsubsection{Direkte Verifikation der Lösungen für $\alpha=\pm\frac12$}
+Tatsächlich führt die Anwendung des Bessel-Operators auf die beiden
+Funktionen auf
+\begin{align*}
+\sqrt{\frac{\pi}2}
+BJ_{\frac12}(x)
+&=
+\sqrt{\frac{\pi}2}
+\biggl(
+x^2J_{\frac12}''(x) + xJ_{\frac12}'(x) + x^2J_{\frac12}(x)
+\biggr)
+\\
+&=
+x^2(x^{-\frac12}\sin x)''
++
+x(x^{-\frac12}\sin x)'
++
+x^2(x^{-\frac12}\sin x)
+\\
+&=
+x^2(
+x^{-{\textstyle\frac12}}\cos x
+-{\textstyle\frac12}x^{-\frac32}\sin x
+)'
++
+x(
+x^{-\frac12}\cos x
+-{\textstyle\frac12}x^{-\frac32}\sin x
+)
++
+x^{\frac32}\sin x
+\\
+&=
+x^2(
+-x^{-\frac12}\sin x
+-{\textstyle\frac12}x^{-\frac32}\cos x
+-{\textstyle\frac12}x^{-\frac32}\cos x
++{\textstyle\frac{3}{4}}x^{-\frac52}\sin x
+)
++
+x^{\frac12}\cos x
++
+x^{-\frac12}(x-{\textstyle\frac12})\sin x
+\\
+&=
+(
+-x^{\frac32}
++{\textstyle\frac34}x^{-\frac12}
++x^{\frac32}
+-{\textstyle\frac12}x^{-\frac12}
+)
+\sin x
+=
+\frac14x^{-\frac12}\sin x
+=
+\frac14
+\sqrt{\frac{\pi}2}
+J_{\frac12}(x)
+\\
+BJ_{\frac12}(x)
+&=
+\biggl(\frac12\biggr)^2 J_{\frac12}(x).
+\end{align*}
+Dies zeigt, dass $J_{\frac12}(x)$ tatsächlich eine Eigenfunktion
+des Bessel-Operators zum Eigenwert $\alpha^2 = \frac14$.
+Analog kann man die Lösung $y_2(x)$ für $-\frac12$ verifizieren.
 
diff --git a/buch/chapters/050-differential/hypergeometrisch.tex b/buch/chapters/050-differential/hypergeometrisch.tex
index 65b3be7..87b9318 100644
--- a/buch/chapters/050-differential/hypergeometrisch.tex
+++ b/buch/chapters/050-differential/hypergeometrisch.tex
@@ -1591,7 +1591,7 @@ x\cdot
 \end{align*}
 als Lösungen.
 Die Differentialgleichung von $\mathstrut_0F_1$ sollte sich in diesem
-Fall also auf die Airy-Differentialgleichung reduzieren lassen.
+Fall also auf die Airy-Dif\-fe\-ren\-tial\-glei\-chung reduzieren lassen.
 
 Bei der Substition der Parameter in die Differentialgleichung
 \eqref{buch:differentialgleichungen:0F1:dgl} beachten wird, dass
diff --git a/buch/chapters/070-orthogonalitaet/bessel.tex b/buch/chapters/070-orthogonalitaet/bessel.tex
index 3e9412a..0ef28fd 100644
--- a/buch/chapters/070-orthogonalitaet/bessel.tex
+++ b/buch/chapters/070-orthogonalitaet/bessel.tex
@@ -1,7 +1,8 @@
 %
 % Besselfunktionen also orthogonale Funktionenfamilie
 %
-\section{Bessel-Funktionen als orthogonale Funktionenfamilie}
+\section{Bessel-Funktionen als orthogonale Funktionenfamilie
+\label{buch:orthogonalitaet:section:bessel}}
 \rhead{Bessel-Funktionen}
 Auch die Besselfunktionen sind eine orthogonale Funktionenfamilie.
 Sie sind Funktionen differenzierbaren Funktionen $f(r)$ für $r>0$
diff --git a/buch/chapters/070-orthogonalitaet/chapter.tex b/buch/chapters/070-orthogonalitaet/chapter.tex
index 4756844..fba1298 100644
--- a/buch/chapters/070-orthogonalitaet/chapter.tex
+++ b/buch/chapters/070-orthogonalitaet/chapter.tex
@@ -8,20 +8,66 @@
 \label{buch:chapter:orthogonalitaet}}
 \lhead{Orthogonalität}
 \rhead{}
+In der linearen Algebra lernt man, dass orthonormierte Basen für die
+Lösung vektorgeometrischer Probleme, bei denen auch das Skalarprodukt
+involviert ist, besonders günstig sind.
+Die Zerlegung eines Vektors in einer Basis verlangt normalerweise nach
+der Lösung eines linearen Gleichungssystems, für orthonormierte
+Basisvektoren beschränkt sie sich auf die Berechnung von Skalarprodukten.
+
+Oft dienen spezielle Funktionen als Basis der Lösungen einer linearen
+partiellen Differentialgleichung (siehe Kapitel~\ref{buch:chapter:pde}).
+Die Randbedingungen müssen dazu in der gewählten Basis von Funktionen
+zerlegt werden.
+Fourier ist es gelungen, die Idee des Skalarproduktes und der Orthogonalität
+auf Funktionen zu verallgemeinern und so zum Beispiel das Wärmeleitungsproblem
+zu lösen.
+
+Der Orthonormalisierungsprozess von Gram-Schmidt wird damit auch auf
+Funktionen anwendbar
+(Abschnitt~\ref{buch:orthogonalitaet:section:orthogonale-funktionen}),
+der Nutzen führt aber noch viel weiter.
+Da $K[x]$ ein Vektorraum ist, führt er von der Basis der Monome
+$\{1,x,x^2,\dots,x^n\}$ 
+auf orthonormierte Polynome.
+Diese haben jedoch eine ganze Reihe weiterer nützlicher Eigenschaften.
+So wird in Abschnitt~\ref{buch:orthogonal:section:drei-term-rekursion}
+gezeigt, dass sich die Werte aller Polynome einer solchen Familie mit
+einer Rekursionsformel effizient berechnen lassen, die höchstens drei
+Terme umfasst.
+In Abschnitt~\ref{buch:orthogonalitaet:section:rodrigues} werden
+die Rodrigues-Formeln vorgeführt, die Polynome durch Anwendung eines
+Differentialoperators hervorbringen.
+In Abschnitt~\ref{buch:orthogonal:section:orthogonale-polynome-und-dgl}
+schliesslich wird gezeigt, dass diese Polynome auch Eigenfunktionen
+eines selbstadjungierten Operators sind.
+Da man in der linearen Algebra auch lernt, dass die Eigenvektoren einer
+symmetrischen Matrix zu verschiedenen Eigenwerten orthogonal sind,
+ist die Orthogonalität plötzlich nicht mehr überraschend.
+
+Die Bessel-Funktionen von
+Abschnitt~\ref{buch:differntialgleichungen:section:bessel}
+sind auch Eigenfunktionen eines Differentialoperators.
+Abschnitt~\ref{buch:orthogonalitaet:section:bessel} findet das zugehörige
+Skalarprodukt, welches andeutet, dass auch für andere Funktionenfamilien
+eine entsprechende Konstruktion möglich ist.
+Das in Abschnitt~\ref{buch:integrale:subsection:sturm-liouville-problem}
+präsentierte Sturm-Liouville-Problem führt sie durch.
+Das Kapitel schliesst mit dem
+Abschnitt~\ref{buch:orthogonal:section:gauss-quadratur}
+über die Gauss-Quadratur, welche die Eigenschaften orthogonaler Polynome
+für einen besonders effizienten numerischen Integrationsalgorithmus
+ausnutzt.
+
 \input{chapters/070-orthogonalitaet/orthogonal.tex}
 \input{chapters/070-orthogonalitaet/rekursion.tex}
 \input{chapters/070-orthogonalitaet/rodrigues.tex}
-%\input{chapters/070-orthogonalitaet/jacobi.tex}
 \input{chapters/070-orthogonalitaet/legendredgl.tex}
 \input{chapters/070-orthogonalitaet/bessel.tex}
 \input{chapters/070-orthogonalitaet/sturm.tex}
 \input{chapters/070-orthogonalitaet/gaussquadratur.tex}
 
-%\section{TODO}
-%\begin{itemize}
-%\end{itemize}
-
-\section*{Übungsaufgaben}
+\section*{Übungsaufgabe}
 \rhead{Übungsaufgaben}
 \aufgabetoplevel{chapters/070-orthogonalitaet/uebungsaufgaben}
 \begin{uebungsaufgaben}
diff --git a/buch/chapters/070-orthogonalitaet/gaussquadratur.tex b/buch/chapters/070-orthogonalitaet/gaussquadratur.tex
index 2e43cec..4a25678 100644
--- a/buch/chapters/070-orthogonalitaet/gaussquadratur.tex
+++ b/buch/chapters/070-orthogonalitaet/gaussquadratur.tex
@@ -1,7 +1,8 @@
 %
 % Anwendung: Gauss-Quadratur
 %
-\section{Anwendung: Gauss-Quadratur}
+\section{Anwendung: Gauss-Quadratur
+\label{buch:orthogonal:section:gauss-quadratur}}
 \rhead{Gauss-Quadratur}
 Orthogonale Polynome haben eine etwas unerwartet Anwendung in einem
 von Gauss erdachten numerischen Integrationsverfahren.
@@ -284,7 +285,7 @@ $p(x)$ sein.
 
 Der Satz~\ref{buch:integral:satz:gaussquadratur} begründet das
 {\em Gausssche Quadraturverfahren}.
-Die in Abschnitt~\ref{buch:integral:section:orthogonale-polynome}
+Die in Abschnitt~\ref{buch:orthogonal:subsection:legendre-polynome}
 bestimmten Legendre-Polynome $P_n$ haben die im Satz
 verlangte Eigenschaft,
 dass sie auf allen Polynomen geringeren Grades orthogonal sind.
diff --git a/buch/chapters/070-orthogonalitaet/legendredgl.tex b/buch/chapters/070-orthogonalitaet/legendredgl.tex
index de8f63f..6401e98 100644
--- a/buch/chapters/070-orthogonalitaet/legendredgl.tex
+++ b/buch/chapters/070-orthogonalitaet/legendredgl.tex
@@ -3,7 +3,8 @@
 %
 % (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule
 %
-\section{Orthogonale Polynome und Differentialgleichungen}
+\section{Orthogonale Polynome und Differentialgleichungen
+\label{buch:orthogonal:section:orthogonale-polynome-und-dgl}}
 \rhead{Differentialgleichungen orthogonaler Polynome}
 Legendre hat einen ganz anderen Zugang zu den nach ihm benannten
 Polynomen gefunden.
@@ -16,6 +17,9 @@ Die Orthogonalität wird dann aus einer Verallgemeinerung der bekannten
 Eingeschaft folgen, dass Eigenvektoren einer symmetrischen Matrix zu 
 verschiedenen Eigenwerten orthogonal sind.
 
+%
+% Legendre-Differentialgleichung
+%
 \subsection{Legendre-Differentialgleichung}
 Die {\em Legendre-Differentialgleichung} ist die Differentialgleichung
 \begin{equation}
@@ -61,7 +65,10 @@ zerlegen, die als Linearkombinationen der beiden Lösungen
 $y(x)$ und $y_s(x)$ ebenfalls Lösungen der Differentialgleichung
 sind.
 
-\subsection{Potenzreihenlösung}
+%
+% Potenzreihenlösungen
+%
+\subsubsection{Potenzreihenlösung}
 Wir suchen eine Lösung in Form einer Potenzreihe um $x=0$ und 
 verwenden dazu den Ansatz
 \[
@@ -170,7 +177,10 @@ eine Polynomlösung $\bar{P}_n(x)$ vom Grad $n$ gibt.
 Dies kann aber nicht erklären, warum die so gefundenen Polynome
 orthogonal sind.
 
-\subsection{Eigenfunktionen}
+%
+% Eigenfunktionen
+% 
+\subsubsection{Eigenfunktionen}
 Die Differentialgleichung
 \eqref{buch:integral:eqn:legendre-differentialgleichung}
 Kann mit dem Differentialoperator
@@ -198,7 +208,10 @@ des Operators $D$ zum Eigenwert $n(n+1)$ sind:
 D\bar{P}_n = -n(n+1) \bar{P}_n.
 \]
 
-\subsection{Orthogonalität von $\bar{P}_n$ als Eigenfunktionen}
+%
+% Orthogonalität von P_n als Eigenfunktionen
+%
+\subsubsection{Orthogonalität von $\bar{P}_n$ als Eigenfunktionen}
 Ein Operator $A$ auf Funktionen heisst {\em selbstadjungiert}, wenn
 für zwei beliebige Funktionen $f$ und $g$ gilt
 \[
@@ -274,7 +287,10 @@ die $\bar{P}_n$ orthogonale Polynome vom Grad $n$ sind, die die
 gleiche Standardierdisierungsbedingung wie die Legendre-Polyonome
 erfüllen, also ist $\bar{P}_n(x)=P_n(x)$.
 
-\subsection{Legendre-Funktionen zweiter Art}
+%
+% Legendre-Funktionen zweiter Art
+%
+\subsubsection{Legendre-Funktionen zweiter Art}
 %Siehe Wikipedia-Artikel \url{https://de.wikipedia.org/wiki/Legendre-Polynom}
 %
 Die Potenzreihenmethode liefert natürlich auch Lösungen der
@@ -368,7 +384,7 @@ Q_1(x) = x \operatorname{artanh}x-1
 verwendet werden.
 
 %
-% 
+%  Laguerre-Differentialgleichung
 %
 \subsection{Laguerre-Differentialgleichung
 \label{buch:orthogonal:subsection:laguerre-differentialgleichung}}
@@ -429,6 +445,9 @@ ein anderer Weg zu einer zweiten Lösung gesucht werden.
 
 XXX TODO: zweite Lösung der Differentialgleichung.
 
+%
+%
+%
 \subsubsection{Die assoziierte Laguerre-Differentialgleichung}
 \index{assoziierte Laguerre-Differentialgleichung}%
 \index{Laguerre-Differentialgleichung, assoziierte}%
diff --git a/buch/chapters/070-orthogonalitaet/orthogonal.tex b/buch/chapters/070-orthogonalitaet/orthogonal.tex
index 677e865..97cd06b 100644
--- a/buch/chapters/070-orthogonalitaet/orthogonal.tex
+++ b/buch/chapters/070-orthogonalitaet/orthogonal.tex
@@ -11,9 +11,13 @@ Funktionenreihen mit Summanden zu bilden, die im Sinne eines
 Skalarproduktes orthogonal sind, welches mit Hilfe eines Integrals
 definiert sind.
 Solche Funktionenfamilien treten jedoch auch als Lösungen von
-Differentialgleichungen.
+Differentialgleichungen auf.
 Besonders interessant wird die Situation, wenn die Funktionen 
 Polynome sind.
+In diesem Abschnitt soll zunächst das Skalarprodukt definiert 
+und an Hand von Beispielen gezeigt werden, wie verschiedenartige
+interessante Familien von orthogonalen Polynomen gewonnen werden
+können.
 
 %
 % Skalarprodukt
@@ -520,7 +524,7 @@ Tabelle~\ref{buch:integral:table:legendre-polynome}.
 Die Graphen sind in Abbildung~\ref{buch:integral:orthogonal:legendregraphen}
 dargestellt.
 Abbildung~\ref{buch:integral:orthogonal:legendreortho} illustriert, 
-dass die die beiden Polynome $P_4(x)$ und $P_7(x)$ orthogonal sind.
+dass die beiden Polynome $P_4(x)$ und $P_7(x)$ orthogonal sind.
 Das Produkt $P_4(x)\cdot P_7(x)$ hat Integral $=0$.
 
 %
diff --git a/buch/chapters/070-orthogonalitaet/rekursion.tex b/buch/chapters/070-orthogonalitaet/rekursion.tex
index dc5531b..c0efc6d 100644
--- a/buch/chapters/070-orthogonalitaet/rekursion.tex
+++ b/buch/chapters/070-orthogonalitaet/rekursion.tex
@@ -31,7 +31,17 @@ für alle $n$, $m$.
 \end{definition}
 
 \subsubsection{Allgemeine Drei-Term-Rekursion für orthogonale Polynome}
-Der folgende Satz besagt, dass $p_n$ eine Rekursionsbeziehung erfüllt.
+Die Multiplikation mit $x$ macht aus einem Polynom vom Grad $n$ ein
+Polynom vom Grad $n+1$.
+Das Polynom $xp_n(x)$ lässt sich daher als Linearkombination der
+Polynome $p_k(x)$ mit $k\le n+1$ schreiben.
+Es muss also eine lineare Beziehung zwischen den Polynomen $p_k(x)$ und
+$xp_n(x)$ geben, die man nach $p_{n+1}(x)$ auflösen kann, um eine lineare
+Darstellung von $p_{n+1}(x)$ durch die $p_k(x)$ und $p_n(x)$ zu
+bekommen.
+A priori muss man damit rechnen, dass sehr viele Summanden nötig sind.
+Der folgende Satz besagt, dass $p_n(x)$ eine Rekursionsbeziehung mit
+nur drei Termen erfüllt.
 
 \begin{satz}
 \label{buch:orthogonal:satz:drei-term-rekursion}
@@ -55,9 +65,13 @@ C_{n+1} = \frac{A_{n+1}}{A_n}\frac{h_{n+1}}{h_n}.
 \end{equation}
 \end{satz}
 
+Die Rekursionsbeziehung~\eqref{buch:orthogonal:eqn:rekursion} bedeutet,
+dass sich die Werte $p_n(x)$ für alle $n$ ausgehend von $p_1(x)$ und
+$p_0(x)$ mit nur $O(n)$ Operationen ermitteln lassen.
+
 \subsubsection{Multiplikationsoperator mit $x$}
-Man kann die Relation auch nach dem Produkt $xp_n(x)$ auflösen, dann
-wird sie
+Man kann die Relation \eqref{buch:orthogonal:eqn:rekursion}
+auch nach dem Produkt $xp_n(x)$ auflösen, dann wird sie
 \begin{equation}
 xp_n(x)
 =
@@ -68,9 +82,12 @@ xp_n(x)
 \frac{C_n}{A_n}p_{n-1}(x).
 \label{buch:orthogonal:eqn:multixrelation}
 \end{equation}
-Die Multiplikation mit $x$ ist eine lineare Abbildung im Raum der Funktionen.
+Die Multiplikation mit $x$ ist eine lineare Abbildung im Raum der Funktionen,
+die wir weiter unten auch $M_x$ abkürzen.
 Die Relation~\eqref{buch:orthogonal:eqn:multixrelation} besagt, dass diese
 Abbildung in der Basis der Polynome $p_k$ tridiagonale Form hat.
+Ein Beispiel dafür ist im nächsten Abschnitt in
+\eqref{buch:orthogonal:eqn:Mx}
 
 \subsubsection{Drei-Term-Rekursion für die Tschebyscheff-Polynome}
 Eine Relation der Form~\eqref{buch:orthogonal:eqn:multixrelation}
@@ -84,6 +101,22 @@ T_{n+1}(x) = 2x\,T_n(x)-T_{n-1}(x),
 \]
 also
 $A_n=2$, $B_n=0$ und $C_n=1$.
+Die Matrixdarstellung des Multiplikationsoperators $M_x$ in der
+Basis der Tschebyscheff-Polynome hat wegen
+\eqref{buch:orthogonal:eqn:multixrelation} die Form
+\begin{equation}
+M_x
+=
+\begin{pmatrix}
+      0&\frac12&      0&      0&      0&\dots  \\
+\frac12&      0&\frac12&      0&      0&\dots  \\
+      0&\frac12&      0&\frac12&      0&\dots  \\
+      0&      0&\frac12&      0&\frac12&\dots  \\
+      0&      0&      0&\frac12&      0&\dots  \\
+ \vdots& \vdots& \vdots& \vdots& \vdots&\ddots 
+\end{pmatrix}.
+\label{buch:orthogonal:eqn:Mx}
+\end{equation}
 
 \subsubsection{Beweis von Satz~\ref{buch:orthogonal:satz:drei-term-rekursion}}
 Die Relation~\eqref{buch:orthogonal:eqn:multixrelation} zeigt auch,
diff --git a/buch/chapters/070-orthogonalitaet/rodrigues.tex b/buch/chapters/070-orthogonalitaet/rodrigues.tex
index 9fded85..9a36bdc 100644
--- a/buch/chapters/070-orthogonalitaet/rodrigues.tex
+++ b/buch/chapters/070-orthogonalitaet/rodrigues.tex
@@ -14,7 +14,8 @@ mit der Ableitung kann man den Grad aber auch senken, man könnte daher
 auch nach einer Rekursionsformel fragen, die bei einem Polynom hohen
 Grades beginnt und mit Hilfe von Ableitungen zu geringeren Graden
 absteigt.
-Solche Formeln heissen Rodrigues-Formeln nach dem Entdecker Olinde
+Solche Formeln heissen {\em Rodrigues-Formeln} nach dem Entdecker Olinde
+\index{Rodriguez, Olinde}%
 Rodrigues, der eine solche Formal als erster für Legendre-Polynome
 gefunden hat.
 
@@ -27,12 +28,17 @@ Die Skalarprodukte sollen
 \]
 sein.
 
+%
+% Pearsonsche Differentialgleichung
+%
 \subsection{Pearsonsche Differentialgleichung}
 Die {\em Pearsonsche Differentialgleichung} ist die Differentialgleichung
 \begin{equation}
 B(x) y' - A(x) y = 0,
 \label{buch:orthogonal:eqn:pearson}
 \end{equation}
+\index{Differentialgleichung!Pearsonsche}%
+\index{Pearsonsche Differentialgleichung}%
 wobei $B(x)$ ein Polynom vom Grad höchstens $2$ ist und $A(x)$ ein
 höchstens lineares Polynom.
 Die Gleichung~\eqref{buch:orthogonal:eqn:pearson}
@@ -45,20 +51,31 @@ Dann kann man die Gleichung umstellen in
 =
 \frac{A(x)}{B(x)}
 \qquad\Rightarrow\qquad
-y = \exp\biggl( \int\frac{A(x)}{B(x)}\biggr)\,dx.
+y
+=
+\exp\biggl(
+\int\frac{A(x)}{B(x)}
+\,dx
+\biggr)
+.
 \]
-Im folgenden nehmen wir zusätzlich an, dass
+Im Folgenden nehmen wir zusätzlich an, dass an den Intervallenden
 \begin{equation}
 \lim_{x\to a+} w(x)B(x) = 0,
 \qquad\text{und}\qquad
-\lim_{x\to b-} w(x)B(x) = 0.
+\lim_{x\to b-} w(x)B(x) = 0
 \end{equation}
+gilt.
+
 Falls $w(x)$ an den Intervallenden einen von $0$ verschiedenen
 Grenzwert hat, bedeutet dies, dass $B(a)=B(b)=0$ sein muss.
 Falls $w(x)$ am Intervallende divergiert, muss $B(x)$ dort eine
 Nullstelle höherer Ordnung haben, was aber für ein Polynom
 zweiten Grades nicht möglich ist.
 
+%
+% Rekursionsformel
+%
 \subsection{Rekursionsformel}
 Multiplikation mit $B(x)$ wird den Grad eines Polynomes typischerweise 
 um $2$ erhöhen, die Ableitung wird ihn wieder um $1$ reduzieren.
@@ -66,12 +83,13 @@ Etwas formeller kann man dies wie folgt formulieren:
 
 \begin{satz}
 Für alle $n\ge 0$ ist
-\[
+\begin{equation}
 q_n(x)
 =
 \frac{1}{w(x)}
 \frac{d^n}{dx^n} B(x)^n w(x)
-\]
+\label{buch:orthogonalitaet:rodrigues:eqn:rekursion}
+\end{equation}
 ein Polynom vom Grad höchstens $n$.
 \end{satz}
 
@@ -85,50 +103,65 @@ r_0(x) B(x)^n w(x)
 \\
 &=
 \frac{d^{n-1}}{dx^{n-1}}
-\bigl(r_0'(x)B(x)+ nB'(x)B(x)^{n-1}w(x) + B(x)^n w'(x) \bigr)
+\bigl(r_0'(x)B(x)+ nr_0(x)B'(x)B(x)^{n-1}w(x) + r_0(x)B(x)^n w'(x) \bigr)
 \\
 &=
 \frac{d^{n-1}}{dx^{n-1}}
-(r_0'(x)B(x)+nB'(x)+A(x)) B(x)^{n-1} w(x)
-=
+(\underbrace{r_0'(x)B(x)+nr_0(x)B'(x)+r_0(x)A(x)}_{\displaystyle = r_1(x)})
+B(x)^{n-1} w(x)
+\\
+&=
 \frac{d^{n-1}}{dx^{n-1}} r_1(x)B^{n-1}(x) w(x).
 \end{align*}
-Für die Funktionen $r_k$ gilt die Rekursionsformel
+Iterativ lässt sich eine Folge von
+Funktionen $r_k(x)$ definieren, für die Rekursionsformel
 \begin{equation}
-r_k(x) = r_{k-1}'(x)B(x) + kB'(x) + A(x).
+r_k(x) = r_{k-1}'(x)B(x) + \bigl((n+1-k)B'(x) + A(x)\bigr)r_{k-1}(x)
 \label{buch:orthogonal:rodrigues:rekursion:beweis1}
 \end{equation}
+gilt.
 Wenn $r_0(x)$ ein Polynom ist, dann sind alle Funktionen $r_k(x)$
 ebenfalls Polynome.
-Durch wiederholte Anwendung dieser Formel kann man schliessen, dass 
+Aus der Konstruktion kann man schliessen, dass
 \[
 \frac{d^n}{dx^n} r_0(x) B(x)^n w(x)
 =
 r_n(x) w(x).
 \]
-Insbesondere folgt für $r_0(x)=1$, dass man durch $w(x)$ dividieren kann
-und dass $r_n(x)=q_n(x)$.
+Insbesondere folgt für $r_0(x)=1$, dass die $n$-te Ableitung den
+Faktor $w(x)$ enthält und dass somit $r_n(x)=q_n(x)$ ein Polynom ist.
+
+Wir müssen auch noch den Grad von $r_k(x)$ bestimmen, wobei wir
+wieder von $r_0(x)=1$ ausgehen.
+Wir behaupten, dass $\deg r_k(x)\le k$ ist, und beweisen dies
+mit vollständiger Induktion.
+Für $k=0$ ist $\deg r_0(x) = 0 \le k$ die Induktionsverankerung.
 
-Wir müssen auch noch den Grad von $r_k(x)$ bestimmen.
-Dazu verwenden wir 
-\eqref{buch:orthogonal:rodrigues:rekursion:beweis1} und berechnen den
-Grad:
+Wir nehmen jetzt also an, dass $\deg r_{k-1}(x)\le k-1$ ist und
+verwenden 
+\eqref{buch:orthogonal:rodrigues:rekursion:beweis1} um den Grad zu berechnen:
 \begin{equation*}
 \deg r_k(x)
 =
 \max \bigl(
-\underbrace{\deg(r_{k-1}'(x) B(x))}_{\displaystyle \deg r_{k-1}(x) -1 + 2}
+\underbrace{\deg(r_{k-1}'(x) B(x))}_{\displaystyle (k-1) -1 + 2}
 ,
-\underbrace{\deg(B'(x))}_{\displaystyle \le 1}
+\underbrace{\deg(r_{k-1}(x)B'(x))}_{\displaystyle \le (k-1)+1}
 ,
-\underbrace{\deg(A(x))}_{\displaystyle \le 1}
+\underbrace{\deg(r_{k-1}(x)A(x))}_{\displaystyle \le (k-1)+1}
 \bigr)
-\le \max r_{k-1}(x) + 1.
+\le k.
 \end{equation*}
-Aus $\deg r_0(x)=0$ kann man jetzt ablesen, dass $\deg r_k(x)\le k$ ist.
-Damit ist gezeigt, dass $\deg q_n(x)\le n$.
+Damit ist der Induktionsschritt und $\deg r_k(x)\le k$ bewiesen.
+Damit ist auch gezeigt, dass $\deg q_n(x)\le n$.
 \end{proof}
 
+Die Rodrigues-Formel~\eqref{buch:orthogonalitaet:rodrigues:eqn:rekursion}
+produziert eine Folge von Polynomen aufsteigenden Grades, es ist aber
+noch nicht klar, dass diese Polynome bezüglich des gewählten Skalarproduktes
+orthogonal sind.
+Dies ist der Inhalt des folgenden Satzes.
+
 \begin{satz}
 Es gibt Konstanten $c_n$ derart, dass
 \[
@@ -140,7 +173,7 @@ gilt.
 \end{satz}
 
 \begin{proof}[Beweis]
-Wir müssen zeigen, dass die Polynome orthogonal sind auf allen Monomen
+Wir zeigen, dass die Polynome orthogonal sind auf allen Monomen
 von geringerem Grad.
 \begin{align*}
 \langle q_n, x^k\rangle_w
@@ -148,15 +181,17 @@ von geringerem Grad.
 \int_a^b q_n(x)x^kw(x)\,dx
 \\
 &=
-\int_a^b \frac{1}{w(x)}\frac{d^n}{dx^n}(B(x)^n w(x)) x^k w(x)\,dx
+\int_a^b \frac{1}{w(x)}
+\biggl(\frac{d^n}{dx^n}\bigl(B(x)^n w(x)\bigr)\biggr)
+x^k w(x)\,dx
 \\
 &=
-\int_a^b \frac{d^n}{dx^n}(B(x)^n w(x)) x^k \,dx
+\int_a^b \frac{d^n}{dx^n}\bigl(B(x)^n w(x)\bigr) x^k \,dx
 \\
 &=
-\biggl[\frac{d^{n-1}}{dx^{n-1}}(B(x)^n w(x)) x^k \biggr]_a^b
+\biggl[\frac{d^{n-1}}{dx^{n-1}}\bigl(B(x)^n w(x)\bigr) x^k \biggr]_a^b
 -
-\int_a^b \frac{d^{n-1}}{dx^{n-1}}(B(x)^n w(x))kx^{k-1}\,dx
+\int_a^b \frac{d^{n-1}}{dx^{n-1}}\bigl(B(x)^n w(x)\bigr)kx^{k-1}\,dx
 \end{align*}
 Durch $n$-fache Iteration wird das Integral auf $0$ reduziert.
 Es bleiben nur die eckigen Klammern stehen, doch wenn man die Produktregel
@@ -164,9 +199,20 @@ auswertet, bleibt immer mindestens ein Produkt $B(x)w(x)$ stehen,
 nach den Voraussetzungen an den Grenzwert dieses Produktes an den
 Intervallenden verschwinden diese Terme alle.
 Damit sind die $q_n(x)$ Polynome, die $w$-orthogonal sind auf allen
-$x^k$ mit $k<n$, also Vielfache der $w$-Orthgonalpolynome.
+$x^k$ mit $k<n$.
+
+Die Polynome $q_k(x)$ mit $k< n$ haben Grad $<n$ und sind daher
+Linearkombinationen von Monomen vom Grad $<n$.
+Soeben wurde gezeigt, dass $q_n(x)$ orthogonal auf diesen Monomen
+ist, also auch auf $q_k(x)$ mit $k<n$.
+Damit ist gezeigt, dass Polynome $q_n(x)$ eine orthogonale Familie
+von Polynomen bilden.
+Durch Normierung müssen sich daraus die Polynome $p_n(x)$ ergeben.
 \end{proof}
 
+%
+% Legendre-Polynome
+%
 \subsubsection{Legendre-Polynome}
 Legendre-Polynome sind orthogonale Polynome zum Standardskalarprodukt
 mit $w(x)=1$.
@@ -195,6 +241,9 @@ P_n(x)
 (x^2-1)^n.
 \]
 
+%
+% Hermite-Polynome
+%
 \subsubsection{Hermite-Polynome}
 Die Hermite-Polynome sind auf ganz $\mathbb{R}$ definiert und verwenden
 die Gewichtsfunktion
@@ -205,13 +254,13 @@ Für jedes beliebige Polynome $B(x)$, auch für höheren Grad als $2$, ist
 \[
 \lim_{x\to-\infty} B(x) w(x)
 =
-\lim_{x\to-\infty} B(x)^e{-x^2}
+\lim_{x\to-\infty} B(x)e^{-x^2}
 =
 0
 \qquad\text{und}\qquad
 \lim_{x\to\infty} B(x) w(x)
 =
-\lim_{x\to\infty} B(x)^e{-x^2}
+\lim_{x\to\infty} B(x)e^{-x^2}
 =
 0,
 \]
@@ -222,7 +271,7 @@ Die Ableitung der Gewichtsfunktion ist
 \[
 w'(x) = -2xe^{-x^2}.
 \]
-Eingsetzt in die Pearsonsche Differentialgleichung findet man
+Eingesetzt in die Pearsonsche Differentialgleichung findet man
 \[
 \frac{w'(x)}{w(x)}
 =
@@ -238,6 +287,8 @@ B(x) = 1.
 \]
 Die Gradbedingung ist also immer erfüllt und es folgt die Rodrigues-Formel
 für die Hermite-Polynome
+\index{Hermite-Polynom}%
+\index{Polynome!Hermite}%
 \begin{equation}
 H_n(x)
 =
@@ -249,13 +300,15 @@ e^{x^2}\frac{d^n}{dx^n} e^{-x^2}.
 \label{buch:orthogonal:eqn:hermite-rodrigues}
 \end{equation}
 
-Die Hermite-Polynome können mit der Rodrigues-Formel berechnen, aber die 
-Form~\eqref{buch:orthogonal:eqn:hermite-rodrigues} ist dazu nicht gut
-geeignet.
-Dazu dient die Berechnung 
+Die Hermite-Polynome können mit der Rodrigues-Formel berechnet werden,
+aber die Form~\eqref{buch:orthogonal:eqn:hermite-rodrigues} ist dazu
+nicht gut geeignet.
+Zur Vereinfachung dient die Berechnung 
 \[
 -\frac{d}{dx}
+\bigl(
 e^{-x^2}f(x)
+\bigr)
 =
 2xe^{-x^2}f(x)
 -
@@ -270,15 +323,15 @@ vertauscht werden kann, wenn er durch die grosse Klammer auf der
 rechten Seite ersetzt wird.
 Die Rodrigues-Formel bekommt daher die Form
 \[
-H_n(x) = \biggl(\frac{d}{dx}-2x\biggr)^n \cdot 1
+H_n(x) = \biggl(2x-\frac{d}{dx}\biggr)^n \cdot 1.
 \]
 
-TODO: Relation zu hypergeometrischen Funktionen $\mathstrut_1F_1$
+%TODO: Relation zu hypergeometrischen Funktionen $\mathstrut_1F_1$
 
 %\url{https://en.wikipedia.org/wiki/Rodrigues%27_formula}
 
 %
-% Jacoib-Gewichtsfunktion
+% Jacobi-Gewichtsfunktion
 %
 \subsubsection{Jacobi-Gewichtsfunktion}
 %(%i1) w: (1-x)^a*(1+x)^b;
@@ -303,6 +356,8 @@ TODO: Relation zu hypergeometrischen Funktionen $\mathstrut_1F_1$
 %                                    x  - 1
 %
 Die Jacobi-Gewichtsfunktion 
+\index{Jacobi-Gewichtsfunktion}%
+\index{Gewichtsfunktion!Jacobi}%
 \[
 w(x)
 =
@@ -357,9 +412,14 @@ Die Konstanten $c_n$ werden durch die Normierung
 % XXX in welchem Abschnitt
 festgelegt.
 
+%
+% Tschebyscheff-Gewichtsfunktion
+%
 \subsubsection{Die Tschebyscheff-Gewichtsfunktion}
 Die Tschebyscheff-Gewichtsfunktion ist der Spezialfall $a=b=-\frac12$
 der Jacobi-Gewichtsfunktion.
+\index{Tschebyscheff-Gewichtsfunktion}%
+\index{Gewichtsfunktion!Tschebyscheff}%
 Die Rodrigues-Formel für die Tschebyscheff-Polynome lautet daher
 \[
 T_n(x)
@@ -373,8 +433,13 @@ c_n\sqrt{1-x^2} \frac{d^n}{dx^n}
 \]
 wobei wir den korrekten Wert von $c_n$ nicht nachgewiesen haben.
 
+%
+% Laguerre Gewichtsfunktion
+%
 \subsubsection{Die Laguerre-Gewichtsfunktion}
 Die Laguerre-Gewichtsfunktion
+\index{Laguerre-Gewichtsfunktion}%
+\index{Gewichtsfunktion!Laguerre}%
 \[
 w_{\text{Laguerre}}(x)
 =
diff --git a/buch/chapters/070-orthogonalitaet/sturm.tex b/buch/chapters/070-orthogonalitaet/sturm.tex
index 35054ab..1ba0ecb 100644
--- a/buch/chapters/070-orthogonalitaet/sturm.tex
+++ b/buch/chapters/070-orthogonalitaet/sturm.tex
@@ -11,6 +11,9 @@ konnte die Orthogonalität der Funktionen dadurch gezeigt werden,
 dass sie als Eigenfunktionen eines bezüglich eines geeigneten
 Skalarproduktes selbstadjungierten Operators erkannt wurden.
 
+%
+% Differentialgleichungen
+%
 \subsection{Differentialgleichung}
 Das klassische Sturm-Liouville-Problem ist das folgende Eigenwertproblem.
 Gesucht sind Lösungen der Differentialgleichung
@@ -30,6 +33,9 @@ erfüllen, wobei $|k_i|^2 + |h_i|^2\ne 0$ mit $i=a,b$.
 Weitere Bedingungen an die Funktionen $p(x)$, $q(x)$, $w(x)$  sowie die
 Lösungsfunktionen $y(x)$ sollen später geklärt werden.
 
+%
+% Das verallgemeinerte Eigenwertproblem für symmetrische Matrizen
+%
 \subsection{Das verallgemeinerte Eigenwertproblem für symmetrische Matrizen}
 Ein zu \eqref{buch:integrale:eqn:sturm-liouville} analoges Eigenwertproblem
 für Matrizen ist das folgende verallgemeinerte Eigenwertproblem.
@@ -175,6 +181,9 @@ ist damit ein gewöhnliches Eigenwertproblem für selbstadjungierte
 Matrizen des Operators $\tilde{A}$ bezüglich des verallgemeinerten
 Skalarproduktes $\langle\,\;,\;\rangle_B$.
 
+%
+% Der Operator L_0 und die Randbedingung
+%
 \subsection{Der Operator $L_0$ und die Randbedingung}
 Die Differentialgleichung kann auch in Operatorform geschrieben werden.
 Dazu schreiben wir
@@ -275,6 +284,9 @@ Ausgeschrieben bedeutet dies, dass die Randbedingung
 \eqref{buch:integrale:sturm:randbedingung}
 erfüllt sein muss.
 
+%
+% Skalarprodukt
+%
 \subsection{Skalarprodukt}
 Das Ziel der folgenden Abschnitte ist, das Sturm-Liouville-Problem als
 Eigenwertproblem für einen selbstadjungierten Operator in einem 
@@ -314,6 +326,9 @@ mit der Gewichtsfunktion $w(x)$ verwendet werden.
 Damit dies ein vernünftiges Skalarprodukt ist, muss $w(x)>0$ im
 Innerend es Intervalls sein.
 
+%
+% Der Vektorraum H
+%
 \subsection{Der Vektorraum $H$}
 Damit können wir jetzt die Eigenschaften der in Frage kommenden
 Funktionen zusammenstellen.
@@ -346,7 +361,10 @@ f\in L^2([a,b],w)\;\bigg|\;
 \biggr\}.
 \]
 
-\subsection{Differentialoperator}
+%
+% Der Sturm-Liouville-Differentialoperator
+%
+\subsection{Der Sturm-Liouville-Differentialoperator}
 Das verallgemeinerte Eigenwertproblem für $A$ und $B$ ist ein
 gewöhnliches Eigenwertproblem für die Operator $\tilde{A}=B^{-1}A$
 bezüglich des modifizierten Skalarproduktes.
@@ -366,12 +384,18 @@ $\lambda$ ist der zu $y(x)$ gehörige Eigenwert.
 Der Operator ist definiert auf Funktionen des im vorangegangenen Abschnitt
 definierten Vektorraumes $H$.
 
+%
+% Beispiele
+%
 \subsection{Beispiele}
 Die meisten der früher vorgestellten Funktionenfamilien stellen sich
 als Lösungen eines geeigneten Sturm-Liouville-Problems heraus.
 Alle Eigenschaften aus der Sturm-Liouville-Theorie gelten daher
 automatisch für diese Funktionenfamilien.
 
+%
+% Trignometrische Funktionen
+%
 \subsubsection{Trigonometrische Funktionen}
 Die trigonometrischen Funktionen sind Eigenfunktionen des Operators
 $d^2/dx^2$, also eines Sturm-Liouville-Operators mit $p(x)=1$, $q(x)=0$
@@ -434,6 +458,9 @@ Dann ist wegen
 die Bedingung~\eqref{buch:integrale:sturm:sabedingung}
 ebenfalls erfüllt, $L_0$ ist in diesem Raum selbstadjungiert.
 
+%
+% Bessel-Funktionen J_n(x)
+%
 \subsubsection{Bessel-Funktionen $J_n(x)$}
 Der Bessel-Operator \eqref{buch:differentialgleichungen:bessel-operator}
 kann wie folgt in die Form eines Sturm-Liouville-Operators gebracht 
@@ -478,6 +505,9 @@ Es folgt damit sofort, dass die Besselfunktionen orthogonale
 Funktionen bezüglich des Skalarproduktes mit der Gewichtsfunktion
 $w(x)=1/x$ sind.
 
+%
+% Bessel-Funktionen J_n(sx)
+%
 \subsubsection{Bessel-Funktionen $J_n(s x)$}
 Das Sturm-Liouville-Problem mit den Funktionen
 \eqref{buch:orthogonal:sturm:bessel:n}
@@ -608,6 +638,9 @@ Damit sind geeignete Randbedingungen für das Sturm-Liouville-Problem
 gefunden.
 \end{proof}
 
+%
+% Laguerre-Polynome
+%
 \subsubsection{Laguerre-Polynome}
 Die Laguerre-Polynome sind orthogonal bezüglich des Skalarprodukts
 mit der Laguerre-Gewichtsfunktion $w(x)=e^{-x}$ und erfüllen die
@@ -646,6 +679,9 @@ also die Laguerre-Differentialgleichung.
 Somit folgt, dass die Laguerre-Polynome orthogonal sind bezüglich
 des Skalarproduktes mit der Laguerre-Gewichtsfunktion.
 
+%
+% Tschebyscheff-Polynome
+%
 \subsubsection{Tschebyscheff-Polynome}
 Die Tschebyscheff-Polynome sind Lösungen der
 Tschebyscheff-Differentialgleichung
@@ -685,10 +721,15 @@ bezüglich des Skalarproduktes
 \langle f,g\rangle = \int_{-1}^1 f(x)g(x)\frac{dx}{\sqrt{1-x^2}}.
 \]
 
+%
+% Jacobi-Polynome
+%
 \subsubsection{Jacobi-Polynome}
 TODO
 
-
+%
+% Hypergeometrische Differentialgleichungen
+%
 \subsubsection{Hypergeometrische Differentialgleichungen}
 %\url{https://encyclopediaofmath.org/wiki/Hypergeometric_equation}
 Auch die Eulersche hypergeometrische Differentialgleichung
diff --git a/buch/chapters/080-funktionentheorie/anwendungen.tex b/buch/chapters/080-funktionentheorie/anwendungen.tex
index 4cdf9be..04c597e 100644
--- a/buch/chapters/080-funktionentheorie/anwendungen.tex
+++ b/buch/chapters/080-funktionentheorie/anwendungen.tex
@@ -5,6 +5,7 @@
 %
 \section{Anwendungen
 \label{buch:funktionentheorie:section:anwendungen}}
+\rhead{Anwendungen}
 
 \input{chapters/080-funktionentheorie/gammareflektion.tex}
 \input{chapters/080-funktionentheorie/carlson.tex}
diff --git a/buch/chapters/080-funktionentheorie/singularitaeten.tex b/buch/chapters/080-funktionentheorie/singularitaeten.tex
index 71d1844..07204ab 100644
--- a/buch/chapters/080-funktionentheorie/singularitaeten.tex
+++ b/buch/chapters/080-funktionentheorie/singularitaeten.tex
@@ -421,7 +421,8 @@ in die ursprüngliche Differentialgleichung ein, verschwindet der
 $\log(z)$-Term und für die verbleibenden Koeffizienten kann die
 bekannte Methode des Koeffizientenvergleichs verwendet werden.
 
-\subsubsection{Bessel-Funktionen zweiter Art}
+\subsubsection{Bessel-Funktionen zweiter Art
+\label{buch:funktionentheorie:subsubsection:bessel2art}}