Let
\(\vx=\langle x_1,\ldots,x_n \rangle\) be near
\(\va\text{,}\) and for
\(i=0,1,\ldots,n\) define
\begin{equation*}
\vx_i
:=
\langle x_1,\ldots,x_i,a_{i+1},\ldots,a_n \rangle,
\end{equation*}
so that
\(\vx_0=\va\) and
\(\vx_n=\vx\text{.}\) Then
\begin{equation*}
f(\vx)-f(\va)
=
\sum_{i=1}^n \bigl(f(\vx_i)-f(\vx_{i-1})\bigr).
\end{equation*}
For each
\(i\text{,}\) consider the one-variable function
\begin{equation*}
g_i(t)
:=
f(\langle x_1,\ldots,x_{i-1},t,a_{i+1},\ldots,a_n \rangle).
\end{equation*}
Since
\(D_if\) exists near
\(\va\text{,}\) the function
\(g_i\) is differentiable, hence continuous. By the one-variable Mean Value Theorem, there exists a number
\(c_i\) between
\(a_i\) and
\(x_i\) such that
\begin{equation*}
f(\vx_i)-f(\vx_{i-1})
=
D_if(\langle x_1,\ldots,x_{i-1},c_i,a_{i+1},\ldots,a_n \rangle)
(x_i-a_i).
\end{equation*}
\begin{equation*}
A_i(\vx)
:=
D_if(\langle x_1,\ldots,x_{i-1},c_i,a_{i+1},\ldots,a_n \rangle)
\end{equation*}
for
\(\vx \ne \va\text{,}\) and set
\(A_i(\va):=D_if(\va)\text{.}\) Then
\begin{equation*}
f(\vx)-f(\va)=\sum_{i=1}^n A_i(\vx)(x_i-a_i).
\end{equation*}
The point
\(\langle x_1,\ldots,x_{i-1},c_i,a_{i+1},\ldots,a_n \rangle\) lies on the segment joining
\(\vx_{i-1}\) and
\(\vx_i\text{,}\) so it tends to
\(\va\) as
\(\vx \to \va\text{.}\) Since
\(D_if\) is continuous at
\(\va\text{,}\) we obtain
\begin{equation*}
A_i(\vx)\to D_if(\va)
\qquad \text{as } \vx \to \va.
\end{equation*}
Thus each
\(A_i\) is continuous at
\(\va\text{.}\) By the scalar-valued form of the Caratheodory criterion in
TheoremΒ 2.3.3, it follows that
\(f\) is differentiable at
\(\va\text{.}\)