Вы можете использовать класс Poly
и метод all_coeffs
. Ссылка доступна здесь: https://docs.sympy.org/latest/modules/polys/reference.html
Это то, что она дала бы в вашем примере, предполагая, что все символы были объявлены:
>>> pol = sp.Poly(((a+x**2) / b + x /c + x*d + k), x); pol
Poly(1/b*x**2 + (c*d + 1)/c*x + (a + b*k)/b, x, domain='ZZ(a,b,c,d,k)')
>>> pol.all_coeffs()
[1/b, (c*d + 1)/c, (a + b*k)/b]
После этого вы можете получить доступ к каждому коэффициенту с помощью егоположение.