朗伯W函數

朗伯W函數（Lambert W function，又稱為歐米加函數或乘積對數），是$f(w)=we^w$的反函數，其中$e^w$是指數函數，$w$是任意複數。對於任何複數$z$，都有：

$$z = W(z)e^{W(z)}.$$

由於函數$f$不是單射，因此函數$W$是多值的（除了0以外）。如果我們把$x$限制為實數，並要求$w$是實數，那麼函數僅對於$x\ge\frac{1}{e}$有定義，在$(-\frac{1}{e},0)$內是多值的；如果加上$w\ge-1$的限制，則定義了一個單值函數$W_0(x)$（見圖）。我們有$W_0(0)=0$，$W_0(-\frac{1}{e}) = -1$。而在$[-\frac{1}{e},0)$內的$w\le-1$分支，則記為$W_{-1}(x)$，從$W_{-1}(-\frac{1}{e}) = -1$遞減為$W_{-1}(0^-) = -\infty$。

朗伯$W$函數不能用初等函數來表示。它在組合數學中有許多用途，例如樹的計算。它可以用來解許多含有指數的方程，也出現在某些微分方程的解中，例如$y(t) = a y(t - 1)$。

<image>

## 微分和積分

朗伯 $W\,$函數的積分形式為

$$W(x)=\frac{x}{\pi}\int_0^{\pi} \frac{\left(1-v\cot v\right)^2+v^2}{x+v\csc v \cdot e^{-v\cot v}} {\rm{d}}v,|\arg\left(x\right)|<\pi\,$$

$$W(x)=\int_{-\infty}^{-\frac{1}{e}}{-\frac{1}{\pi}}\Im \left[\frac{{\rm{d}}}{{\rm{d}}x}W(x)\right]\ln \left(1-\frac{z}{x}\right){\rm{d}}x\,$$

若 $x\not\in\left[-\frac{1}{e},0\right],k\in{\mathbb{Z}}\,$ ,若 $x\in\left(-\frac{1}{e},0\right),k=1,\pm2,\pm3,...\,$

$$W_k(x)=1+\left(\ln x-1+2k\pi {{\rm{i}}}\right)e^{\frac{{\rm{i}}}{2\pi}\int_0^{\infty}\ln \frac{t-\ln t+\ln x+(2k+1)\pi{\rm{i}}}{t-\ln t+\ln x+(2k-1)\pi{\rm{i}}}\cdot\frac{{\rm{d}}t}{t+1}}=1+\left(\ln x-1+2k\pi {{\rm{i}}}\right)e^{\frac{{\rm{i}}}{2\pi}\int_0^{\infty}\ln \frac{\left(t-\ln t+\ln x\right)^2+\left(4k^2-1\right)\pi^2+2\pi\left(t-\ln t+\ln x\right){\rm{i}}}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}}\,$$

把被積函數的實部和虛部分離出來：

$$W_k(x)=1+\left(\ln x-1+2k\pi {{\rm{i}}}\right)e^{\frac{{\rm{i}}}{2\pi}\int_0^\infty\left[\frac{1}{2}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}+{\rm{i}}\arctan\frac{2\pi\left(t-\ln t+\ln x\right)}{\left(t-\ln t+\ln x\right)^2+\left(4k^2-1\right)\pi^2}\right]\cdot\frac{{\rm{d}}t}{t+1}}$$

$${}_{W_k(x)=1+\frac{\left(\ln x-1\right)\cos\frac{1}{4\pi}\int_0^{\infty}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}-2k\pi\sin\frac{1}{4\pi}\int_0^{\infty}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}+{\rm{i}}\left[\left(\ln x-1\right)\sin\frac{1}{4\pi}\int_0^{\infty}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}+2k\pi\cos\frac{1}{4\pi}\int_0^{\infty}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}\right]}{e^{\frac{1}{2\pi}\int_0^\infty\arctan\frac{2\pi\left(t-\ln t+\ln x\right)}{\left(t-\ln t+\ln x\right)^2+\left(4k^2-1\right)\pi^2}\cdot\frac{\rm{d}t}{t+1}}}}$$

設 $W_k(x)=u+v{\rm{i}},x=t+s{\rm{i}}$，則有 $\left(u+v{\rm{i}}\right)e^{u+v{\rm{i}}}=t+s{\rm{i}}$，展開分離出實部和虛部。

$e^u\left(u\cos v-v\sin v\right)=t,e^u\left(u\sin v+v\cos v\right)=s$,當$s=0$時，易知 $u=-v\cot v$

