Legendresche Identität

Die Legendresche Identität oder auch Legendresche Relation ist eine mathematische Identität aus der Infinitesimalrechnung. Sie handelt von vollständigen Elliptischen Integralen erster und zweiter Art. Diese Identität wurde vom französischen Mathematiker Adrien-Marie Legendre entdeckt und nach diesem benannt.

Definition

Folgende Formel definiert und beschreibt die Legendresche Identität:

$K(\varepsilon )E({\sqrt {1-\varepsilon ^{2}}})+E(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}})-K(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}})={\tfrac {1}{2}}\pi$

Diese Formel ist für alle reellen Werte 0 < ε < 1 gültig. Sie stellt die vollständigen elliptischen Integrale erster und zweiter Art von einem elliptischen Modul ε und dessen Pythagoräischen Gegenstück zueinander in Beziehung. In leicht abgewandelter Form kann die Legendresche Identität für denselben Definitionsbereich von ε auch in Bezug auf tangentielle Gegenstücke von elliptischen Modulen formuliert werden:

$(1+\varepsilon )K(\varepsilon )E({\tfrac {1-\varepsilon }{1+\varepsilon }})+{\tfrac {2}{1+\varepsilon }}E(\varepsilon )K({\tfrac {1-\varepsilon }{1+\varepsilon }})-2K(\varepsilon )K({\tfrac {1-\varepsilon }{1+\varepsilon }})={\tfrac {1}{2}}\pi$

Die vollständigen elliptischen Integrale selbst sind so definiert:

$K(\varepsilon )=\int _{0}^{\pi /2}{\frac {1}{\sqrt {1-\varepsilon ^{2}\sin(\varphi )^{2}}}}\,\mathrm {d} \varphi =2\int _{0}^{1}{\frac {1}{\sqrt {(x^{2}+1)^{2}-4\varepsilon ^{2}x^{2}}}}\,\mathrm {d} x$

$E(\varepsilon )=\int _{0}^{\pi /2}{\sqrt {1-\varepsilon ^{2}\sin(\varphi )^{2}}}\,\mathrm {d} \varphi =2\int _{0}^{1}{\frac {\sqrt {(x^{2}+1)^{2}-4\varepsilon ^{2}x^{2}}}{(x^{2}+1)^{2}}}\,\mathrm {d} x$

Nach einer exemplarischen Ausführung der obersten Formel über die Pythagoräischen Gegenstücke gilt somit beispielsweise:

$K({\tfrac {3}{5}})E({\tfrac {4}{5}})+E({\tfrac {3}{5}})K({\tfrac {4}{5}})-K({\tfrac {3}{5}})K({\tfrac {4}{5}})={\tfrac {1}{2}}\pi$

Und nach einer exemplarischen Ausführung der zweitobersten Formel über die tangentiellen Gegenstücke gilt zum Beispiel:

${\tfrac {4}{3}}K({\tfrac {1}{3}})E({\tfrac {1}{2}})+{\tfrac {3}{2}}E({\tfrac {1}{3}})K({\tfrac {1}{2}})-2K({\tfrac {1}{3}})K({\tfrac {1}{2}})={\tfrac {1}{2}}\pi$

Geschichte

Der Mathematiker Adrien-Marie Legendre schrieb in seinem Werk Exercices de calcul intégral sur divers ordres de transcendantes et sur les quadratures aus dem Jahre 1811 den in der soeben genannten Definition beschriebenen Zusammenhang nieder. In diesem Werk gründete er die sogenannte Legendresche Normalform. Darin führte er auch die Aufteilung der elliptischen Integrale in drei Kategorien ein, nämlich in die erster Art, zweiter Art und dritter Art. Zu dieser Zeit gehörte Legendre der Académie des sciences in Paris an. In einem weiteren Werk, Traité des fonctions elliptiques et des intégrales eulériennes aus dem Jahre 1825, leitete er seine Identität noch ausführlicher her. In dem Werk analysierte er vor allem die Additionstheoreme der elliptischen Funktionen.

Beweisführung

Spezielle Legendresche Identität für den lemniskatischen Fall

Für den lemniskatischen Fall ist das elliptische Modul beziehungsweise die spezifische Exzentrizität ε gleich der Hälfte der Quadratwurzel aus zwei. Die elliptischen Integrale erster Art handeln von den Bogenmaßen der Lemniskate von Bernoulli und die elliptischen Integrale zweiter Art von den Bogenmaßen einer Ellipse mit der Quadratwurzel aus zwei als zugehöriges Halbachsenverhältnis. Die Legendresche Identität für den lemniskatischen Fall kann so bewiesen werden:

