A0 15*a0_const "Simulation domain radius" kC 1/(4*pi*epsilon0_const) "Coulomb constant" h0 0.2*a0_const "Mesh size at origin" hg 1/(2*a0_const) "Mesh growth scaling" RH me_const*e_const^4/(8*h_const^2*epsilon0_const^2)/(1[eV]) "Rydberg constant in units of eV" E1 -RH "n=1 eigenenergy (eV)" E2 -RH/4 "n=2 eigenenergy (eV)" Eext 0.01[V]/a0_const "External electric field" E2StarkLow E2-3*e_const*Eext*a0_const/1[eV] "Stark splitting low eigenenergy for n=2 (eV)" E2StarkHigh E2+3*e_const*Eext*a0_const/1[eV] "Stark splitting high eigenenergy for n=2 (eV)"