Taylor Series for functions $f:R^n\rightarrow R^n$ I've been told that we can write a taylor series for functions $f:R^n\rightarrow R$ but we can't write one for $f:R^n\rightarrow R^n$.  I'm not quite sure why this not possible, but I suspect it have something to do with the mean value theorem.  Could anyone shed some light on this?
 A: If $f : \mathbb{R}^n \to \mathbb{R}^n$, then you can split $f = (f_1, f_2, \ldots, f_n)$, and (under the usual smoothness assumptions) each component can be expanded in a Taylor series.
A: For  each $f_j:\mathbb{R}^n\to\mathbb{R}$ in $f=(f_1,\ldots,f_j,\ldots,f_n)$ use a  notation for multi index $\alpha=(\alpha_1,\ldots,\alpha_i,\ldots,\alpha_n)$ in Higher-order partial derivatives, norms, factorials $!$ and coordinates of vectors $x,v\in\mathbb{R}^n$. Then for $f:\mathbb{R}^n\to\mathbb{R}^n$ we have Taylor series of $f(x+v)=$
\begin{align}
\bigg(\sum_{k=0}^{\infty}\sum_{|\alpha|=k}\frac{1}{\alpha!}\partial^\alpha f_1(x)\cdot v^\alpha,\ldots,\sum_{k=0}^{\infty}\sum_{|\alpha|=k}\frac{1}{\alpha!}\partial^\alpha f_j(x)\cdot v^\alpha,\ldots, \sum_{k=0}^{\infty}\sum_{|\alpha|=k}\frac{1}{\alpha!}\partial^\alpha f_n(x)\cdot v^\alpha\bigg)
\end{align}
A: The same formula that applies for scalar fields applies for vector fields.
$$f(x) = f(a) + [(x-a) \cdot \nabla] f \Big|_a + \frac{1}{2!} [(x-a) \cdot \nabla]^2 f \Big |_a + \ldots$$
$f$ can be scalar or vector.