Nach der Kettenregel gelten diese vier Ableitungen:

${\frac {\mathrm {d} }{\mathrm {d} y}}\left[K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-F{\biggl [}\arccos(xy);{\frac {1}{2}}{\sqrt {2}}{\biggr ]}\right]={\frac {{\sqrt {2}}\,x}{\sqrt {1-x^{4}y^{4}}}}$

${\frac {\mathrm {d} }{\mathrm {d} y}}\left[2E{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-2E{\biggl [}\arccos(xy);{\frac {1}{2}}{\sqrt {2}}{\biggr ]}+F{\biggl [}\arccos(xy);{\frac {1}{2}}{\sqrt {2}}{\biggr ]}\right]={\frac {{\sqrt {2}}\,x^{3}y^{2}}{\sqrt {1-x^{4}y^{4}}}}$

${\frac {\mathrm {d} }{\mathrm {d} x}}\left[{\frac {1}{y^{2}}}(y^{2}+1){\biggl [}{\text{artanh}}{\bigl (}y^{2}{\bigr )}-{\text{artanh}}{\biggl (}{\frac {{\sqrt {1-x^{4}}}\,y^{2}}{\sqrt {1-x^{4}y^{4}}}}{\biggr )}{\biggr ]}\right]={\frac {2\,x^{3}(y^{2}+1)}{\sqrt {(1-x^{4})(1-x^{4}y^{4})}}}$

${\frac {\mathrm {d} }{\mathrm {d} y}}{\biggl [}2\arctan(y)-{\frac {1}{y}}(1-y^{2})\,{\text{artanh}}(y^{2}){\biggr ]}={\frac {1}{y^{2}}}(y^{2}+1)\operatorname {artanh} (y^{2})$

Dann ist folgende Formel gültig:

${\frac {\sqrt {2}}{\sqrt {1-x^{4}}}}{\biggl \{}2E{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-2E{\biggl [}\arccos(x);{\frac {1}{2}}{\sqrt {2}}{\biggr ]}+F{\biggl [}\arccos(x);{\frac {1}{2}}{\sqrt {2}}{\biggr ]}{\biggr \}}\,+$

$+\,{\frac {{\sqrt {2}}\,x^{2}}{\sqrt {1-x^{4}}}}{\biggl \{}K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-F{\biggl [}\arccos(x);{\frac {1}{2}}{\sqrt {2}}{\biggr ]}{\biggr \}}=\int _{0}^{1}{\frac {2\,x^{3}(y^{2}+1)}{\sqrt {(1-x^{4})(1-x^{4}\,y^{4})}}}\,\mathrm {d} y$

Durch die Bildung der Ursprungsstammfunktion bezüglich x von der nun gezeigten Funktion entsteht diese Formel:

${\biggl \{}K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-F{\biggl [}\arccos(x);{\frac {1}{2}}{\sqrt {2}}{\biggr ]}{\biggr \}}{\biggl \{}2E{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-2E{\biggl [}\arccos(x);{\frac {1}{2}}{\sqrt {2}}{\biggr ]}+F{\biggl [}\arccos(x);{\frac {1}{2}}{\sqrt {2}}{\biggr ]}{\biggr \}}=$

$=\int _{0}^{1}{\frac {1}{y^{2}}}(y^{2}+1){\biggl [}{\text{artanh}}(y^{2})-{\text{artanh}}{\bigl (}{\frac {{\sqrt {1-x^{4}}}\,y^{2}}{\sqrt {1-x^{4}y^{4}}}}{\bigr )}{\biggr ]}\mathrm {d} y$

Nach der Regel von de L’Hospital gilt:

$\lim _{y\rightarrow 0}{\frac {1}{y}}(1-y^{2})\,{\text{artanh}}(y^{2})=0$

$\lim _{y\rightarrow 1}{\frac {1}{y}}(1-y^{2})\,{\text{artanh}}(y^{2})=0$

Wenn der Wert $x=1$ in die zuletzt genannte Integralidentität eingesetzt wird, dann entsteht folgende Identität:

$K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}{\biggl [}2\,E{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}{\biggr ]}=\int _{0}^{1}{\frac {1}{y^{2}}}(y^{2}+1)\,{\text{artanh}}(y^{2})\,\mathrm {d} y=$

$={\biggl [}2\arctan(y)-{\frac {1}{y}}(1-y^{2})\,{\text{artanh}}(y^{2}){\biggr ]}_{y=0}^{y=1}=2\arctan(1)={\frac {\pi }{2}}$

So kommt dann dieser Ausschnitt aus der Legendreschen Identität hervor:

$2E{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}^{2}={\frac {\pi }{2}}$

Verallgemeinerung für den nicht lemniskatischen Gesamtfall

Nach der soeben durchgeführten Herleitung gilt das genannte Resultat:

$2E{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}^{2}={\frac {\pi }{2}}$

Nun soll im Folgenden der moduläre Allgemeinfall bewiesen werden. Hierfür werden die Ableitungen der vollständigen elliptischen Integrale hergeleitet. Und im Anschluss wird die Ableitung der Legendreschen Identitätsbilanz ermittelt.

Beweis für die Ableitung des elliptischen Integrals erster Art:

${\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}K(\varepsilon )={\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}\int _{0}^{1}{\frac {1}{\sqrt {(1-x^{2})(1-\varepsilon ^{2}x^{2})}}}\mathrm {d} x=\int _{0}^{1}{\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}{\frac {1}{\sqrt {(1-x^{2})(1-\varepsilon ^{2}x^{2})}}}\mathrm {d} x=\int _{0}^{1}{\frac {\varepsilon x^{2}}{\sqrt {(1-x^{2})(1-\varepsilon ^{2}x^{2})^{3}}}}\mathrm {d} x=$

$=\int _{0}^{1}{\frac {\sqrt {1-\varepsilon ^{2}x^{2}}}{\varepsilon (1-\varepsilon ^{2}){\sqrt {(1-x^{2})}}}}\mathrm {d} x-\int _{0}^{1}{\frac {1}{\varepsilon {\sqrt {(1-x^{2})(1-\varepsilon ^{2}x^{2})}}}}\mathrm {d} x-\int _{0}^{1}{\frac {\varepsilon (1-2x^{2}+\varepsilon ^{2}x^{4})}{(1-\varepsilon ^{2}){\sqrt {(1-x^{2})(1-\varepsilon ^{2}x^{2})^{3}}}}}\mathrm {d} x=$

$={\frac {1}{\varepsilon (1-\varepsilon ^{2})}}E(\varepsilon )-{\frac {1}{\varepsilon }}K(\varepsilon )-\int _{0}^{1}{\frac {\mathrm {d} }{\mathrm {d} x}}{\frac {\varepsilon x{\sqrt {1-x^{2}}}}{(1-\varepsilon ^{2}){\sqrt {1-\varepsilon ^{2}x^{2}}}}}\mathrm {d} x={\frac {1}{\varepsilon (1-\varepsilon ^{2})}}{\bigl [}E(\varepsilon )-(1-\varepsilon ^{2})K(\varepsilon ){\bigr ]}$

Beweis für die Ableitung des elliptischen Integrals zweiter Art:

${\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}E(\varepsilon )={\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}\int _{0}^{1}{\frac {\sqrt {1-\varepsilon ^{2}x^{2}}}{\sqrt {1-x^{2}}}}\mathrm {d} x=\int _{0}^{1}{\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}{\frac {\sqrt {1-\varepsilon ^{2}x^{2}}}{\sqrt {1-x^{2}}}}\mathrm {d} x=\int _{0}^{1}{\frac {-\varepsilon x^{2}}{\sqrt {(1-x^{2})(1-\varepsilon ^{2}x^{2})}}}\mathrm {d} x=$

$=-\int _{0}^{1}{\frac {1}{\varepsilon {\sqrt {(1-x^{2})(1-\varepsilon ^{2}x^{2})}}}}\mathrm {d} x+\int _{0}^{1}{\frac {\sqrt {1-\varepsilon ^{2}x^{2}}}{\varepsilon {\sqrt {(1-x^{2})}}}}\mathrm {d} x=-{\frac {1}{\varepsilon }}{\bigl [}K(\varepsilon )-E(\varepsilon ){\bigr ]}$

Für die Pythagoräischen Gegenmodule gilt mit der Kettenregel dann:

${\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}K({\sqrt {1-\varepsilon ^{2}}})={\frac {1}{\varepsilon (1-\varepsilon ^{2})}}{\bigl [}\varepsilon ^{2}K({\sqrt {1-\varepsilon ^{2}}})-E({\sqrt {1-\varepsilon ^{2}}}){\bigr ]}$

${\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}E({\sqrt {1-\varepsilon ^{2}}})={\frac {\varepsilon }{1-\varepsilon ^{2}}}{\bigl [}K({\sqrt {1-\varepsilon ^{2}}})-E({\sqrt {1-\varepsilon ^{2}}}){\bigr ]}$

Denn die Ableitung der Kreisfunktion ist das negative Produkt aus der identischer Abbildungsfunktion und dem Kehrwert der Kreisfunktion. Die Legendresche Identität beinhaltet Produkte von jeweils zwei vollständigen elliptischen Integralen. Für die Ableitung der Funktionsseite von der Gleichungswaage der Legendreschen Identität wird die Produktregel im nun Folgenden angewendet:

${\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}K(\varepsilon )E({\sqrt {1-\varepsilon ^{2}}})={\frac {1}{\varepsilon (1-\varepsilon ^{2})}}{\bigl [}E(\varepsilon )E({\sqrt {1-\varepsilon ^{2}}})-K(\varepsilon )E({\sqrt {1-\varepsilon ^{2}}})+\varepsilon ^{2}K(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}}){\bigr ]}$

