integrate (exp (x)/(1+exp (2*x)))
-I*log(exp(x) - I)/2 + I*log(exp(x) + I)/2
integrate(exp(x)/(exp(2*x) + 1), x)
plot(RootSum(4*_z**2 + 1, Lambda(_i, _i*log(2*_i + exp(x)))))
solve(RootSum(4*_z**2 + 1, Lambda(_i, _i*log(2*_i + exp(x)))), x)
diff(RootSum(4*_z**2 + 1, Lambda(_i, _i*log(2*_i + exp(x)))), x)
series(RootSum(4*_z**2 + 1, Lambda(_i, _i*log(2*_i + exp(x)))), x, 0, 10)
Experiment with SymPy at SymPy Live.