function poisson(lambda,k)
for i=1:length(k)
    P(i) = lambda.^k(i)*exp(-lambda)/factorial(k(i));
end
ss=stem(k, P); set(ss(1),'markersize',1.5)
box off;title(sprintf('lambda=%.3f', lambda));
