aboutsummaryrefslogtreecommitdiffstats
path: root/buch/chapters/110-elliptisch/ellintegral.tex
diff options
context:
space:
mode:
Diffstat (limited to '')
-rw-r--r--buch/chapters/110-elliptisch/ellintegral.tex703
1 files changed, 688 insertions, 15 deletions
diff --git a/buch/chapters/110-elliptisch/ellintegral.tex b/buch/chapters/110-elliptisch/ellintegral.tex
index 46659cd..466aeb7 100644
--- a/buch/chapters/110-elliptisch/ellintegral.tex
+++ b/buch/chapters/110-elliptisch/ellintegral.tex
@@ -7,7 +7,7 @@
\label{buch:elliptisch:section:integral}}
\rhead{Elliptisches Integral}
Bei der Berechnung des Ellipsenbogens in
-Abschnitt~\ref{buch:geometrie:subsection:hyperbeln-und-ellipsen}
+Abschnitt~\ref{buch:geometrie:subsection:kegelschnitte}
sind wir auf ein Integral gestossen, welches sich nicht in geschlossener
Form ausdrücken liess.
Um solche Integrale in den Griff zu bekommen, ist es nötig, sie als
@@ -172,11 +172,193 @@ die {\em Jacobi-Normalform} heisst.
\subsubsection{Vollständige elliptische Integrale als hypergeometrische
Funktionen}
-XXX Als hypergeometrische Funktionen \url{https://www.youtube.com/watch?v=j0t1yWrvKmE} \\
-
+%XXX Als hypergeometrische Funktionen \url{https://www.youtube.com/watch?v=j0t1yWrvKmE} \\
+Das vollständige elliptische Integral $K(k)$ kann mit Hilfe der
+Binomialreihe umgeformt werden in eine hypergeometrische Reihe.
+Da im Integral nur $k^2$ auftaucht, wird sich $K(k)$ als
+hypergeometrische Funktion von $k^2$ ausdrücken lassen.
+\begin{satz}
+\index{Satz!vollständiges elliptisches Integral als hypergeometrische Funktion}%
+\label{buch:elliptisch:satz:hyperK}
+Das vollständige elliptische Integral $K(k)$ lässt sich durch die
+hypergeometrische Funktion $\mathstrut_2F_1$ als
+\[
+K(k)
+=
+\frac{\pi}2
+\cdot
+\mathstrut_2F_1\biggl(
+\begin{matrix}\frac12,\frac12\\1\end{matrix};1;k^2
+\biggr)
+\]
+ausdrücken.
+\end{satz}
+\begin{proof}[Beweis]
+Zunächst ist das vollständige elliptische Integral in der Legendre-Form
+\begin{align}
+K(k)
+&=
+\int_0^{\frac{\pi}2}
+\frac{d\vartheta}{\sqrt{1-k^2\sin^2\vartheta}}
+%\notag
+%\\
+%&
+=
+\int_0^{\frac{\pi}2}
+\bigl(
+1-(k\sin\vartheta)^2
+\bigr)^{-\frac12}\,d\vartheta.
+\notag
+\intertext{Die Wurzel im letzten Integral kann mit Hilfe der binomischen
+Reihe vereinfacht werden zu}
+&=
+\sum_{n=0}^\infty
+(-1)^n k^2\binom{-\frac12}{n}
+\int_0^{\frac{\pi}2}
+\sin^{2n}\vartheta
+\,d\vartheta.
+\label{buch:elliptisch:beweis:ellharm2}
+\end{align}
+Der verallgemeinerte Binomialkoeffizient lässt sich nach
+\begin{align*}
+\binom{-\frac12}{n}
+&=
+\frac{(-\frac12)(-\frac32)(-\frac52)\cdot\ldots\cdot(-\frac12-n+1)}{n!}
+=
+(-1)^n
+\cdot
+\frac{1}{n!}
+\cdot
+\frac12\cdot\frac32\cdot\frac52\cdot\ldots\cdot\biggl(\frac12+n-1\biggr)
+=
+(-1)^n\frac{(\frac12)_n}{n!}
+\end{align*}
+vereinfachen.
+Setzt man dies in \eqref{buch:elliptisch:beweis:ellharm2} ein, erhält
+man
+\begin{align*}
+K(k)
+&=
+\sum_{n=0}^\infty
+(-1)^n k^{2n}
+\cdot
+(-1)^n
+\frac{(\frac12)_n}{n!}
+\cdot
+\int_0^{\frac{\pi}2} \sin^{2n}\vartheta\,d\vartheta
+=
+\sum_{n=0}^\infty
+\frac{(\frac12)_n}{n!}
+\int_0^{\frac{\pi}2} \sin^{2n}\vartheta\,d\vartheta
+\cdot (k^2)^n.
+\end{align*}
+Es muss jetzt also nur noch das Integral von $\sin^{2n}\vartheta$
+berechnet werden.
+Mit partieller Integration kann man
+\begin{align*}
+\int \sin^m\vartheta\,d\vartheta
+&=
+\int
+\underbrace{\sin \vartheta}_{\uparrow}
+\underbrace{\sin^{m-1}\vartheta}_{\downarrow}
+\,d\vartheta
+\\
+&=
+-\cos\vartheta\sin^{m-1}\vartheta
++
+\int \cos^2\vartheta (m-1)\sin^{m-2}\vartheta\,d\vartheta
+\\
+&=
+-\cos\vartheta \sin^{m-1}\vartheta
++
+(m-1)
+\int
+(1-\sin^2\vartheta)
+\sin^{m-2}\vartheta\,d\vartheta.
+\end{align*}
+Wegen $\sin 0=0$ und
+$\cos\frac{\pi}2=0$ verschwindet der erste Term im bestimmten Integral
+und der zweite wird
+\begin{align*}
+\int_0^{\frac{\pi}2}
+\sin^{m} \vartheta
+\,d\vartheta
+&=
+(m-1)
+\int_0^{\frac{\pi}2}
+\sin^{m-2}\vartheta\,d\vartheta
+-
+(m-1)
+\int_0^{\frac{\pi}2}
+\sin^m \vartheta\,d\vartheta
+\\
+m
+\int_0^{\frac{\pi}2}
+\sin^{m} \vartheta\,d\vartheta
+&=
+(m-1)
+\int_0^{\frac{\pi}2}
+\sin^{m-2} \vartheta\,d\vartheta
+\\
+\int_0^{\frac{\pi}2}
+\sin^{m} \vartheta\,d\vartheta
+&=
+\frac{m-1}{m}
+\int_0^{\frac{\pi}2}
+\sin^{m-2} \vartheta\,d\vartheta.
+\end{align*}
+Mit dieser Rekursionsformel kann jetzt das Integral berechnet werden.
+Es folgt
+\begin{align*}
+\int_0^{\frac{\pi}2}
+\sin^{2n}\vartheta\,d\vartheta
+&=
+\frac{2n-1}{2n}
+\int_0^{\frac{\pi}2}
+\sin^{2n-2}\vartheta\,d\vartheta
+\\
+&=
+\frac{2n-1}{2n}
+\frac{2n-3}{2n-2}
+\frac{2n-5}{2n-4}
+\cdots
+\frac{2n-(2n-1)}{2(n-1)}
+\int_0^{\frac{\pi}2}
+\sin^{2n-4}\vartheta\,d\vartheta
+\\
+&=
+\frac{
+(n-\frac12)(n-\frac32)(n-\frac52)\cdot\ldots\cdot\frac32\cdot\frac12
+}{
+n!
+}
+\int_0^{\frac{\pi}2} 1\,d\vartheta
+\\
+&=
+\frac{(\frac12)_n}{n!}
+\cdot
+\frac{\pi}2.
+\end{align*}
+Damit wird die Reihenentwicklung für $K(k)$ jetzt zu
+\[
+K(k)
+=
+\frac{\pi}2
+\sum_{n=0}^\infty
+\frac{(\frac12)_n(\frac12)_n}{n!} \cdot \frac{(k^2)^n}{n!}
+=
+\frac{\pi}2
+\cdot
+\mathstrut_2F_1\biggl(\begin{matrix}\frac12,\frac12\\1\end{matrix};k^2\biggr),
+\]
+dies beweist die Behauptung.
+\end{proof}
+%
+% Umfang einer Ellipse
+%
\subsubsection{Umfang einer Ellipse}
\begin{figure}
\centering
@@ -247,13 +429,354 @@ Für den extremen Wert $\varepsilon=0$ entsteht der Umfang einer Ellipse,
also $E(0)=\frac{\pi}2$.
Für $\varepsilon=1$ ist $a=0$, es entsteht eine Strecke mit Länge $E(1)=1$.
-\subsubsection{Komplementäre Integrale}
+\begin{satz}
+\label{buch:elliptisch:satz:hyperE}
+Das vollständige elliptische Integral $E(k)$ ist
+\[
+E(k)
+=
+\int_0^{\frac{\pi}2} \sqrt{1-k^2\sin^2\vartheta}\,d\vartheta
+=
+\frac{\pi}2
+\cdot
+\mathstrut_2F_1\biggl(
+\begin{matrix}-\frac12,\frac12\\1\end{matrix};
+k^2
+\biggr).
+\]
+\end{satz}
+
+\begin{proof}[Beweis]
+Die Identität kann wie im Satz~\ref{buch:elliptisch:satz:hyperK} mit
+Hilfe einer Entwicklung der Wurzel mit der Binomialreihe gefunden
+werden.
+\end{proof}
+
+Die Darstellung von $E(k)$ als hypergeometrische Reihe ermöglicht
+jetzt zum Beispiel auch die Berechnung der Ableitung nach dem
+Parameter $k$ mit der Ableitungsformel für die Funktion $\mathstrut_2F_1$.
+
+
+%
+% Berechnung mit dem arithmetisch-geometrischen Mittel
+%
+\subsection{Berechnung mit dem arithmetisch-geometrischen Mittel
+\label{buch:elliptisch:subsection:agm}}
+Die numerische Berechnung von elliptischer Integrale mit gewöhnlichen
+numerischen Integrationsroutinen ist nicht sehr effizient.
+Das in diesem Abschnitt vorgestellte arithmetisch-geometrische Mittel
+\index{arithmetisch-geometrisches Mittel}%
+liefert einen Algorithmus mit sehr viel besserer Konvergenz.
+Die Methode lässt sich auch auf die unvollständigen elliptischen
+Integrale von Abschnitt~\eqref{buch:elliptisch:subsection:unvollstintegral}
+verallgemeinern.
+Sie ist ein Speziallfall der sogenannten Landen-Transformation,
+\index{Landen-Transformation}%
+welche ausser für die elliptischen Integrale auch für die
+Jacobischen elliptischen Funktionen formuliert werden kann und
+für letztere ebenfalls sehr schnelle numerische Algorithmen liefert
+(siehe dazu auch die
+Aufgaben~\ref{buch:elliptisch:aufgabe:2}--\ref{buch:elliptisch:aufgabe:4}).
+Sie kann auch verwendet werden, um die Werte der Jacobischen elliptischen
+Funktionen für komplexe Argument zu berechnen.
+Eine weiter Anwendung ist die Berechnung einer grossen Zahl von
+Stellen der Kreiszahl $\pi$, siehe Aufgaben~\ref{buch:elliptisch:aufgabe:5}.
+
+%
+% Das arithmetisch-geometrische Mittel
+%
+\subsubsection{Das arithmetisch-geometrische Mittel}
+Seien $a$ und $b$ zwei nichtnegative reelle Zahlen.
+Aus $a$ und $b$ werden jetzt zwei Folgen konstruiert, deren Glieder
+durch
+\begin{align*}
+a_0&=a &&\text{und}& a_{n+1} &= \frac{a_n+b_n}2 &&\text{arithmetisches Mittel}
+\\
+b_0&=b &&\text{und}& b_{n+1} &= \sqrt{a_nb_n} &&\text{geometrisches Mittel}
+\end{align*}
+definiert sind.
+
+\begin{satz}
+\index{Satz!arithmetisch-geometrisches Mittel}%
+Falls $a>b>0$ ist, nimmt die Folge $(a_k)_{k\ge 0}$ monoton ab und
+$(b_k)_{k\ge 0}$ nimmt monoton zu.
+Beide konvergieren quadratisch gegen einen gemeinsamen Grenzwert.
+\end{satz}
+
+\begin{definition}
+Der gemeinsame Grenzwert von $a_n$ und $b_n$ heisst das
+{\em arithmetisch-geometrische Mittel} und wird mit
+\[
+M(a,b)
+=
+\lim_{n\to\infty} a_n
+=
+\lim_{n\to\infty} b_n
+\]
+bezeichnet.
+\index{arithmetisch-geometrisches Mittel}%
+\end{definition}
+
+\begin{proof}[Beweis]
+Zunächst ist zu zeigen, dass die Folgen monoton sind.
+Dies folgt sofort aus der Definition der Folgen:
+\begin{align*}
+a_{n+1} &= \frac{a_n+b_n}{2} \ge \frac{a_n+a_n}{2} = a_n
+\\
+b_{n+1} &= \sqrt{a_nb_n} \ge \sqrt{b_nb_n} = b_n.
+\end{align*}
+Die Konvergenz folgt aus
+\[
+a_{n+1}-b_{n+1}
+\le
+a_{n+1}-b_n
+=
+\frac{a_n+b_n}{2}-b_n
+=
+\frac{a_n-b_n}2
+\le
+\frac{a-b}{2^{n+1}}.
+\]
+Dies zeigt jedoch nur, dass die Konvergenz mindestens ein
+Bit in jeder Iteration ist.
+Aus
+\[
+a_{n+1}^2 - b_{n+1}^2
+=
+\frac{(a_n+b_n)^2}{4} - a_nb_n
+=
+\frac{a_n^2 -2a_nb_n+b_n^2}{4}
+=
+\frac{(a_n-b_n)^2}{4}
+\]
+folgt
+\[
+a_{n+1}-b_{n+1}
+=
+\frac{(a_n-b_n)^2}{2(a_{n+1}+b_{n+1})}.
+\]
+Da der Nenner gegen $2M(a,b)$ konvergiert, wird der Fehler für in
+jeder Iteration quadriert, die Zahl korrekter Stellen verdoppelt sich
+in jeder Iteration, es liegt also quadratische Konvergenz vor.
+\end{proof}
-\subsubsection{Ableitung}
-XXX Ableitung \\
-XXX Stammfunktion \\
+%
+% Transformation des elliptischen Integrals
+%
+\subsubsection{Transformation des elliptischen Integrals}
+In diesem Abschnitt soll das Integral
+\[
+I(a,b)
+=
+\int_0^{\frac{\pi}2}
+\frac{dt}{\sqrt{a^2\cos^2 t + b^2\sin^2t}}
+\]
+berechnet werden.
+Es ist klar, dass
+\[
+I(sa,sb)
+=
+\frac{1}{s} I(a,b).
+\]
-\subsection{Unvollständige elliptische Integrale}
+Gauss hat gefunden, dass die Substitution
+\begin{equation}
+\sin t
+=
+\frac{2a\sin t_1}{a+b+(a-b)\sin^2 t_1}
+\label{buch:elliptisch:agm:subst}
+\end{equation}
+zu
+\begin{equation}
+\frac{dt}{\sqrt{a^2_{\phantom{1}}\cos^2 t + b^2_{\phantom{1}} \sin^2 t}}
+=
+\frac{dt_1}{\sqrt{a_1^2\cos^2 t_1 + b_1^2 \sin^2 t_1}}
+\label{buch:elliptisch:agm:dtdt1}
+\end{equation}
+führt.
+Um dies nachzuprüfen, muss man zunächst
+\eqref{buch:elliptisch:agm:subst}
+nach $t_1$ ableiten, was
+\[
+\frac{d}{dt_1}\sin t
+=
+\cos t
+\frac{dt}{dt_1}
+\qquad\Rightarrow\qquad
+\biggl(
+\frac{d}{dt_1}\sin t
+\biggr)^2
+=
+(1-\sin^2t)\biggl(\frac{dt}{dt_1}\biggr)^2
+\]
+ergibt.
+Die Ableitung von $t$ nach $t_1$ kann auch aus
+\eqref{buch:elliptisch:agm:dtdt1}
+ableiten, es ist
+\[
+\biggl(
+\frac{dt}{dt_1}
+\biggr)^2
+=
+\frac{a^2_{\phantom{1}} \cos^2 t + b^2_{\phantom{1}} \sin^2 t}{a_1^2 \cos^2 t_1 + b_1^2 \sin^2 t_1}.
+\]
+Man muss also nachprüfen, dass
+\begin{equation}
+\frac{1}{1-\sin^2 t}
+\frac{d}{dt_1}\sin t
+=
+\frac{a^2 \cos^2 t + b^2 \sin^2 t}{a_1^2 \cos^2 t_1 + b_1^2 \sin^2 t_1}.
+\label{buch:elliptisch:agm:deq}
+\end{equation}
+Dazu muss man zunächst $a_1=(a+b)/2$, $b_1=\!\sqrt{ab}$ setzen.
+Ausserdem muss man $\cos^2 t$ durch $1-\sin^2t$ ersetzen und
+$\sin t$ durch \eqref{buch:elliptisch:agm:subst}.
+Auch $\cos^2 t_1$ muss man durch $1-\sin^2t_1$ ersetzt werden.
+Dann kann man nach einer langwierigen Rechnung, die sich am leichtesten
+mit einem Computer-Algebra-System ausführen lässt finden, dass
+\eqref{buch:elliptisch:agm:deq}
+tatsächlich korrekt ist.
+
+\begin{satz}
+\index{Satz!Gauss-Integrale}%
+\label{buch:elliptisch:agm:integrale}
+Für $a_1=(a+b)/2$ und $b_1=\sqrt{ab}$ gilt
+\[
+\int_0^{\frac{\pi}2}
+\frac{dt}{a^2\cos^2 t + b^2 \sin^2 t}
+=
+\int_0^{\frac{\pi}2}
+\frac{dt_1}{a_1^2\cos^2 t_1 + b_1^2 \sin^2 t_1}.
+\]
+\end{satz}
+
+Der Satz~\ref{buch:elliptisch:agm:integrale} zeigt, dass die Ersetzung
+von $a$ und $b$ durch $a_1$ und $b_1$ das Integral $I(a,b)$ nicht ändert.
+Dies gilt natürlich für alle Glieder der Folge zur Bestimmung des
+arithmetisch-geometrischen Mittels.
+
+\begin{satz}
+\index{Satz!Iab@$I(a,b)$ und arithmetisch geometrisches Mittel}%
+Für $a\ge b>0$ gilt
+\begin{equation}
+I(a,b)
+=
+\int_0^{\frac{\pi}2}
+\frac{dt}{a^2\cos^2 t + b^2\sin^2t}
+=
+\frac{\pi}{2M(a,b)}
+\end{equation}
+\end{satz}
+
+\begin{proof}[Beweis]
+Zunächst folgt aus Satz~\ref{buch:elliptisch:agm:integrale}, dass
+\[
+I(a,b)
+=
+I(a_1,b_1)
+=
+\dots
+=
+I(a_n,b_n).
+\]
+Ausserdem ist $a_n\to M(a,b)$ und $b_n\to M(a,b)$,
+damit wird
+\[
+I(a,b)
+=
+\frac{1}{M(a,b)}
+\int_0^{\frac{\pi}2}
+\frac{dt}{\sqrt{\cos^2 t + \sin^2 t}}
+=
+\frac{\pi}{2M(a,b)}.
+\qedhere
+\]
+\end{proof}
+
+%
+% Berechnung des elliptischen Integrals
+%
+\subsubsection{Berechnung des elliptischen Integrals}
+Das elliptische Integral erster Art hat eine Form, die dem Integral
+$I(a,b)$ bereits sehr ähnlich ist.
+Im die Verbindung herzustellen, berechnen wir
+\begin{align*}
+I(a,b)
+&=
+\int_0^{\frac{\pi}2}
+\frac{dt}{\sqrt{a^2\cos^2 t + b^2 \sin^2 t}}
+\\
+&=
+\frac{1}{a}
+\int_0^{\frac{\pi}2}
+\frac{dt}{\sqrt{1-\sin^2 t + \frac{b^2}{a^2} \sin^2 t}}
+\\
+&=
+\frac{1}{a}
+\int_0^{\frac{\pi}2}
+\frac{dt}{\sqrt{1-(1-\frac{b^2}{a^2})\sin^2 t}}
+=
+K(k)
+\qquad\text{mit}\qquad
+k'=\frac{b^2}{a^2},\;
+k=\sqrt{1-k^{\prime 2}}
+\end{align*}
+
+\begin{satz}
+\index{Satz!vollständige elliptische Integrale und arithmetisch-geometrisches Mittel}%
+\label{buch:elliptisch:agm:satz:Ek}
+Für $0<k\le 1$ ist
+\[
+K(k) = I(1,\sqrt{1-k^2}) = \frac{\pi}{2M(1,\sqrt{1-k^2})}
+\]
+\end{satz}
+
+%
+% Numerisches Beispiel
+%
+\subsubsection{Numerisches Beispiel}
+\begin{table}
+\centering
+\begin{tabular}{|>{$}c<{$}|>{$}c<{$}|>{$}c<{$}|>{$}c<{$}|}
+\hline
+n& a_n & b_n & \pi/2a_n \mathstrut\text{\vrule height12pt depth6pt width0pt}\\
+\hline
+\text{\vrule height12pt depth0pt width0pt}%
+0 & 1.0000000000000000000 & 0.7071067811865475243 & 1.5707963267948965579 \\
+1 & 0.8535533905932737621 & 0.8408964152537145430 & 1.\underline{8}403023690212201581 \\
+2 & 0.8472249029234941526 & 0.8472012667468914603 & 1.\underline{8540}488143993356315 \\
+3 & 0.8472130848351928064 & 0.8472130847527653666 & 1.\underline{854074677}2111781089 \\
+4 & 0.8472130847939790865 & 0.8472130847939790865 & 1.\underline{854074677301371}8463 \\
+\infty& & & 1.8540746773013719184%
+\text{\vrule height12pt depth6pt width0pt}\\
+\hline
+\end{tabular}
+\caption{Die Berechnung des arithmetisch-geometrischen Mittels für
+$a=1$ und $b=\sqrt{2}/2$ zeigt die sehr rasche Konvergenz.
+\label{buch:elliptisch:agm:numerisch}}
+\end{table}
+In diesem Abschnitt soll als Zahlenbeispiel $E(k)$ für $k=\sqrt{2}/2$
+berechnet werden.
+In diesem speziellen Fall ist $k'=k$.
+Tabelle~\ref{buch:elliptisch:agm:numerisch} zeigt die sehr rasche
+Konvergenz der Berechnung des arithmetisch-geometrischen Mittels
+von $1$ und $\sqrt{2}/2$.
+Mit Satz~\ref{buch:elliptisch:agm:satz:Ek} folgt jetzt
+\[
+K(\!\sqrt{2}/2)
+=
+\frac{\pi}{2M(1,\!\sqrt{2}/2)}
+=
+1.854074677301372.
+\]
+Die Berechnung hat nur 4 Mittelwerte, 4 Produkte, 4 Quadratwurzeln und
+eine Division erfordert.
+
+%
+% Unvollständige elliptische Integrale
+%
+\subsection{Unvollständige elliptische Integrale
+\label{buch:elliptisch:subsection:unvollstintegral}}
Die Funktionen $K(k)$ und $E(k)$ sind als bestimmte Integrale über ein
festes Intervall definiert.
Die {\em unvollständigen elliptischen Integrale} entstehen, indem die
@@ -318,12 +841,18 @@ Die Abbildung~\ref{buch:elliptisch:fig:unvollstaendigeintegrale}
zeigt Graphen der unvollständigen elliptischen Integrale für verschiedene
Werte des Parameters.
+%
+% Symmetrieeigenschaften
+%
\subsubsection{Symmetrieeigenschaften}
Die Integranden aller drei unvollständigen elliptischen Integrale
sind gerade Funktionen der reellen Variablen $t$.
Die Funktionen $F(x,k)$, $E(x,k)$ und $\Pi(n,x,k)$ sind daher
ungeraden Funktionen von $x$.
+%
+% Elliptische Integrale als komplexe Funktionen
+%
\subsubsection{Elliptische Integrale als komplexe Funktionen}
Die unvollständigen elliptischen Integrale $F(x,k)$, $F(x,k)$ und $\Pi(n,x,k)$
in Jacobi-Form lassen sich auch für komplexe Argumente interpretieren.
@@ -334,10 +863,14 @@ Die Faktoren, die in den Integranden der unvollständigen elliptischen
Integrale vorkommen, haben Nullstellen bei $\pm1$, $\pm1/k$ und
$\pm 1/\sqrt{n}$
-XXX Additionstheoreme \\
-XXX Parameterkonventionen \\
+% XXX Additionstheoreme \\
+% XXX Parameterkonventionen \\
+%
+% Wertebereich
+%
\subsubsection{Wertebereich}
+\label{buch:elliptische:subsubsection:wertebereich}
Die unvollständigen elliptischen Integrale betrachtet als reelle Funktionen
haben nur positive relle Werte.
Zum Beispiel nimmt das unvollständige elliptische Integral erster Art
@@ -427,6 +960,9 @@ l({\textstyle\frac{1}{k}})=\int_1^{\frac1{k}}
\end{equation}
ausgewertet werden.
+%
+% Komplementärmodul
+%
\subsubsection{Komplementärmodul}
Im vorangegangen Abschnitt wurde gezeigt, dass der Wertebereicht des
unvollständigen elliptischen Integrals der ersten Art als komplexe
@@ -447,7 +983,7 @@ werden, dass $1-k'^2=k^2$ ist.
\begin{definition}
Ist $0\le k\le 1$ der Modul eines elliptischen Integrals, dann heisst
-$k' = \sqrt{1-k^2}$ er {\em Komplementärmodul} oder {\em Komplement
+$k' = \sqrt{1-k^2}$ der {\em Komplementärmodul} oder {\em Komplement
des Moduls}. Es ist $k^2+k'^2=1$.
\end{definition}
@@ -530,6 +1066,9 @@ in das blaue.
\label{buch:elliptisch:fig:rechteck}}
\end{figure}
+%
+% Reelle Argument > 1/k
+%
\subsubsection{Reelle Argument $> 1/k$}
Für Argument $x> 1/k$ sind beide Faktoren im Integranden des
unvollständigen elliptischen Integrals negativ, das Integral kann
@@ -576,7 +1115,141 @@ F(x,k) = iK(k') - F\biggl(\frac1{kx},k\biggr)
für die Werte des elliptischen Integrals erster Art für grosse Argumentwerte
fest.
-\subsection{Potenzreihe}
-XXX Potenzreihen \\
-XXX Als hypergeometrische Funktionen \url{https://www.youtube.com/watch?v=j0t1yWrvKmE} \\
-XXX Berechnung mit der Landen-Transformation https://en.wikipedia.org/wiki/Landen%27s_transformation
+%
+% AGM und Berechnung von F(x,k)
+%
+\subsubsection{Berechnung von $F(x,k)$ mit dem arithmetisch-geometrischen
+Mittel\label{buch:elliptisch:subsubection:berechnung-fxk-agm}}
+Wie das vollständige elliptische Integral $K(k)$ kann auch das
+unvollständige elliptische Integral
+\begin{align*}
+F(x,k)
+&=
+\int_0^x \frac{d\xi}{\sqrt{(1-\xi^2)(1-k^{\prime 2}\xi^2)}}
+=
+\int_0^{\varphi}
+\frac{dt}{\sqrt{1-k^2 \sin^2 t}}
+&&\text{mit $x=\sin\varphi$}
+\\
+&=
+a
+\int_0^{\varphi} \frac{dt}{a^2 \cos^2 t + b^2 \sin^2 t}
+&&\text{mit $k=b/a$}
+\end{align*}
+mit dem arithmetisch-geometrischen Mittel berechnet werden.
+Dazu muss die Substitution
+\eqref{buch:elliptisch:agm:subst}
+verwendet werden, um auch den Winkel $\varphi_1$ zu berechnen.
+Zunächst wird \eqref{buch:elliptisch:agm:subst} nach $x_1=\sin t_1$
+aufgelöst.
+Durch Multiplikation mit dem Nenner erhält man mit der Abkürzung
+$x=\sin t$ %und $x_1=\sin t_1$
+die quadratische Gleichung
+\[
+(a-b)x x_1^2
+-
+2ax_1
++
+(a+b)x
+=
+0,
+\]
+mit der Lösung
+\begin{equation}
+x_1
+=
+\frac{a-\sqrt{a^2-(a^2-b^2)x^2}}{(a-b)x}.
+\label{buch:elliptisch:unvollstagm:xrek}
+\end{equation}
+Der Algorithmus zur Berechnung des arithmetisch-geometrischen Mittels
+muss daher verallgemeinert werden zu
+\begin{equation}
+\left.
+\begin{aligned}
+a_{n+1} &= \frac{a_n+b_n}2, &\qquad a_0 &= a
+\\
+b_{n+1} &= \sqrt{a_nb_n}, & b_0 &= b
+\\
+x_{n+1} &= \frac{a_n-\sqrt{a_n^2-(a_n^2-b_n^2)x_n^2}}{(a_n-b_n)x_n}, & x_0 &= x
+\end{aligned}
+\quad
+\right\}
+\label{buch:elliptisch:unvollstagm:rek}
+\end{equation}
+Die Folge $x_n$ konvergiert gegen einen Wert $x_{\infty} = \lim_{n\to\infty} x_n$.
+Der Wert des unvollständigen elliptischen Integrals ist dann der Grenzwert
+\[
+F(x,k)
+=
+\lim_{n\to\infty}
+\frac{\arcsin x_n}{M(a_n,b_n)}
+=
+\frac{\arcsin x_{\infty}}{M(1,\sqrt{1-k^2})}.
+\]
+
+In dieser Form ist die Berechnung allerdings nicht praktisch durchführbar.
+Das Problem ist, dass die Differenz $a_n-b_n$, die in
+\eqref{buch:elliptisch:unvollstagm:rek}
+im Nenner vorkommt, sehr schnell gegen Null geht.
+Ausserdem ist die Quadratwurzel im Zähler fast gleich gross wie
+$a_n$, was zu Auslöschung und damit ungenauen Resultaten führt.
+\label{buch:elliptisch:agm:ellintegral-stabilitaet}
+
+Eine Möglichkeit, das Problem zu entschärfen, ist, die Rekursionsformel
+nach $\varepsilon = a-b$ zu entwickeln.
+Mit $a+b=2a+\varepsilon$ kann man $b$ aus der Formel elimineren und erhält
+mit Hilfe der binomischen Reihe
+\begin{align*}
+x_1
+&=
+\frac{a}{x\varepsilon}
+\left(1-\sqrt{1-\varepsilon(2a-\varepsilon)x^2/a^2}\right)
+\\
+&=
+\frac{a}{x\varepsilon}
+\biggl(
+1-\sum_{k=0}^\infty
+(-1)^k
+\frac{(\frac12)_k}{k!} \varepsilon^k(2a-\varepsilon)^k\frac{x^{2k}}{a^{2k}}
+\biggr)
+\\
+&=
+\sum_{k=1}^\infty
+(-1)^{k-1}
+\frac{(\frac12)_k}{k!} \varepsilon^{k-1}(2a-\varepsilon)^k\frac{x^{2k-1}}{a^{2k-1}}
+\\
+&=
+\frac{\frac12}{1!}(2a-\varepsilon)\frac{x}{a}
+-
+\frac{\frac12\cdot(\frac12-1)}{2!}\varepsilon(2a-\varepsilon)^2\frac{x^3}{a^3}
++
+\frac{\frac12\cdot(\frac12-1)(\frac12-2)}{3!}\varepsilon^2(2a-\varepsilon)^3\frac{x^5}{a^5}
+-
+\dots
+\\
+&=
+x\biggl(1-\frac{\varepsilon}{2a}\biggr)
+\biggl(
+1
+-
+\frac{\frac12-1}{2!}\varepsilon(2a-\varepsilon)\frac{x^2}{a^2}
++
+\frac{(\frac12-1)(\frac12-2)}{3!}\varepsilon^2(2a-\varepsilon)^2\frac{x^4}{a^4}
+-
+\dots
+\biggr)
+\\
+&=
+x\biggl(1-\frac{\varepsilon}{2a}\biggr)
+\cdot
+\mathstrut_2F_1\biggl(
+\begin{matrix}-\frac12,1\\2\end{matrix};-\varepsilon(2a-\varepsilon)\frac{x^2}{a^2}
+\biggr).
+\end{align*}
+Diese Form ist wesentlich besser, aber leider kann es bei der numerischen
+Rechnung passieren, dass $\varepsilon < 0$ wird.
+
+%\subsection{Potenzreihe}
+%XXX Potenzreihen \\
+%XXX Als hypergeometrische Funktionen \url{https://www.youtube.com/watch?v=j0t1yWrvKmE} \\
+%XXX Berechnung mit der Landen-Transformation https://en.wikipedia.org/wiki/Landen%27s_transformation