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()