Sinon, si tu veut reprendre les notations correspondant à :
\ =\ f(X_o)\ +\ d_f(X_o)(H)\ +\ \frac{1}{2!}d^2_f(X_o)(H,H)\ +\ \frac{1}{3!}d^3_f(X_o)(H,H,H)+...)
(je préfère cette formule à celle donnée par Arkhnor du fait qu'elle montre mieux qu'il n'y a aucune différence avec le cas des fonctions d'une seule variable)
alors, en considérant que les coordonnées des

sont
)
, on a :
(avec

termes dans la somme !!!)
Qui est une aplication linéaire en chacun des

et qui est en fait symétrique en les

(Théorème de Schwarz)
C'est justement cette symétrie qui permet de regrouper des termes lorsque l'on calcule
(H,H,...,H))
dans la formule de Taylors
et qui fait apparaitre les fameux coefficients multinomiaux :
Ils correspondent trés précisément aux nombres de p-uplets
\in\{1..n\}^p)
tels qu'il y ait

d'entre eux égaux à 1,

d'entre eux égaux à 2... ,

d'entre eux égaux à n (avec évidement

)
EDIT : Il y a bien des fractions dans le résultat final qui proviennent en fait des

avant le

dans la formule de Taylors "classique" (i.e. celle en haut de ce post)
Par contre, si tu utilise les regroupement sus-mentionnée, il apparrait en plus les coeff multinomiaux

(avec

) qui, évidement sont entiers (c'est
le nombre de p-uplet tels que...)