$${}_{W_k(x)=\frac{\left(1-\ln x\right)\sin\frac{1}{4\pi}\int_0^{\infty}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}-2k\pi\cos\frac{1}{4\pi}\int_0^{\infty}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}}{e^{\frac{1}{2\pi}\int_0^\infty\arctan\frac{2\pi\left(t-\ln t+\ln x\right)}{\left(t-\ln t+\ln x\right)^2+\left(4k^2-1\right)\pi^2}\cdot\frac{\rm{d}t}{t+1}}}\cot\frac{\left(\ln x-1\right)\sin\frac{1}{4\pi}\int_0^{\infty}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}+2k\pi\cos\frac{1}{4\pi}\int_0^{\infty}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}}{e^{\frac{1}{2\pi}\int_0^\infty\arctan\frac{2\pi\left(t-\ln t+\ln x\right)}{\left(t-\ln t+\ln x\right)^2+\left(4k^2-1\right)\pi^2}\cdot\frac{\rm{d}t}{t+1}}} +\frac{\left(\ln x-1\right)\sin\frac{1}{4\pi}\int_0^{\infty}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}+2k\pi\cos\frac{1}{4\pi}\int_0^{\infty}\ln\frac{\left(t-\ln t+\ln x\right)^2+\left(2k+1\right)^2\pi^2}{\left(t-\ln t+\ln x\right)^2+\left(2k-1\right)^2\pi^2}\cdot\frac{{\rm{d}}t}{t+1}}{e^{\frac{1}{2\pi}\int_0^\infty\arctan\frac{2\pi\left(t-\ln t+\ln x\right)}{\left(t-\ln t+\ln x\right)^2+\left(4k^2-1\right)\pi^2}\cdot\frac{\rm{d}t}{t+1}}}{\rm{i}},}$$

$$W_0(x)=1+\left(\ln x-1\right)e^{-\frac{1}{\pi}\int_0^\infty\arg\left(t-\ln t+\ln x+\pi{\rm{i}}\right)\cdot\frac{\rm{d}t}{t+1}},x>0$$

若 $x>\frac{1}{e}$，上式還可化為$W_0(x)=1+\left(\ln x-1\right)e^{-\frac{1}{\pi}\int_0^\infty\arctan\frac{\pi}{t-\ln t+\ln x}\cdot\frac{\rm{d}t}{t+1}}$

由隱函數的求導法則，朗伯$W\,$函數滿足以下的微分方程：

$z\left[1+W(z)\right]\frac{{\rm{d}}}{{\rm{d}}z}W(z)=W(z)$，$z\neq -\frac{1}{e}\,,$

因此：

$\frac{{\rm{d}}}{{\rm{d}}z}W(z)=\frac{ W(z) }{z\left[1 + W(z)\right]}$，$z\neq -\frac{1}{e} \,.$

函數$W(x)\,$，以及許多含有$W(x)\,$的表達式，都可以用$w=W(x)\,$的變量代換來積分，也就是說$x=we^w\,$

$$\int W(x) {\rm{d}}x = x \left[ W(x)+ \frac{1}{W (x) }-1 \right] + C$$

$$\int_0^1 W(x) {\rm{d}}x = \Omega+\frac{1}{\Omega} -2\approx 0.330366$$

其中$\Omega$為歐米加常數。

## 性質

$1\,$、$z^{z^{z^{z^{z^{.^{.^{.}}}}}}}=\lim_{n \to \infty}(z \upuparrows n) = - \frac{W(-\ln z)}{\ln z}$。

其中$\upuparrows$是高德納箭號表示法。

$2\,$、若$z>0 \,$，則$\ln W(z)=\ln z-W(z)\,$

## 泰勒級數

$W_0 \,$在$x=0 \,$的泰勒級數如下：

$$W_0 (x) = \sum_{n=1}^\infty \frac{(-n)^{n-1}}{n!}\ x^n = x - x^2 + \frac{3}{2}x^3 - \frac{8}{3}x^4 + \frac{125}{24}x^5 - \cdots$$

收斂半徑為 $\frac{1}{e}\,$。

## 加法定理

$$W(x)+W(y)=W\left[\frac{xy}{W(x)}+\frac{xy}{W(y)}\right]\,$$

$$x>0,y>0\,$$

## 複數值

實部

$\Re\left[W(x+y{\rm{i}})\right]=\sum_{k=1}^{\infty}\frac{(-k)^{k-1}}{k!}\sqrt{(x^2+y^2)^k}\cos \left(k\arctan\frac{x}{y}\right)\,$ , $x^2+y^2<\frac{1}{e^2}\,$

