% calculate poisson probabilities
clear all
t = 4
for j = 1:8
    f(j)=t.^(j-1).*exp(-t)/factorial(j-1);
end
f(1)
1-sum(f)