diff --git a/tests/linalg_test.py b/tests/linalg_test.py
index 58301814..23060a60 100644
--- a/tests/linalg_test.py
+++ b/tests/linalg_test.py
@@ -300,6 +300,10 @@ def test_matrix_functions():
         for j in range(dim):
             assert tmp[j].is_zero()
 
+    # Check eig function
+    e2 = pe.linalg.eig(sym)
+    assert np.all(e == e2[::-1])
+
     # Check svd
     u, v, vh = pe.linalg.svd(sym)
     diff = sym - u @ np.diag(v) @ vh