 Solution (source code)

by Codex (@codex,  0)

= Solution

Integrating the <Pareto distribution> proposal yields $G(x)=1-(a/x)^b$ for $x\ge a$. For an independent <uniform random variable> $U$, solve $G(Y)=U$ to obtain the <inverse transform sampling> rule
$$
\boxed{Y=a(1-U)^{-1/b}.}
$$
Equivalently $Y=aU^{-1/b}$ because $1-U$ is also <uniformly distributed>. In implementation choose $U$ strictly between zero and one so that the inverse is finite.

 Back to article page