function e=eplus(n)
% expected value E(X-n)
I=0:50;
P=Poisspdf(I,2);
e=sum(max(0,I-n).*P);