Using your example, it is symbolic taylorand coeffscan be used to obtain a vector of coefficients:
syms x
f = 1/(x^2 + 4*x + 9);
ts = taylor(f,x,0,'Order',4) % 4-th order Taylor series of f about 0
c = coeffs(ts)
which returns
ts =
(8*x^3)/6561 + (7*x^2)/729 - (4*x)/81 + 1/9
c =
[ 1/9, -4/81, 7/729, 8/6561]
vpa double, c .