function x = Pois(L)
%POIS generates a Poisson random number
n = rand;
l = 0;
u = exp(-L);
i=0;

while n > u
    i = i+1;
    l = u;
    u = u + (L^i)/factorial(i)*exp(-L);
end

x=i;