diff --git a/src/solver/newton/newton.hpp b/src/solver/newton/newton.hpp
index c0fedd3af..8cff8bd3d 100644
--- a/src/solver/newton/newton.hpp
+++ b/src/solver/newton/newton.hpp
@@ -35,7 +35,7 @@ namespace newton {
  */
 
 static constexpr int MAX_ITER = 50;
-static constexpr double EPS = 1e-12;
+static constexpr double EPS = 1e-13;
 
 template <int N>
 bool is_converged(const Eigen::Matrix<double, N, 1>& X,