This time with the list CCed... Sorry I keep screwing up on the CC.
I think there is an issue with function *gsl_ran_poisson_pdf* (unsigned int
k, double mu):
gsl_cdf_poisson_P(-1,5) returns 1.0000
while
gsl_cdf_poisson_P(1,5) returns 0.040428
The second number is fine, but the first one should be either equal to the
second one, or if the sign makes a difference, it should be zero.
If you turn on warnings, does the compiler shout at you for providing
a negative number (-1) as an unsigned argument (unsigned int)?
- Rhys