Hermitesche Matrix

Eine hermitesche Matrix ist in der Mathematik eine komplexe quadratische Matrix, die gleich ihrer adjungierten Matrix ist. Die Einträge einer hermiteschen Matrix oberhalb der Hauptdiagonale ergeben sich demnach durch Spiegelung der Einträge unterhalb der Diagonale und nachfolgender komplexer Konjugation; die Einträge auf der Hauptdiagonale selbst sind alle reell. Hermitesche Matrizen sind nach dem Mathematiker Charles Hermite benannt.

Hermitesche Matrizen weisen eine Reihe besonderer Eigenschaften auf. Die Summe zweier hermitescher Matrizen ist stets wieder hermitesch. Jede komplexe quadratische Matrix lässt sich eindeutig als Summe einer hermiteschen und einer schiefhermiteschen Matrix schreiben. Das Produkt zweier hermitescher Matrizen ist wiederum hermitesch, sofern die beiden Matrizen kommutieren. Eine hermitesche Matrix ist stets normal und selbstadjungiert, sie besitzt nur reelle Eigenwerte und sie ist stets unitär diagonalisierbar. Eine wichtige Klasse hermitescher Matrizen sind positiv definite Matrizen, bei denen alle Eigenwerte positiv sind. Eine hermitesche Matrix mit reellen Einträgen ist symmetrisch.

In der linearen Algebra werden hermitesche Matrizen zur Beschreibung hermitescher Sesquilinearformen verwendet. Die Darstellungsmatrix einer komplexen selbstadjungierten Abbildung bezüglich einer Orthonormalbasis ist ebenfalls stets hermitesch. Lineare Gleichungssysteme mit hermitescher Koeffizientenmatrix lassen sich effizient und numerisch stabil lösen. Weiterhin werden hermitesche Matrizen bei Orthogonalprojektionen und bei der Polarzerlegung von Matrizen verwendet. Hermitesche Matrizen besitzen Anwendungen unter anderem in der Quantenmechanik.

Definition

Eine komplexe quadratische Matrix $A=(a_{jk})\in {ℂ}^{n\times n}$ heißt hermitesch, wenn für ihre Einträge

$a_{jk}=\overline{a_{kj}}$

für $j,k=1,\dots ,n$ gilt. Eine hermitesche Matrix stimmt daher mit ihrer adjungierten Matrix $A^H$ überein, das heißt, es gilt

$A=A^H$.

Äquivalent dazu ist eine Matrix genau dann hermitesch, wenn ihre transponierte Matrix $A^T$ gleich ihrer konjugierten Matrix $\overline{A}$ ist, also

$A^T=\bar{A}$

gilt. Eine hermitesche Matrix ist also bis auf komplexe Konjugation aller Einträge spiegelsymmetrisch bezüglich ihrer Hauptdiagonale.

Beispiele

Beispiele für hermitesche Matrizen sind ($\mathrm{i}$ stellt die imaginäre Einheit dar):

$\left(\begin{array}{c}2\end{array}\right),\left(\begin{array}{cc}1& \mathrm{i}\\ -\mathrm{i}& 1\end{array}\right),\left(\begin{array}{ccc}1& 3-\mathrm{i}& 4\\ 3+\mathrm{i}& -2& -6+\mathrm{i}\\ 4& -6-\mathrm{i}& 5\end{array}\right)$.

Allgemein haben hermitesche Matrizen der Größe $2\times 2$, $3\times 3$ und $4\times 4$ die Struktur

$\left(\begin{array}{cc}a& \overline{b}\\ b& c\end{array}\right),\left(\begin{array}{ccc}a& \overline{b}& \overline{c}\\ b& d& \overline{e}\\ c& e& f\end{array}\right),\left(\begin{array}{cccc}a& \overline{b}& \overline{c}& \overline{d}\\ b& e& \overline{f}& \overline{g}\\ c& f& h& \overline{i}\\ d& g& i& j\end{array}\right)$

mit reellen Zahlen auf der Hauptdiagonale.

Algebraische Eigenschaften

Einträge

Die Diagonaleinträge einer hermiteschen Matrix sind aufgrund von

$a_{jj}=\overline{a_{jj}}$

stets reell. Die Matrix aus den Realteilen einer hermiteschen Matrix ist stets symmetrisch, denn

