Because the difference between h αβ and the gradient (g α = h αα ) vectors is the element of the differentiated symmetrized effective Hamiltonian, the algorithm for computing h αβ is almost equal to that for the gradient vector, and the algorithm in the previous section is used without additional complications.