diff --git a/ho_basis.jl b/ho_basis.jl
index 44c8427..b2dc916 100644
--- a/ho_basis.jl
+++ b/ho_basis.jl
@@ -35,7 +35,7 @@ function get_2p_basis(E_max)
     l2s = Int[]
 
     # E = 2*n1 + l1 + 2*n2 + l2
-    for E in 0 : 2 * E_max
+    for E in 0 : E_max
         for n1 in 0 : E ÷ 2
             for n2 in 0 : (E - 2*n1) ÷ 2
                 for l1 in 0 : (E - 2*n1 - 2*n2)
diff --git a/ho_basis_3body.jl b/ho_basis_3body.jl
index ae7f221..7ea84e6 100644
--- a/ho_basis_3body.jl
+++ b/ho_basis_3body.jl
@@ -2,7 +2,7 @@ using Arpack, SparseArrays
 include("ho_basis.jl")
 include("p_space.jl")
 
-E_max = 20
+E_max = 40
 μ1ω1_gen = 1/2
 μ2ω2_gen = 2
 Λ = 0
diff --git a/test/moshinsky.jl b/test/moshinsky.jl
index 1406188..772b6cd 100644
--- a/test/moshinsky.jl
+++ b/test/moshinsky.jl
@@ -4,7 +4,7 @@ using LinearAlgebra
 
 include("../ho_basis.jl")
 
-E_max = 15
+E_max = 30
 Λ = 0
 
 println("No of threads = ", Threads.nthreads())