k=0:10;
y=poisson3(0.1,10);
stem(k,y)
axis([-1 10 0 1])
