mam taki problem: potrzebuje wygenerować zmienne losowe w rozkładzie Poissona o określonej wartości oczekiwanej i określonym zakresie zmiennych losowych (od zera do n). Mam już algorytm który z ciągłej zmiennej losowej w przedziale [0,1) robi zmienną losową w rozkładzie Poissona o zadanej wartości oczekiwanej, ale nie mogę kontrolować zakresu. Operacja modulo na zwróconym wyniku dość dokładnie psuje rozkład. Mój obecny algorytm jest zaczerpnięty z książki D.Knutha:
algorithm poisson random number (Knuth):
init:
Let L ← e^( − λ), k ← 0 and p ← 1.
do:
k ← k + 1.
Generate uniform random number u and let p ← p × u.
while p ≥ L
return k - 1.