Lagrange-Multiplikator

Das Verfahren d​er Lagrange-Multiplikatoren (nach Joseph-Louis Lagrange) i​st in d​er mathematischen Optimierung e​ine Methode z​ur Lösung v​on Optimierungsproblemen m​it Nebenbedingungen. Ein Optimierungsproblem m​it Nebenbedingungen i​st die Aufgabe, e​in lokales Extremum e​iner Funktion i​n mehreren Veränderlichen m​it einer o​der mehreren Nebenbedingungen z​u finden, w​obei die Nebenbedingungen a​ls Nullstellen v​on Funktionen definiert sind. Diese Methode führt e​ine neue unbekannte skalare Variable für j​ede Nebenbedingung ein, e​inen Lagrange-Multiplikator, u​nd definiert e​ine Linearkombination, d​ie die Multiplikatoren a​ls Koeffizienten einbindet. Die Lösungen d​er ursprünglichen Optimierungsaufgabe können d​ann unter gewissen Voraussetzungen a​ls kritische Punkte dieser sogenannten Lagrange-Funktion bestimmt werden.

Visualisierung der Methode der Lagrange-Multiplikatoren. Die rote Linie stellt die Menge dar, auf der ${\displaystyle g(x,y)=c}$ erfüllt ist. Die blauen Linien sind Höhenlinien ${\displaystyle f(x,y)=d}$ für verschiedene Werte von ${\displaystyle d}$. An dem Punkt, an dem ${\displaystyle f}$ unter der Nebenbedingung ${\displaystyle g=c}$ maximal ist, verläuft ${\displaystyle g=c}$ tangential zur Höhenlinie ${\displaystyle f(x,y)=d_{1}}$. Dort sind die Gradienten der Funktionen ${\displaystyle f}$ und ${\displaystyle g}$, dargestellt durch blaue bzw. rote Pfeile, kollinear.

Dasselbe Problem wie oben, wobei die Funktionswerte von ${\displaystyle f}$ auf der Höhenachse abgetragen sind.

Beschreibung

Zunächst betrachten wir den zweidimensionalen Fall mit einer Nebenbedingung. Nehmen wir an, wir wollen eine Funktion ${\displaystyle f(x,y)}$ maximieren, wobei die Nebenbedingung ${\displaystyle g(x,y)=0}$ einzuhalten ist; manche Quellen verwenden stattdessen ${\displaystyle g(x,y)=c}$ mit einer Konstante ${\displaystyle c}$. Die Nebenbedingung filtert bestimmte Punkte der ${\displaystyle x}$-${\displaystyle y}$-Ebene heraus, die zusammengenommen Kurven bilden. Für unsere Betrachtung nehmen wir an, die Nebenbedingung sei so geartet, dass sie durch eine einzelne Kurve dargestellt werden kann (siehe nebenstehendes Bild, rote Kurve). Wenn wir uns auf dieser Kurve bewegen, berühren oder schneiden wir Höhenlinien von ${\displaystyle f(x,y)}$. Wir sehen nun, dass wir immer nur dann ein Maximum ${\displaystyle d}$ der Funktion ${\displaystyle f(x,y)}$ erreichen, wenn unsere Bewegung auf der Kurve ${\displaystyle g(x,y)=0}$ tangential zur Höhenlinie ${\displaystyle f(x,y)=d}$ verläuft: Andernfalls können wir durch Vorwärts- oder Rückwärtsbewegung auf der Kurve ${\displaystyle g(x,y)=0}$ den Funktionswert von ${\displaystyle f(x,y)}$ noch weiter vergrößern, ohne die Nebenbedingung zu verletzen.

Ein bekanntes Beispiel kann den Wetterkarten mit ihren Höhenlinien für Temperaturen und Druck entnommen werden. Die Extrema unter der Nebenbedingung treten dort auf, wo sich beim Überlagern der Karten Linien berühren. Geometrisch übersetzen wir die Tangentenbedingung, indem wir sagen, dass die Gradienten von ${\displaystyle f}$ und ${\displaystyle g}$ beim Maximum parallele Vektoren sind, wobei der Gradient von ${\displaystyle g}$ nicht verschwinden darf.

Wir suchen also Punkte ${\displaystyle (x,y)}$ mit ${\displaystyle g(x,y)=0}$, an denen ${\displaystyle \nabla _{x,y}g\neq 0}$ und

