Reason: First published version. Restricts to an affine path, forms the auxiliary function whose value at 1 is the Taylor remainder, applies the mean value theorem twice, and estimates using continuity of the second partial derivatives and symmetry of the Hessian.
Step 1 (a path interval). Since x∈U and U is open, there is r with 0<r such that every point at Euclidean distance less than r from x lies in U. Suppose h∈Rn satisfies ∥h∥<r⋅2−1. For τ∈R with ∣τ∣<2, the point x+τh is at distance ∣τ∣∥h∥ from x, since the squared distance is ∑i(τhi)2=τ2∥h∥2=(∣τ∣∥h∥)2 and both quantities are nonnegative. By claim 10 this is less than 2⋅(r⋅2−1)=r when ∥h∥=0, and equals 0<r otherwise. Hence x+τh∈U for every τ in the open interval(−2,2), which contains 0 and 1.
Define on (−2,2) the functions F(τ)=f(x+τh) and Φi(τ)=∂if(x+τh). By Chain Rule Along an Affine Path, applied to f and to each ∂if, these are differentiable at every τ∈(−2,2), with
Step 3 (an auxiliary function). Put Q=∑i=1n∑j=1n∂ij2f(x)hihj and define G:(−2,2)→R by
G(τ)=F(τ)−F(0)−τF′(0)−τ2Q⋅2−1.
The map τ↦F(0)+τF′(0)+τ2Q⋅2−1 is differentiable at each τ0 with derivative F′(0)+τ0Q: its difference quotient at τ0 for an increment k=0 equals F′(0)+(2τ0+k)Q⋅2−1, which differs from F′(0)+τ0Q by kQ⋅2−1, and given ε with 0<ε one may take δ to be ε if Q=0 and otherwise the value ε(∣Q∣⋅2−1)−1, positive by claims 5 and 7, so that ∣kQ⋅2−1∣=∣k∣∣Q∣⋅2−1<ε by claims 4 of Properties of the Absolute Value in an Ordered Field and 10. By Derivative of a Sum and of a Difference, G is therefore differentiable at every τ0∈(−2,2) with
G′(τ0)=F′(τ0)−F′(0)−τ0Q.
Also G(0)=0, and G(1) is exactly the quantity whose absolute value is to be estimated, since F(1)=f(x+h), F(0)=f(x) and F′(0)=∑i∂if(x)hi.
Step 4 (two applications of the mean value theorem). By Mean Value Theorem on an Open Interval applied to G on (−2,2) with the points 0<1, there is θ∈(0,1) with G(1)−G(0)=G′(θ), so G(1)=G′(θ).
For each i, apply Mean Value Theorem on an Open Interval to Φi on (−2,2) with the points 0<θ: there is θi∈(0,θ) with Φi(θ)−Φi(0)=Φi′(θi)θ. Hence, using Step 2 and F′(θ)−F′(0)=∑i(Φi(θ)−Φi(0))hi,
Step 5 (the estimate). Let ε∈R with 0<ε. Since n≥1, the element n⋅n is positive (a sum of copies of 1, positive by claim 6 and claim 3), so ε′=ε(n⋅n)−1 is positive by claims 7 and 5. Each ∂ji2f is continuous at x, by condition 2 of C^2 Real-Valued Map on an Open Subset of Euclidean Space; taking the least of the finitely many radii by repeated use of claim 9, there is θ0 with 0<θ0 such that every z∈U at distance less than θ0 from x satisfies ∣∂ji2f(z)−∂ji2f(x)∣<ε′ for all i,j.
Let δ be the least, by claim 9, of r⋅2−1 and θ0, so 0<δ. Suppose ∥h∥<δ. Then x+h∈U by Step 1 with τ=1, and Steps 1 to 4 apply. For each i we have 0<θi<θ<1, so the distance from x+θih to x is θi∥h∥≤∥h∥<θ0 by claim 10 and claim 2; hence ∣Rji(θi)∣<ε′ for all i,j.