Problem:
Find the first order expansion about $x = 5$ of the composition $f(g(x))$ for $f(x) = x^5$ and $g(x) = $x^5$
from sympy import *
from sympy.abc import x

f = x**5
g = x**5
series = f.subs(x, g).series(x, 5, 1)
print(series)
