s1rout3b.m

Industrial Mathematics

% Routine 1.3b
lT = 100;
for m=90:105
 q = 0;
 p = 1;
 for k=1:m
  p = p*lT/k;
  q = q + p;
  end;
  q = q*exp(-lT);
  dE(m) = .25 - .75*q;
 m,dE(m)
end;