integrate (exp (x)/(1+exp (2*x)))

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.