diff --git a/folie/analysis/overdamped_1d.py b/folie/analysis/overdamped_1d.py
index 68d7f37..274cab0 100644
--- a/folie/analysis/overdamped_1d.py
+++ b/folie/analysis/overdamped_1d.py
@@ -21,10 +21,10 @@ def grad_V(x, _):
         diff_prime_val = model.diffusion.grad_x(x).ravel()
         drift_val = model.drift(x).ravel()
         diff_val = model.diffusion(x).ravel()
-
+        # print(np.min(diff_val), x[np.argmin(diff_val)])
         return (-drift_val + diff_prime_val) / diff_val
 
-    sol = solve_ivp(grad_V, [x.min() - 1e-10, x.max() + 1e10], np.array([0.0]), t_eval=x)  # Add some epsilon to range to ensure inclusion of x
+    sol = solve_ivp(grad_V, [x.min() - 1e-10, x.max() + 1e-10], np.array([0.0]), t_eval=x)  # Add some epsilon to range to ensure inclusion of x
 
     V = sol.y.ravel()