Die Lambert-W-Funktion

Die Lambert-W-Funktion, 1758 von Johann Heinrich Lambert entdeckt und von Leonhard Euler erforscht, ist die Umkehrfunktion von f(w) = w*e^w. Formell: wenn W(x) = w, dann w*e^w = x. Sie ist keine elementare Funktion, erscheint aber in unzähligen Problemen der Physik, chemischen Kinetik, Biologie und Algorithmenanalyse.

Äste der Funktion

Da f(w) = w*e^w nicht monoton ist (Minimum bei w = -1 mit f(-1) = -1/e), ist die Umkehrung mehrwertig. Auf den reellen Zahlen gibt es zwei Äste:

• W0(x) (Hauptast): definiert für x >= -1/e. W0(-1/e) = -1, W0(0) = 0, W0(1) = Omega ~= 0.5671, W0(e) = 1.
• W-1(x) (Nebenast): definiert für -1/e <= x 0-.

Die Omega-Konstante

Omega = W0(1) = 0.5671432904... ist die eindeutige Lösung von Omega*e^Omega = 1. Sie ist transzendent (Seleshi 2015) und erscheint in analogen Kontexten wie e oder pi.

Halley-Verfahren

Der Rechner verwendet Halley-Iteration mit kubischer Konvergenz (verdreifacht korrekte Stellen pro Iteration). Gegeben f(w) = w*e^w - x:

w_{n+1} = w_n - f(w_n) / [f'(w_n) - f''(w_n) * f(w_n) / (2 * f'(w_n))]

Typischerweise reichen 4-7 Iterationen für 14 Dezimalstellen.

Fünf Modi

1. W0(x): Hauptast, x >= -1/e.
2. W-1(x): Nebenast, -1/e <= x < 0.
3. Löse x*e^x = y: x = W0(y).
4. Löse x*ln(x) = y: x = e^W(y).
5. Unendlicher Potenzturm x^x^x^...: konvergiert für x in [e^-e, e^(1/e)].

Reihen

Für kleines |x|: W(x) = x - x^2 + 3x^3/2 - 8x^4/3 + ..., Radius 1/e. Für großes x: W(x) = L1 - L2 + L2/L1 + ..., mit L1 = ln(x), L2 = ln(ln(x)).

Physikalische Anwendungen

• Statistische Mechanik: Boltzmann-Verteilungen mit nichtlinearen Nebenbedingungen.
• Enzymkinetik: exakte Lösung Michaelis-Menten.
• Wiensches Gesetz: Maximum des Planck-Spektrums.
• Dioden: explizite Shockley-Gleichung für I-V.
• Kosmologie: Radius der Hubble-Sphäre in LCDM-Modellen.

Geschichte der Lambert-W-Funktion

Die Funktion ist nach Johann Heinrich Lambert benannt, einem Mathematiker und Physiker des 18. Jahrhunderts, der 1758 mit dem Problem verbundene transzendente Gleichungen untersuchte. Leonhard Euler vertiefte ihre Eigenschaften wenige Jahre später. Lange blieb die Funktion ohne Standardnamen und einheitliche Notation und erschien vereinzelt in Problemen der Physik und Technik. Erst 1996 gab ihr eine berühmte Arbeit von Robert Corless, Gaston Gonnet, David Hare, David Jeffrey und Donald Knuth den endgültigen Namen "Lambert-W-Funktion" und die Notation W und erkannte ihre Bedeutung. Ihre Wiederentdeckung wurde durch das Aufkommen symbolischer Rechensysteme wie Maple und Mathematica vorangetrieben, die ihren Gebrauch praktisch machten: aus einer mathematischen Kuriosität wurde ein Standardwerkzeug, heute in allen wichtigen wissenschaftlichen Bibliotheken implementiert.

Warum sie keine elementare Funktion ist

Die W-Funktion ist als Umkehrung von f(w) = w·eʷ definiert: W·e^W = x zu lösen bedeutet, den unbekannten Exponenten zu finden, der zugleich als multiplikative Basis und im Exponential erscheint. Diese Art von Gleichung, bei der die Unbekannte gleichzeitig innerhalb und außerhalb eines Exponentials auftritt, besitzt keine Lösung in elementaren Funktionen (Polynome, Wurzeln, Exponential- und Logarithmusfunktionen und deren endliche Kombinationen). Genau deshalb verdient die W-Funktion einen eigenen Namen: so wie die Quadratwurzel x² = a und der Logarithmus eˣ = a löst, löst die Lambert-W eine Klasse sonst nicht auflösbarer Gleichungen. Ihre Berechnung erfordert iterative numerische Verfahren wie das von Halley, die dank guter, je Ast verschiedener Anfangsschätzungen rasch konvergieren.

Mathematische und informatische Anwendungen

Die W-Funktion erscheint in zahlreichen scheinbar unverbundenen Problemen. In der Mathematik löst sie Verzögerungsdifferentialgleichungen, liefert die geschlossene Form des unendlichen Potenzturms x^x^x^... und zählt beschriftete Bäume in der Kombinatorik (die Baumfunktion ist −W(−x)). In der Physik beschreibt sie das Wiensche Verschiebungsgesetz für die Schwarzkörperstrahlung, die Enzymkinetik und Potenziale in der Quantenmechanik. In der Informatik tritt sie in der Algorithmenanalyse auf: die Iterationszahl gewisser Verfahren, die erwartete Tiefe zufälliger Datenstrukturen und die Lösung von Rekurrenzen der Form T(n) im Zusammenhang mit n·log(n) lassen sich natürlich durch W ausdrücken. Ihre Allgegenwart macht sie zu einer der nützlichsten speziellen Funktionen der modernen angewandten Mathematik.

Aktualisiert im Juli 2026