$\mathrm{R}\mathrm{e}(a_{jk})=\mathrm{R}\mathrm{e}(a_{kj})$,

und die Matrix aus den Imaginärteilen einer hermiteschen Matrix stets schiefsymmetrisch, denn

$\mathrm{I}\mathrm{m}(a_{jk})=-\mathrm{I}\mathrm{m}(a_{kj})$.

Daher wird eine hermitesche Matrix durch

$n+\frac{n(n-1)}{2}+\frac{n(n-1)}{2}=n^2$

reelle Zahlen eindeutig charakterisiert. Im Vergleich dazu wird eine allgemeine komplexe $(n\times n)$-Matrix durch $2n^2$ reelle Zahlen beschrieben, also gerade doppelt so viele.

Summe

Die Summe $A+B$ zweier hermitescher Matrizen $A,B\in {ℂ}^{n\times n}$ ist stets wieder hermitesch, denn

$(A+B{)}^H=A^H+B^H=A+B$.

Zudem lässt sich jede komplexe quadratische Matrix $M\in {ℂ}^{n\times n}$ eindeutig als Summe $M=A+B$ einer hermiteschen Matrix $A$ und einer schiefhermiteschen Matrix $B$ schreiben, indem

$A=\frac{1}{2}(M+M^H)$   und   $B=\frac{1}{2}(M-M^H)$

gewählt werden.

Skalarmultiplikation

Das Produkt $cA$ einer hermiteschen Matrix $A\in {ℂ}^{n\times n}$ mit einem Skalar $c\in ℂ$ ist nur wieder hermitesch, wenn $c$ reell ist, denn dann gilt

$(cA{)}^H=\overline{c}{A}^H=cA$.

Wenn $c$ rein imaginär ist, dann ist das Produkt $cA$ schiefhermitesch. Die hermiteschen Matrizen bilden demnach keinen Untervektorraum im $ℂ$-Vektorraum der komplexen quadratischen Matrizen, sondern lediglich einen Untervektorraum im $ℝ$-Vektorraum der komplexen quadratischen Matrizen. Dieser Untervektorraum hat die Dimension $n^2$, wobei die Standardmatrizen $E_{jj}$, $1\le j\le n$, $E_{jk}+E_{kj}$ und $\mathrm{i}(E_{jk}-E_{kj})$, $1\le j<k\le n$, darin eine Basis bilden. Im Raum der hermiteschen Matrizen bilden wiederum die reellen symmetrischen Matrizen einen Untervektorraum.

Produkt

Das Produkt $AB$ zweier hermitescher Matrizen $A,B\in {ℂ}^{n\times n}$ ist im Allgemeinen nicht wieder hermitesch. Das Produkt hermitescher Matrizen ist genau dann hermitesch, wenn $A$ und $B$ kommutieren, also wenn $AB=BA$ gilt, denn dann ergibt sich

$(AB{)}^H=B^H{A}^H=BA=AB$.

Insbesondere sind damit für eine hermitesche Matrix $A$ auch alle ihre Potenzen $A^k$ mit $k\in \mathbb{N}$ und daher auch ihr Matrixexponential $e^A$ wieder hermitesch. Für eine beliebige komplexe Matrix $M\in {ℂ}^{m\times n}$ sind sowohl die $m\times m$-Matrix $M{M}^H$ als auch die $n\times n$-Matrix $M^{H}M$ stets hermitesch.

Normalität

Eine hermitesche Matrix $A\in {ℂ}^{n\times n}$ ist stets normal, denn es gilt

$A^{H}A=AA=A{A}^H$.

Jede hermitesche Matrix kommutiert also mit ihrer Adjungierten. Es gibt allerdings auch normale Matrizen, die nicht hermitesch sind, beispielsweise schiefhermitesche Matrizen.

Kongruenz

Jede komplexe Matrix $B\in {ℂ}^{n\times n}$, die kongruent zu einer hermiteschen Matrix $A\in {ℂ}^{n\times n}$ ist, ist ebenfalls hermitesch, denn es gilt

$B^H=(S^{H}AS{)}^H=S^H{A}^{H}S=S^{H}AS=B$,

wobei $S\in {ℂ}^{n\times n}$ die zugehörige Transformationsmatrix ist. Matrizen, die ähnlich zu einer hermiteschen Matrix sind, müssen jedoch nicht notwendigerweise ebenfalls hermitesch sein.

