diff --git a/src/Rod2d.cpp b/src/Rod2d.cpp
index cf29906..a73b636 100644
--- a/src/Rod2d.cpp
+++ b/src/Rod2d.cpp
@@ -19,11 +19,11 @@ Rod2d::Rod2d(double L) : m_pos({0, 0}), m_e({1, 0}) {
         const double D0 = 3.0 * M_Pl / L;
         const double p = L;
         auto D_para =
-            D0 / (M_Pl * 2.0) * (log(p) - 0.207 + 0.980 / p - 0.133 / (p * p));
+            D0 / (M_Pl * 2.0) * (log(p) - 0.1401 + 1.034 / p - 0.228 / (p * p));
         auto D_ortho =
-            D0 / (M_Pl * 2.0) * (log(p) + 0.839 + 0.185 / p + 0.233 / (p * p));
+            D0 / (M_Pl * 4.0) * (log(p) + 0.8369 + 0.5551 / p - 0.06066 / (p * p));
         m_D_rot = 3 * D0 / (M_Pl * L * L) *
-                  (log(p) - 0.662 + 0.917 / p - 0.050 / (p * p));
+                  (log(p) - 0.3512 + 0.7804 / p - 0.09801 / (p * p));
         m_Diff << D_para, 0, 0, D_ortho;
     }
     m_Diff_sqrt = m_Diff;