From a3ff2834c46b9cea3217d8515d7dc081f5d927c8 Mon Sep 17 00:00:00 2001
From: Hezy Amiel <hezy.amiel@gmail.com>
Date: Tue, 16 Jun 2020 15:02:52 +0300
Subject: [PATCH] 	modified:   vdP.py

---
 vdP.py | 31 ++++++++++++++++++-------------
 1 file changed, 18 insertions(+), 13 deletions(-)

diff --git a/vdP.py b/vdP.py
index 046f760..63799fc 100644
--- a/vdP.py
+++ b/vdP.py
@@ -4,31 +4,36 @@
 @author: Hezy Amiel
 """
 
+import numpy as np
+import matplotlib.pyplot as plt
+from scipy import optimize
+
+
 """
 find the sheet resistance Rs from R1 = R_AB,CD and R2 = R_BC,AD
-the solution is from the equation in a paper L. J. van der Pauw "A method of measuring the resistivity and hall coefficients of Lamellae of arbitrary shape", Philips Technical Review, Vol 26, 220.
+the solution is from the equation in a paper L. J. van der Pauw
+"A method of measuring the resistivity and hall coefficients of
+Lamellae of arbitrary shape", Philips Technical Review, Vol 26, 220.
 """
 
-import numpy as np
-import matplotlib.pyplot as plt
-from scipy import optimize
+
+# define the equation
+def fun(Rs):
+    return np.exp(-np.pi*R1/Rs) + np.exp(-np.pi*R2/Rs) - 1
+
 
 # these are the resistance R_AB,CD and R_BC,AD respectivly
 R1 = 65
 R2 = 89
 
-# define the equation
-def fun(Rs):
-    return np.exp(-np.pi*R1/Rs) + np.exp(-np.pi*R2/Rs) - 1
-     
 # plot the fuction
-Rs = np.linspace (0, 10*np.sqrt(R1*R2), 200)
-# Rs = np.linspace (0, 10000, 200)
-plt.plot(Rs,fun(Rs))
-plt.axhline(color = "r")
+Rs = np.linspace(0, 10*np.sqrt(R1*R2), 200)
+# Rs = np.linspace(0, 10000, 200)
+plt.plot(Rs, fun(Rs))
+plt.axhline(color="r")
 plt.show()
 
 # find the solution
 root = optimize.root(fun, [1])
 
-print (root)
+print(root)