Since we are looking for a number of components able to describe
completely the observed variable, one additional assumption
is that the EOFs φk are orthogonal and the PCs Ak vary independently.
Based on these assumptions the minimization process
reduces to solving an eigenvalue/eigenvector problem (secular
equation)