diff --git a/GHA_triaxial/panou.py b/GHA_triaxial/panou.py
index 9f7c6ca..5b265ff 100644
--- a/GHA_triaxial/panou.py
+++ b/GHA_triaxial/panou.py
@@ -65,7 +65,7 @@ def gha1_num(ell, point, alpha0, s, num):
     werte = rk.rk_chat(F, v_init, s, num)
     x1, _, y1, _, z1, _ = werte[-1]
 
-    return x1, y1, z1
+    return x1, y1, z1, werte
 
 
 def checkLiouville(ell: ellipsoide.EllipsoidTriaxial, points):
@@ -181,11 +181,14 @@ if __name__ == "__main__":
     re = ellipsoide.EllipsoidBiaxial.init_name("Bessel")
 
     # Panou 2013, 7, Table 1, beta0=60°
-    beta0, lamb0, beta1, lamb1, c, alpha0, alpha1, s = nep.get_example(table=1, example=5)
+    beta0, lamb0, beta1, lamb1, c, alpha0, alpha1, s = nep.get_example(table=1, example=3)
     P0 = ell.ell2cart(beta0, lamb0)
     P1 = ell.ell2cart(beta1, lamb1)
 
     # P1_num = gha1_num(ell, P0, alpha0, s, 1000)
     P1_num = gha1_num(ell, P0, alpha0, s, 10000)
     P1_ana = gha1_ana(ell, P0, alpha0, s, 30)
+
+    beta, lamb = ellipsoide.EllipsoidTriaxial.cart2ell(ell, P1_num)
+
     pass
\ No newline at end of file