%poisson.m
%distribuzione di Poisson
function dp=poisson(n,t)
dp=(t^n)*exp(-t)/factorial(n);