Write β₯β
β₯ for the Euclidean norm and β£β
β£ for the absolute value. Claims 1 and 2 of Elementary Properties of the Euclidean Norm on Rn give β₯zβ₯2=zβ
z=βi=1nβziβziβ and β₯zβ₯=dEβ(z,0Rnβ) for the Euclidean distance and the origin 0Rnβ, and claim 5 gives β₯ΞΌzβ₯=β£ΞΌβ£β₯zβ₯. Claims 2 and 4 of Properties of the Absolute Value in an Ordered Field are used as symmetry and multiplicativity, claims 2 and 3 of Properties of Finite Sums as additivity and homogeneity of finite sums, and claims 1 and 2 of Monotonicity of Squaring on the Nonnegative Elements of an Ordered Field as the strict and weak comparison of squares of nonnegative elements. Order arithmetic is taken from Elementary Order Arithmetic in an Ordered Field (claim 5 product of positive elements, claim 7 inverse of a positive element, claim 8 halving, claim 10 strict compatibility with multiplication by a positive element). The field axioms of the field R and the vector space axioms of Euclidean Space Rn is a Real Vector Space are used for rearrangement; in particular
p+(t0β+s)h=x0β+sh(sβR).
Claim 1. Put L=βi=1nββxiββfβ(x0β)hiβ. By Gradient of a Real-Valued Function on a Euclidean Open Set the ith coordinate of Df(x0β) is βxiββfβ(x0β), so L=hβ
Df(x0β) by Difference, Dot Product, and Orthogonality in Rn and commutativity of multiplication. We show that L has the property required of gβ²(t0β).
Suppose first that h is the origin. Then p+th=p for every t, so g is constant and every difference quotient of g vanishes; moreover hiβ=0 for every i, so homogeneity of finite sums with the factor 0 gives L=0. Hence for any Ξ΅ with 0<Ξ΅ the choice Ξ΄=1 works.
Assume now that h is not the origin, so that 0<β₯hβ₯ by the axioms of Metric Space, the Euclidean distance being a metric by Euclidean Distance is a Metric on Rn. By C^1 Maps on Euclidean Open Sets are Differentiable the function f, regarded as a map into Rm with m=1, is differentiable at x0β. For kβRn with x0β+kβU write
R(k)=f(x0β+k)βf(x0β)βi=1βnββxiββfβ(x0β)kiβ.
The defining condition of differentiability reads: for every Ξ΅β² with 0<Ξ΅β² there is Ξ΄0β with 0<Ξ΄0β such that 0<βi=1nβki2β<Ξ΄02β and x0β+kβU imply R(k)2β€Ξ΅β²2βi=1nβki2β. Since βi=1nβki2β=β₯kβ₯2, since R(k)2=β£R(k)β£2 by multiplicativity of the absolute value, and since Ξ΅β²2β₯kβ₯2=(Ξ΅β²β₯kβ₯)2, the weak comparison of squares of nonnegative elements turns that conclusion into
β£R(k)β£β€Ξ΅β²β₯kβ₯.
Let Ξ΅ with 0<Ξ΅ be given. Put Ξ΅β²=21βΞ΅β₯hβ₯β1, which is positive, take Ξ΄0β as above for this Ξ΅β², and put Ξ΄=Ξ΄0ββ₯hβ₯β1, which is positive. Let sβR satisfy 0<β£sβ£<Ξ΄ and t0β+sβJ, and put k=sh. Then β₯kβ₯=β£sβ£β₯hβ₯ is positive, and multiplying β£sβ£<Ξ΄ by the positive element β₯hβ₯ gives β₯kβ₯<Ξ΄0β; the strict comparison of squares then gives 0<βi=1nβki2β<Ξ΄02β. Moreover x0β+k=p+(t0β+s)hβU because t0β+sβJ. Homogeneity of finite sums gives
i=1βnββxiββfβ(x0β)kiβ=sL,
so g(t0β+s)βg(t0β)=f(x0β+k)βf(x0β)=R(k)+sL and therefore
sg(t0β+s)βg(t0β)ββL=sR(k)β.
Using multiplicativity of the absolute value, β£R(k)/sβ£=β£R(k)β£β£sβ£β1β€Ξ΅β²β₯kβ₯β£sβ£β1=Ξ΅β²β₯hβ₯=21βΞ΅<Ξ΅, the last step by halving. Hence g is differentiable at t0β with gβ²(t0β)=L.
Claim 2. Let jβ{1,β¦,n}. By clause 2 of C^2 Real-Valued Map on an Open Subset of Euclidean Space the function βf/βxjβ is of class C1 on U, so claim 1 applied to it in place of f shows that the function Οjβ:JβR with Οjβ(t)=βxjββfβ(p+th) is differentiable at t0β with
Οjβ²β(t0β)=i=1βnββxiββxjββ2fβ(x0β)hiβ=i=1βnβ(D2f(x0β))ijβhiβ,
the second equality by Hessian Matrix of a C^2 Function.
By commutativity of multiplication, g1β(t)=βj=1nβhjβΟjβ(t) for every tβJ, so Derivative of a Finite Linear Combination of Real Functions, applied with the constants hjβ, shows that g1β is differentiable at t0β with
g1β²β(t0β)=j=1βnβhjβΟjβ²β(t0β)=j=1βnβΒ i=1βnβ(D2f(x0β))ijβhiβhjβ,
the inner factor hjβ being moved inside the inner sum by homogeneity of finite sums. Interchanging the two summations by Interchange of a Finite Double Sum and applying claim 4 of Linearity of the Matrix-Vector Product and the Quadratic Form as a Double Sum with M=D2f(x0β) and w=z=h gives
g1β²β(t0β)=i=1βnβΒ j=1βnβ(D2f(x0β))ijβhiβhjβ=hβ
(D2f(x0β)h).