import sympy as sym
x, a = sym.symbols('x, a')
f = sym.sqrt(a**2+x**2)
print(sym.series(f, x, 0, 4))
print(sym.series(sym.sinh(x)*sym.sin(x),x, 0, 12))