Sampling from the product of two betas
from pymc3 import Beta, Deterministic, Model, sample with Model(): a = Beta('a', 3, 1) b = Beta('b', 2, 5) a_times_b = Deterministic('a_times_b', a*b) trace = sample(2000)