Sea \( f:\mathbb{R}^n \to R \) entonces:
Editado
\( \displaystyle E(x) = f(x_0 + x) - f(x) - \sum_{i=1}^n D_if(x_0)x_i \).
Tenemos que por existir las derivadas parciales en una bola de radio r y ser continuas en \( x_0 \) obtenemos que dado \( \epsilon > 0 \) existe \( r >\delta > 0 \) tal que si:
\( h \in B(x_0,\delta) \) entonces \( |D_if(h) - D_if(x_o)| < \epsilon \)
Ahora viene el punto dificil:
Editado
Sea \( x \in B(\color{red}0_n\color{black},\delta) \) tomemos los puntos definidos por \( u_o = x_0 \) y \( u_i =x_0 + \displaystyle \sum_{k=1}^i x_{\color{red}k\color{black}} \cdot e_{\color{red}k\color{black}} \) donde \( i \in \{1,2, \cdots , n\} \)
Luego tenemos:
\( f(x_o+x) - f(x_0) =\displaystyle \sum_{i=1}^{n} (f(u_i) - f(u_{i-1})) = \sum_{i=1}^{n} (f(u_{i-1} + x_i \cdot e_i) - f(u_{i-1}) = \sum_{i=1}^{n} D_i\color{red}f\color{black}(u_{i-1} + \varepsilon_i e_i) \cdot x_i \)
Como \( u_{i-1} + \varepsilon_i e_i \in B(x_0 , \delta) \) tenemos que \( |D_i\color{red}f\color{black}(u_{i-1} + \varepsilon_i e_i) -D_i(x_0)| < \epsilon \).
Entonces:
\( |E(x)| = |f(x_0 + x) - f(x) - \sum_{i=1}^n D_if(x_0)x_i| = \displaystyle |\sum_{k=1}^n (D_i\color{red}f\color{black}(u_{i-1} + \varepsilon_i e_i) -D_i\color{red}f\color{black}(x_0))x_i |< \epsilon \cdot \sum_{k=1}^n |x_i| < \epsilon \cdot n \cdot ||x|| \)