cdf = cdf / cdf[-1] * n / (n + 1.5e-8)