powered by
Calculation of Log multinomial Beta value.
multinomial_lbeta_function(x)
sum(lgamma(x[j])) - lgamma(sum(x))
a vector