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