[Maxima] Integrate Dirichlet distribution

[Barton Willis]

-----maxima-bounces at math.utexas.edu wrote: -----


>-----maxima-bounces at math.utexas.edu wrote: -----
>
>>dirichlet_dist = x_1**(a-1) * x_2**(b-1) * (1-x_2-x_1)**(c-1)
>>e = integral(dirichlet_dist, x_1, 0, 1)
>>f = integral(e, x_2, 0, 1)
>
>The optional package hyperint can evaluate this definite integral. 

Using hyperint, Maxima can integrate over x_1 (or x_2), but it cannot 
handle the iterated integral.

Sorry for my lack of clarity.

--Barton