series(RootSum(4*_z**2 + 1, Lambda(_i, _i*log(2*_i + exp(x)))), x, 0, 10)

See what
Wolfram|Alpha has to say.

Want to compute something more complicated? Try a full Python/SymPy console at SymPy Live.