 That is, we want to construct

where the are small constants that are to be determined.Substituting (4) into (3) and rearranging gives

Because the eigenvectors are -orthogonal when is positive definite, we can remove the summations by left-multiplying by formula_20