# Fastest way to calculate Taylor expansion coefficients where the coefficients are functions of another variable that is evaluated at several poinst

I have an expression say V which is a function of 18 variables, 17 of which will be Taylor expanded to some order and the last one which will be explicitly evaluated (it is an angle) at several points as the angle is rotated. Therefore, I have an expression t defined by

t = ReplaceAll[D[V, Delete[powers[[i]], 0]]/Product[Factorial[powers[[i]][[k]][[2]]], {k, 1, Length[powers[[i]]]}], EqCoordsSubs];


which essentially is derivative of V with the number for each coordinate defined in powers[[i]] and the EqCoordSub is just the substitution of the 17 coordinates to their equilibrium values. t is then evaluated at several angles (coordinate 18). However, I'm now aware of SeriesCoefficient which got me thinking of potentially more efficient ways of doing this. Would this (or perhaps another function) be more suitable for what I am trying to do?