Se começarmos com , o somatório do telescópio fornece para o CDF (modificado). Inverter isso e cuidar do caso especial , fornece o seguinte algoritmo (codificado , receio, mas você pode tomá-lo como pseudocódigo para uma implementação em Python):1 - 1 / k k = 1k=21−1/kk=1R
rsoliton <- function(n.values, n=2) {
x <- runif(n.values) # Uniform values in [0,1)
i <- ceiling(1/x) # Modified soliton distribution
i[i > n] <- 1 # Convert extreme values to 1
i
}
Como exemplo de seu uso (e um teste), vamos desenhar valores para : N = 10105N=10
n.trials <- 10^5
i <- rsoliton(n.trials, n=10)
freq <- table(i) / n.trials # Tabulate frequencies
plot(freq, type="h", lwd=6)