clc; clearvars;
lambda = 0:0.01:10;
k = 0:5;
lk = length(k);
ll = length(lambda);
prob = zeros(ll,1);
dt = 5;
for n=1:lk
    for l = 1:ll
        prob(l,n)  = llegadasPoisson( lambda(l),k(n),dt );
    end   
end
figure();
plot((lambda'.*dt), prob)
grid on