朗伯函数

Question

我在python上使用了Lambert函数，但没有得到正确的结果。
我的实际代码是

E = np.linspace(0,1,100)
y = []
for f in E:
   W = scipy.special.lambertw((f-1)/np.exp(1), k=0, tol=1e-10)
   y.append(W)
   

在上面的代码中，理论上，值应该介于-1和0之间。但这里的情况似乎并非如此

我注意到，就主要分支机构而言，结果略有差异。Wo（-1/e）=-1但当我在python中使用以下代码时，答案是错误的

scipy.special.lambertw(-1/np.exp(1), k=0, tol=1e-8)
>>> (nan+nanj)

有人能帮我吗
谢谢大家!