N1s = (1/(Sqrt[π]))
N2p = (1/(4 (Sqrt[2 π])))
V0p = (1/3) N1s N2p 2 3 4 ((((2/3)))^5) 2 π 2
EL = (1/10)
dE = -(1/8) + (1/2)
cpz = -((EL V0p)/dE)
s = N1s (E^(-(Sqrt[(x^2) + (y^2)])))
p = 
      N2p x (E^(-(1/2) (Sqrt[(x^2) + (y^2)])))
f = s + cpz p
ContourPlot[s, {x, -2, 2}, {y, -2, 2}, 
      PlotPoints -> 60]
ContourPlot[f, {x, -2, 2}, {y, -2, 2}, 
      PlotPoints -> 60]