Write A=(Aijβ) for iβ{1,β¦,m}, jβ{1,β¦,n}, and write Inβ=(Ξ΄ijβ) as in Identity Matrix, where Ξ΄ijβ=1 if i=j and Ξ΄ijβ=0 otherwise.
By Product of Real Matrices, for all iβ{1,β¦,m} and jβ{1,β¦,n},
(AInβ)ijβ=l=1βnβAilβΞ΄ljβ.
Every term with lξ =j vanishes, and the term with l=j equals Aijβ; hence (AInβ)ijβ=Aijβ for all i,j, so AInβ=A.
Likewise,
(ImβA)ijβ=l=1βmβΞ΄ilβAljβ=Aijβ,
so ImβA=A. The square case is the case m=n. β