We discuss the HFB equations in coordinate representation,a suitable method for handling the full effects of the continuous quasiparticle spectrum. We show how the continuum HFB equations can be solved with the correct asymptotic conditions instead of the discretization conditions which are commonly used in the literature. The continuum HFB method is illustrated with a model where the mean field and pairing field have simple forms. The relationship with the continuum Hartree-Fock-BCS (HF-BCS) approximation is also discussed. Realistic HFB and HF-BCS calculations based on Skyrme interactions are compared for the case of a neutron-rich nucleus.