${\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}E(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}})={\frac {1}{\varepsilon (1-\varepsilon ^{2})}}{\bigl [}-E(\varepsilon )E({\sqrt {1-\varepsilon ^{2}}})+E(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}})-(1-\varepsilon ^{2})K(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}}){\bigr ]}$

${\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}K(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}})={\frac {1}{\varepsilon (1-\varepsilon ^{2})}}{\bigl [}E(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}})-K(\varepsilon )E({\sqrt {1-\varepsilon ^{2}}})-(1-2\varepsilon ^{2})K(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}}){\bigr ]}$

Wenn von diesen drei Gleichungen die beiden oberen Gleichungen addiert werden und die unterste Gleichung subtrahiert wird, dann entsteht dieses Resultat:

${\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}{\bigl [}K(\varepsilon )E({\sqrt {1-\varepsilon ^{2}}})+E(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}})-K(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}}){\bigr ]}=0$

Bezüglich ε ergibt die Bilanz konstant den Wert Null.

Für den Modul ε = 1/sqrt(2) gilt das zuvor ermittelte Resultat:

$2E{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}^{2}={\frac {\pi }{2}}$

Die Kombination der beiden zuletzt genannten Formeln ruft folgendes Ergebnis hervor:

$K(\varepsilon )E({\sqrt {1-\varepsilon ^{2}}})+E(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}})-K(\varepsilon )K({\sqrt {1-\varepsilon ^{2}}})={\tfrac {1}{2}}\pi$

