This second-order differential equation is solved numerically for varying electrode potentials. Numerical solutions provide the electric potential, electric field, dielectric constant, and ion distribution as functions of the distance from the electrode. The model's results depend on the parameters defined in the free energy expression, with detailed parameter information provided in the supplementary information.