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