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)