Lambert W function

From Citizendium
Revision as of 18:27, 26 April 2007 by imported>Fredrik Johansson (→‎Calculus)
Jump to navigation Jump to search

The Lambert W function is used in mathematics to solve equations in which the unknown appears both as a linear term and in the argument of an exponential function or a logarithm, such as Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 3x+2=e^x} or Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x=\ln(4x)} . Such equations cannot otherwise, except in special cases, be solved explicitly in terms of algebraic operations, exponentials and logarithms.

Definition

The Lambert W function is defined as the multivalued function Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle W} that satisfies

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z = W(z)e^{W(z)}\;}

for any complex number Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle z} . Equivalently, it may be defined as the inverse function of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(w) = we^w} . An equation can be solved by rewriting it in the canonical form Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle C = we^w} , after which the solution is given by Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle w = W(C)} . The multivaluedness of the Lambert W function means that such an equation may have more than one solution. The graph of the Lambert W function in the real numbers looks as follows:

Lambertw.png

The function has two real branches in the interval which join at Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x = -1/e} . Concretely, this means that the equation Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x = w e^w} has two real solutions if Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle -1/e < x < 0} . For example, if Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x = -0.15} (which is about half way between Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle -1/e} and 0), there is one solution Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle w_0 \approx -0.17} that lies on the blue graph and another solution Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle w_1 \approx -3.0} that lies on the dashed red graph.

The single-valued function corresponding to the blue graph for Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle -1/e\, \le \, x \,< \,\infty} is called the principal branch of the Lambert W function and is denoted by Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle W_0} . The single-valued function corresponding to the dashed red graph for Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle -1/e \, \le \, x \, < \, 0} is called the negative branch, denoted by Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle W_{-1}} . The negative branch asymptotically approaches Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle -\infty} as Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle x \to 0} while the principal branch grows slowly but unboundedly (asymptotically like Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle \ln x} ) as Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle x \to \infty} .

Besides Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle W(-1/e) = -1} , the Lambert W function has the special values Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle W_0(0) = 0} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle W_0(e) = 1} . The value Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \scriptstyle W_0(1) = 0.567143\ldots} is called the omega constant.

Calculus

The principal branch of the Lambert W function has the Taylor series expansion

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle W_0(z) = \sum_{n=1}^\infty \frac{(-n)^{n-1}}{n!}z^n = z - z^2 + \begin{matrix}\frac{3}{2}\end{matrix}z^3 - \begin{matrix}\frac{8}{3}\end{matrix}z^4 + \begin{matrix}\frac{125}{24}\end{matrix}z^5 - \cdots}

around Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z=0} , which can be obtained the via Lagrange inversion theorem. Due to the singularity at Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z=-1/e} , the series converges for Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |z|<1/e} . The Lambert W function has the derivative

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle W'(x) = \frac{1}{(1+W(x))\exp(W(x))}}

and the antiderivative

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \int W(x) \,dx = x\left(W(x)-1+\frac{1}{W(x)}\right)+C.}

It is also possible to find antiderivatives of more complicated expressions containing the Lambert W function. This matter is discussed further in Corless et al.[1]

History and application

Equations of the kind that can be solved analytically with the Lambert W function are common in mathematics and science, yet the utility of such a function was not realized until recently. The Lambert W function was introduced in the 1980's as a function in the Maple computer algebra system, whose interface required an explicit notation for solutions of equations. The function's history highlights the importance of good mathematical notation: due to previously not being recognized as a function in its own right, it had not been studied systematically, despite its most important properties requiring only elementary complex analysis. An account of the function and its history that helped popularize it is given in a 1996 paper by R. M. Corless et al. (with Donald Knuth a notable co-author).[1]

The basic theory behind the Lambert W function was investigated in 1779 by Leonhard Euler.[2] The Maple developers chose the name of Johann Heinrich Lambert instead of Euler's since Euler had referenced work by Lambert in his paper, and possibly because "naming yet another function after Euler would not be useful".[3]

Since its introduction, the Lambert W function has been applied to problems ranging from quantum physics to the complexity of algorithms. Cranmer[4] discusses the application of the Lambert W function in solar wind physics and writes in the conclusion: "The Lambert W function used in these solutions was defined and publicized only about a decade ago, but it has rapidly become a convenient tool for mathematical physicists. The elegance of explicit solutions to equations thought previously to be expressible only implicitly is clear, but there also are many practical benefits to having explicit solutions as well."

See also

References

  1. 1.0 1.1 Corless, R. M., Gonnet, G. H., Hare, D. E. G., Jeffrey, D. J. & Knuth, D. E. (1996). "On the Lambert W function". Adv. Computational Maths. 5, 329–359
  2. Euler, L. (1779). "De serie Lambertina plurimisque eius insignibus proprietatibus". Originally published in Acta Academiae Scientarum Imperialis Petropolitinae 1779, 1783, pp. 29–51. Also in Opera Omnia: Series 1, Volume 6, pp. 350 - 369. See E532 in The Euler Archive for a scanned copy.
  3. Corless, R., Jeffrey, D. & Knuth, D. E. (1997). "A Sequence of Series for the Lambert W Function", ISSAC: Proceedings of the ACM SIGSAM International Symposium on Symbolic and Algebraic Computation (formerly SYMSAM, SYMSAC, EUROSAM, EUROCAL) (also sometimes in cooperation with the Symbolic and Algebraic Manipulation Groupe in Europe (SAME))
  4. Cranmer, S. R. (2004). "New views of the solar wind with the Lambert W function". Am. J. Phys. 72, 1397