Statement

Let f∈C∞((−ε,ε)×U)f\in C^\infty(( -\varepsilon,\varepsilon)\times U), with U⊆RnU\subseteq\mathbb R^n open. If m≥1m\ge1 and

∂tjf(0,x)=0(0≤j<m),\partial_t^j f(0,x)=0\qquad(0\le j<m),

then f(t,x)=tmg(t,x)f(t,x)=t^m g(t,x) for a smooth function gg, given by

g(t,x)=1(m−1)!∫01(1−s)m−1∂tmf(st,x) ds.g(t,x)=\frac1{(m-1)!}\int_0^1(1-s)^{m-1}\partial_t^m f(st,x)\,ds.

This is smooth division by the normal coordinate.

Proof and boundary value

Repeated use of the gives the integral remainder formula. Smoothness follows by on the compact interval [0,1][0,1]. In particular,

g(0,x)=∂tmf(0,x)m!.g(0,x)=\frac{\partial_t^m f(0,x)}{m!}.
Limitation

This divides by a finite power of tt. Flatness of two functions at zero does not by itself make their quotient smooth or bounded; for example e−1/t2/e−2/t2=e1/t2e^{-1/t^2}/e^{-2/t^2}=e^{1/t^2} on t>0t>0.