1
{ ******************************************************************
3
****************************************************************** }
12
function PPoisson(Mu : Float; K : Integer) : Float;
13
{ Probability of Poisson distrib. }
17
function PPoisson(Mu : Float; K : Integer) : Float;
22
if (Mu <= 0.0) or (K < 0) then
23
PPoisson := DefaultVal(FDomain, 0.0)
24
else if (- Mu) < MinLog then
25
PPoisson := DefaultVal(FUnderflow, 0.0)
27
PPoisson := DefaultVal(FOk, Exp(- Mu))
31
for I := 2 to K do { P = Mu^K / K! }
33
PPoisson := Exp(- Mu) * P;
b'\\ No newline at end of file'