
function y=poisson3(lambda,n)
for k=0:n
    y(k+1)=lambda^k/factorial(k)*exp(-lambda);    
end