虛部

$\Im\left[W(x+y{\rm{i}})\right]=\sum_{k=1}^{\infty}\frac{(-k)^{k-1}}{k!}\sqrt{(x^2+y^2)^k}\sin \left(k\arctan\frac{x}{y}\right)\,$, $x^2+y^2<\frac{1}{e^2}\,$

模長

$$|W(x+y{\rm{i}})|=W(\sqrt{x+y})\,$$

模角

$\arg\left[W(x+y{\rm{i}})\right]=\sum_{k=1}^{\infty}\frac{(-k)^{k-1}}{k!}\arctan\left[\cot(k\arctan\frac{x}{y})\right]\,$, $x^2+y^2<\frac{1}{e^2}\,$

共軛值

$\overline{W(x+y{\rm{i}})}=\sum_{k=1}^{\infty}\frac{(-k)^{k-1}}{k!}\sqrt{(x^2+y^2)^k}\left[\cos \left(k\arctan\frac{x}{y}\right)-{\rm{i}}\sin \left(k\arctan\frac{x}{y}\right)\right]\,$, $x^2+y^2<\frac{1}{e^2}\,$

## 特殊值

$$W\left(-\frac{\pi}{2}\right) = \frac{\pi}{2}i$$

$$W\left(-\frac{\ln 2}{2}\right)= -\ln 2$$

$$W\left(-{1\over e}\right) = -1$$

$W\left(1\right) = \Omega=\frac{1}{\int_{-\infty}^{\infty}\frac{{\rm{d}}x}{(e^x-x)^2+\pi^2}}-1\approx 0.56714329\dots \,$（歐米加常數）

$$W(e) = 1\,$$

$$W(e^{e+1}) = e\,$$

$$W\left(\frac{1}{e^{1- \frac{1}{e}}}\right)= \frac{1}{e}$$

$$W({\pi}e^{\pi})=\pi$$

$W(k{\ln k})={\ln k}$ $(k>0)$

$$W({\rm{i}}\pi)=-{\rm{i}}\pi$$

$$W(-{\rm{i}}\pi)={\rm{i}}\pi$$

$$W({\rm{i}}\cos1-\sin1)={\rm{i}}$$

$$W(-\frac{3}{2}{\pi})=-\frac{3}{2}{\pi}{\rm{i}}$$

$$W(-\frac{\sqrt[7]{8}}{7}{\ln 2})=-\frac{32}{7}{\ln 2}$$

$$W(-\frac{\sqrt{3}}{54}{\ln 3})=-\frac{9}{2}{\ln 3}$$

$$W(-\frac{\ln 2}{4})=-4{\ln 2}$$

$$W\left(-1\right)=\frac{e^{\frac{1}{2\pi}\int_0^\infty{1\over t+1}\arctan{2\pi\over t-\ln t}{\rm{d}}t}-\cos\left[\frac{1}{4\pi}\int_0^\infty{1\over t+1}\ln{\left(t-\ln t\right)^2\over 4\pi^2+\left(t-\ln t\right)^2}{\rm{d}}t\right]+\pi\sin\left[\frac{1}{4\pi}\int_0^\infty{1\over t+1}\ln{\left(t-\ln t\right)^2\over 4\pi^2+\left(t-\ln t\right)^2}{\rm{d}}t\right]-{\rm{i}}\left\{\pi\cos\left[\frac{1}{4\pi}\int_0^\infty{1\over t+1}\ln{\left(t-\ln t\right)^2\over 4\pi^2+\left(t-\ln t\right)^2}{\rm{d}}t\right]+\sin\left[\frac{1}{4\pi}\int_0^\infty{1\over t+1}\ln{\left(t-\ln t\right)^2\over 4\pi^2+\left(t-\ln t\right)^2}{\rm{d}}t\right]\right\}}{e^{\frac{1}{2\pi}\int_0^\infty{1\over t+1}\arctan{2\pi\over t-\ln t}{\rm{d}}t}}\approx -0.31813-1.33723{\rm{i}}$$

$$W(-\frac{\ln k}{k})=-\ln k$$

$$W\left[-\frac{\ln (x+1)}{x(x+1)^{\frac{1}{x}}}\right]=-\frac{x+1}{x}\ln (x+1)>,-1<x<0$$

## 應用

許多含有指數的方程都可以用$W\,$函數來解出。一般的方法是把未知數都移到方程的一側，並設法化為$Y= Xe^X \,$的形式。

