  We present a formulation of the Hartree-Fock-Bogoliubov (HFB) equations which
solves the problem directly in the basis of natural orbitals. This provides a
very efficient scheme which is particularly suited for large scale calculations
on coordinate-space grids.