Inverse

Ist eine hermitesche Matrix $A\in {ℂ}^{n\times n}$ invertierbar, dann ist auch ihre Inverse $A^{-1}$ wieder hermitesch, denn es gilt

$(A^{-1}{)}^H=(A^H{)}^{-1}=A^{-1}$.

Für eine reguläre hermitesche Matrix $A$ sind demnach auch alle Potenzen $A^{-k}$ mit $k\in \mathbb{N}$ wieder hermitesch.

Spektrale Eigenschaften

Selbstadjungiertheit

Eine hermitesche Matrix $A\in {ℂ}^{n\times n}$ ist stets selbstadjungiert, denn es gilt mit dem komplexen Standardskalarprodukt $⟨\cdot ,\cdot ⟩$

$⟨Ax,y⟩=(Ax{)}^{H}y=x^H{A}^{H}y=x^{H}Ay=⟨x,Ay⟩$

für alle Vektoren $x,y\in {ℂ}^n$. Es gilt auch die Umkehrung und jede komplexe selbstadjungierte Matrix ist hermitesch.

Eigenwerte

Die Eigenwerte einer hermiteschen Matrix $A\in {ℂ}^{n\times n}$, das heißt die Lösungen der Eigenwertgleichung $Ax=\lambda x$, sind stets reell. Ist nämlich $\lambda \in ℂ$ ein komplexer Eigenwert von $A$ mit zugehörigem Eigenvektor $x\in {ℂ}^n$, $x\ne 0$, dann gilt mit der Selbstadjungiertheit von $A$

$\lambda ⟨x,x⟩=⟨x,\lambda x⟩=⟨x,Ax⟩=⟨Ax,x⟩=⟨\lambda x,x⟩=\overline{\lambda }⟨x,x⟩$.

Nachdem $⟨x,x⟩\ne 0$ für $x\ne 0$ ist, muss $\lambda =\overline{\lambda }$ gelten und der Eigenwert $\lambda$ damit reell sein.

Vielfachheiten

Bei jeder hermiteschen Matrix $A\in {ℂ}^{n\times n}$ stimmen die algebraischen und die geometrischen Vielfachheiten aller Eigenwerte überein. Ist nämlich $\lambda$ ein Eigenwert von $A$ mit geometrischer Vielfachheit $m$, dann existiert eine Orthonormalbasis $\{x_1,\dots ,x_m\}$ des Eigenraums von $\lambda$, welche durch $\{x_{m+1},\dots ,x_n\}$ zu einer Orthonormalbasis des Gesamtraums ${ℂ}^n$ ergänzt werden kann. Mit der unitären Basistransformationsmatrix $S=(x_1\mid \cdots \mid x_n)$ ergibt sich damit die transformierte Matrix

$C=S^{-1}AS=S^{H}AS=\left(\begin{array}{cc}\lambda I& 0\\ 0& X\end{array}\right)$

als Blockdiagonalmatrix mit den Blöcken $\lambda I\in {ℂ}^{m\times m}$ und $X\in {ℂ}^{(n-m)\times (n-m)}$. Für die Einträge $c_{jk}$ von $C$ mit $\min \{j,k\}\le m$ gilt nämlich mit der Selbstadjungiertheit von $A$ und der Orthonormalität der Basisvektoren $x_1,\dots ,x_n$

$c_{jk}=⟨x_j,A{x}_k⟩=⟨A{x}_j,x_k⟩=\lambda ⟨x_j,x_k⟩=\lambda {\delta }_{jk}$,

wobei ${\delta }_{jk}$ das Kronecker-Delta darstellt. Da $x_{m+1},\dots ,x_n$ nach Voraussetzung keine Eigenvektoren zum Eigenwert $\lambda$ von $A$ sind, kann $\lambda$ kein Eigenwert von $X$ sein. Die Matrix $C$ besitzt daher nach der Determinantenformel für Blockmatrizen den Eigenwert $\lambda$ genau mit algebraischer Vielfachheit $m$ und aufgrund der Ähnlichkeit der beiden Matrizen damit auch $A$.

Diagonalisierbarkeit