Denn wenn die Ableitung einer kontinuierlichen Funktion konstant den Wert Null annimmt, dann ist die betroffene Funktion eine konstante Funktion. Das bedeutet, dass diese Funktion für jeden Abszissenwert ε den gleichen Funktionswert ergibt und der zugehörige Funktionsgraph somit eine waagrechte Gerade ist.

Anwendung

Reihe für den Kehrwert der Kreiszahl

Gültig sind diese Maclaurinschen Reihen für alle reellen Werte −1 < ε < 1:

$K(\varepsilon )={\frac {\pi }{2}}\sum _{k=0}^{\infty }{\frac {1}{16^{k}}}{\binom {2k}{k}}^{2}\varepsilon ^{2k}$

$E(\varepsilon )={\frac {\pi }{2}}\sum _{k=0}^{\infty }{\frac {1}{16^{k}(1-2k)}}{\binom {2k}{k}}^{2}\varepsilon ^{2k}$

Deswegen gilt auch jenes Formelpaar:

$K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}={\frac {\pi }{2}}\sum _{k=0}^{\infty }{\frac {1}{32^{k}}}{\binom {2k}{k}}^{2}$

$E{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}={\frac {\pi }{2}}\sum _{k=0}^{\infty }{\frac {1}{32^{k}(1-2k)}}{\binom {2k}{k}}^{2}$

Diese beiden Formeln können in jene Formel eingesetzt werden:

$2E{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}-K{\bigl (}{\frac {1}{2}}{\sqrt {2}}{\bigr )}^{2}={\frac {\pi }{2}}$

Dann kann folgende Reihenentwicklung synthetisiert werden:

