sqrt(-2-4 lambertw(-1,-1/(2 sqrt(e))))