diff --git a/ho_basis_3body.jl b/ho_basis_3body.jl
index 70f9fc7..0eb3d4e 100644
--- a/ho_basis_3body.jl
+++ b/ho_basis_3body.jl
@@ -29,9 +29,9 @@ println("Constructing KE matrices")
 println("Constructing PE matrices")
 V_l(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω) + Vb * V_Gaussian(Rb, 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" Moshinsky_mat = Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ)
-@time "V2 via transform" V2 = Moshinsky_mat * V_relative
+@time "V relative" V_relative = V1 + sp_V_matrix(V_l, n2s, l2s)
+@time "Moshinsky brackets" U = Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ)
+@time "V2" V2 =  U' * V_relative * U
 
 println("Calculating spectrum")
 @time "H" H = T1 + T2 + V1 + V2