poisson-random* ( mean rnd -- n )

Vocabulary

random

Inputs

mean	an object
rnd	an object

Outputs

n

an object

Definition

USING: kernel math math.functions ;

IN: random

:: poisson-random* ( mean rnd -- n )
    -1 0 mean
    [ 2dup < ]
    [ [ 1 + ] 2dip [ rnd random-unit* log neg + ] dip ] while
    2drop ;