We present the first set of results of solving the Hartree-Fock-Bogoliubov equations, which describe the self-consistent mean field theory with pairing interaction. Calculations for even-even nuclei are carried out on a two-dimensional axially symmetric lattice, in coordinate space. An important aspect of our method is the proper representation of the quasi-particle continuum wavefunctions, which are considered for energies up to 60 MeV. This stage is essential for a proper description of nuclei near the drip lines, due to the strong coupling between weakly bound states and the particle continuum for such nuclei. High accuracy is achieved by representing the operators and wavefunctions using the technique of basis-splines. Calculations for Sn isotopes are presented to demonstrate the reliability of the method.