beta random draws
[rnd]=beta_rnd(n,a,b)
n = size of the vector of draws
a = beta distribution parameter, a = scalar
b = beta distribution parameter b = scalar
rnd = n-vector of random draws from the beta(a,b) distribution
x=beta_rnd(100,10,5)