From 44df47ce7756a9ca03b2c447138d2b1c9d1f30c4 Mon Sep 17 00:00:00 2001 From: imran-marwat Date: Thu, 16 Mar 2017 13:12:41 +0000 Subject: [PATCH] Made it R^2 instead of R --- Main | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/Main b/Main index 6b1e4ea..1eb9b07 100644 --- a/Main +++ b/Main @@ -176,8 +176,8 @@ while t < maxTime: iposition_info = (str(particles[i].label) + " " + str(particles[i].position[0]) + " " + str(particles[i].position[1]) + " " + str(particles[i].position[2])) position_info.append(iposition_info) updatePos(particles[i],oldForce[i],dt,t) - R = np.linalg.norm(particles[i].position) - SDTemp.append(R) + R^2 = (np.linalg.norm(particles[i].position))**2 + SDTemp.append(R^2) if l%5== 0: SDContain.append(SDTemp)