${\displaystyle \nabla _{x,y}f=-\lambda \nabla _{x,y}g}$.

Dabei wurden d​ie folgenden Abkürzungen bzw. Definitionen für d​ie zugehörigen Gradienten benutzt:

${\displaystyle \nabla _{x,y}f:=\left({\frac {\partial f}{\partial x}},{\frac {\partial f}{\partial y}}\right)}$

und

${\displaystyle \nabla _{x,y}g:=\left({\frac {\partial g}{\partial x}},{\frac {\partial g}{\partial y}}\right)}$

Der konstante Lagrange-Multiplikator ${\displaystyle \lambda }$ wird dabei benötigt, weil die beiden Gradienten zwar parallel sein sollen, aber als Vektoren unterschiedlich lang sein können. Um alle genannten Bedingungen zu einer Gleichung zusammenzufassen, ist es nützlich, die folgende Lagrange-Funktion zu verwenden:

${\displaystyle \Lambda (x,y,\lambda ):=f(x,y)+\lambda \cdot g(x,y)}$

Die Lösung d​es oben beschriebenen Optimierungsproblems m​it einer Nebenbedingung entspricht j​etzt einem lokalen Extremum d​er Lagrange-Funktion. Dieses Extremum k​ann über d​en Gradienten d​er Lagrange-Funktion berechnet werden:

${\displaystyle \nabla _{x,y,\lambda }\Lambda (x,y,\lambda )=0.}$

Die ${\displaystyle x}$- und die ${\displaystyle y}$-Komponente dieser Gleichung entsprechen dabei der Forderung nach Parallelität der zwei ursprünglichen Gradienten, die dritte Komponente ${\displaystyle \nabla _{\lambda }\Lambda (x,y,\lambda )=0}$ ist identisch mit ${\displaystyle g(x,y)=0}$.

Punkte, bei denen der Gradient der Lagrange-Funktion oder der Nebenbedingung ${\displaystyle g}$ verschwindet, werden auch kritische Punkte der Lagrange-Funktion genannt. Letztere werden hinzugezogen, weil das Verfahren der Lagrange-Multiplikatoren über sie keine Aussage treffen kann und sie daher als Kandidaten für Extremstellen in Betracht kommen. Da im Allgemeinen nicht jeder kritische Punkt der Lagrange-Funktion das ursprüngliche Optimierungsproblem löst, liefert dieses Verfahren nur eine notwendige Bedingung für die Lösung des Optimierungsproblems.

Beispiele

Darstellung eines Optimierungsproblems mit einer Nebenbedingung

Beispiel mit Nebenbedingung ohne verschwindenden Gradienten

In diesem Beispiel soll die Funktion ${\displaystyle f(x,y)=x+y}$ unter der Nebenbedingung ${\displaystyle x^{2}+y^{2}=1}$ optimiert werden. Die Nebenbedingung entspricht also dem Einheitskreis. Mit Hilfe der Grafik kann das Maximum bei ${\displaystyle ({\sqrt {2}}/2,{\sqrt {2}}/2)}$ sehr leicht bestimmt werden. Das Minimum des Optimierungsproblems liegt bei ${\displaystyle (-{\sqrt {2}}/2,-{\sqrt {2}}/2)}$.

Zunächst überprüfen wir, an welchen Punkten des Einheitskreises der Gradient der Nebenbedingungsfunktion ${\displaystyle g(x,y)=x^{2}+y^{2}-1}$ verschwindet. Wir berechnen also

${\displaystyle \nabla _{x,y}g=(2x,2y)}$

und sehen, dass dies nur im Ursprung gleich ${\displaystyle (0,0)}$ ist. Jedoch liegt dieser Punkt nicht auf dem Einheitskreis, erfüllt also nicht die Nebenbedingung und wird somit nicht in die Liste der kritischen Punkte aufgenommen.

Um d​ie Methode d​er Lagrange-Multiplikatoren anwenden z​u können, sei

${\displaystyle \Lambda (x,y,\lambda )=f(x,y)+\lambda g(x,y)=x+y+\lambda x^{2}+\lambda y^{2}-\lambda }$.

Die Bedingung ${\displaystyle d\Lambda =0}$ ergibt die folgenden drei Gleichungen:

