From ead73809198598919871e51a011f1f097e6f7e36 Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Tue, 12 Mar 2024 20:00:36 -0400
Subject: [PATCH] New 3 body system

---
 ho_basis_3body.jl | 8 +++-----
 1 file changed, 3 insertions(+), 5 deletions(-)

diff --git a/ho_basis_3body.jl b/ho_basis_3body.jl
index 0eb3d4e..38c84c5 100644
--- a/ho_basis_3body.jl
+++ b/ho_basis_3body.jl
@@ -7,10 +7,8 @@ E_max = 20
 Λ = 0
 m = 1.0
 
-Va = -5
-Ra = sqrt(3)
-Vb = 2
-Rb = sqrt(10)
+Va = -2
+Ra = 2
 
 μ1 = m * 1/2
 μ2 = m * 2/3
@@ -27,7 +25,7 @@ println("Constructing KE matrices")
 @time "T2" T2 = sp_T_matrix(n2s, l2s; ω=ω, μ=c2^2 * μ2)
 
 println("Constructing PE matrices")
-V_l(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω) + Vb * V_Gaussian(Rb, l, n1, n2; ω=ω)
+V_l(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω)
 @time "V1" V1 = sp_V_matrix(V_l, n1s, l1s)
 @time "V relative" V_relative = V1 + sp_V_matrix(V_l, n2s, l2s)
 @time "Moshinsky brackets" U = Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ)