Aim: Compute complementary cumulative Poisson distribution over some range. compCumulPoisson(k) = Sum(k to infinity) P(k)
The routine declaration is:
subroutine compCumulPoisson(mu, kmin, kmax, cvf) ! Input: ! mu : average value ! kmin : minimum number of counts ! kmax : maximum number of counts ! Output: ! cvf(1:kmax-kmin+1) : complementary cumulative Poisson distribution ! from kmin to kmax real(double), intent(in) :: mu integer, intent(in) :: kmin, kmax real(double), dimension(:), intent(out) :: cvf