Nachdem bei einer hermiteschen Matrix $A\in {ℂ}^{n\times n}$ algebraische und geometrische Vielfachheiten aller Eigenwerte übereinstimmen und da Eigenvektoren zu verschiedenen Eigenwerten stets linear unabhängig sind, kann aus Eigenvektoren von $A$ eine Basis des ${ℂ}^n$ gebildet werden. Daher ist eine hermitesche Matrix stets diagonalisierbar, das heißt, es gibt eine reguläre Matrix $S\in {ℂ}^{n\times n}$ und eine Diagonalmatrix $D\in {ℂ}^{n\times n}$ (sogar $D\in {ℝ}^{n\times n}$), sodass

$S^{-1}AS=D$

gilt. Die Matrix $S=(x_1\mid \cdots \mid x_n)$ hat dabei die Eigenvektoren $x_1,\dots ,x_n$ als Spalten und die Matrix $D=\mathrm{diag}({\lambda }_1,\dots ,{\lambda }_n)$ hat die zu diesen Eigenvektoren jeweils zugehörigen Eigenwerte ${\lambda }_1,\dots ,{\lambda }_n$ auf der Diagonale. Durch eine Permutation der Eigenvektoren kann dabei die Reihenfolge der Diagonaleinträge von $D$ beliebig gewählt werden. Daher sind zwei hermitesche Matrizen genau dann zueinander ähnlich, wenn sie die gleichen Eigenwerte besitzen. Weiterhin sind zwei hermitesche Matrizen genau dann simultan diagonalisierbar, wenn sie kommutieren.

Unitäre Diagonalisierbarkeit

Die Eigenvektoren $x_j,x_k$ zu zwei verschiedenen Eigenwerten ${\lambda }_j\ne {\lambda }_k$ einer hermiteschen Matrix $A\in {ℂ}^{n\times n}$ sind stets orthogonal. Es gilt nämlich wiederum mit der Selbstadjungiertheit von $A$

${\lambda }_j⟨x_j,x_k⟩=⟨{\lambda }_j{x}_j,x_k⟩=⟨A{x}_j,x_k⟩=⟨x_j,A{x}_k⟩=⟨x_j,{\lambda }_k{x}_k⟩={\lambda }_k⟨x_j,x_k⟩$.

Da ${\lambda }_j$ und ${\lambda }_k$ als verschieden angenommen wurden, folgt daraus dann $⟨x_j,x_k⟩=0$. Daher kann aus Eigenvektoren von $A$ eine Orthonormalbasis des ${ℂ}^n$ gebildet werden. Damit ist eine hermitesche Matrix sogar unitär diagonalisierbar, das heißt, es gibt eine unitäre Matrix $S$, mit der

$S^{H}AS=D$

gilt. Diese Darstellung bildet die Grundlage für die Hauptachsentransformation und ist die einfachste Version des Spektralsatzes.

Kenngrößen

Aufgrund der Diagonalisierbarkeit einer hermiteschen Matrix $A\in {ℂ}^{n\times n}$ gilt für ihre Spur

$\mathrm{spur}(A)={\lambda }_1+\dots +{\lambda }_n$

und für ihre Determinante entsprechend

$\det (A)={\lambda }_1\cdot \dots \cdot {\lambda }_n$.

Spur und Determinante einer hermiteschen Matrix sind demnach stets reell. Der Rang einer hermiteschen Matrix ist gleich der Anzahl der Eigenwerte ungleich Null, also mit dem Kronecker-Delta

$\mathrm{rang}(A)=n-({\delta }_{{\lambda }_1,0}+\dots +{\delta }_{{\lambda }_n,0})$.

Eine hermitesche Matrix ist genau dann invertierbar, wenn keiner ihrer Eigenwerte Null ist. Die Spektralnorm einer hermiteschen Matrix ist

$\Vert A{\Vert }_2=\max \{|{\lambda }_1|,\dots ,|{\lambda }_n|\}$

und damit gleich dem Spektralradius der Matrix. Die Frobeniusnorm ergibt sich aufgrund der Normalität entsprechend zu

$\Vert A{\Vert }_F=\sqrt{{\lambda }_1^2+\dots +{\lambda }_n^2}$.

Abschätzungen

Nach dem Satz von Courant-Fischer liefert der Rayleigh-Quotient Abschätzungen für den kleinsten und den größten Eigenwert einer hermiteschen Matrix $A\in {ℂ}^{n\times n}$ der Form

