Let $f: \mathbb{R}^{n^{2}} \rightarrow \mathbb{R}^{n^{2}}$ defined by $f(X)=X^6$, where $\mathbb{R}^{n^{2}}$ is the space of square $n\times n$ matrices.
Show that $f$ is a differential application.
I try to prove that using the definition of a differential application in this way:
Let $X \in \mathbb{R}^{n^{2}}$ and $Y\neq0 \in \mathbb{R}^{n^{2}}$.
$$f(X+Y) - f(X) = f'(X)\cdot Y + r(Y) $$ $$\Leftrightarrow r(Y)=f(X+Y) - f(X) - f'(X)\cdot Y $$ $$\Leftrightarrow \frac{r(Y)}{|Y|}=\frac{f(X+Y) - f(X) - f'(X)\cdot Y}{|Y|}.$$ Then $$\lim_{|Y|\rightarrow0}\frac{r(Y)}{|Y|}=\lim_{|Y|\rightarrow0}\frac{f(X+Y) - f(X) - f'(X)\cdot Y}{|Y|} $$ $$\Leftrightarrow\lim_{|Y|\rightarrow0}\frac{r(Y)}{|Y|}=\lim_{|Y|\rightarrow0}\frac{(X+Y)^6 - X^6 - f'(X)\cdot Y}{|Y|}.$$
I want to prove that $\lim_{|Y|\rightarrow0}\frac{r(Y)}{|Y|}=0$. If I could get this, the result follows by the differential application definition. The problem is that the algebraic manipulations aren't so simple, and I don't have any additional properties for these matrices.
I try something else: There is a result that says every function $f:\mathbb{R^n}\rightarrow \mathbb{R^n}$ such that $f \in C^1$, then $f$ is differentiable. I thought to calculate the partial differentials of $f$ in each variable $x_{i} \in \mathbb{R}$ and proof that derivative is continuous in each $x_i$. But I don't know how to calculate these derivatives because they are matrices. Could someone help me, please?