This solves the Gouy-Chapman-Stern model self-consistently: the compact Stern layer is a
parallel-plate capacitor, the diffuse layer obeys the full nonlinear Poisson-Boltzmann (Grahame)
relation, and the potential at the Stern/diffuse boundary φ_S is the value that makes both layers
carry the same charge magnitude — solved here by bisection on every slider change.
C_Stern (φ0 − φ_S) = √(8 εr ε0 kT c NA) · sinh(zeφ_S / 2kT)
tanh(zeφ(x)/4kT) = tanh(zeφ_S/4kT) · exp(−κ(x−x_S))