Some new GP features lambertw lambertw now supports complex arguments and all branches. ? l=lambertw(-1/4) %66 = -0.3574029561813889031 ? l*exp(l) %67 = -0.2500000000000000000 ? l=lambertw(-1) %68 = -0.3181315052047641353+1.337235701430689409*I ? l*exp(l) %69 = -1.000000000000000000+5.421010862427522170E-20 ? l=lambertw(-1,1) %70 = -2.062277729598283885+7.588631178472512623*I ? l*exp(l) %71 = -1.000000000000000000+1.355252715606880543E-19