m:9;

P(x):=m**x/x!*exp(-m);

wxplot2d([P(x)], [x,0,18], [y,0,0.2])$

makelist(floor(P(x)*1000+0.5)/1000.0,x,0,18);

transpose(%);