The beta density function with parameters a and b at the values of x
f=x.^(a-1) .* (1-x).^(b-1) ./ beta(a,b)1_{0<x<1}