Barton Willis 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.
Just to make completely clear: You're saying Maxima can't do the
iterated integral *and* there's no easy way to extend it to do so?
John.