With all lower states A i known, this eigenvalue equation can be solved iteratively together with the expectation constraint, the normalization constraint, and the orthogonality constraints, which will give solutions of f n , g ni , |A n ⟩, and Ẽn . When all constraints are satisfied, the total energy will serve as the effective PES for the nth excited constrained minimized energy state: