y'' = lambertW(y')