diff --git a/cython/core/triangulation.pyx b/cython/core/triangulation.pyx
index d643cd20..bdd12321 100644
--- a/cython/core/triangulation.pyx
+++ b/cython/core/triangulation.pyx
@@ -3040,6 +3040,7 @@ cdef class Triangulation():
         """
         Extend the Neumann-Zagier Matrix to one which is symplectic (up to factors of 2)
         using oscillating curves. Verify parameter explicitly tests if the resulting matrix is symplectic.
+        Only accepts triangulations with 1 cusp.
 
         >>> M = Manifold("4_1")
         >>> M.symplectic_basis()                    # doctest:
@@ -3081,6 +3082,10 @@ cdef class Triangulation():
         if self.c_triangulation is NULL:
             raise ValueError('The Triangulation is empty.')
 
+        # current get_symplectic_eqns() implementation requires 1 cusp
+        if self.num_cusps() > 1:
+            raise ValueError('Triangulation contains {self.num_cusps()} cusps, only accepts triangulations with 1 cusp')
+
         eqns = []
 
         peripheral_curves(self.c_triangulation)
diff --git a/dev/symplectic_basis/symplectic_basis_main.c b/dev/symplectic_basis/symplectic_basis_main.c
index 3723307b..6b4cb4b3 100644
--- a/dev/symplectic_basis/symplectic_basis_main.c
+++ b/dev/symplectic_basis/symplectic_basis_main.c
@@ -20,8 +20,8 @@ int main(void) {
     int fromFile = 1;
 
     int count = 1;
-    int numTet[] = {6};
-    int index[] = {286};
+    int numTet[] = {5};
+    int index[] = {6};
 
     char *error[] = {
             "CuspedCensusData/link-60819.tri",  /* curve holonomy */
diff --git a/dev/symplectic_basis/symplectic_basis_test.c b/dev/symplectic_basis/symplectic_basis_test.c
index 0e83009a..13471340 100644
--- a/dev/symplectic_basis/symplectic_basis_test.c
+++ b/dev/symplectic_basis/symplectic_basis_test.c
@@ -44,7 +44,7 @@ void testDual(void) {
             if (get_orientability(theTriangulation) == nonorientable_manifold)
                 continue;
 
-//            printf("Num Tet: %d Index: %d \n", index[i][0], j);
+            printf("Num Tet: %d Index: %d \n", index[i][0], j);
 
             basis = get_symplectic_basis(theTriangulation, &dual_rows, &dual_cols, 0);