### 例子

例子1：$$2^t = 5 t\,$$

$$\Rightarrow 1 = \frac{5 t}{2^t}\,$$

$$\Rightarrow 1 = 5 t \, e^{-t \ln 2}\,$$

$$\Rightarrow \frac{1}{5} = t \, e^{-t \ln 2}\,$$

$$\Rightarrow -\frac{\ln 2}{5} = ( - \, t \, \ln 2 ) \, e^{-t \ln 2}\,$$

$$\Rightarrow -t \ln 2 = W_k \left (-\frac{\ln 2}{5} \right )\,$$

$$\Rightarrow t = -\frac{ W_k \left ( -\frac{\ln 2}{5} \right )}{\ln 2}\,$$

更一般地，以下的方程

$$Q^{a x + b} = c x + d \,$$

其中

$$Q > 0 \land Q \neq 1\land c \neq 0$$

兩邊同乘: $\frac{a}{c}$。

得到：$\frac{a}{c} Q^{ax+b} = ax + \frac{ad}{c} \,$

同除以：$Q^{ax} \,$。

得到：$\frac{a}{c} Q^{b} = \left(ax + \frac{ad}{c} \right)Q^{-ax} \,$

同除：$Q^{\frac{ad}{c}} \,$。

$$\frac{a}{c} Q^{b-\frac{ad}{c}}= \left(ax + \frac{ad}{c}\right)Q^{-\left(ax+\frac{ad}{c}\right)} \,$$

可以用變量代換

令$t = a x + \frac{a d}{c}$

化為

$$t Q^{-t} = \frac{a}{c} Q^{b-\frac{a d}{c}}$$

即：$t \left(e^{\ln Q}\right) ^{-t} = \frac{a}{c} Q^{b-\frac{a d}{c}}$

同乘：${\ln Q} \,$

得出

$$t{\ln Q} \cdot e^{-t \ln Q}={\ln Q} \cdot \frac{a}{c} Q^{b-\frac{a d}{c}}$$

故$t{\ln Q}=-W_k\left(-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}}\right)$

帶入$t= a x + \frac{a d}{c}$

為

$$\left(ax+\frac{ad}{c}\right){\ln Q}=-W_k\left(-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}}\right)$$

因此最終的解為

$$x = -\frac{W_k\left(-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}}\right)}{a\ln Q} - \frac{d}{c}$$

若輔助方程：$xe^x=-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}}$中。

$-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}} \in \left(-\infty ,-\frac{1}{e} \right)$,

輔助方程無實數解，原方程亦無實解。

若：$-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}} \in \left\{-\frac{1}{e}\right\} \cup\mathbf [0,+\infty )$,

輔助方程有一實數解，原方程有一實解：

$$x = -\frac{W_k\left(-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}}\right)}{a\ln Q} - \frac{d}{c}$$

若: $-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}} \in \left(-\frac{1}{e},0 \right)$,

輔助方程有二實解，設為$W\left(-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}}\right)$。

${\rm{W}}_{-1}\left(-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}}\right)$。

為

$$x_1=-\frac{W\left(-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}}\right)}{a\ln Q} - \frac{d}{c}$$

$$x_2=-\frac{{\rm{W}}_{-1}\left(-\frac{a\ln Q}{c}\,Q^{b-\frac{a d}{c}}\right)}{a\ln Q} - \frac{d}{c}$$

例子2

用類似的方法，可知以下方程的解

$$x^x={\mathrm{t}}\, ,$$

為

$$x=\frac{\ln{\rm{t}}}{W(\ln {\rm{t}})}\,$$

或

$$x=\exp\left(W_k\left[\ln({\rm{t}})\right]\right).$$

例子3

以下方程的解

$$x \log_b {x} = a \,$$

具有形式

$$x = \frac{a {\ln b}}{W_k\left(a {\ln b}\right)}$$

例子4：$$x^a-b^x=0\,$$

$a > 0 \,$ : $b > 0 \,$ : $x > 0 \,$

取對數。

$$a \ln x=x \ln b \,$$

$$\frac{\ln x}{x}=\frac{\ln b}{a}\,$$

$$e^{\frac{\ln x}{x}}=e^{\frac{\ln b}{a}} \,$$

$$x^{\frac{1}{x}}=b^{\frac{1}{a}}\,$$

取倒數。

$$\left(\frac{1}{x} \right)^{\frac{1}{x}}=b^{-\frac{1}{a}}\,$$

