diff --git a/test/3body_Berggren_orthogonality.jl b/test/3body_Berggren_orthogonality.jl
new file mode 100644
index 0000000..82aca8a
--- /dev/null
+++ b/test/3body_Berggren_orthogonality.jl
@@ -0,0 +1,8 @@
+include("../berggren_3body_resonance.jl")
+
+@time "Eigenvectors" evals, evecs = eigs(H, sigma=target, maxiter=5000, tol=1e-5, ritzvec=true, check=1)
+
+weights_mat = spdiagm(repeat(kron(ws, ws), jmax + 1))
+N = transpose(evecs) * weights_mat * evecs
+
+display(abs.(N))
\ No newline at end of file