plot 2+LambertW(-x*exp(-x)) for x=1 to 2