diff --git a/dune/gfe/hyperbolichalfspacepoint.hh b/dune/gfe/hyperbolichalfspacepoint.hh
index 307ad7f20a932ef7f49672e42df35e274864b3bb..192236b0bad7ed72d2ca8f8b598f8084bda70c96 100644
--- a/dune/gfe/hyperbolichalfspacepoint.hh
+++ b/dune/gfe/hyperbolichalfspacepoint.hh
@@ -107,7 +107,9 @@ class HyperbolicHalfspacePoint
                     
                 } else if (i==N-1 and j==N-1) {
                     
-                    dFdpdq[i][j] = -1/(p[N-1]*p[N-1]*q[N-1]) - (p[N-1]-q[N-1]) / (p[N-1]*q[N-1]*q[N-1]) + diffNormSquared / (p[N-1]*p[N-1]*q[N-1]*q[N-1]);
+                    dFdpdq[i][j] = -1/(p[N-1]*q[N-1]) 
+                                   - (p[N-1]-q[N-1]) / (p[N-1]*q[N-1]*q[N-1]) 
+                                   + (p[N-1]-q[N-1]) / (p[N-1]*p[N-1]*q[N-1]) + diffNormSquared / (2*p[N-1]*p[N-1]*q[N-1]*q[N-1]);
                 
                 }