I happened to encounter lots of scientific/business scenarios where a Zipf/Pareto/powerlaw describes well my data. However, whenever the mean of the distribution is large enough, the fact that these distributions have infinite variance turns into unphysical results. I then found articles mentioning "hooked powerlaw" that I guess is meant to address this. What's a mathematical formulation of it? Any python/numpy code for that?