Note that to evaluate the Taylor polynomial as an approximation of f, e.g. with polyval(), 
the coefficients must be reversed, 

and the point of the Taylor expansion must be subtracted from the argument

( http://docs.sympy.org/dev/modules/mpmath/calculus/approximation.html#taylor-series-taylor )
from sympy import Symbol, exp, series
from sympy.mpmath import polyval,taylor,mp
mp.dps = 15; mp.pretty = True

def getpolyval(f,p,n=6,center=0):
    return polyval(taylor(f,center,n)[::-1],p)

x = Symbol('x')

def f(x):
    return 1/(exp(x)+1)

print series(f(x),x,)
print getpolyval(f,x,n=6).expand()

print 1./48,-1./480