1D radial proton-electron nucleus model. Inner electron domain [0, Me] (charge magnitude Z, homogenised into a single wavefunction with reduced self-repulsion (Z−1)/Z) symmetrically treated with four outer proton shells P1..P4 (each with its own reduced self-repulsion). All five domains have free Bernoulli boundaries that adjust to enforce u continuity. The attractive potential between electron and each proton shell is computed self-consistently from the converged densities; no static kernel approximation.