diff --git a/dune/gfe/unitvector.hh b/dune/gfe/unitvector.hh
index dedc20d1f9f9c2c2b8b25a1c64b1b724fa9c9f3d..de6df3f9e5d209207204778359adca342f71a7f1 100644
--- a/dune/gfe/unitvector.hh
+++ b/dune/gfe/unitvector.hh
@@ -15,7 +15,7 @@ class UnitVector
 {
     /** \brief Computes sin(x/2) / x without getting unstable for small x */
     static double sinc(const double& x) {
-        return (x < 1e-4) ? 1 + (x*x/6) : std::sin(x)/x;
+        return (x < 1e-4) ? 1 - (x*x/6) : std::sin(x)/x;
     }
 
     /** \brief Compute the derivative of arccos^2 without getting unstable for x close to 1 */