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