Claim 1. Since J is positive definite and z is nonzero, zβ
(Jz)>0 by the definition of positive definiteness.
Two elementary identities. Let M be a real matrix with k rows and k columns and let u,vβRk. By the index formulas for the matrix-vector product and the dot product,
uβ
(Mv)=Ξ³=1βkβΞ΄=1βkβMΞ³Ξ΄βuΞ³vΞ΄,
from which two facts follow. First, the expression is linear in u for fixed v and linear in v for fixed u. Second, if M is symmetric, so that Mδγβ=MΞ³Ξ΄β by the definition of the transpose, then exchanging the names of the two summation indices gives uβ
(Mv)=vβ
(Mu).
Claim 2. Fix xβRk. By claim 1 the real number ΞΈ=(xβ
z)/(zβ
(Jz)) is well defined; put y=xβΞΈJz. The matrix Jβ1 is symmetric positive definite, hence positive semidefinite, so yβ
(Jβ1y)β₯0. Expanding by the bilinearity just recorded,
yβ
(Jβ1y)=xβ
(Jβ1x)βΞΈxβ
(Jβ1(Jz))βΞΈ(Jz)β
(Jβ1x)+ΞΈ2(Jz)β
(Jβ1(Jz)).
By the defining property of the inverse, Jβ1J=Ikβ with Ikβ the identity matrix, and the matrix product satisfies (Jβ1J)z=Jβ1(Jz) by the associativity of the index sums; hence Jβ1(Jz)=z. The second term is therefore βΞΈ(xβ
z) and the fourth is ΞΈ2(Jz)β
z=ΞΈ2zβ
(Jz), the dot product being symmetric. For the third, the second identity above applied to the symmetric matrix Jβ1 gives (Jz)β
(Jβ1x)=xβ
(Jβ1(Jz))=xβ
z, so the third term is βΞΈ(xβ
z) as well. Consequently
0Β β€Β yβ
(Jβ1y)Β =Β xβ
(Jβ1x)β2ΞΈ(xβ
z)+ΞΈ2zβ
(Jz)Β =Β xβ
(Jβ1x)βzβ
(Jz)(xβ
z)2β,
the last equality by the choice of ΞΈ. Multiplying by zβ
(Jz)>0 gives claim 2.
Claim 3. Put
M=Jβ1βzβ
(Jz)1β(zβz),
a difference of symmetric matrices formed entrywise, hence symmetric. For xβRk the index formulas give
xβ
((zβz)x)=Ξ³=1βkβΞ΄=1βkβzΞ³zΞ΄xΞ³xΞ΄=(xβ
z)2,
so that, by the linearity in the matrix argument evident from the same formula and by claim 2,
xβ
(Mx)=xβ
(Jβ1x)βzβ
(Jz)(xβ
z)2βΒ β₯Β 0.
Thus M is positive semidefinite, which by the definition of the semidefinite order is exactly claim 3.