diff --git a/test/Pb208_nucleon_dens.jl b/test/Pb208_nucleon_dens.jl
index f869a79..a95c7cb 100644
--- a/test/Pb208_nucleon_dens.jl
+++ b/test/Pb208_nucleon_dens.jl
@@ -18,7 +18,7 @@ A_interp = linear_interpolation(xs, As)
 N_p = 82
 N_n = 126
 r_max = maximum(xs)
-E_min = 880
+E_min = 860
 E_max = 939
 divs = 400
 
diff --git a/test/Pb208_nucleon_spectrum.jl b/test/Pb208_nucleon_spectrum.jl
index 095f494..0489162 100644
--- a/test/Pb208_nucleon_spectrum.jl
+++ b/test/Pb208_nucleon_spectrum.jl
@@ -16,9 +16,9 @@ R_interp = linear_interpolation(xs, Rs)
 A_interp = linear_interpolation(xs, As)
 
 p = true
-N = 82
+N = p ? 82 : 126
 r_max = maximum(xs)
-E_min = 880
+E_min = 860
 E_max = 939
 
 κs, Es = findAllOrbitals(p, S_interp, V_interp, R_interp, A_interp, r_max, E_min, E_max)