Description

The LambertW function satisfies: LambertW(x) * exp(LambertW(x)) = x.

As the equation y× exp(y) = x has an infinite number of solutions y for each (non-zero) value of x, LambertW has an infinite number of branches. Exactly one of these branches is analytic at 0. In Maple this branch is referred to as the principal branch of LambertW, and is denoted by LambertW(x). The other branches all have a branch point at 0, and these branches are denoted in Maple by LambertW(k, x), where k is any non-zero integer. (The principal branch can also be referred to as LambertW(0, x)).

The principal branch and the pair of branches LambertW(-1, x) and LambertW(1, x) share an order 2 branch point at -exp(-1). The branch cut dividing these branches is the subset of the real line from -infinity to -exp(-1), and the values of the branches of LambertW on this branch cut are assigned using the rule of counter-clockwise continuity around the branch point. This means that Lambert W(x) is real-valued for x in the range -exp(-1).. infinity, while the image of -infinity.. -exp(-1) under LambertW(x) is the curve -y cot(y) + y*I, for y in 0.. Pi.

For all the branches other than the principal branch, the branch cut dividing them is the negative real axis. The branches are numbered up and down from the real axis (this is very similar to the way the branches of the logarithm are indexed by the multiple of 2*Pi*I which must be subtracted from the imaginary part to recover the principal branch). Again, the values of the branches of LambertW along the branch cut are determined by the rule of counter-clockwise continuity around the branch point at 0. Thus, the image of the negative real axis under the branch LambertW(k, x) is the curve -y cot(y) + y*I, for y in 2*k*Pi.. (2*k + 1)*Pi if k > 0 and y in (2*k + 1)*Pi. (2*k + 2)*Pi if k < -1. These curves, therefore, bound the ranges of the branches of LambertW, and in each case, the upper boundary of the region is included in the range of the corresponding branch.

The asymptotic behavior of LambertW at complex infinity and at 0 (for the non-principal branches) is given by

LambertW(k, x)~log(k, x) -log(log(k, x) +sum(c(m, n)*log(log(k, x))^(m+1)/log(k, x)^(m+n+1),

0 < = m, n < = infinity)

where log(x) denotes the principal branch of the logarithm, log(k, x) = log(x) + 2*k*Pi*I and the c(m, n) are constants independent of k. The expansion for LambertW(-1, x) is not valid for x tending to 0 along the negative real axis (the effect of the branch point at -exp(-1) must be considered), but holds otherwise.

Examples

> LambertW(0); LambertW(-exp(-1));

> LambertW(1.5+2.5*I);

> plot(LambertW(x), x=-1..5);

График 10.3. Функция Ламберта.