Statement

Fix c>0c>0, an integer j0j\ge0, and a smooth function b(y,η)b(y,\eta) for y<ε|y|<\varepsilon, ηURd\eta\in U\subseteq\mathbb R^d. For 0<y<ε0<y<\varepsilon, define

H(y,η)=ec/y20yec/u2ujb(u,η)du.H(y,\eta)=e^{c/y^2}\int_0^y e^{-c/u^2}u^j b(u,\eta)\,du.

Then H(y,η)=yj+3B(y,η)H(y,\eta)=y^{j+3}B(y,\eta), where BB extends smoothly through y=0y=0 and

B(0,η)=b(0,η)2c.B(0,\eta)=\frac{b(0,\eta)}{2c}.
Fixed-domain integral proof

Use the s=u2y2s=u^{-2}-y^{-2}. It gives

B(y,η)=120ecs(1+y2s)(j+3)/2b ⁣(y1+y2s,η)ds.B(y,\eta)=\frac12\int_0^\infty e^{-cs} (1+y^2s)^{-(j+3)/2} b\!\left(\frac{y}{\sqrt{1+y^2s}},\eta\right)\,ds.

The right side is defined for both signs of small yy. On compact parameter sets, every derivative of its integrand is bounded by C(1+s)NecsC(1+s)^N e^{-cs}, for some C,NC,N depending on that derivative. This is integrable, so proves smoothness. Set y=0y=0 to obtain the stated boundary value.

Positive smooth multipliers

If a(y)=ec/y2g(y)a(y)=e^{-c/y^2}g(y) for y>0y>0, with gg smooth and g(0)>0g(0)>0, then

1a(y)0ya(u)ujb(u,η)du\frac1{a(y)}\int_0^y a(u)u^j b(u,\eta)\,du

also equals yj+3y^{j+3} times a smooth coefficient near zero. Apply the formula with gbgb and divide the resulting coefficient by the nonvanishing function g(y)g(y). This conclusion uses the specific exponential profile, not just flatness.