Как отмечает @Severin, гамма выглядит разумно. например:
import matplotlib.pyplot as plt
import numpy as np
import scipy.stats as sps
x = np.linspace(75)
plt.plot(x, sps.gamma.pdf(x, 4) '.-')
plt.show()
если бы им действительно нужно было сложить 1, а не интегрировать, я бы использовал cdf
, а затем np.diff
для результата