I'm trying to show the statement in the title is true for Taylor series centered at any point $\textbf{a}$ in $\mathbb{R}^n$.
It's easy for me to show this for when $\textbf{a}=\vec{0}$ using the fact that the multivariate Taylor expansion for a function $f$ is:
$$f(x) = \sum_{|\alpha|\leq n}\frac{\partial ^{\alpha}f(\textbf{a})}{\alpha!}(\textbf{x} - \textbf{a})^{\alpha} + r_{n,\textbf{a}}(x)$$
where $\alpha$ is the multi-index, and $r$ being the remainder term. Since the multivariate polynomial's variable portion is exactly like $(\textbf{x} - \textbf{a})^{\alpha}$ when $\textbf{a} = \vec{0}$, and the fractional term reduces to the same coefficient as in the original polynomial, we get back the original polynomial. But I am completely lost as to how to proceed in the case where $\textbf{a} = \vec{0}$, a hint would be greatly appreciated it.