Aiming at microscopically understanding the shape-coexistence phenomen
a, a new algorithm for obtaining many self-consistent Hartree-Fock sta
tes is developed. In contrast with the conventional numerical method o
f the constrained Hartree-Fock equation, which gives the most energeti
cally favorable state under a given constrained condition, it can find
many high-lying Hartree-Fock states as well as many continuous constr
ained Hartree-Fock solutions by dictating their configurations through
some reference state. Numerical calculation is performed by using the
Skyrme III.