$\min \{{\lambda }_1,\dots ,{\lambda }_n\}\le \frac{⟨x,Ax⟩}{⟨x,x⟩}\le \max \{{\lambda }_1,\dots ,{\lambda }_n\}$

für alle $x\in {ℂ}^n$ mit $x\ne 0$. Gleichheit gilt dabei jeweils genau dann, wenn $x$ ein Eigenvektor zum jeweiligen Eigenwert ist. Der kleinste und der größte Eigenwert einer hermiteschen Matrix kann demnach durch Minimierung beziehungsweise Maximierung des Rayleigh-Quotienten ermittelt werden. Eine weitere Möglichkeit zur Eigenwertabschätzung bieten die Gerschgorin-Kreise, die für hermitesche Matrizen die Form von Intervallen haben.

Definitheit

Ist $A\in {ℂ}^{n\times n}$ eine hermitesche Matrix, dann wird der Ausdruck

$Q_A(x)=x^{H}Ax=⟨x,Ax⟩$

mit $x\in {ℂ}^n$ quadratische Form von $A$ genannt. Je nachdem ob $Q_A(x)$ größer als, größer gleich, kleiner als oder kleiner gleich null für alle $x\ne 0$ ist, heißt die Matrix $A$ positiv definit, positiv semidefinit, negativ definit oder negativ semidefinit. Kann $Q_A(x)$ sowohl positive, als auch negative Vorzeichen annehmen, so heißt $A$ indefinit. Die Definitheit einer hermiteschen Matrix kann anhand der Vorzeichen ihrer Eigenwerte ermittelt werden. Sind alle Eigenwerte positiv, ist die Matrix positiv definit, sind sie alle negativ, ist die Matrix negativ definit und so weiter. Das Tripel bestehend aus den Anzahlen der positiven, negativen und Null-Eigenwerte einer hermiteschen Matrix wird Signatur der Matrix genannt. Nach dem Trägheitssatz von Sylvester bleibt die Signatur einer hermiteschen Matrix unter Kongruenztransformationen erhalten.

Verwendung

Hermitesche Sesquilinearformen

Ist $V$ ein $n$-dimensionaler komplexer Vektorraum, dann lässt sich jede Sesquilinearform $b:V\times V\to ℂ$ nach Wahl einer Basis $\{v_1,\dots ,v_n\}$ für $V$ durch die Darstellungsmatrix

$A_b=(b(v_j,v_k))\in {ℂ}^{n\times n}$

beschreiben. Ist die Sesquilinearform hermitesch, gilt also $b(v,w)=\overline{b(w,v)}$ für alle $v,w\in V$, dann ist auch die Darstellungsmatrix $A_b$ hermitesch. Umgekehrt definiert jede hermitesche Matrix $A\in {ℂ}^{n\times n}$ mittels

$b_A(x,y)=x^{H}Ay$

eine hermitesche Sesquilinearform $b_A:{ℂ}^n\times {ℂ}^n\to ℂ$. Ist eine hermitesche Matrix $A\in {ℂ}^{n\times n}$ zudem positiv definit, dann stellt $b_A$ ein Skalarprodukt im unitären Raum ${ℂ}^n$ dar.

Selbstadjungierte Abbildungen

Ist $(V,⟨\cdot ,\cdot ⟩)$ ein $n$-dimensionaler komplexer Skalarproduktraum, dann lässt sich jede lineare Abbildung $f:V\to V$ nach Wahl einer Orthonormalbasis $\{e_1,\dots ,e_n\}$ für $V$ durch die Abbildungsmatrix

$A_f=(a_{jk})\in {ℂ}^{n\times n}$

darstellen, wobei $f(e_k)=a_{1k}{e}_1+\dots +a_{nk}{e}_n$ für $k=1,\dots ,n$ ist. Die Abbildungsmatrix $A_f$ ist nun genau dann hermitesch, wenn die Abbildung $f$ selbstadjungiert ist. Dies folgt aus

$⟨f(v),w⟩=(A_{f}x{)}^{H}y=x^H{A}_f^{H}y=x^H{A}_{f}y=x^H(A_{f}y)=⟨v,f(w)⟩$,

wobei $v=x_1{e}_1+\dots +x_n{e}_n$ und $w=y_1{e}_1+\dots +y_n{e}_n$ sind.

Projektionen und Spiegelungen

