diff --git a/ComplexScaling-FV-bound.jl b/ComplexScaling-FV-bound.jl
index 2208c61..f852c44 100644
--- a/ComplexScaling-FV-bound.jl
+++ b/ComplexScaling-FV-bound.jl
@@ -9,12 +9,13 @@ d = 3
 n = 2
 N = 30
 L = 6
+n_imag = 1
 
 open("ComplexScaling-FV-bound.dat", "w") do f
     for ϕ = range(0.0, 0.5, length=11)
         println("Calculating ϕ=", ϕ)
         s = system{T}(d, n, N, L)
-        H = Hamiltonian{T}(s, V_gauss, ϕ, 0, mode)
+        H = Hamiltonian{T}(s, V_gauss, ϕ, n_imag, mode)
         @time evals, _, info = eig(H, 10, resonances = false)
 
         dataline = vcat([ϕ], hcat(real.(evals), imag.(evals))'[:])
diff --git a/ComplexScaling-FV-res.jl b/ComplexScaling-FV-res.jl
index e3b1cfe..2d7f869 100644
--- a/ComplexScaling-FV-res.jl
+++ b/ComplexScaling-FV-res.jl
@@ -9,12 +9,13 @@ d = 3
 n = 2
 N = 96
 ϕ = pi/6
+n_imag = 1
 
 open("ComplexScaling-FV-res.dat", "w") do f
     for L = range(20, 35, length=16)
         println("Calculating L=", L)
         s = system{T}(d, n, N, L)
-        H = Hamiltonian{T}(s, V_gauss, ϕ, 0, mode)
+        H = Hamiltonian{T}(s, V_gauss, ϕ, n_imag, mode)
         @time evals, _, info = eig(H, 40)
 
         dataline = vcat([L], hcat(real.(evals), imag.(evals))'[:])