一般分布的乘积

Question

假设在[0,1]之间的随机变量有两个非规范化的非参数分布，例如：

unnormalized_pdf_A = abs(sin(linspace(1,10,100)))

以及

^{pr2}$

我想获得变量C = A*B的pdf，以及它的一些基本统计数据（例如平均值、方差和/或贝叶斯区间）

用statsmodels、PyMC或scipy.stats实现这一点吗？我怎样才能开始解决这个问题？在