from scipy.stats import beta import numpy as np import matplotlib.pyplot as plt a=0.5 b=0.5 x=np.arange(0.01,1,0.01) y=beta.pdf(x,a,b) plt.plot(x,y) plt.title('Beta') plt.xlabel('x') plt.ylabel('density') plt.show()