diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh
index e0ac46af0995e35a4e98f430c7a6718201130cb3..28198eff3ecfb75a72a6c9f54c43e77f335e4f25 100644
--- a/dune/gfe/rotation.hh
+++ b/dune/gfe/rotation.hh
@@ -143,7 +143,7 @@ class Rotation<T,3> : public Quaternion<T>
 
     /** \brief Computes sin(x/2) / x without getting unstable for small x */
     static T sincHalf(const T& x) {
-        return (x < 1e-4) ? 0.5 - (x*x/48) : std::sin(x/2)/x;
+        return (x < 1e-4) ? 0.5 - (x*x/48)  + (x*x*x*x)/3840 : std::sin(x/2)/x;
     }
 
 public: