The first-order perturbation, c', of an eigenvector c of a symmetric matrix H satisfies the inhomogeneous equation (H-yI)c' = -(H'-y'I)c where I is the unit matrix, y the eigenvalue belonging to c, and the prime denotes the appropriate derivative. The derivative c' can be determined from y, y', c, and a particular solution of the inhomogeneous equation; a convenient computational procedure is demonstrated. In contrast, the conventional perturbation theory requires knowledge of a complete set of eigenvectors of H. Another advantage of the present method is that arithmetic checks are easily applied in numerical work. The method lends itself to the calculation of polarizabilities in the simple molecular-orbital theory. In this connection a complication associated with partly occupied degenerate orbitals is pointed out, the present method is convenient for dealing with this complication.
A more rapid method of reducing a matrix by use of symmetry than that commonly used is presented. It is useful in the present perturbation problem but it should prove equally convenient in other problems.
The determination of c' in certain cases of degeneracy is described and an error in the literature pointed out.
The analysis throughout is general, and may be of interest outside molecular-orbital theory.