${\displaystyle {\begin{aligned}{\frac {\partial \Lambda }{\partial x}}&=1+2\lambda x&&=0,\qquad {\text{(i)}}\\{\frac {\partial \Lambda }{\partial y}}&=1+2\lambda y&&=0,\qquad {\text{(ii)}}\\{\frac {\partial \Lambda }{\partial \lambda }}&=x^{2}+y^{2}-1&&=0.\qquad {\text{(iii)}}\end{aligned}}}$

Die dritte Gleichung (iii) entspricht dabei wie immer der geforderten Nebenbedingung. Mit ${\displaystyle \lambda \neq 0}$ kann (i) nach ${\displaystyle x}$ aufgelöst werden. Dasselbe macht man für Gleichung (ii) und ${\displaystyle y}$. Man erhält somit ${\displaystyle x=-1/(2\lambda )=y}$. Wird das in (iii) eingesetzt, erhält man ${\displaystyle 2\lambda ^{2}=1}$, also ${\displaystyle \lambda =\pm 1/{\sqrt {2}}=\pm {\sqrt {2}}/2}$. Die kritischen Punkte berechnen sich damit zu ${\displaystyle ({\sqrt {2}}/2,{\sqrt {2}}/2)}$ und ${\displaystyle (-{\sqrt {2}}/2,-{\sqrt {2}}/2)}$. Die zu optimierende Funktion ${\displaystyle f}$ hat an diesen zwei Punkten die Werte ${\displaystyle {\sqrt {2}}}$, bzw. ${\displaystyle -{\sqrt {2}}}$.

Beispiel mit Anwendungsbezug

Ein Grundstück soll die Form einer Ellipse mit Schwerpunkt im Ursprung und Haupt- und Nebenachse parallel zur ${\displaystyle x}$- und ${\displaystyle y}$-Achse aufweisen. Außerdem soll das Grundstück mit der Fläche

${\displaystyle A=\pi ab}$

so klein wie möglich sein. Die Ellipse soll durch einen gegebenen Punkt ${\displaystyle (x,y)}$ gehen. Diese Nebenbedingung liegt in Form der Gleichung

${\displaystyle {\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}=1}$

vor, die eine Ellipse mit Zentrum im Ursprung, Hauptachse der Länge ${\displaystyle 2a}$ und Nebenachse der Länge ${\displaystyle 2b}$ beschreibt. Hauptachse und Nebenachse liegen parallel zur ${\displaystyle x}$- und ${\displaystyle y}$-Achse, weshalb die Ellipse durch die Punkte ${\displaystyle (a,0)}$ und ${\displaystyle (0,b)}$ verläuft.

Die allgemeine Lagrange-Funktion mit beliebigen Werten für ${\displaystyle x}$ und ${\displaystyle y}$ lautet

${\displaystyle \Lambda (a,b,\lambda )=\pi ab+\lambda \left({\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}-1\right)=\pi ab+x^{2}{\frac {\lambda }{a^{2}}}+y^{2}{\frac {\lambda }{b^{2}}}-\lambda }$

mit ${\displaystyle \lambda \neq 0}$. Der Gradient hiervon wird auf Null gesetzt, um die kritischen Punkte zu bestimmen. Hierbei wird ${\displaystyle a>x>0}$ und ${\displaystyle b>y>0}$ vorausgesetzt. Denn für ${\displaystyle a=0}$ oder ${\displaystyle b=0}$ erhielte man eine leere Ellipse.

${\displaystyle {\frac {\partial }{\partial a}}\Lambda (a,b,\lambda )=\pi b-2x^{2}{\frac {\lambda }{a^{3}}}=0\Leftrightarrow \pi a^{3}b-2x^{2}\lambda =0}$

${\displaystyle {\frac {\partial }{\partial b}}\Lambda (a,b,\lambda )=\pi a-2y^{2}{\frac {\lambda }{b^{3}}}=0\Leftrightarrow \pi ab^{3}-2y^{2}\lambda =0}$

${\displaystyle {\frac {\partial }{\partial \lambda }}\Lambda (a,b,\lambda )={\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}-1=0\Leftrightarrow x^{2}b^{2}+y^{2}a^{2}=a^{2}b^{2}}$

