Как умножить 2 бета-функции в julia0.6.2?

Question

У меня есть две случайные величины:

B1 ~ бета (α1, β1)

B2 ~ бета (α2,β2)

Я хотел бы вычислить распределение Z = B1*B2 (произведение случайных величин)

С помощью JuliaКак я могу получить PDF одной бета-версии с:

x = collect(0:0.001:1);
b1 = Beta(11,41);
b2 = Beta(41,11);
posterior = pdf.(b1.*b2,x);

Я получаю эту ошибку MethodError: no method matching *(::Distributions.Beta{Float64}, ::Distributions.Beta{Float64})