Empirically determined values of the nuclear volume and surface symmetry energy coefficients from nuclear masses are expressed in terms of density distributions of nucleons in heavy nuclei in the local density approximation. This is then used to extract the value of the symmetry energy slope parameter L. The density distributions in both spherical and well deformed nuclei calculated within microscopic framework with different energy density functionals give L = 59.0 ± 13.0 MeV. Application of the method also helps in a precision determination of the neutron skin thickness of nuclei that are difficult to measure accurately. The nuclear symmetry energy measures the energy transfer in converting symmetric nuclear matter to the asymmetric one. The density dependence gives information on the isospin-dependent part of the equation of state (EOS) of asymmetric nuclear matter. The density content of symmetry energy is mostly encoded in the symmetry energy coefficient C v , the symmetry slope parameter L and the symmetry incompressibility K sym , all evaluated at the nuclear saturation density ρ 0 . Here C v (ρ) is the volume symmetry energy per nucleon of homogeneous nuclear matter at density ρ,andIn Eq. (1), e is the energy per nucleon and X = (ρ n − ρ p )/(ρ n +ρ p ) is the isospin asymmetry of the system. The parameters C v (ρ 0 ), L and K sym deem to be of fundamental importance in both nuclear physics and astrophysics. The nuclear binding energies, the position of the nuclear drip lines, the neutron skin thickness or the neutron density distribution in neutron-rich nuclei, − all of these are known to have affectations from [1-8] the symmetry parameters so mentioned. Many astrophysical phenomena also depend sensitively on the symmetry slope parameter L. Most notable among them are the dynamical evolution of the core collapse of a massive star and the associated explosive nucleosynthesis [9, 10], the cooling of proto-neutron stars through neutrino convection [11] or the radii of cold neutron stars [11]. The nature and stability of phases within a neutron star, its crustal com-[12] also seem to be strongly influenced by the symmetry energy and its density dependence. A glimmer of their import could further be seen in relation to some issues of new physics beyond the standard model [13,14]. Attempts on estimates of L have been made in the last few years from analyses of diverse experimental data, uncertainties still linger, however. Pygmy dipole resonance [15] in 68 N i and 132 Sn predicts a weighted average in the range L=64.8 ± 15.7 MeV, but giant dipole resonance in 208 P b [16] points to a value of L ∼ 52 ±7 MeV. Nucleon emission ratios from heavy ion collisions [17] favor a value close to it, L ∼ 55 MeV, but isoscaling gives L ∼ 65 MeV [18] and isospin diffusion shifts the value further up, L = 88 ± 25 MeV [19,20]. Recently, analyzing nuclear energies within the standard Skyrme-Hartree-Fock approach, Chen [21] brought down L to 52.5 ± 20 MeV. However, from the fit [3] of the experimental nuclear masses to the calc...