# Principal branch of the Lambert W function W(z) cplot(lambertw, [-1,1], [-1,1], points=50000)