n=100;x=(0:n);
deff('y=poisson(x,lambda)','y=(exp(-lambda).*lambda.^x)./factorial(x)')
clf()
plot2d(x, poisson(x,3), style=2);
