import scipy.stats x=[] for u in range(100): x1=0.01*u x.append(x1) rv = scipy.stats.beta(2,2) y = rv.pdf(x) print(y)