function possion(lambda) r=poissrnd(lambda,10000,1); mean(r) var(r) rmin=min(r); rmax=max(r); x=linspace(rmin,rmax,rmax-rmin+1); yy=hist(r,x); yy=yy/length(r); bar(x,yy) end