From b7f276a8dd206fd4dbf73997ce5b3558c6b9ba47 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Daniel=20Kn=C3=BCttel?= <daniel.knuettel@daknuett.eu>
Date: Wed, 31 Jul 2019 12:01:10 +0200
Subject: [PATCH] removed r=0 singularity

---
 force.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/force.py b/force.py
index 941ac0d..79cae25 100644
--- a/force.py
+++ b/force.py
@@ -10,7 +10,7 @@ force_function = UFuncWrapper(0, c)
 potential_function = UFuncWrapper(2, c)
 
 # Plot the force and potential.
-r = np.arange(0, 100, 0.02, dtype=np.float16)
+r = np.arange(0.01, 100, 0.02, dtype=np.float16)
 f, = plt.plot(r, force_function(r), label="force")
 p, = plt.plot(r, potential_function(r), label="potential")
 plt.legend(handles=[f, p])