${\biggl [}\sum _{k=0}^{\infty }{\frac {1}{32^{k}}}{\binom {2k}{k}}^{2}{\biggr ]}{\biggl [}\sum _{k=0}^{\infty }{\frac {1+2k}{32^{k}(1-2k)}}{\binom {2k}{k}}^{2}{\biggr ]}={\frac {2}{\pi }}$

Die Konvergenzgeschwindigkeit für diese Reihenformel verhält sich bezüglich der Nachkommastellen linear:

Obergrenze vom Index	Wert der Summe	Dezimale Nachkommastellen
0	1	1
1	45/64	0,70312500
2	43065/65536	0,65711975
3	2701125/4194304	0,64399838
4	43945661025/68719476736	0,63949353
5	2805051005757/4398046511104	0,63779475

Die Nachkommastellenresultate wurden durch Abrundung hervorgerufen.

Der Bruch 2/π hat die folgenden ersten dezimalen Nachkommastellen:

${\frac {2}{\pi }}\approx 0{,}636619772367581343$

Ableitung vom elliptischen Nomen

Das elliptische Nomen ist so definiert:

$q(x)=\exp {\bigl [}-\pi \,K({\sqrt {1-x^{2}}})\,K(x)^{-1}{\bigr ]}$

Für die Ableitung des vollständigen elliptischen Integrals erster Art gilt diese Formel:

${\frac {\mathrm {d} }{\mathrm {d} x}}K(x)={\frac {1}{x(1-x^{2})}}{\bigl [}E(x)-(1-x^{2})K(x){\bigr ]}$

Mit der Kettenregel und der Quotientenregel kann dann die Ableitung des elliptischen Nomens ermittelt werden:

${\frac {\mathrm {d} }{\mathrm {d} x}}q(x)={\frac {\mathrm {d} }{\mathrm {d} x}}\exp {\bigl [}-\pi K({\sqrt {1-x^{2}}})K(x)^{-1}{\bigr ]}={\biggl [}{\frac {\mathrm {d} }{\mathrm {d} x}}\,{\frac {-\pi K({\sqrt {1-x^{2}}})}{K(x)}}{\biggr ]}\exp {\bigl [}-\pi K({\sqrt {1-x^{2}}})K(x)^{-1}{\bigr ]}=$

$={\biggl [}{\frac {\mathrm {d} }{\mathrm {d} x}}\,{\frac {-\pi K({\sqrt {1-x^{2}}})}{K(x)}}{\biggr ]}q(x)={\frac {\pi }{x(1-x^{2})K(x)^{2}}}{\bigl [}K(x)E({\sqrt {1-x^{2}}})+E(x)K({\sqrt {1-x^{2}}})-K(x)K({\sqrt {1-x^{2}}}){\bigr ]}q(x)=$

$={\frac {\pi ^{2}}{2x(1-x^{2})K(x)^{2}}}q(x)$

Das elliptische Nomen stellt die Beziehung zwischen der Jacobischen Thetafunktion und dem vollständigen elliptischen Integral erster Art her:

$\vartheta _{00}[q(x)]=\prod _{n=1}^{\infty }{\bigl [}1-q(x)^{2n}{\bigr ]}{\bigl [}1+q(x)^{2n-1}{\bigr ]}^{2}={\sqrt {2\pi ^{-1}K(x)}}$

$\vartheta _{01}[q(x)]=\prod _{n=1}^{\infty }{\bigl [}1-q(x)^{2n}{\bigr ]}{\bigl [}1-q(x)^{2n-1}{\bigr ]}^{2}={\sqrt[{4}]{1-x^{2}}}{\sqrt {2\pi ^{-1}K(x)}}$

Literatur

• Duren, Peter (1991), "The Legendre relation for elliptic integrals", in Ewing, John H.; Gehring, F. W. (eds.), Paul Halmos. Celebrating 50 years of mathematics, New York: Springer-Verlag, pp. 305-315, doi:10.1007/978-1-4612-0967-6_32, ISBN 0-387-97509-8, MR 1113282

• Karatsuba, E. A.; Vuorinen, M. (2001), "On hypergeometric functions and generalizations of Legendre's relation", J. Math. Anal. Appl., 260 (2): 623–640, MR 1845572

• Legendre, A.M. (1811), Exercices de calcul intégral sur divers ordres de transcendantes et sur les quadratures, vol. I, Paris

• Legendre, A.M. (1825), Traité des fonctions elliptiques et des intégrales eulériennes, vol. I, Paris