  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.