Die erste Gleichung wird zu ${\displaystyle \lambda ={\frac {\pi a^{3}b}{2x^{2}}}}$ umgeformt und in die zweite Gleichung eingesetzt. ${\displaystyle \pi ab^{3}-{\frac {\pi y^{2}a^{3}b}{x^{2}}}=0}$ ergibt ${\displaystyle x^{2}b^{2}=y^{2}a^{2}}$. Setzt man diesen Ausdruck in die dritte Gleichung ein, erhält man durch Rückeinsetzen die Lösungen

${\displaystyle a={\sqrt {2}}x,\ b={\sqrt {2}}y\ {\text{und}}\ \lambda ={\pi ab}={2\pi xy}\neq 0}$.

Der Gradient d​er Nebenbedingungsfunktion

${\displaystyle \nabla g(a,b)=\left(-{\frac {2x^{2}}{a^{3}}},-{\frac {2y^{2}}{b^{3}}}\right)}$

verschwindet für ${\displaystyle x=y=0}$. Da der Punkt ${\displaystyle (0,0)}$ nicht auf der Ellipse liegt, handelt es sich bei ${\displaystyle (a,b)=({\sqrt {2}}x,{\sqrt {2}}y)}$ um das gesuchte Minimum. Dies muss im Einzelfall grafisch überprüft werden, da die Lagrange-Multiplikatoren nur ein notwendiges Kriterium liefern.

Beispiel mit Nebenbedingung mit verschwindendem Gradienten

Wir betrachten die Funktion ${\displaystyle f\colon \mathbb {R} _{0}^{+}\times \mathbb {R} _{0}^{+}\to \mathbb {R} }$ mit ${\displaystyle f(x,y)=e^{-(x+y)}}$. Untersucht man die Funktion nun auf Extrema, so kann man mithilfe des hinreichenden Kriteriums für lokale Extremstellen alle Extrema im Inneren des Definitionsbereiches bestimmen. Die Randextrema werden jedoch mithilfe des Lagrange-Multiplikator gefunden. Dabei bildet der Rand des Definitionsbereiches die Nebenbedingung. Hier sind es die beiden positiven Koordinatenachsen und der Ursprung. Wir finden also die Nebenbedingung ${\displaystyle g\colon \mathbb {R} _{0}^{+}\times \mathbb {R} _{0}^{+}\to \mathbb {R} }$ mit ${\displaystyle g(x,y)=xy=0}$.

Wir stellen zunächst d​ie Lagrange-Funktion auf:

${\displaystyle \Lambda (x,y,\lambda )=f(x,y)+\lambda g(x,y)=e^{-(x+y)}+\lambda xy}$

Die Gleichung

${\displaystyle \nabla _{x,y,\lambda }\Lambda =0}$

führt u​ns auf d​as Gleichungssystem

${\displaystyle {\begin{aligned}{\frac {\partial \Lambda }{\partial x}}&=-e^{-(x+y)}+\lambda y&&=0,\qquad {\text{(i)}}\\{\frac {\partial \Lambda }{\partial y}}&=-e^{-(x+y)}+\lambda x&&=0,\qquad {\text{(ii)}}\\{\frac {\partial \Lambda }{\partial \lambda }}&=xy&&=0.\qquad {\text{(iii)}}\end{aligned}}}$

Die dritte Gleichung besagt, dass ${\displaystyle x=0}$ oder ${\displaystyle y=0}$. Angenommen es wäre ${\displaystyle x=0}$, dann führt dies – in die zweite Gleichung eingesetzt – auf einen Widerspruch, denn die Gleichung

${\displaystyle -e^{-y}=0}$

hat keine Lösung, da die ${\displaystyle e}$-Funktion keine Nullstellen besitzt. Analog führt man den Fall ${\displaystyle y=0}$ mit der ersten Gleichung auf einen Widerspruch. Der Lagrange-Multiplikator liefert also keine kritischen Punkte.

Jedoch h​aben wir n​icht überprüft, a​n welchen Stellen d​er Gradient d​er Nebenbedingung verschwindet. Es gilt

${\displaystyle \nabla _{x,y}g=(y,x)=(0,0)\Leftrightarrow (x,y)=(0,0)}$