$$\frac{1}{x} =-\frac{\ln b}{aW\left(-\frac{1}{a} \ln b\right)}\,$$

最終解為 : $x=-\frac{a}{\ln b}W_k\left(-\frac{\ln b}{a}\right)\,$

例子5：$$(ax+b)^n=u^{cx+d} \,$$

兩邊開$n \,$次方並除以$a \,$得

$$x+\frac{b}{a}=\frac{u^{\frac{c}{n}x+\frac{d}{n}}}{a}\left(\cos\frac{2k\pi}{n}+{\rm{i}}\sin\frac{2k\pi}{n}\right)\,$$

令$u=e^{\ln u}\,$。

化為

$$x+\frac{b}{a}=\frac{e^{\frac{c\ln u}{n}x+\frac{d\ln u}{n}}}{a}\left(\cos\frac{2k\pi}{n}+{\rm{i}}\sin\frac{2k\pi}{n}\right)\,$$

兩邊同乘

$-\frac{c\ln u}{n}u^{-\frac{c}{n}x-\frac{cb}{na}}\,$。

$$\left(-\frac{c\ln u}{n}x-\frac{cb\ln u}{na}\right)e^{-\frac{c\ln u}{n}x-\frac{cb\ln u}{na}}=-\frac{c\ln u}{na}u^{\frac{d}{n}-\frac{cb}{na}}\left(\cos\frac{2k\pi}{n}+{\rm{i}}\sin\frac{2k\pi}{n}\right)\,$$

最終得

$$x_k=-\frac{n}{c\ln u}W_k\left[-\frac{c\ln u}{na}u^{\frac{d}{n}-\frac{cb}{na}}\left(\cos\frac{2k\pi}{n}+{\rm{i}}\sin\frac{2k\pi}{n}\right)\right]-\frac{b}{a}\,$$

$$k\in{\mathbb{Z}}\,$$

## 一般化

標準的 Lambert W 函數可用來表示以下超越代數方程式的解：

$$e^{-c x} = a_o (x-r) ~~\quad\qquad\qquad\qquad\qquad(1)$$

其中 a0, c 與 r 為實常數。

其解為$x = r +\tfrac{W\left(\frac{ce^{-c r}}{a_o}\right)}{c}$

Lambert W 函數之一般化包括:

- 一項在低維空間內廣義相對論與量子力學的應用（量子引力），實際上一種以前未知的 連結 於此二區域中，如 「Journal of Classical and Quantum Gravity」所示其 (1) 的右邊式現為二維多項式 x：

$$e^{-c x} = a_o (x-r_1 ) (x-r_2 ) ~~\qquad\qquad(2)$$

其中 r1 和 r2 是不同實常數，為二維多項式的根。於此函數解有單一引數 x 但 ri 和 ao 為函數的參數。如此一來，此一般式類似於 「hypergeometric」（超幾何分布）函數與 「Meijer G「，但屬於不同類函數。當 r1 = r2，(2)的兩方可分解為 (1) 因此其解簡化為標準 W 函數。(2)式代表著 「dilaton」（軸子）場的方程，可據此推導線性，雙體重力問題 1+1 維（一空間維與一時間維）當兩不等（靜止）質量，以及，量子力學的特徵能Delta位勢阱給不等電位於一維空間。

- 量子力學的一特例特徵能的分析解三體問題，亦即（三維）氫分子離子。於此 (1)（或 (2)）的右手邊現為無限級數多項式之比於 x：

$$e^{-c x} = a_o \frac{\prod_{i=1}^{\infty} (x-r_i )}{ \prod_{i=1}^{\infty} (x-s_i)} \qquad \qquad\qquad(3)$$

其中 ri 與 si 是相異實常數而 x 是特徵能和內核距離R之函數。式 (3) 與其特例表示於 (1) 和 (2) 是與一更大類型延遲微分方程。由於哈代的「虛假導數」概念，多根的特殊情況得以解決。

Lambert "W" 函數於基礎物理問題之應用並未完全即使標準情況如 (1) 最近在原子，分子，與光學物理領域可見以及黎曼假設的 Keiper-Li 準則

## 圖象

<image>

- z = Re(W0(x + i y))

<image>

- z = Im(W0(x + i y))

<image>

- W0(x + i y)

<image>

## 計算

W函數可以用以下的遞推關係算出：

$$w_{j+1}=w_j-\frac{w_j e^{w_j}-z}{e^{w_j}(w_j+1)-\frac{(w_j+2)(w_je^{w_j}-z)} {2w_j+2}}$$