Gammadist

Gamma probability density function
import numpy as np import scipy.stats as stats from matplotlib import pyplot as plt x = np.linspace(0,20,500) y1 = stats.gamma.pdf(x, a=5, scale=1) y2 = stats.gamma.pdf(x, a=7.5, scale=1) plt.plot(x,y1) plt.plot(x,y2,color='red') plt.show