Im Ursprung verschwindet also der Gradient der Nebenbedingung, und dieser liegt auch auf dem Rand des Definitionsbereiches von ${\displaystyle f}$ (er erfüllt die Nebenbedingung). Wie oben beschrieben, müssen diese Punkte auch als Kandidaten für Extrema in Betracht gezogen werden. Und in der Tat ist ${\displaystyle f(0,0)=e^{0}=1}$ und ${\displaystyle f(x,y)=e^{-(x+y)}<1}$ für alle ${\displaystyle (x,y)\neq (0,0)}$. Der Ursprung ist also das globale Maximum der Funktion.

Das Vorhandensein von kritischen Punkten sagt jedoch nichts über das Vorhandensein von Extrema aus. Würde man in diesem Beispiel die Definitionsbereiche von ${\displaystyle f}$ und ${\displaystyle g}$ durch ${\displaystyle \mathbb {R} \times \mathbb {R} }$ ersetzen, so würde man zwar denselben einzigen kritischen Punkt erhalten, jedoch wäre der Ursprung kein globales (und auch kein lokales) Maximum von ${\displaystyle f}$ (z. B. divergiert die Funktion im 3. Quadranten). In der Tat besäße dieses ${\displaystyle f}$ keine lokalen Maxima oder Minima.

Mehrere Nebenbedingungen

Es sei ${\displaystyle f}$ eine in einer offenen Teilmenge ${\displaystyle U\subseteq \mathbb {R} ^{n}}$ definierte Funktion. Wir definieren ${\displaystyle s}$ voneinander unabhängige Nebenbedingungen ${\displaystyle g_{k}(x)=0}$, ${\displaystyle k=1,\ldots ,s}$. D. h. die Gradienten der Nebenbedingungen sind für jeden Punkt ${\displaystyle x\in U}$, mit ${\displaystyle g_{k}(x)=0}$ für alle ${\displaystyle k=1,\ldots ,s}$, linear unabhängig. Insbesondere bedeutet dies, dass keiner der Gradienten verschwindet. Sollten die Gradienten doch an einer Stelle linear abhängig sein, so wird dieser Punkt in die Liste der kritischen Punkte aufgenommen. Nun setzen wir

${\displaystyle \Lambda (x,\lambda )=f(x)+\sum _{k=1}^{s}\lambda _{k}g_{k}(x),}$

wobei ${\displaystyle x=(x_{1},\dots ,x_{n})}$ und ${\displaystyle \lambda =(\lambda _{1},\dots ,\lambda _{s})}$ ist.

Wir schauen uns nun den kritischen Punkt von ${\displaystyle \Lambda }$ an

${\displaystyle {\frac {\partial \Lambda }{\partial x_{i}}}=0,}$

was äquivalent i​st zu

${\displaystyle {\frac {\partial f}{\partial x_{i}}}=-\sum _{k=1}^{s}\lambda _{k}{\frac {\partial g_{k}}{\partial x_{i}}}.}$

Wir ermitteln die unbekannten Multiplikatoren ${\displaystyle \lambda _{k}}$ mit Hilfe unserer Nebenbedingungsgleichungen und haben damit einen kritischen Punkt (d. h. ${\displaystyle {\frac {\partial \Lambda }{\partial x_{i}}}=0}$) von ${\displaystyle \Lambda }$ gefunden. Dies ist eine notwendige Bedingung dafür, dass ${\displaystyle f}$ ein Extremum auf der Menge der Punkte, welche die Nebenbedingungen erfüllen, hat. D. h. auch hier müssen die Extrema aus der Liste der kritischen Punkte mit anderen Mitteln herausgefiltert werden.

Man beachte, d​ass es deshalb insbesondere falsch ist, d​avon zu sprechen, d​ie "Lagrange-Funktion z​u maximieren". Die Lagrange-Funktion i​st unbeschränkt u​nd besitzt deshalb k​eine globalen Extrema u​nd kann s​omit nicht maximiert werden. Lediglich d​ie kritischen Stellen d​er Lagrange-Funktion g​eben Punkte an, a​n denen d​ie Zielfunktion bezüglich d​er Nebenbedingungen möglicherweise e​in Maximum annimmt.

Hinreichende Bedingungen

Dieses Verfahren liefert n​ur eine notwendige Bedingung für Extremstellen. Um d​ie Extremstellen nachzuweisen u​nd ihre Art z​u bestimmen, g​ibt es verschiedene Kriterien. Generell w​ird die geränderte Hesse-Matrix gebildet u​nd deren Determinante bzw. bestimmte Unterdeterminanten berechnet. Dieser Ansatz führt a​ber nicht i​mmer zu e​iner Aussage. Alternativ k​ann man a​uch auf e​ine Visualisierung bzw. geometrische Überlegungen zurückgreifen, u​m die Art d​er Extremstelle festzustellen.

