Recently, I have come across the Taylor's formula for functions between Banach spaces. To better understand the machinery, I try to see how higher-order derivatives look like in a simple example, i.e., $f:\mathbb R \to \mathbb R, x \mapsto x^3$.
Let $E,F$ be Banach spaces and $f:E \to F$ be $m$-times continuously differentiable. Here we use Fréchet derivative, i.e., $\partial^m f : E \to \mathcal L^m(E, F)$ with $\mathcal L^m(E, F)$ the space of continuous multilinear maps from $E^m$ to $F$. For clarity, I will enclose the arguments of multilinear maps in square brackets.
Let $E = F := \mathbb R$. We take $f: x \mapsto x^3$ as an example. It's easy to see that $\partial f : x \mapsto (h_3 \mapsto 3x^2h_3)$. I pretend that $\partial f$ has an untrue form $x \mapsto 3x^2h_3$, which is much simpler than its true form $\partial f:E \to \mathcal L(E, F)$. I differentiate this untrue form and get an untrue $\partial^2 f: x \mapsto (h_2 \mapsto 6xh_3h_2)$. I guess the correct form is $\partial^2 f : x \mapsto \big ( (h_2, h_3) \mapsto 6xh_3h_2 \big)$.
To make sure that above process is valid. I aim at proving
$$ \partial^m f(x)[h_1, \ldots, h_m] = \partial g(x)[h_1] \text{ with }g:y \mapsto \partial^{m-1}f(y)[h_2,\ldots, h_m] $$ for all $x,h_1,\ldots, h_m \in E$.
Assume we know $\partial^{m-1} f$, this result reduces the computation of $\partial^m f:E \to \mathcal L^m(E, F)$ to that of $g:E \to F$.
Could you have a check on my attempt?
I posted my proof separately so that I can accept my own answer and thus remove my question from unanswered list. If other people post answers, I will happily accept theirs.