from sympy import init_printing, symbols, Function
init_printing()
x, h = symbols('x,h')
f = Function('f')
pprint(f(x).series(x, x0=h, n=3))
Как сделать символическое разложение Тейлора неизвестной функции $f(x)$ с помощью sympy