Hello all, I submitted a PR to add the upper bounded Pareto (power-law) distribution to scipy, in the line of the existing trunc- norm, expon and weibull_min. It has a bounded support and pops up regularly when explaining scientific observations. I overloaded most methods from rv_continuous, and parameter fitting works somehow (not close to target parameters, but qualitatively good agreement between e.g. the fitted PDF and a histogram). Cheers, Nico PR: https://github.com/scipy/scipy/pull/16317 Wikipedia: https://en.wikipedia.org/wiki/Pareto_distribution#Bounded_Pareto_distributio...