f(x)我想计算多元函数wrt 到dirichlet 分布的期望值。我的问题是“五项式”（即 5 个变量），因此计算期望值的显式形式似乎不合理。有没有办法对它进行数值积分？是否有一些包（例如在 R 中）可以让我这样做？我意识到我可以使用第一原理来解决它，但我希望计算速度非常快。

f(x) = \sum_{0,4}(q_i*log(n/q_i))