import scipy.stats import matplotlib.pyplot as plt x=[] for u in range(100): x1=0.01*u x.append(x1) rv = scipy.stats.beta(2,2) y = rv.pdf(x) plt.plot(x, y) plt.show()