diff --git a/dune/gfe/hyperbolichalfspacepoint.hh b/dune/gfe/hyperbolichalfspacepoint.hh
index 96ef67f4b14b7fab2c57dd2a6fac86ccf38e3541..aa77c3ce2a5226a000f75b9bd78af8871960bbfa 100644
--- a/dune/gfe/hyperbolichalfspacepoint.hh
+++ b/dune/gfe/hyperbolichalfspacepoint.hh
@@ -387,10 +387,7 @@ public:
                 
                     } else if (i==N-1 and j==N-1 and k==N-1) {
                     
-                        dFdqdqdq[i][j][k] = -2.0/Dune::Power<3>::eval(q[N-1]) 
-                                          - (2*p[N-1]*p[N-1]*q[N-1] - p[N-1]*q[N-1]*q[N-1]) / (p[N-1]*p[N-1]*Dune::Power<4>::eval(q[N-1]))
-                                          + 2 * (p[N-1]-q[N-1]) / (p[N-1]*Dune::Power<3>::eval(q[N-1])) 
-                                          - 3 * diffNormSquared / (p[N-1]*Dune::Power<4>::eval(q[N-1]));
+                        dFdqdqdq[i][j][k] = (2.0*p[N-1] -1 -4*p[N-1] -3*diffNormSquared)/(p[N-1]*Dune::Power<3>::eval(q[N-1])); 
                 
                     }