The theorem states that if A, B are matrices of size p × n and n × p respectively, then
where Ia is the identity matrix of order a.
This theorem is useful in developing a Bayes estimator for multivariate Gaussian distributions.
Sylvester (1857) stated this theorem without proof.