https://operator.pythonprogramminghelp.com/think-you-know-how-to-negative-binomialsampling-distribution/