  In order to study structure of the crust in neutron stars, we develop a
finite-temperature Skyrme-Hartree-Fock method in the full three-dimensional
coordinate space using the Fermion operator expansion method. It provides us
with a possible order-N approach to non-uniform neutron-star matters.
