Lambert W-Function
The Lambert π W
-function,
also called the omega function, is the inverse function
of
| π f(W)=We^W. |
(1)
|
The plot above shows the function along the real axis. The principal value of the Lambert π W
-function is implemented in the Wolfram
Language as [z].
Different branches of the function are available in the Wolfram
Language as [k,
z], where π k
is any integer and π k=0
corresponds to the principal value. Although undocumented, [k,
z] autoevaluates to [k,
z] in the Wolfram Language.
Lambert (1758) considered the solution to
now known as Lambert's transcendental equation. Euler heard about Lambert's paper in 1764 when Lambert moved from Zurich
to Berlin. After some private disputes about the priorities of some related series
expansions in 1770/1771, Euler (1783) wrote a paper about Lambert's
transcendental equation in which he introduced a special case which reduces to
π wa^w=lx
,
which is nearly the definition of π W(x)
, although Euler proposed defining a function more like
π -W(-x)
.
Euler considered series solutions in this paper and, in the first paragraph, explicitly
quotes Lambert as the one who first considered this equation.
Eisenstein (1844) considered the series of the infinite power tower
which can be expressed in closed form as
PΓ³lya and SzegΓΆ (1925) were the first to use the symbol π W
for the Lambert function.
Banwell and Jayakumar (2000) showed that a π W
-function describes the relation between voltage, current,
and resistance in a diode, and Packel and Yuen (2004) applied the π W
-function to a ballistic projectile in the presence of air
resistance. Other applications have been discovered in statistical mechanics, quantum
chemistry, combinatorics, enzyme kinetics, the physiology of vision, the engineering
of thin films, hydrology, and the analysis of algorithms (Hayes 2005).
The Lambert π W
-function
is illustrated above in the complex plane.
The real (left) and imaginary (right) parts of the analytic continuation of π W(z)
over the complex plane are illustrated
above (M. Trott, pers. comm.).
π W(x)
is real for π x>=-1/e
.
It has the special values
| π W(-1/2pi) | π = | π 1/2ipi |
(5)
|
| π W(-e^(-1)) | π = | π -1 |
(6)
|
| π W(0) | π = | π 0 |
(7)
|
| π W(1) | π = | π 0.567143.... |
(8)
|
π W(1)=0.567143...
(OEIS A030178) is called the omega
constant and can be considered a sort of "golden
ratio" of exponentials since
giving
The Lambert π W
-function
obeys the identity
(pers. comm. from R. Corless to O. Marichev, Sep. 29, 2015).
The function π W(ze^z)/z
has a very complicated structure in the complex plane, but is simply equal to 1 for
π R[z]>=1
and a slightly extended region above and below the real axis.
The Lambert π W
-function
has the series expansion
| π W(x) | π = | π sum_(n=1)^(infty)((-1)^(n-1)n^(n-2))/((n-1)!)x^n |
(12)
|
| π Image | π = | π x-x^2+3/2x^3-8/3x^4+(125)/(24)x^5-(54)/5x^6+(16807)/(720)x^7+.... |
(13)
|
The Lagrange inversion theorem gives the equivalent series expansion
where π n!
is a factorial. However, this series oscillates between
ever larger positive and negative
values for real π z>~0.4
, and so cannot be used for practical numerical computation.
An asymptotic formula which yields reasonably accurate results for π z>~3
is
where
(Corless et al. 1996), correcting a typographical error in de Bruijn (1981). Another expansion due to Gosper (pers. comm., July 22, 1996) is the double series
| π W(x)=a+sum_(n=0)^infty{sum_(k=0)^n(S_1(n,k))/([ln(x/a)-a]^(k-1)(n-k+1)!)}[1-(ln(x/a))/a]^n, |
(19)
|
where π S_1
is a nonnegative Stirling number
of the first kind and π a
is a first approximation which can be used to select between
branches. The Lambert π W
-function is two-valued for π -1/e<=x<0
. For π W(x)>=-1
, the function is denoted π W_0(x)
or simply π W(x)
, and this is called the principal
branch. For π W(x)<=-1
,
the function is denoted π W_(-1)(x)
. The derivative of
π W
is
for π x!=0
.
For the principal branch when π z>0
,
The π n
th
derivatives of the Lambert π W
-function are given by
where π a_(kn)
is the number triangle
(OEIS A042977). This has exponential generating function
See also
Abel Polynomial, Digit-Shifting Constants, Lambert's Transcendental Equation, Omega Constant, Power TowerRelated Wolfram sites
http://functions.wolfram.com/ElementaryFunctions/ProductLog/Explore with Wolfram|Alpha
More things to try:
References
Banwell, T. C. and Jayakumar, A. "Exact Analytical Solution for Current Flow Through Diode with Series Resistance." Electronics Lett. 36, 291-292, 2000.Barry, D. J., Culligen-Hensley, P. J.; and Barry, S. J. "Real Values of the π WFunction." ACM Trans. Math. Software 21, 161-171, 1995.Borwein, J. M. and Corless, R. M. "Emerging Tools for Experimental Mathematics." Amer. Math. Monthly 106, 899-909, 1999.Briggs, K. "π W
-ology, or, Some Exactly Solvable Growth Models." http://keithbriggs.info/W-ology.html.Briggs, K. "Graph Theory and Lambert's π W
function." http://keithbriggs.info/graph_theory_and_W.html.Corless, R. M. "Material on the Lambert π W
Function." (Ed. W. W. KΓΌchlin). New York: ACM, pp. 197-204, 1997. http://www.apmaths.uwo.ca/~rcorless/frames/PAPERS/LambertW/.Corless, R. M.; Gonnet, G. H.; Hare, D. E. G.; Jeffrey, D. J.; and Knuth, D. E. "On the Lambert π W
Function." Adv. Comput. Math. 5, 329-359, 1996.Corless, R. M.; Gonnet, G. H.; Hare, D. E. G.; and Jeffrey, D. J. "Lambert's π W
Function in Maple." Maple Technical Newsletter 9, 12-22, Spring 1993.Corless, R. M. and Jeffrey, D. J. "The Wright π omega
Function." In Artificial Intelligence, Automated Reasoning, and Symbolic Computation (Ed. J. Calmet, B. Benhamou, O. Caprotti, L. Henocque and V. Sorge). Berlin: Springer-Verlag, pp. 76-89, 2002.Corless, R. M.; Jeffrey, D. J.; and Knuth, D. E. "A Sequence of Series for the Lambert π W
Function." In Proceedings of the 1997 International Symposium on Symbolic and Algebraic Computation, Maui, Hawaii. New York: ACM Press, pp. 197-204, 1997.de Bruijn, N. G. Asymptotic Methods in Analysis. New York: Dover, pp. 27-28, 1981.Eisenstein, G. "Entwicklung von π alpha^(alpha^(alpha^...))
." J. reine angew. Math. 28, 49-52, 1844.Euler, L. "De serie Lambertina Plurimisque eius insignibus proprietatibus." Acta Acad. Scient. Petropol. 2, 29-51, 1783. Reprinted in Euler, L. Opera Omnia, Series Prima, Vol. 6: Commentationes Algebraicae. Leipzig, Germany: Teubner, pp. 350-369, 1921.Fritsch, F. N.; Shafer, R. E.; and Crowley, W. P. "Algorithm 443: Solution of the Transcendental Equation π we^w=x
." Comm. ACM 16, 123-124, 1973.Gosper, R. W. Jr. "The Solutions of π ye^(y^2)=x
and π ye^y=x
." ACM SIGSAM Bull. 32, 8-10, 1998.Gosper, R. W. "Re: 2nd Order Eulerians." posting, July 22, 1996.Gray, J. J. and Tilling, L. "Johann Heinrich Lambert, Mathematician and Scientist 1728-1777." Historia Math. 5, 13-41, 1978.Hayes, B. "Why π W
?" Amer. Sci. 93, 104-108, 2005.Jeffrey, D. J.; Hare, D. E. G.; and Corless, R. M. "Unwinding the Branches of the Lambert π W
Function." Math. Scientist 21, 1-7, 1996.Jeffrey, D. J.; Hare, D. E. G.; and Corless, R. M. "Exact Rational Solutions of a Transcendental Equation." C. R. Math. Acad. Sci. Canada 20, 71-76, 1998.Jeffrey, D. J.; Corless, R. M.; Hare, D. E. G.; and Knuth, D. E. "Sur l'inversion de π y^ae^y
au moyen des nombres de Stirling associes." Comptes Rendus Acad. Sci. Paris 320, 1449-1452, 1995.Kalman, D. "A Generalized Logarithm for Exponential-Linear Equations." College Math. J. 32, 2-14, 2001.Lambert, J. H. "Observations variae in Mathesin Puram." Acta Helvitica, physico-mathematico-anatomico-botanico-medica 3, 128-168, 1758.ORCCA. "The Lambert W Function." http://www.orcca.on.ca/LambertW.Packel, E. and Yuen, D. "Projectile Motion with Resistance and the Lambert π W
Function." College Math. J. 35, 337-350, 2004.PΓ³lya, G. and SzegΓΆ, G. Aufgaben und LehrsΓ€tze der Analysis. Berlin, 1925. Reprinted as Problems and Theorems in Analysis I. Berlin: Springer-Verlag, 1998.Sloane, N. J. A. Sequences A030178 and A042977 in "The On-Line Encyclopedia of Integer Sequences.""Time for a New Elementary Function?" FOCUS: Newsletter Math. Assoc. Amer. 20, 2, Feb. 2000.Valluri, S. R.; Jeffrey, D. J.; and Corless, R. M. "Some Applications of the Lambert π W
Function to Physics." Canad. J. Phys. 78, 823-831, 2000.Wright, E. M. "Solution of the Equation π ze^z=a
." Bull. Amer. Math. Soc. 65, 89-93, 1959.
Referenced on Wolfram|Alpha
Lambert W-FunctionCite this as:
Weisstein, Eric W. "Lambert W-Function." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/LambertW-Function.html
