the meson fields are obtained, all the quantities that depend on them, as the energy, pressure densities, chemical potentials, baryonic densities, etc are also computed. The solution is a droplet with a certain proton fraction surrounded by a gas of neutrons. If stable nuclei are calculated, the gas vanishes because the energy of the system lies below the neutron drip line and the finite nuclei properties are easily calculated. This is the general method, but the results depend strongly on the model used because of the reasons discussed in Section VI.

  --------------------------------------------------------------------------------------------------------------------------------- --
   ![Binodal section for the NL3, TW and NL$\omega\rho$ parametrizations.[]{data-label="binodais"}](fig1.eps "fig:"){width="8.cm"}  
  --------------------------------------------------------------------------------------------------------------------------------- --

Numerical result within a harmonic oscillator basis
===================================================

Here a different prescription for solving the equations of motion and the thermodynamical quantities within the Thomas-Fermi approximation is used. According to [@ring], meson field equations of motion of the Klein-Gordon type with sources can be carried out by an expansion in a complete set of basis states. The harmonic oscillator functions with orbital angular momentum equal to zero are then chosen. The oscillator length is given by $$b_B=\frac{b_0}{\sqrt 2}, \quad b_0=\sqrt{\frac{\hbar}{M \omega_0}},$$ where $M$ is the nucleon mass and $\omega_0$ is the oscillator frequency. The meson 