Perturbation theory is derived in such a way that the wave function renormalization terms appear automatically. The derivation is based on the Feynman theorem for the rate of change of an eigenvalue with respect to a parameter, and a corresponding theorem for the rate of change of an eigenfunction.