(28/48/50) Lambert W Function

[Albert Chan]

(01-30-2024 12:18 AM)Albert Chan Wrote:  With limited domain (A ≈ -1/e), fixed loops is simple to implement.

For my Python code (~ 15 digits precision), it never get over 5 loops.
Perhaps you can use fixed Newton loops instead of testing convergence.

Also, you may use this to implement complex log1p(x)

y=1+x --> log1p(x) ≈ ln(y) - (y-1-x)/y

https://www.hpmuseum.org/forum/thread-55...#pid111682