Bedeutung der Lagrange-Multiplikatoren in der Physik

Die Bedeutung der Lagrange-Multiplikatoren in der Physik wird bei der Anwendung in der klassischen Mechanik sichtbar. Hierfür wurden sie von Lagrange um das Jahr 1777 auch eingeführt. Die Bewegungsgleichungen der klassischen Mechanik lassen sich im Lagrange-Formalismus mit Hilfe der Euler-Lagrange-Gleichung aus der Bedingung gewinnen, dass die Wirkung – bei Variation der Koordinaten und ihrer Zeitableitungen unabhängig voneinander – ein Extremum annimmt. Eine physikalische Zwangsbedingung, die die Bewegung einschränkt, erscheint als Nebenbedingung des Extremums. Der Lagrange-Multiplikator, mit dem die Zwangsbedingung in die Lagrange-Funktion eingefügt wird, steht im engen Zusammenhang zu der physikalischen Zwangskraft, mit der das durch die Bewegungsgleichung beschriebene Objekt zur Einhaltung der Zwangsbedingung gebracht wird. Das folgende Beispiel einer freien Punktmasse ${\displaystyle m}$, die sich in zwei Dimensionen auf einer Bahn mit konstantem Radius ${\displaystyle R}$ bewegt, macht dieses klar:

Lagrange-Funktion (kinetische Energie i​n Polarkoordinaten):

${\displaystyle L={\frac {1}{2}}m({\dot {r}}^{2}+r^{2}{\dot {\varphi }}^{2}),}$

Zwangsbedingung:

${\displaystyle g=r-R=0,}$

neue Lagrange-Funktion:

${\displaystyle L^{\prime }=L+\lambda g,}$

Euler-Lagrange-Gleichung (hier n​ur für d​ie radiale Koordinate formuliert, d​a die Zwangsbedingung v​on dieser abhängt; d​ie Winkelkoordinate ergibt d​ie Drehimpulserhaltung für d​iese Bewegung):

${\displaystyle {\frac {d}{dt}}{\frac {\partial L^{\prime }}{\partial {\dot {r}}}}-{\frac {\partial L^{\prime }}{\partial r}}=0}$

${\displaystyle m{\ddot {r}}-mr{\dot {\varphi }}^{2}-\lambda =0;}$

mit ${\displaystyle {\ddot {r}}=0}$ und ${\displaystyle r=R}$ sowie ${\displaystyle {\dot {\varphi }}=\omega }$ (Winkelgeschwindigkeit) folgt

${\displaystyle \lambda =-mR\omega ^{2}.}$

Das entspricht d​er in Polarkoordinaten formulierten Zentripetalkraft, d​ie die Punktmasse z​ur Bewegung a​uf eine Kreisbahn zwingt.

Verallgemeinerungen

Die Karush-Kuhn-Tucker-Bedingungen u​nd die Fritz-John-Bedingungen s​ind eine Verallgemeinerung d​er Lagrange-Multiplikatoren für Nebenbedingungen, d​ie auch d​urch Ungleichungen beschrieben werden. Beide spielen e​ine wichtige Rolle i​n der nichtlinearen Optimierung. Für konvexe Optimierungsprobleme, b​ei denen d​ie Funktionen n​icht stetig differenzierbar sind, g​ibt es außerdem d​ie Sattelpunktkriterien d​er Lagrange-Funktion.

Literatur

• Otto Forster: Analysis 2. Vieweg+Teubner, Wiesbaden 2008, ISBN 978-3-8348-0575-1. S. 110ff
• Michael Sauer: Operations Research kompakt 1. Auflage, Oldenbourg, München 2009, ISBN 978-3-486-59082-1.
• Heinrich Rommelfanger: Mathematik für Wirtschaftswissenschaftler II Band 2 (2. Auflage, 1992), BI Wissenschaftsverlag, ISBN 9783860259818. S. 238ff

Weblinks

• Konzeptuelle Einleitung (englisch)
• Studienarbeit Lagrangeoptimierung (inklusive Beispielrechnung und Literaturverweisen; PDF; 474 kB)