Ist wieder $(V,⟨\cdot ,\cdot ⟩)$ ein $n$-dimensionaler komplexer Skalarproduktraum und ist $U$ ein $m$-dimensionaler Untervektorraum von $V$, wobei $x_1,\dots ,x_m$ die Koordinatenvektoren einer Orthonormalbasis für $U$ sind, dann ist die Orthogonalprojektionsmatrix auf diesen Untervektorraum

$A_U=x_1{x}_1^H+\dots +x_m{x}_m^H\in {ℂ}^{n\times n}$

als Summe hermitescher Rang-Eins-Matrizen ebenfalls hermitesch. Auch die Orthogonalprojektionsmatrix auf den Komplementärraum $U^{\mathrm{\perp }}$ ist aufgrund der Darstellung $A_{U^{\mathrm{\perp }}}=I-A_U$ stets hermitesch. Mit Hilfe der Projektionsmatrizen $A_U$ und $A_{U^{\perp }}$ lässt sich jeder Vektor $v\in V$ in zueinander orthogonale Vektoren $u\in U$ und $u^{\perp }\in U^{\perp }$ zerlegen. Auch die Spiegelungsmatrix $I-2A_U$ an einem Untervektorraum $U$ ist stets hermitesch.

Lineare Gleichungssysteme

Das Auffinden der Lösung eines linearen Gleichungssystems $Ax=b$ mit hermitescher Koeffizientenmatrix $A$ vereinfacht sich, wenn man die Hermitizität der Koeffizientenmatrix ausnutzt. Auf Grund der Hermitizität lässt sich die Koeffizientenmatrix $A$ als Produkt

$A=LD{L}^H$

mit einer unteren Dreiecksmatrix $L$ mit lauter Einsen auf der Diagonale und einer Diagonalmatrix $D$ schreiben. Diese Zerlegung wird beispielsweise bei der Cholesky-Zerlegung positiv definiter hermitescher Matrizen verwendet, um die Lösung des Gleichungssystems zu berechnen. Beispiele moderner Verfahren zur numerischen Lösung großer linearer Gleichungssysteme mit dünnbesetzter hermitescher Koeffizientenmatrix sind das CG-Verfahren und das MINRES-Verfahren.

Polarzerlegung

Jede quadratische Matrix $A\in {ℂ}^{n\times n}$ kann mittels der Polarzerlegung auch als Produkt

$A=QP$

einer unitären Matrix $Q\in {ℂ}^{n\times n}$ und einer positiv semidefiniten hermiteschen Matrix $P\in {ℂ}^{n\times n}$ faktorisiert werden. Die Matrix $P$ ergibt sich dabei als die Quadratwurzel von $A^{H}A$. Ist $A$ regulär, so ist $P$ positiv definit und die Polarzerlegung eindeutig mit $Q=A{P}^{-1}$.

Quantenmechanik

Die in der Quantenmechanik verwendeten Pauli-Matrizen

${\sigma }_1={\sigma }_x=\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right),{\sigma }_2={\sigma }_y=\left(\begin{array}{cc}0& -\mathrm{i}\\ \mathrm{i}& 0\end{array}\right),{\sigma }_3={\sigma }_z=\left(\begin{array}{cc}1& 0\\ 0& -1\end{array}\right)$

sind hermitesch und spurfrei. Die Pauli-Matrizen werden unter anderem zur Beschreibung von Isospin-Symmetrien verwendet. Die Gell-Mann-Matrizen sind hermitesche $(3\times 3)$-Matrizen, die in der Quantenchromodynamik eingesetzt werden.

Siehe auch

• Hermitescher Operator, eine Verallgemeinerung hermitescher Matrizen auf unendlichdimensionale Räume
• Selbstadjungierter Operator, eine weitere Verallgemeinerung hermitescher Matrizen auf unendlichdimensionale Räume

Literatur

• Gerd Fischer: Lineare Algebra. (Eine Einführung für Studienanfänger). 13. durchgesehene Auflage. Vieweg, Braunschweig u. a. 2002, ISBN 3-528-97217-3.
• Vorlage:Literatur
• Hans-Rudolf Schwarz, Norbert Köckler: Numerische Mathematik. 5. überarbeitete Auflage. Teubner, Stuttgart u. a. 2004, ISBN 3-519-42960-8.

Weblinks

• Vorlage:EoM
• Vorlage:MathWorld
• Vorlage:PlanetMath