Theorem: Derivative of Composite Function/Jacobians

Let $k, m, n \in \N$.
Let $U$ be an open set of $\R^n$.
Let $V$ be an open set of $\R^m$.
Let $\mathbf f = \paren {f_1, f_2, \ldots, f_m}^\intercal: U \to V$ be a vector valued function, differentiable at $\mathbf x = \paren {x_1, x_2, \ldots, x_n}^\intercal \in U$.
Let $\mathbf g = \paren {g_1, g_2, \ldots, g_k}^\intercal: V \to \R^k$ be a vector valued function, differentiable at $\map {\mathbf f} {\mathbf x} = \paren {\map {f_1} {\mathbf x}, \map {f_2} {\mathbf x}, \ldots, \map {f_m} {\mathbf x} }^\intercal \in V$.
The Jacobian matrix of a composite function $\mathbf g \circ \mathbf f$ is obtained by multiplying the Jacobian matrices of $\mathbf f$ and $\mathbf g$:
:$\map {\mathbf J_{\mathbf g \circ \mathbf f} } {\mathbf x} = \map {\mathbf J_{\mathbf g} } {\map {\mathbf f} {\mathbf x} } \map {\mathbf J_{\mathbf f} } {\mathbf x}$

Proof:
{{MissingLinks}}
The Jacobian matrix of $\mathbf g \circ \mathbf f$ at $\mathbf x$ is defined to be the $k \times n$ matrix:
:<nowiki>$\map {\mathbf J_{\mathbf g \circ \mathbf f} } {\mathbf x} = \begin{pmatrix}
 \map {\dfrac {\map \partial {g_1 \circ \mathbf f} } {\partial x_1} } {\mathbf x} & \cdots & \map {\dfrac {\map \partial {g_1 \circ \mathbf f} } {\partial x_n} } {\mathbf x} \\
 \vdots & \ddots & \vdots \\
 \map {\dfrac {\map \partial {g_k \circ \mathbf f} } {\partial x_1} } {\mathbf x} & \cdots & \map {\dfrac {\map \partial {g_k \circ \mathbf f} } {\partial x_n} } {\mathbf x}
\end{pmatrix}$</nowiki>
The Jacobian matrix of $\mathbf g$ at $\map {\mathbf f} {\mathbf x}$ is defined to be the $k \times m$ matrix:
:<nowiki>$\map {\mathbf J_{\mathbf g} } {\map {\mathbf f} {\mathbf x} } := \begin{pmatrix}
 \map {\dfrac {\partial g_1} {\partial x_1} } {\map {\mathbf f} {\mathbf x} } & \cdots & \map {\dfrac {\partial g_1} {\partial x_m} } {\map {\mathbf f} {\mathbf x} } \\
 \vdots & \ddots & \vdots \\
 \map {\dfrac {\partial g_k} {\partial x_1} } {\map {\mathbf f} {\mathbf x} } & \cdots & \map {\dfrac {\partial g_k} {\partial x_m} } {\map {\mathbf f} {\mathbf x} }
\end{pmatrix}$</nowiki>
The Jacobian matrix of $\mathbf f$ at $\mathbf x$ is defined to be the $m \times n$ matrix:
:<nowiki>$\map {\mathbf J_{\mathbf f} } {\mathbf x} = \begin{pmatrix}
 \map {\dfrac {\partial f_1} {\partial x_1} } {\mathbf x} & \cdots & \map {\dfrac {\partial f_1} {\partial x_n} } {\mathbf x} \\
 \vdots & \ddots & \vdots \\
 \map {\dfrac {\partial f_m} {\partial x_1} } {\mathbf x} & \cdots & \map {\dfrac {\partial f_m} {\partial x_n} } {\mathbf x}
\end{pmatrix}$</nowiki>
The entry in row $r$, column $c$ of $\map {\mathbf J_{\mathbf g \circ \mathbf f} } {\mathbf x}$ is:
:$\ds \map {\dfrac {\map \partial {g_r \circ \mathbf f} } {\partial x_c} } {\mathbf x}$
The entry in row $r$, column $c$ of the matrix product $\map {\mathbf J_{\mathbf g} } {\map {\mathbf f} {\mathbf x} } \map {\mathbf J_{\mathbf f} } {\mathbf x}$ is:
:$\ds \sum_{i \mathop = 1}^m \map {\dfrac {\partial g_r} {\partial x_i} } {\map {\mathbf f} {\mathbf x} } \map {\dfrac {\partial f_i} {\partial x_c} } {\mathbf x}$
for each $r \in \set {1, 2, \ldots, k}, c \in \set{1, 2, \ldots, n}$.
By Derivative of Composite Function:
:$\ds \map {\dfrac {\map \partial {g_r \circ \mathbf f} } {\partial x_c} } {\mathbf x} = \sum_{i \mathop = 1}^m \map {\dfrac {\partial g_r} {\partial x_i} } {\map {\mathbf f} {\mathbf x} } \map {\dfrac {\partial f_i} {\partial x_c} } {\mathbf x}$
Hence the result.
{{qed}}