diff --git a/test_jacobi.jl b/test_jacobi.jl
new file mode 100644
index 0000000..ab8e320
--- /dev/null
+++ b/test_jacobi.jl
@@ -0,0 +1,21 @@
+include("Hamiltonian.jl")
+
+T=Float32
+
+function V_test(r2)
+    return -4*exp(-r2/4)
+end
+
+for (n,N) in [(2,16), (3,8)]
+    println("\n$n-body system with N=$N")
+    n_image=0
+    for L::T in 5.0:9.0
+        print("L=$L: ")
+        s=system{T}(3,n,N,L)
+        print("Constructing H...")
+        H=Hamiltonian{T}(s,V_test,0.0,n_image,cpu_tensor)
+        print("Diagonalizing...")
+        evals,_,_ = eig(H,5)
+        println(real.(evals))
+    end
+end