Poisson cumulative distribution function
Returns the sum of the first k terms of the Poisson distribution: sum(exp(-m) * m**j / j!, j=0..k) = gammaincc( k+1, m). Arguments must both be positive and k an integer.
scipy.special.nrdtrisd
scipy.special.pdtrc