for each element of x, F=Prob(X<x) where X is a random variable with beta density
x.^(a-1) .* (1-x).^(b-1) ./ beta(a,b)1_{0<x<1}