https://functions.pythonprogramminghelp.com/how-to-unlock-gammasampling-distribution/