Ext Erlang C P(A) fix

A-Herzog · May 19, 2024 · 0aa129e · 0aa129e
1 parent e05e0de
commit 0aa129e
Show file tree

Hide file tree

Showing 2 changed files with 84 additions and 83 deletions.
diff --git a/example_analytic_erlang_c_ext.ipynb b/example_analytic_erlang_c_ext.ipynb
diff --git a/queuesim/analytic/erlang_c_ext.py b/queuesim/analytic/erlang_c_ext.py
@@ -168,7 +168,8 @@ def PA(self) -> float:
             float: Cancelation probability P(A)
         """
         p0 = self.pn(0)
-        return sum([self.__nu / self.__l * (n - self.__c) * p0 * self.Cn(n) for n in range(self.__c + 1, self.__K + 1)])
+        input_reject = self.pn(self.__K)
+        return sum([self.__nu / (self.__l * (1 - input_reject)) * (n - self.__c) * p0 * self.Cn(n) for n in range(self.__c + 1, self.__K + 1)])
 
     def Pt(self, t: float) -> float:
         """Probability that clients have to wait t or less seconds ( = extended Erlang C formula)