modified distribution.c

In documentation , p is the possibility of success.
so scale θ = p/(1 − p) rather than (1-p)/p
cite: http://en.wikipedia.org/wiki/Negative_binomial_distribution