adds B/2 term to Paasch element, fixes ECSHackWeek#251

mdmurbach · Feb 20, 2023 · 5ed337b · 5ed337b
1 parent e754f8b
commit 5ed337b
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/impedance/models/circuits/elements.py b/impedance/models/circuits/elements.py
@@ -344,7 +344,7 @@ def T(p, f):
     -----
     .. math::
 
-        Z = A\\frac{\\coth{\\beta}}{\\beta} + B\\frac{1}{\\beta\\sinh{\\beta}}
+        Z = A\\frac{\\coth{\\beta}}{\\beta} + B\\frac{1}{\\beta\\sinh{\\beta}} + \\frac{B}{2}
 
     where
 
@@ -377,7 +377,7 @@ def T(p, f):
         else:
             sinh.append(1e10)
 
-    Z = A / (beta * np.tanh(beta)) + B / (beta * np.array(sinh))
+    Z = A / (beta * np.tanh(beta)) + B / (beta * np.array(sinh)) + B / 2
     return Z