marginal_distribution function

Fit marginal distribution