Numerical analysis was executed with the Gaussian 16 code. The B3LYP method has been chosen, following its former successful tests in the electron density calculations in atoms. The aug-cc-pvqz basis set has been routinely applied, except for the Ca atom (cc-pvqz). Using the pVTZ basis set recommended by Sadlej was necessary for the potassium atom. The ground states of atoms have been identified by the specification of atomic electronic terms. The numerical DFT method has been used to calculate the electron density and electron density gradient to overview the local softness of atoms and ions. The effect of degeneracy of frontier orbitals has been avoided by using the integral electron density for every atom and ion.