diff --git a/test/Linear_nucleon_spectrum.jl b/test/Linear_nucleon_spectrum.jl
index e1b52fb..5714e9c 100644
--- a/test/Linear_nucleon_spectrum.jl
+++ b/test/Linear_nucleon_spectrum.jl
@@ -10,20 +10,20 @@ Vs = test_data[:, 3]
 Rs = test_data[:, 4]
 As = test_data[:, 5]
 
-S_interp = linear_interpolation(xs, Ss)
-V_interp = linear_interpolation(xs, Vs)
-R_interp = linear_interpolation(xs, Rs)
-A_interp = linear_interpolation(xs, As)
-
 p = false
 r_max = maximum(xs)
 divs = length(xs) - 1
 s = system(8, 8, r_max, divs)
 
+s.Φ0 = Ss
+s.W0 = Vs
+s.B0 = Rs
+s.A0 = As
+
 E_min = 860
 E_max = 939
 
-κs, Es = findAllOrbitals(p, S_interp, V_interp, R_interp, A_interp, s, E_min, E_max)
+κs, Es = findAllOrbitals(p, s, E_min, E_max)
 Ebinds = (p ? M_p : M_n) .- Es
 occ = fillNucleons(Z_or_N(s, p), κs, Es)