Let $f\in C^r(A\subset \mathbb R^n,\mathbb R^m)$,
$Df:A\subset\mathbb R^n\to L(\mathbb R^n,\mathbb R^m)$ so that $Df(x):\mathbb R^n\to\mathbb R^m$ is $f$'s total derivative,
(abusing notation) $D^2f(x): \mathbb R^n \times \mathbb R^n\to\mathbb R^m$ be the bilinear map defined by $$[D^2f(x)](x_1,x_2)=[[D(Df)(x)] (x_2)](x_1),$$ (still abusing notation) $D^3f(x): \mathbb R^n \times \mathbb R^n\times \mathbb R^n\to\mathbb R^m$ be the trilinear map defined by $$[D^3f(x)](x_1,x_2,x_3)=[[D(D^2f)(x)] (x_3)](x_1,x_2),$$ and so on.
For $ r\geq 3:$
How do I show that $D^rf(x)$ is symmetric (i.e. returns the same value for every permutation of inputs)?
And why is $D^rf(x)$ being symmetric equivalent to the fact that for each component, $f$'s $r$th partials can be taken in any order?
P.S. I have read the proof for the case of second-order derivatives.
Reference: Marsden's Elementary Classical Analysis