  Motivated by applications in Bayesian analysis we introduce a
multidimensional beta distribution in an ordered simplex. We study properties
of this distribution and connect them with the generalized incomplete beta
function. This function is crucial in applications of multidimensional beta
distribution, thus we present two efficient numerical algorithms for computing
the generalized incomplete beta function, one based on Taylor series expansion
and another based on Chebyshev polynomials.
