By Transpose of a Real Matrix, Aβ€ is the real nΓm matrix with (Aβ€)jiβ=Aijβ for iβ[m] and jβ[n], and by Real Matrix and the Set of Real Matrices two matrices of the same size are equal exactly when all of their entries agree. Sums below are the finite sums of the field of real numbers.
Claim 1. For iβ[m] and jβ[n], ((Aβ€)β€)ijβ=(Aβ€)jiβ=Aijβ.
Claim 2. By Sum of Real Matrices, ((A+B)β€)jiβ=(A+B)ijβ=Aijβ+Bijβ=(Aβ€)jiβ+(Bβ€)jiβ=(Aβ€+Bβ€)jiβ. By Scalar Multiple of a Real Matrix, ((ΞΌA)β€)jiβ=(ΞΌA)ijβ=ΞΌAijβ=ΞΌ(Aβ€)jiβ=(ΞΌAβ€)jiβ.
Claim 3. By Identity Matrix the entry (Inβ)ijβ is 1 if i=j and 0 otherwise, a condition unchanged when i and j are interchanged. Hence (Inβ€β)jiβ=(Inβ)ijβ=(Inβ)jiβ.
Claim 4. Let iβ[m] and jβ[p]. By Product of Real Matrices,
((AC)β€)jiβ=(AC)ijβ=k=1βnβAikβCkjβ,(Cβ€Aβ€)jiβ=k=1βnβ(Cβ€)jkβ(Aβ€)kiβ=k=1βnβCkjβAikβ.
The two families of summands agree by commutativity of multiplication in R, so the two sums agree.
Claim 5. By Difference, Dot Product, and Orthogonality in Rn and Matrix-Vector Product, and then by claim 3 of Properties of Finite Sums applied to each inner sum,
wβ
(Av)=i=1βmβwiβ(Av)iβ=i=1βmβwiβj=1βnβAijβvjβ=i=1βmβj=1βnβAijβwiβvjβ,
using commutativity and associativity of multiplication to arrange each summand. In the same way
(Aβ€w)β
v=j=1βnβ(Aβ€w)jβvjβ=j=1βnβi=1βmβAijβwiβvjβ.
Interchanging the order of summation in the second double sum by Interchange of a Finite Double Sum shows that the two are equal.
Claim 6. By Inverse Matrix and Invertible Real Square Matrix we have AAβ1=Inβ and Aβ1A=Inβ. Applying the transpose to each identity and using claims 4 and 3,
(Aβ1)β€Aβ€=Inβ€β=Inβ,Aβ€(Aβ1)β€=Inβ€β=Inβ.
Hence Aβ€ is invertible with inverse (Aβ1)β€, and this inverse is the only one by Uniqueness of the Matrix Inverse.