# W(x) = LambertW(x), is inverse of exp(x)*x;

W(b_*exp(b_)):=if(b>= -1,b);
W(exp(1)):=1;
W(0):=0;