The output contains functions with complex values, but the end result is real:
int = Integrate[(1/(1 + Exp[-x]))*Log[1/(1 + Exp[-x])], x]
int2 = FullSimplify[ComplexExpand[int], Element[x, Reals]]
FullSimplify[ComplexExpand[int] == int2, Element[x, Reals]]
Plot[int - int2, {x, -1, 1}]
Unfortunately we cannot check the derivative of int2
, because int2
contains Re
, which is not differentiable in the complex sense.