What's wrong with my code? (poisson distribution)

I tried to get a bar graph of poisson pdf without poisspdf function. For example, x=0:20; p=3.^x / factorial(x) * exp(-3); bar(x,p)

but it shows wrong graph. What is the problem?